Predicting gene expression changes upon epigenomic drug treatment

Predicting gene expression changes upon... | F1000Research "use strict";function _typeof(t){return(_typeof="function"==typeof Symbol&&"symbol"==typeof Symbol.iterator?function(t){return typeof t}:function(t){return t&&"function"==typeof Symbol&&t.constructor===Symbol&&t!==Symbol.prototype?"symbol":typeof t})(t)}!function(){var t=function(){var t,e,o=[],n=window,r=n;for(;r;){try{if(r.frames.__tcfapiLocator){t=r;break}}catch(t){}if(r===n.top)break;r=r.parent}t||(!function t(){var e=n.document,o=!!n.frames.__tcfapiLocator;if(!o)if(e.body){var r=e.createElement("iframe");r.style.cssText="display:none",r.name="__tcfapiLocator",e.body.appendChild(r)}else setTimeout(t,5);return!o}(),n.__tcfapi=function(){for(var t=arguments.length,n=new Array(t),r=0;r 3&&2===parseInt(n[1],10)&&"boolean"==typeof n[3]&&(e=n[3],"function"==typeof n[2]&&n[2]("set",!0)):"ping"===n[0]?"function"==typeof n[2]&&n[2]({gdprApplies:e,cmpLoaded:!1,cmpStatus:"stub"}):o.push(n)},n.addEventListener("message",(function(t){var e="string"==typeof t.data,o={};if(e)try{o=JSON.parse(t.data)}catch(t){}else o=t.data;var n="object"===_typeof(o)&&null!==o?o.__tcfapiCall:null;n&&window.__tcfapi(n.command,n.version,(function(o,r){var a={__tcfapiReturn:{returnValue:o,success:r,callId:n.callId}};t&&t.source&&t.source.postMessage&&t.source.postMessage(e?JSON.stringify(a):a,"*")}),n.parameter)}),!1))};"undefined"!=typeof module?module.exports=t:t()}(); dataLayer = dataLayer || []; // Standard GTM initialization - Google Consent Mode handles consent automatically (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], j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src= 'https://www.googletagmanager.com/gtm.js?id='+i+dl+ '>m_auth=hzk0Vc3qFsQYhCrIoHz68A>m_preview=env-1>m_cookies_win=x';f.parentNode.insertBefore(j,f); })(window,document,'script','dataLayer','GTM-MWFK8L5J'); ;window.NREUM||(NREUM={});NREUM.init={distributed_tracing:{enabled:true},privacy:{cookies_enabled:true},ajax:{deny_list:["bam.nr-data.net"]}}; ;NREUM.loader_config={accountID:"438030",trustKey:"438030",agentID:"772317073",licenseKey:"97f8f67f26",applicationID:"772317073"} ;NREUM.info={beacon:"bam.nr-data.net",errorBeacon:"bam.nr-data.net",licenseKey:"97f8f67f26",applicationID:"772317073",sa:1} ;/*! For license information please see nr-loader-spa-1.236.0.min.js.LICENSE.txt */ (()=>{"use strict";var e,t,r={5763:(e,t,r)=>{r.d(t,{P_:()=>l,Mt:()=>g,C5:()=>s,DL:()=>v,OP:()=>T,lF:()=>D,Yu:()=>y,Dg:()=>h,CX:()=>c,GE:()=>b,sU:()=>_});var n=r(8632),i=r(9567);const o={beacon:n.ce.beacon,errorBeacon:n.ce.errorBeacon,licenseKey:void 0,applicationID:void 0,sa:void 0,queueTime:void 0,applicationTime:void 0,ttGuid:void 0,user:void 0,account:void 0,product:void 0,extra:void 0,jsAttributes:{},userAttributes:void 0,atts:void 0,transactionName:void 0,tNamePlain:void 0},a={};function s(e){if(!e)throw new Error("All info objects require an agent identifier!");if(!a[e])throw new Error("Info for ".concat(e," was never set"));return a[e]}function c(e,t){if(!e)throw new Error("All info objects require an agent identifier!");a[e]=(0,i.D)(t,o),(0,n.Qy)(e,a[e],"info")}var u=r(7056);const d=()=>{const e={blockSelector:"[data-nr-block]",maskInputOptions:{password:!0}};return{allow_bfcache:!0,privacy:{cookies_enabled:!0},ajax:{deny_list:void 0,enabled:!0,harvestTimeSeconds:10},distributed_tracing:{enabled:void 0,exclude_newrelic_header:void 0,cors_use_newrelic_header:void 0,cors_use_tracecontext_headers:void 0,allowed_origins:void 0},session:{domain:void 0,expiresMs:u.oD,inactiveMs:u.Hb},ssl:void 0,obfuscate:void 0,jserrors:{enabled:!0,harvestTimeSeconds:10},metrics:{enabled:!0},page_action:{enabled:!0,harvestTimeSeconds:30},page_view_event:{enabled:!0},page_view_timing:{enabled:!0,harvestTimeSeconds:30,long_task:!1},session_trace:{enabled:!0,harvestTimeSeconds:10},harvest:{tooManyRequestsDelay:60},session_replay:{enabled:!1,harvestTimeSeconds:60,sampleRate:.1,errorSampleRate:.1,maskTextSelector:"*",maskAllInputs:!0,get blockClass(){return"nr-block"},get ignoreClass(){return"nr-ignore"},get maskTextClass(){return"nr-mask"},get blockSelector(){return e.blockSelector},set blockSelector(t){e.blockSelector+=",".concat(t)},get maskInputOptions(){return e.maskInputOptions},set maskInputOptions(t){e.maskInputOptions={...t,password:!0}}},spa:{enabled:!0,harvestTimeSeconds:10}}},f={};function l(e){if(!e)throw new Error("All configuration objects require an agent identifier!");if(!f[e])throw new Error("Configuration for ".concat(e," was never set"));return f[e]}function h(e,t){if(!e)throw new Error("All configuration objects require an agent identifier!");f[e]=(0,i.D)(t,d()),(0,n.Qy)(e,f[e],"config")}function g(e,t){if(!e)throw new Error("All configuration objects require an agent identifier!");var r=l(e);if(r){for(var n=t.split("."),i=0;i {r.d(t,{D:()=>i});var n=r(50);function i(e,t){try{if(!e||"object"!=typeof e)return(0,n.Z)("Setting a Configurable requires an object as input");if(!t||"object"!=typeof t)return(0,n.Z)("Setting a Configurable requires a model to set its initial properties");const r=Object.create(Object.getPrototypeOf(t),Object.getOwnPropertyDescriptors(t)),o=0===Object.keys(r).length?e:r;for(let a in o)if(void 0!==e[a])try{"object"==typeof e[a]&&"object"==typeof t[a]?r[a]=i(e[a],t[a]):r[a]=e[a]}catch(e){(0,n.Z)("An error occurred while setting a property of a Configurable",e)}return r}catch(e){(0,n.Z)("An error occured while setting a Configurable",e)}}},6818:(e,t,r)=>{r.d(t,{Re:()=>i,gF:()=>o,q4:()=>n});const n="1.236.0",i="PROD",o="CDN"},385:(e,t,r)=>{r.d(t,{FN:()=>a,IF:()=>u,Nk:()=>f,Tt:()=>s,_A:()=>o,il:()=>n,pL:()=>c,v6:()=>i,w1:()=>d});const n="undefined"!=typeof window&&!!window.document,i="undefined"!=typeof WorkerGlobalScope&&("undefined"!=typeof self&&self instanceof WorkerGlobalScope&&self.navigator instanceof WorkerNavigator||"undefined"!=typeof globalThis&&globalThis instanceof WorkerGlobalScope&&globalThis.navigator instanceof WorkerNavigator),o=n?window:"undefined"!=typeof WorkerGlobalScope&&("undefined"!=typeof self&&self instanceof WorkerGlobalScope&&self||"undefined"!=typeof globalThis&&globalThis instanceof WorkerGlobalScope&&globalThis),a=""+o?.location,s=/iPad|iPhone|iPod/.test(navigator.userAgent),c=s&&"undefined"==typeof SharedWorker,u=(()=>{const e=navigator.userAgent.match(/Firefox[/\s](\d+\.\d+)/);return Array.isArray(e)&&e.length>=2?+e[1]:0})(),d=Boolean(n&&window.document.documentMode),f=!!navigator.sendBeacon},1117:(e,t,r)=>{r.d(t,{w:()=>o});var n=r(50);const i={agentIdentifier:"",ee:void 0};class o{constructor(e){try{if("object"!=typeof e)return(0,n.Z)("shared context requires an object as input");this.sharedContext={},Object.assign(this.sharedContext,i),Object.entries(e).forEach((e=>{let[t,r]=e;Object.keys(i).includes(t)&&(this.sharedContext[t]=r)}))}catch(e){(0,n.Z)("An error occured while setting SharedContext",e)}}}},8e3:(e,t,r)=>{r.d(t,{L:()=>d,R:()=>c});var n=r(2177),i=r(1284),o=r(4322),a=r(3325);const s={};function c(e,t){const r={staged:!1,priority:a.p[t]||0};u(e),s[e].get(t)||s[e].set(t,r)}function u(e){e&&(s[e]||(s[e]=new Map))}function d(){let e=arguments.length>0&&void 0!==arguments[0]?arguments[0]:"",t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:"feature";if(u(e),!e||!s[e].get(t))return a(t);s[e].get(t).staged=!0;const r=[...s[e]];function a(t){const r=e?n.ee.get(e):n.ee,a=o.X.handlers;if(r.backlog&&a){var s=r.backlog[t],c=a[t];if(c){for(var u=0;s&&u {let[t,r]=e;return r.staged}))&&(r.sort(((e,t)=>e[1].priority-t[1].priority)),r.forEach((e=>{let[t]=e;a(t)})))}function f(e,t){var r=e[1];(0,i.D)(t[r],(function(t,r){var n=e[0];if(r[0]===n){var i=r[1],o=e[3],a=e[2];i.apply(o,a)}}))}},2177:(e,t,r)=>{r.d(t,{c:()=>f,ee:()=>u});var n=r(8632),i=r(2210),o=r(1284),a=r(5763),s="nr@context";let c=(0,n.fP)();var u;function d(){}function f(e){return(0,i.X)(e,s,l)}function l(){return new d}function h(){u.aborted=!0,u.backlog={}}c.ee?u=c.ee:(u=function e(t,r){var n={},c={},f={},g=!1;try{g=16===r.length&&(0,a.OP)(r).isolatedBacklog}catch(e){}var p={on:b,addEventListener:b,removeEventListener:y,emit:v,get:x,listeners:w,context:m,buffer:A,abort:h,aborted:!1,isBuffering:E,debugId:r,backlog:g?{}:t&&"object"==typeof t.backlog?t.backlog:{}};return p;function m(e){return e&&e instanceof d?e:e?(0,i.X)(e,s,l):l()}function v(e,r,n,i,o){if(!1!==o&&(o=!0),!u.aborted||i){t&&o&&t.emit(e,r,n);for(var a=m(n),s=w(e),d=s.length,f=0;fn,p:()=>i});var n=r(2177).ee.get("handle");function i(e,t,r,i,o){o?(o.buffer([e],i),o.emit(e,t,r)):(n.buffer([e],i),n.emit(e,t,r))}},4322:(e,t,r)=>{r.d(t,{X:()=>o});var n=r(5546);o.on=a;var i=o.handlers={};function o(e,t,r,o){a(o||n.E,i,e,t,r)}function a(e,t,r,i,o){o||(o="feature"),e||(e=n.E);var a=t[o]=t[o]||{};(a[r]=a[r]||[]).push([e,i])}},3239:(e,t,r)=>{r.d(t,{bP:()=>s,iz:()=>c,m$:()=>a});var n=r(385);let i=!1,o=!1;try{const e={get passive(){return i=!0,!1},get signal(){return o=!0,!1}};n._A.addEventListener("test",null,e),n._A.removeEventListener("test",null,e)}catch(e){}function a(e,t){return i||o?{capture:!!e,passive:i,signal:t}:!!e}function s(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2],n=arguments.length>3?arguments[3]:void 0;window.addEventListener(e,t,a(r,n))}function c(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2],n=arguments.length>3?arguments[3]:void 0;document.addEventListener(e,t,a(r,n))}},4402:(e,t,r)=>{r.d(t,{Ht:()=>u,M:()=>c,Rl:()=>a,ky:()=>s});var n=r(385);const i="xxxxxxxx-xxxx-4xxx-yxxx-xxxxxxxxxxxx";function o(e,t){return e?15&e[t]:16*Math.random()|0}function a(){const e=n._A?.crypto||n._A?.msCrypto;let t,r=0;return e&&e.getRandomValues&&(t=e.getRandomValues(new Uint8Array(31))),i.split("").map((e=>"x"===e?o(t,++r).toString(16):"y"===e?(3&o()|8).toString(16):e)).join("")}function s(e){const t=n._A?.crypto||n._A?.msCrypto;let r,i=0;t&&t.getRandomValues&&(r=t.getRandomValues(new Uint8Array(31)));const a=[];for(var s=0;s {r.d(t,{Bq:()=>n,Hb:()=>o,oD:()=>i});const n="NRBA",i=144e5,o=18e5},7894:(e,t,r)=>{function n(){return Math.round(performance.now())}r.d(t,{z:()=>n})},7243:(e,t,r)=>{r.d(t,{e:()=>o});var n=r(385),i={};function o(e){if(e in i)return i[e];if(0===(e||"").indexOf("data:"))return{protocol:"data"};let t;var r=n._A?.location,o={};if(n.il)t=document.createElement("a"),t.href=e;else try{t=new URL(e,r.href)}catch(e){return o}o.port=t.port;var a=t.href.split("://");!o.port&&a[1]&&(o.port=a[1].split("/")[0].split("@").pop().split(":")[1]),o.port&&"0"!==o.port||(o.port="https"===a[0]?"443":"80"),o.hostname=t.hostname||r.hostname,o.pathname=t.pathname,o.protocol=a[0],"/"!==o.pathname.charAt(0)&&(o.pathname="/"+o.pathname);var s=!t.protocol||":"===t.protocol||t.protocol===r.protocol,c=t.hostname===r.hostname&&t.port===r.port;return o.sameOrigin=s&&(!t.hostname||c),"/"===o.pathname&&(i[e]=o),o}},50:(e,t,r)=>{function n(e,t){"function"==typeof console.warn&&(console.warn("New Relic: ".concat(e)),t&&console.warn(t))}r.d(t,{Z:()=>n})},2587:(e,t,r)=>{r.d(t,{N:()=>c,T:()=>u});var n=r(2177),i=r(5546),o=r(8e3),a=r(3325);const s={stn:[a.D.sessionTrace],err:[a.D.jserrors,a.D.metrics],ins:[a.D.pageAction],spa:[a.D.spa],sr:[a.D.sessionReplay,a.D.sessionTrace]};function c(e,t){const r=n.ee.get(t);e&&"object"==typeof e&&(Object.entries(e).forEach((e=>{let[t,n]=e;void 0===u[t]&&(s[t]?s[t].forEach((e=>{n?(0,i.p)("feat-"+t,[],void 0,e,r):(0,i.p)("block-"+t,[],void 0,e,r),(0,i.p)("rumresp-"+t,[Boolean(n)],void 0,e,r)})):n&&(0,i.p)("feat-"+t,[],void 0,void 0,r),u[t]=Boolean(n))})),Object.keys(s).forEach((e=>{void 0===u[e]&&(s[e]?.forEach((t=>(0,i.p)("rumresp-"+e,[!1],void 0,t,r))),u[e]=!1)})),(0,o.L)(t,a.D.pageViewEvent))}const u={}},2210:(e,t,r)=>{r.d(t,{X:()=>i});var n=Object.prototype.hasOwnProperty;function i(e,t,r){if(n.call(e,t))return e[t];var i=r();if(Object.defineProperty&&Object.keys)try{return Object.defineProperty(e,t,{value:i,writable:!0,enumerable:!1}),i}catch(e){}return e[t]=i,i}},1284:(e,t,r)=>{r.d(t,{D:()=>n});const n=(e,t)=>Object.entries(e||{}).map((e=>{let[r,n]=e;return t(r,n)}))},4351:(e,t,r)=>{r.d(t,{P:()=>o});var n=r(2177);const i=()=>{const e=new WeakSet;return(t,r)=>{if("object"==typeof r&&null!==r){if(e.has(r))return;e.add(r)}return r}};function o(e){try{return JSON.stringify(e,i())}catch(e){try{n.ee.emit("internal-error",[e])}catch(e){}}}},3960:(e,t,r)=>{r.d(t,{K:()=>a,b:()=>o});var n=r(3239);function i(){return"undefined"==typeof document||"complete"===document.readyState}function o(e,t){if(i())return e();(0,n.bP)("load",e,t)}function a(e){if(i())return e();(0,n.iz)("DOMContentLoaded",e)}},8632:(e,t,r)=>{r.d(t,{EZ:()=>u,Qy:()=>c,ce:()=>o,fP:()=>a,gG:()=>d,mF:()=>s});var n=r(7894),i=r(385);const o={beacon:"bam.nr-data.net",errorBeacon:"bam.nr-data.net"};function a(){return i._A.NREUM||(i._A.NREUM={}),void 0===i._A.newrelic&&(i._A.newrelic=i._A.NREUM),i._A.NREUM}function s(){let e=a();return e.o||(e.o={ST:i._A.setTimeout,SI:i._A.setImmediate,CT:i._A.clearTimeout,XHR:i._A.XMLHttpRequest,REQ:i._A.Request,EV:i._A.Event,PR:i._A.Promise,MO:i._A.MutationObserver,FETCH:i._A.fetch}),e}function c(e,t,r){let i=a();const o=i.initializedAgents||{},s=o[e]||{};return Object.keys(s).length||(s.initializedAt={ms:(0,n.z)(),date:new Date}),i.initializedAgents={...o,[e]:{...s,[r]:t}},i}function u(e,t){a()[e]=t}function d(){return function(){let e=a();const t=e.info||{};e.info={beacon:o.beacon,errorBeacon:o.errorBeacon,...t}}(),function(){let e=a();const t=e.init||{};e.init={...t}}(),s(),function(){let e=a();const t=e.loader_config||{};e.loader_config={...t}}(),a()}},7956:(e,t,r)=>{r.d(t,{N:()=>i});var n=r(3239);function i(e){let t=arguments.length>1&&void 0!==arguments[1]&&arguments[1],r=arguments.length>2?arguments[2]:void 0,i=arguments.length>3?arguments[3]:void 0;return void(0,n.iz)("visibilitychange",(function(){if(t)return void("hidden"==document.visibilityState&&e());e(document.visibilityState)}),r,i)}},1214:(e,t,r)=>{r.d(t,{em:()=>v,u5:()=>N,QU:()=>S,_L:()=>I,Gm:()=>L,Lg:()=>M,gy:()=>U,BV:()=>Q,Kf:()=>ee});var n=r(2177);const i="nr@original";var o=Object.prototype.hasOwnProperty,a=!1;function s(e,t){return e||(e=n.ee),r.inPlace=function(e,t,n,i,o){n||(n="");var a,s,c,u="-"===n.charAt(0);for(c=0;c 2?n-2:0),o=2;o {r(A[T],e,w),r(E[T],e,w)})),r(l._A,"fetch",y),t.on(y+"end",(function(e,r){var n=this;if(r){var i=r.headers.get("content-length");null!==i&&(n.rxSize=i),t.emit(y+"done",[null,r],n)}else t.emit(y+"done",[e],n)})),t}const O={},j=["pushState","replaceState"];function S(e){const t=function(e){return(e||n.ee).get("history")}(e);return!l.il||O[t.debugId]++||(O[t.debugId]=1,s(t).inPlace(window.history,j,"-")),t}var P=r(3239);const C={},R=["appendChild","insertBefore","replaceChild"];function I(e){const t=function(e){return(e||n.ee).get("jsonp")}(e);if(!l.il||C[t.debugId])return t;C[t.debugId]=!0;var r=s(t),i=/[?&](?:callback|cb)=([^&#]+)/,o=/(.*)\.([^.]+)/,a=/^(\w+)(\.|$)(.*)$/;function c(e,t){var r=e.match(a),n=r[1],i=r[3];return i?c(i,t[n]):t[n]}return r.inPlace(Node.prototype,R,"dom-"),t.on("dom-start",(function(e){!function(e){if(!e||"string"!=typeof e.nodeName||"script"!==e.nodeName.toLowerCase())return;if("function"!=typeof e.addEventListener)return;var n=(a=e.src,s=a.match(i),s?s[1]:null);var a,s;if(!n)return;var u=function(e){var t=e.match(o);if(t&&t.length>=3)return{key:t[2],parent:c(t[1],window)};return{key:e,parent:window}}(n);if("function"!=typeof u.parent[u.key])return;var d={};function f(){t.emit("jsonp-end",[],d),e.removeEventListener("load",f,(0,P.m$)(!1)),e.removeEventListener("error",l,(0,P.m$)(!1))}function l(){t.emit("jsonp-error",[],d),t.emit("jsonp-end",[],d),e.removeEventListener("load",f,(0,P.m$)(!1)),e.removeEventListener("error",l,(0,P.m$)(!1))}r.inPlace(u.parent,[u.key],"cb-",d),e.addEventListener("load",f,(0,P.m$)(!1)),e.addEventListener("error",l,(0,P.m$)(!1)),t.emit("new-jsonp",[e.src],d)}(e[0])})),t}var k=r(5763);const H={};function L(e){const t=function(e){return(e||n.ee).get("mutation")}(e);if(!l.il||H[t.debugId])return t;H[t.debugId]=!0;var r=s(t),i=k.Yu.MO;return i&&(window.MutationObserver=function(e){return this instanceof i?new i(r(e,"fn-")):i.apply(this,arguments)},MutationObserver.prototype=i.prototype),t}const z={};function M(e){const t=function(e){return(e||n.ee).get("promise")}(e);if(z[t.debugId])return t;z[t.debugId]=!0;var r=n.c,o=s(t),a=k.Yu.PR;return a&&function(){function e(r){var n=t.context(),i=o(r,"executor-",n,null,!1);const s=Reflect.construct(a,[i],e);return t.context(s).getCtx=function(){return n},s}l._A.Promise=e,Object.defineProperty(e,"name",{value:"Promise"}),e.toString=function(){return a.toString()},Object.setPrototypeOf(e,a),["all","race"].forEach((function(r){const n=a[r];e[r]=function(e){let i=!1;[...e||[]].forEach((e=>{this.resolve(e).then(a("all"===r),a(!1))}));const o=n.apply(this,arguments);return o;function a(e){return function(){t.emit("propagate",[null,!i],o,!1,!1),i=i||!e}}}})),["resolve","reject"].forEach((function(r){const n=a[r];e[r]=function(e){const r=n.apply(this,arguments);return e!==r&&t.emit("propagate",[e,!0],r,!1,!1),r}})),e.prototype=a.prototype;const n=a.prototype.then;a.prototype.then=function(){var e=this,i=r(e);i.promise=e;for(var a=arguments.length,s=new Array(a),c=0;c e())),t};function m(e,t){i.inPlace(t,["onreadystatechange"],"fn-",E)}function b(){var e=this,t=r.context(e);e.readyState>3&&!t.resolved&&(t.resolved=!0,r.emit("xhr-resolved",[],e)),i.inPlace(e,f,"fn-",E)}if(function(e,t){for(var r in e)t[r]=e[r]}(o,p),p.prototype=o.prototype,i.inPlace(p.prototype,J,"-xhr-",E),r.on("send-xhr-start",(function(e,t){m(e,t),function(e){h.push(e),a&&(y?y.then(A):u?u(A):(w=-w,x.data=w))}(t)})),r.on("open-xhr-start",m),a){var y=c&&c.resolve();if(!u&&!c){var w=1,x=document.createTextNode(w);new a(A).observe(x,{characterData:!0})}}else t.on("fn-end",(function(e){e[0]&&e[0].type===d||A()}));function A(){for(var e=0;e {r.d(t,{t:()=>n});const n=r(3325).D.ajax},6660:(e,t,r)=>{r.d(t,{A:()=>i,t:()=>n});const n=r(3325).D.jserrors,i="nr@seenError"},3081:(e,t,r)=>{r.d(t,{gF:()=>o,mY:()=>i,t9:()=>n,vz:()=>s,xS:()=>a});const n=r(3325).D.metrics,i="sm",o="cm",a="storeSupportabilityMetrics",s="storeEventMetrics"},4649:(e,t,r)=>{r.d(t,{t:()=>n});const n=r(3325).D.pageAction},7633:(e,t,r)=>{r.d(t,{Dz:()=>i,OJ:()=>a,qw:()=>o,t9:()=>n});const n=r(3325).D.pageViewEvent,i="firstbyte",o="domcontent",a="windowload"},9251:(e,t,r)=>{r.d(t,{t:()=>n});const n=r(3325).D.pageViewTiming},3614:(e,t,r)=>{r.d(t,{BST_RESOURCE:()=>i,END:()=>s,FEATURE_NAME:()=>n,FN_END:()=>u,FN_START:()=>c,PUSH_STATE:()=>d,RESOURCE:()=>o,START:()=>a});const n=r(3325).D.sessionTrace,i="bstResource",o="resource",a="-start",s="-end",c="fn"+a,u="fn"+s,d="pushState"},7836:(e,t,r)=>{r.d(t,{BODY:()=>A,CB_END:()=>E,CB_START:()=>u,END:()=>x,FEATURE_NAME:()=>i,FETCH:()=>_,FETCH_BODY:()=>v,FETCH_DONE:()=>m,FETCH_START:()=>p,FN_END:()=>c,FN_START:()=>s,INTERACTION:()=>l,INTERACTION_API:()=>d,INTERACTION_EVENTS:()=>o,JSONP_END:()=>b,JSONP_NODE:()=>g,JS_TIME:()=>T,MAX_TIMER_BUDGET:()=>a,REMAINING:()=>f,SPA_NODE:()=>h,START:()=>w,originalSetTimeout:()=>y});var n=r(5763);const i=r(3325).D.spa,o=["click","submit","keypress","keydown","keyup","change"],a=999,s="fn-start",c="fn-end",u="cb-start",d="api-ixn-",f="remaining",l="interaction",h="spaNode",g="jsonpNode",p="fetch-start",m="fetch-done",v="fetch-body-",b="jsonp-end",y=n.Yu.ST,w="-start",x="-end",A="-body",E="cb"+x,T="jsTime",_="fetch"},5938:(e,t,r)=>{r.d(t,{W:()=>o});var n=r(5763),i=r(2177);class o{constructor(e,t,r){this.agentIdentifier=e,this.aggregator=t,this.ee=i.ee.get(e,(0,n.OP)(this.agentIdentifier).isolatedBacklog),this.featureName=r,this.blocked=!1}}},9144:(e,t,r)=>{r.d(t,{j:()=>m});var n=r(3325),i=r(5763),o=r(5546),a=r(2177),s=r(7894),c=r(8e3),u=r(3960),d=r(385),f=r(50),l=r(3081),h=r(8632);function g(){const e=(0,h.gG)();["setErrorHandler","finished","addToTrace","inlineHit","addRelease","addPageAction","setCurrentRouteName","setPageViewName","setCustomAttribute","interaction","noticeError","setUserId"].forEach((t=>{e[t]=function(){for(var r=arguments.length,n=new Array(r),i=0;i 1?r-1:0),i=1;i {e.exposed&&e.api[t]&&o.push(e.api[t](...n))})),o.length>1?o:o[0]}(t,...n)}}))}var p=r(2587);function m(e){let t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:{},m=arguments.length>2?arguments[2]:void 0,v=arguments.length>3?arguments[3]:void 0,{init:b,info:y,loader_config:w,runtime:x={loaderType:m},exposed:A=!0}=t;const E=(0,h.gG)();y||(b=E.init,y=E.info,w=E.loader_config),(0,i.Dg)(e,b||{}),(0,i.GE)(e,w||{}),(0,i.sU)(e,x),y.jsAttributes??={},d.v6&&(y.jsAttributes.isWorker=!0),(0,i.CX)(e,y),g();const T=function(e,t){t||(0,c.R)(e,"api");const h={};var g=a.ee.get(e),p=g.get("tracer"),m="api-",v=m+"ixn-";function b(t,r,n,o){const a=(0,i.C5)(e);return null===r?delete a.jsAttributes[t]:(0,i.CX)(e,{...a,jsAttributes:{...a.jsAttributes,[t]:r}}),x(m,n,!0,o||null===r?"session":void 0)(t,r)}function y(){}["setErrorHandler","finished","addToTrace","inlineHit","addRelease"].forEach((e=>h[e]=x(m,e,!0,"api"))),h.addPageAction=x(m,"addPageAction",!0,n.D.pageAction),h.setCurrentRouteName=x(m,"routeName",!0,n.D.spa),h.setPageViewName=function(t,r){if("string"==typeof t)return"/"!==t.charAt(0)&&(t="/"+t),(0,i.OP)(e).customTransaction=(r||"http://custom.transaction")+t,x(m,"setPageViewName",!0)()},h.setCustomAttribute=function(e,t){let r=arguments.length>2&&void 0!==arguments[2]&&arguments[2];if("string"==typeof e){if(["string","number"].includes(typeof t)||null===t)return b(e,t,"setCustomAttribute",r);(0,f.Z)("Failed to execute setCustomAttribute.\nNon-null value must be a string or number type, but a type of was provided."))}else(0,f.Z)("Failed to execute setCustomAttribute.\nName must be a string type, but a type of was provided."))},h.setUserId=function(e){if("string"==typeof e||null===e)return b("enduser.id",e,"setUserId",!0);(0,f.Z)("Failed to execute setUserId.\nNon-null value must be a string type, but a type of was provided."))},h.interaction=function(){return(new y).get()};var w=y.prototype={createTracer:function(e,t){var r={},i=this,a="function"==typeof t;return(0,o.p)(v+"tracer",[(0,s.z)(),e,r],i,n.D.spa,g),function(){if(p.emit((a?"":"no-")+"fn-start",[(0,s.z)(),i,a],r),a)try{return t.apply(this,arguments)}catch(e){throw p.emit("fn-err",[arguments,this,"string"==typeof e?new Error(e):e],r),e}finally{p.emit("fn-end",[(0,s.z)()],r)}}}};function x(e,t,r,i){return function(){return(0,o.p)(l.xS,["API/"+t+"/called"],void 0,n.D.metrics,g),i&&(0,o.p)(e+t,[(0,s.z)(),...arguments],r?null:this,i,g),r?void 0:this}}function A(){r.e(439).then(r.bind(r,7438)).then((t=>{let{setAPI:r}=t;r(e),(0,c.L)(e,"api")})).catch((()=>(0,f.Z)("Downloading runtime APIs failed...")))}return["actionText","setName","setAttribute","save","ignore","onEnd","getContext","end","get"].forEach((e=>{w[e]=x(v,e,void 0,n.D.spa)})),h.noticeError=function(e,t){"string"==typeof e&&(e=new Error(e)),(0,o.p)(l.xS,["API/noticeError/called"],void 0,n.D.metrics,g),(0,o.p)("err",[e,(0,s.z)(),!1,t],void 0,n.D.jserrors,g)},d.il?(0,u.b)((()=>A()),!0):A(),h}(e,v);return(0,h.Qy)(e,T,"api"),(0,h.Qy)(e,A,"exposed"),(0,h.EZ)("activatedFeatures",p.T),T}},3325:(e,t,r)=>{r.d(t,{D:()=>n,p:()=>i});const n={ajax:"ajax",jserrors:"jserrors",metrics:"metrics",pageAction:"page_action",pageViewEvent:"page_view_event",pageViewTiming:"page_view_timing",sessionReplay:"session_replay",sessionTrace:"session_trace",spa:"spa"},i={[n.pageViewEvent]:1,[n.pageViewTiming]:2,[n.metrics]:3,[n.jserrors]:4,[n.ajax]:5,[n.sessionTrace]:6,[n.pageAction]:7,[n.spa]:8,[n.sessionReplay]:9}}},n={};function i(e){var t=n[e];if(void 0!==t)return t.exports;var o=n[e]={exports:{}};return r[e](o,o.exports,i),o.exports}i.m=r,i.d=(e,t)=>{for(var r in t)i.o(t,r)&&!i.o(e,r)&&Object.defineProperty(e,r,{enumerable:!0,get:t[r]})},i.f={},i.e=e=>Promise.all(Object.keys(i.f).reduce(((t,r)=>(i.f[r](e,t),t)),[])),i.u=e=>(({78:"page_action-aggregate",147:"metrics-aggregate",242:"session-manager",317:"jserrors-aggregate",348:"page_view_timing-aggregate",412:"lazy-feature-loader",439:"async-api",538:"recorder",590:"session_replay-aggregate",675:"compressor",733:"session_trace-aggregate",786:"page_view_event-aggregate",873:"spa-aggregate",898:"ajax-aggregate"}[e]||e)+"."+{78:"ac76d497",147:"3dc53903",148:"1a20d5fe",242:"2a64278a",317:"49e41428",348:"bd6de33a",412:"2f55ce66",439:"30bd804e",538:"1b18459f",590:"cf0efb30",675:"ae9f91a8",733:"83105561",786:"06482edd",860:"03a8b7a5",873:"e6b09d52",898:"998ef92b"}[e]+"-1.236.0.min.js"),i.o=(e,t)=>Object.prototype.hasOwnProperty.call(e,t),e={},t="NRBA:",i.l=(r,n,o,a)=>{if(e[r])e[r].push(n);else{var s,c;if(void 0!==o)for(var u=document.getElementsByTagName("script"),d=0;d {s.onerror=s.onload=null,clearTimeout(h);var i=e[r];if(delete e[r],s.parentNode&&s.parentNode.removeChild(s),i&&i.forEach((e=>e(n))),t)return t(n)},h=setTimeout(l.bind(null,void 0,{type:"timeout",target:s}),12e4);s.onerror=l.bind(null,s.onerror),s.onload=l.bind(null,s.onload),c&&document.head.appendChild(s)}},i.r=e=>{"undefined"!=typeof Symbol&&Symbol.toStringTag&&Object.defineProperty(e,Symbol.toStringTag,{value:"Module"}),Object.defineProperty(e,"__esModule",{value:!0})},i.j=364,i.p="https://js-agent.newrelic.com/",(()=>{var e={364:0,953:0};i.f.j=(t,r)=>{var n=i.o(e,t)?e[t]:void 0;if(0!==n)if(n)r.push(n[2]);else{var o=new Promise(((r,i)=>n=e[t]=[r,i]));r.push(n[2]=o);var a=i.p+i.u(t),s=new Error;i.l(a,(r=>{if(i.o(e,t)&&(0!==(n=e[t])&&(e[t]=void 0),n)){var o=r&&("load"===r.type?"missing":r.type),a=r&&r.target&&r.target.src;s.message="Loading chunk "+t+" failed.\n("+o+": "+a+")",s.name="ChunkLoadError",s.type=o,s.request=a,n[1](s)}}),"chunk-"+t,t)}};var t=(t,r)=>{var n,o,[a,s,c]=r,u=0;if(a.some((t=>0!==e[t]))){for(n in s)i.o(s,n)&&(i.m[n]=s[n]);if(c)c(i)}for(t&&t(r);u {i.r(o);var e=i(3325),t=i(5763);const r=Object.values(e.D);function n(e){const n={};return r.forEach((r=>{n[r]=function(e,r){return!1!==(0,t.Mt)(r,"".concat(e,".enabled"))}(r,e)})),n}var a=i(9144);var s=i(5546),c=i(385),u=i(8e3),d=i(5938),f=i(3960),l=i(50);class h extends d.W{constructor(e,t,r){let n=!(arguments.length>3&&void 0!==arguments[3])||arguments[3];super(e,t,r),this.auto=n,this.abortHandler,this.featAggregate,this.onAggregateImported,n&&(0,u.R)(e,r)}importAggregator(){let e=arguments.length>0&&void 0!==arguments[0]?arguments[0]:{};if(this.featAggregate||!this.auto)return;const r=c.il&&!0===(0,t.Mt)(this.agentIdentifier,"privacy.cookies_enabled");let n;this.onAggregateImported=new Promise((e=>{n=e}));const o=async()=>{let t;try{if(r){const{setupAgentSession:e}=await Promise.all([i.e(860),i.e(242)]).then(i.bind(i,3228));t=e(this.agentIdentifier)}}catch(e){(0,l.Z)("A problem occurred when starting up session manager. This page will not start or extend any session.",e)}try{if(!this.shouldImportAgg(this.featureName,t))return void(0,u.L)(this.agentIdentifier,this.featureName);const{lazyFeatureLoader:r}=await i.e(412).then(i.bind(i,8582)),{Aggregate:o}=await r(this.featureName,"aggregate");this.featAggregate=new o(this.agentIdentifier,this.aggregator,e),n(!0)}catch(e){(0,l.Z)("Downloading and initializing ".concat(this.featureName," failed..."),e),this.abortHandler?.(),n(!1)}};c.il?(0,f.b)((()=>o()),!0):o()}shouldImportAgg(r,n){return r!==e.D.sessionReplay||!1!==(0,t.Mt)(this.agentIdentifier,"session_trace.enabled")&&(!!n?.isNew||!!n?.state.sessionReplay)}}var g=i(7633),p=i(7894);class m extends h{static featureName=g.t9;constructor(r,n){let i=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];if(super(r,n,g.t9,i),("undefined"==typeof PerformanceNavigationTiming||c.Tt)&&"undefined"!=typeof PerformanceTiming){const n=(0,t.OP)(r);n[g.Dz]=Math.max(Date.now()-n.offset,0),(0,f.K)((()=>n[g.qw]=Math.max((0,p.z)()-n[g.Dz],0))),(0,f.b)((()=>{const t=(0,p.z)();n[g.OJ]=Math.max(t-n[g.Dz],0),(0,s.p)("timing",["load",t],void 0,e.D.pageViewTiming,this.ee)}))}this.importAggregator()}}var v=i(1117),b=i(1284);class y extends v.w{constructor(e){super(e),this.aggregatedData={}}store(e,t,r,n,i){var o=this.getBucket(e,t,r,i);return o.metrics=function(e,t){t||(t={count:0});return t.count+=1,(0,b.D)(e,(function(e,r){t[e]=w(r,t[e])})),t}(n,o.metrics),o}merge(e,t,r,n,i){var o=this.getBucket(e,t,n,i);if(o.metrics){var a=o.metrics;a.count+=r.count,(0,b.D)(r,(function(e,t){if("count"!==e){var n=a[e],i=r[e];i&&!i.c?a[e]=w(i.t,n):a[e]=function(e,t){if(!t)return e;t.c||(t=x(t.t));return t.min=Math.min(e.min,t.min),t.max=Math.max(e.max,t.max),t.t+=e.t,t.sos+=e.sos,t.c+=e.c,t}(i,a[e])}}))}else o.metrics=r}storeMetric(e,t,r,n){var i=this.getBucket(e,t,r);return i.stats=w(n,i.stats),i}getBucket(e,t,r,n){this.aggregatedData[e]||(this.aggregatedData[e]={});var i=this.aggregatedData[e][t];return i||(i=this.aggregatedData[e][t]={params:r||{}},n&&(i.custom=n)),i}get(e,t){return t?this.aggregatedData[e]&&this.aggregatedData[e][t]:this.aggregatedData[e]}take(e){for(var t={},r="",n=!1,i=0;i t.max&&(t.max=e),e 2&&void 0!==arguments[2])||arguments[2];super(e,r,j.t,n),c.il&&((0,t.OP)(e).initHidden=Boolean("hidden"===document.visibilityState),(0,N.N)((()=>(0,s.p)("docHidden",[(0,p.z)()],void 0,j.t,this.ee)),!0),(0,O.bP)("pagehide",(()=>(0,s.p)("winPagehide",[(0,p.z)()],void 0,j.t,this.ee))),this.importAggregator())}}var P=i(3081);class C extends h{static featureName=P.t9;constructor(e,t){let r=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];super(e,t,P.t9,r),this.importAggregator()}}var R,I=i(2210),k=i(1214),H=i(2177),L={};try{R=localStorage.getItem("__nr_flags").split(","),console&&"function"==typeof console.log&&(L.console=!0,-1!==R.indexOf("dev")&&(L.dev=!0),-1!==R.indexOf("nr_dev")&&(L.nrDev=!0))}catch(e){}function z(e){try{L.console&&z(e)}catch(e){}}L.nrDev&&H.ee.on("internal-error",(function(e){z(e.stack)})),L.dev&&H.ee.on("fn-err",(function(e,t,r){z(r.stack)})),L.dev&&(z("NR AGENT IN DEVELOPMENT MODE"),z("flags: "+(0,b.D)(L,(function(e,t){return e})).join(", ")));var M=i(6660);class B extends h{static featureName=M.t;constructor(r,n){let i=!(arguments.length>2&&void 0!==arguments[2])||arguments[2];super(r,n,M.t,i),this.skipNext=0;try{this.removeOnAbort=new AbortController}catch(e){}const o=this;o.ee.on("fn-start",(function(e,t,r){o.abortHandler&&(o.skipNext+=1)})),o.ee.on("fn-err",(function(t,r,n){o.abortHandler&&!n[M.A]&&((0,I.X)(n,M.A,(function(){return!0})),this.thrown=!0,(0,s.p)("err",[n,(0,p.z)()],void 0,e.D.jserrors,o.ee))})),o.ee.on("fn-end",(function(){o.abortHandler&&!this.thrown&&o.skipNext>0&&(o.skipNext-=1)})),o.ee.on("internal-error",(function(t){(0,s.p)("ierr",[t,(0,p.z)(),!0],void 0,e.D.jserrors,o.ee)})),this.origOnerror=c._A.onerror,c._A.onerror=this.onerrorHandler.bind(this),c._A.addEventListener("unhandledrejection",(t=>{const r=function(e){let t="Unhandled Promise Rejection: ";if(e instanceof Error)try{return e.message=t+e.message,e}catch(t){return e}if(void 0===e)return new Error(t);try{return new Error(t+(0,D.P)(e))}catch(e){return new Error(t)}}(t.reason);(0,s.p)("err",[r,(0,p.z)(),!1,{unhandledPromiseRejection:1}],void 0,e.D.jserrors,this.ee)}),(0,O.m$)(!1,this.removeOnAbort?.signal)),(0,k.gy)(this.ee),(0,k.BV)(this.ee),(0,k.em)(this.ee),(0,t.OP)(r).xhrWrappable&&(0,k.Kf)(this.ee),this.abortHandler=this.#e,this.importAggregator()}#e(){this.removeOnAbort?.abort(),this.abortHandler=void 0}onerrorHandler(t,r,n,i,o){"function"==typeof this.origOnerror&&this.origOnerror(...arguments);try{this.skipNext?this.skipNext-=1:(0,s.p)("err",[o||new F(t,r,n),(0,p.z)()],void 0,e.D.jserrors,this.ee)}catch(t){try{(0,s.p)("ierr",[t,(0,p.z)(),!0],void 0,e.D.jserrors,this.ee)}catch(e){}}return!1}}function F(e,t,r){this.message=e||"Uncaught error with no additional information",this.sourceURL=t,this.line=r}let U=1;const q="nr@id";function G(e){const t=typeof e;return!e||"object"!==t&&"function"!==t?-1:e===c._A?0:(0,I.X)(e,q,(function(){return U++}))}function V(e){if("string"==typeof e&&e.length)return e.length;if("object"==typeof e){if("undefined"!=typeof ArrayBuffer&&e instanceof ArrayBuffer&&e.byteLength)return e.byteLength;if("undefined"!=typeof Blob&&e instanceof Blob&&e.size)return e.size;if(!("undefined"!=typeof FormData&&e instanceof FormData))try{return(0,D.P)(e).length}catch(e){return}}}var X=i(7243);class W{constructor(e){this.agentIdentifier=e,this.generateTracePayload=this.generateTracePayload.bind(this),this.shouldGenerateTrace=this.shouldGenerateTrace.bind(this)}generateTracePayload(e){if(!this.shouldGenerateTrace(e))return null;var r=(0,t.DL)(this.agentIdentifier);if(!r)return null;var n=(r.accountID||"").toString()||null,i=(r.agentID||"").toString()||null,o=(r.trustKey||"").toString()||null;if(!n||!i)return null;var a=(0,_.M)(),s=(0,_.Ht)(),c=Date.now(),u={spanId:a,traceId:s,timestamp:c};return(e.sameOrigin||this.isAllowedOrigin(e)&&this.useTraceContextHeadersForCors())&&(u.traceContextParentHeader=this.generateTraceContextParentHeader(a,s),u.traceContextStateHeader=this.generateTraceContextStateHeader(a,c,n,i,o)),(e.sameOrigin&&!this.excludeNewrelicHeader()||!e.sameOrigin&&this.isAllowedOrigin(e)&&this.useNewrelicHeaderForCors())&&(u.newrelicHeader=this.generateTraceHeader(a,s,c,n,i,o)),u}generateTraceContextParentHeader(e,t){return"00-"+t+"-"+e+"-01"}generateTraceContextStateHeader(e,t,r,n,i){return i+"@nr=0-1-"+r+"-"+n+"-"+e+"----"+t}generateTraceHeader(e,t,r,n,i,o){if(!("function"==typeof c._A?.btoa))return null;var a={v:[0,1],d:{ty:"Browser",ac:n,ap:i,id:e,tr:t,ti:r}};return o&&n!==o&&(a.d.tk=o),btoa((0,D.P)(a))}shouldGenerateTrace(e){return this.isDtEnabled()&&this.isAllowedOrigin(e)}isAllowedOrigin(e){var r=!1,n={};if((0,t.Mt)(this.agentIdentifier,"distributed_tracing")&&(n=(0,t.P_)(this.agentIdentifier).distributed_tracing),e.sameOrigin)r=!0;else if(n.allowed_origins instanceof Array)for(var i=0;i 2&&void 0!==arguments[2])||arguments[2];super(r,n,Z.t,i),(0,t.OP)(r).xhrWrappable&&(this.dt=new W(r),this.handler=(e,t,r,n)=>(0,s.p)(e,t,r,n,this.ee),(0,k.u5)(this.ee),(0,k.Kf)(this.ee),function(r,n,i,o){function a(e){var t=this;t.totalCbs=0,t.called=0,t.cbTime=0,t.end=E,t.ended=!1,t.xhrGuids={},t.lastSize=null,t.loadCaptureCalled=!1,t.params=this.params||{},t.metrics=this.metrics||{},e.addEventListener("load",(function(r){_(t,e)}),(0,O.m$)(!1)),c.IF||e.addEventListener("progress",(function(e){t.lastSize=e.loaded}),(0,O.m$)(!1))}function s(e){this.params={method:e[0]},T(this,e[1]),this.metrics={}}function u(e,n){var i=(0,t.DL)(r);i.xpid&&this.sameOrigin&&n.setRequestHeader("X-NewRelic-ID",i.xpid);var a=o.generateTracePayload(this.parsedOrigin);if(a){var s=!1;a.newrelicHeader&&(n.setRequestHeader("newrelic",a.newrelicHeader),s=!0),a.traceContextParentHeader&&(n.setRequestHeader("traceparent",a.traceContextParentHeader),a.traceContextStateHeader&&n.setRequestHeader("tracestate",a.traceContextStateHeader),s=!0),s&&(this.dt=a)}}function d(e,t){var r=this.metrics,i=e[0],o=this;if(r&&i){var a=V(i);a&&(r.txSize=a)}this.startTime=(0,p.z)(),this.listener=function(e){try{"abort"!==e.type||o.loadCaptureCalled||(o.params.aborted=!0),("load"!==e.type||o.called===o.totalCbs&&(o.onloadCalled||"function"!=typeof t.onload)&&"function"==typeof o.end)&&o.end(t)}catch(e){try{n.emit("internal-error",[e])}catch(e){}}};for(var s=0;s 1?e[1]=i:e.push(i)}else e[0]&&e[0].headers&&s(e[0].headers,n)&&(this.dt=n);function s(e,t){var r=!1;return t.newrelicHeader&&(e.set("newrelic",t.newrelicHeader),r=!0),t.traceContextParentHeader&&(e.set("traceparent",t.traceContextParentHeader),t.traceContextStateHeader&&e.set("tracestate",t.traceContextStateHeader),r=!0),r}}function x(e,t){this.params={},this.metrics={},this.startTime=(0,p.z)(),this.dt=t,e.length>=1&&(this.target=e[0]),e.length>=2&&(this.opts=e[1]);var r,n=this.opts||{},i=this.target;"string"==typeof i?r=i:"object"==typeof i&&i instanceof Y?r=i.url:c._A?.URL&&"object"==typeof i&&i instanceof URL&&(r=i.href),T(this,r);var o=(""+(i&&i instanceof Y&&i.method||n.method||"GET")).toUpperCase();this.params.method=o,this.txSize=V(n.body)||0}function A(t,r){var n;this.endTime=(0,p.z)(),this.params||(this.params={}),this.params.status=r?r.status:0,"string"==typeof this.rxSize&&this.rxSize.length>0&&(n=+this.rxSize);var o={txSize:this.txSize,rxSize:n,duration:(0,p.z)()-this.startTime};i("xhr",[this.params,o,this.startTime,this.endTime,"fetch"],this,e.D.ajax)}function E(t){var r=this.params,n=this.metrics;if(!this.ended){this.ended=!0;for(var o=0;o 2&&void 0!==arguments[2])||arguments[2];super(e,t,we.t,r),this.importAggregator()}}new class{constructor(e){let t=arguments.length>1&&void 0!==arguments[1]?arguments[1]:(0,_.ky)(16);c._A?(this.agentIdentifier=t,this.sharedAggregator=new y({agentIdentifier:this.agentIdentifier}),this.features={},this.desiredFeatures=new Set(e.features||[]),this.desiredFeatures.add(m),Object.assign(this,(0,a.j)(this.agentIdentifier,e,e.loaderType||"agent")),this.start()):(0,l.Z)("Failed to initial the agent. Could not determine the runtime environment.")}get config(){return{info:(0,t.C5)(this.agentIdentifier),init:(0,t.P_)(this.agentIdentifier),loader_config:(0,t.DL)(this.agentIdentifier),runtime:(0,t.OP)(this.agentIdentifier)}}start(){const t="features";try{const r=n(this.agentIdentifier),i=[...this.desiredFeatures];i.sort(((t,r)=>e.p[t.featureName]-e.p[r.featureName])),i.forEach((t=>{if(r[t.featureName]||t.featureName===e.D.pageViewEvent){const n=function(t){switch(t){case e.D.ajax:return[e.D.jserrors];case e.D.sessionTrace:return[e.D.ajax,e.D.pageViewEvent];case e.D.sessionReplay:return[e.D.sessionTrace];case e.D.pageViewTiming:return[e.D.pageViewEvent];default:return[]}}(t.featureName);n.every((e=>r[e]))||(0,l.Z)("".concat(t.featureName," is enabled but one or more dependent features has been disabled (").concat((0,D.P)(n),"). This may cause unintended consequences or missing data...")),this.features[t.featureName]=new t(this.agentIdentifier,this.sharedAggregator)}})),(0,T.Qy)(this.agentIdentifier,this.features,t)}catch(e){(0,l.Z)("Failed to initialize all enabled instrument classes (agent aborted) -",e);for(const e in this.features)this.features[e].abortHandler?.();const r=(0,T.fP)();return delete r.initializedAgents[this.agentIdentifier]?.api,delete r.initializedAgents[this.agentIdentifier]?.[t],delete this.sharedAggregator,r.ee?.abort(),delete r.ee?.get(this.agentIdentifier),!1}}}({features:[J,m,S,class extends h{static featureName=oe;constructor(t,r){if(super(t,r,oe,!(arguments.length>2&&void 0!==arguments[2])||arguments[2]),!c.il)return;const n=this.ee;let i;(0,k.QU)(n),this.eventsEE=(0,k.em)(n),this.eventsEE.on(se,(function(e,t){this.bstStart=(0,p.z)()})),this.eventsEE.on(ae,(function(t,r){(0,s.p)("bst",[t[0],r,this.bstStart,(0,p.z)()],void 0,e.D.sessionTrace,n)})),n.on(ce+ne,(function(e){this.time=(0,p.z)(),this.startPath=location.pathname+location.hash})),n.on(ce+ie,(function(t){(0,s.p)("bstHist",[location.pathname+location.hash,this.startPath,this.time],void 0,e.D.sessionTrace,n)}));try{i=new PerformanceObserver((t=>{const r=t.getEntries();(0,s.p)(te,[r],void 0,e.D.sessionTrace,n)})),i.observe({type:re,buffered:!0})}catch(e){}this.importAggregator({resourceObserver:i})}},C,xe,B,class extends h{static featureName=de;constructor(e,r){if(super(e,r,de,!(arguments.length>2&&void 0!==arguments[2])||arguments[2]),!c.il)return;if(!(0,t.OP)(e).xhrWrappable)return;try{this.removeOnAbort=new AbortController}catch(e){}let n,i=0;const o=this.ee.get("tracer"),a=(0,k._L)(this.ee),s=(0,k.Lg)(this.ee),u=(0,k.BV)(this.ee),d=(0,k.Kf)(this.ee),f=this.ee.get("events"),l=(0,k.u5)(this.ee),h=(0,k.QU)(this.ee),g=(0,k.Gm)(this.ee);function m(e,t){h.emit("newURL",[""+window.location,t])}function v(){i++,n=window.location.hash,this[ve]=(0,p.z)()}function b(){i--,window.location.hash!==n&&m(0,!0);var e=(0,p.z)();this[pe]=~~this[pe]+e-this[ve],this[ye]=e}function y(e,t){e.on(t,(function(){this[t]=(0,p.z)()}))}this.ee.on(ve,v),s.on(be,v),a.on(be,v),this.ee.on(ye,b),s.on(ge,b),a.on(ge,b),this.ee.buffer([ve,ye,"xhr-resolved"],this.featureName),f.buffer([ve],this.featureName),u.buffer(["setTimeout"+le,"clearTimeout"+fe,ve],this.featureName),d.buffer([ve,"new-xhr","send-xhr"+fe],this.featureName),l.buffer([me+fe,me+"-done",me+he+fe,me+he+le],this.featureName),h.buffer(["newURL"],this.featureName),g.buffer([ve],this.featureName),s.buffer(["propagate",be,ge,"executor-err","resolve"+fe],this.featureName),o.buffer([ve,"no-"+ve],this.featureName),a.buffer(["new-jsonp","cb-start","jsonp-error","jsonp-end"],this.featureName),y(l,me+fe),y(l,me+"-done"),y(a,"new-jsonp"),y(a,"jsonp-end"),y(a,"cb-start"),h.on("pushState-end",m),h.on("replaceState-end",m),window.addEventListener("hashchange",m,(0,O.m$)(!0,this.removeOnAbort?.signal)),window.addEventListener("load",m,(0,O.m$)(!0,this.removeOnAbort?.signal)),window.addEventListener("popstate",(function(){m(0,i>1)}),(0,O.m$)(!0,this.removeOnAbort?.signal)),this.abortHandler=this.#e,this.importAggregator()}#e(){this.removeOnAbort?.abort(),this.abortHandler=void 0}}],loaderType:"spa"})})(),window.NRBA=o})(); window.jQuery || document.write(' ') CKEDITOR_BASEPATH='https://f1000research.com/js/vendor/ckeditor/' window.reactTheme = 'research'; window.MathJax = { CommonHTML: { linebreaks: { automatic: true } }, 'HTML-CSS': { linebreaks: { automatic: true } }, SVG: { linebreaks: { automatic: true } }, AuthorInit: function() { MathJax.Hub.Register.MessageHook('End Process', function () { let timeout = false; // holder for timeout id const delay = 250; // delay after event is "complete" to run callback const reflowMath = function() { const dispFormulas = document.querySelectorAll('.disp-formula.panel'); if (!dispFormulas) { return; } for (const dispFormula of dispFormulas) { const child = dispFormula.querySelector('.MathJax_Preview').nextSibling.firstChild; const isMultiline = MathJax.Hub.getAllJax(dispFormula)[0].root.isMultiline; if (dispFormula.offsetWidth < child.offsetWidth || isMultiline) { MathJax.Hub.Queue(['Rerender', MathJax.Hub, dispFormula]); } } }; window.addEventListener('resize', function() { clearTimeout(timeout); // clear the timeout timeout = setTimeout(reflowMath, delay); // start timing for event "completion" }); }); }, }; if (window.location.hash == '#_=_'){ window.location = window.location.href.split('#')[0] } !function(f,b,e,v,n,t,s){if(f.fbq)return;n=f.fbq=function() {n.callMethod? n.callMethod.apply(n,arguments):n.queue.push(arguments)} ;if(!f._fbq)f._fbq=n; n.push=n;n.loaded=!0;n.version='2.0';n.queue=[];t=b.createElement(e);t.async=!0; t.src=v;s=b.getElementsByTagName(e)[0];s.parentNode.insertBefore(t,s)}(window, document,'script','https://connect.facebook.net/en_US/fbevents.js'); fbq('init', '1641728616063202'); fbq('track', "PixelInitialized", {}); (function(h,o,t,j,a,r){ h.hj=h.hj||function(){(h.hj.q=h.hj.q||[]).push(arguments)}; h._hjSettings={hjid:2318163,hjsv:6}; a=o.getElementsByTagName('head')[0]; r=o.createElement('script');r.async=1; r.src=t+h._hjSettings.hjid+j+h._hjSettings.hjsv; a.appendChild(r); })(window,document,'https://static.hotjar.com/c/hotjar-','.js?sv='); search file_upload Submit your research search menu close search Browse Gateways & Collections How to Publish Submit your Research My Submissions Article Guidelines Article Guidelines (New Versions) Open Data, Software and Code Guidelines Open Data and Accessible Source Materials Guidelines (HSS) Open Data, Software and Code Guidelines (PSE) Prepublication Checks Production Process Posters and Slides Guidelines Document Guidelines Article Processing Charges Peer Review Finding Article Reviewers About How it Works For Reviewers Our Advisors Policies Glossary FAQs For Developers Newsroom Contact My Research Submissions Content and Tracking Alerts My Details Sign In file_upload Submit your research { "@context": "https://schema.org", "@type": "ScholarlyArticle", "mainEntityOfPage": { "@type": "WebPage", "@id": "https://f1000research.com/articles/12-1089" }, "headline": "Predicting gene expression changes upon epigenomic drug treatment", "datePublished": "2023-09-01T13:02:01", "dateModified": "2025-05-02T10:35:13", "author": [ { "@type": "Person", "name": "Piyush Agrawal" }, { "@type": "Person", "name": "Vishaka Gopalan" }, { "@type": "Person", "name": "Monjura Afrin Rumi" }, { "@type": "Person", "name": "Sridhar Hannenhalli" } ], "publisher": { "@type": "Organization", "name": "F1000Research", "logo": { "@type": "ImageObject", "url": "https://f1000research.com/img/AMP/F1000Research_image.png", "height": 480, "width": 60 } }, "image": { "@type": "ImageObject", "url": "https://f1000research.com/img/AMP/F1000Research_image.png", "height": 1200, "width": 150 }, "description": " Background Tumors are characterized by global changes in epigenetic modifications such as DNA methylation and histone modifications that are functionally linked to tumor progression. Accordingly, several drugs targeting the epigenome have been proposed for cancer therapy, notably, histone deacetylase inhibitors (HDACi) such as vorinostat and DNA methyltransferase inhibitors (DNMTi) such as zebularine. However, a fundamental challenge with such approaches is the lack of genomic specificity, i.e., the transcriptional changes at different genomic loci can be highly variable, thus making it difficult to predict the consequences on the global transcriptome and drug response. For instance, treatment with DNMTi may upregulate the expression of not only a tumor suppressor but also an oncogene, leading to unintended adverse effect. Methods Given the pre-treatment transcriptome and epigenomic profile of a sample, we assessed the extent of predictability of locus-specific changes in gene expression upon treatment with HDACi using machine learning. Results We found that in two cell lines (HCT116 treated with Largazole at eight doses and RH4 treated with Entinostat at 1µM) where the appropriate data (pre-treatment transcriptome and epigenome as well as post-treatment transcriptome) is available, our model distinguished the post-treatment up versus downregulated genes with high accuracy (up to ROC of 0.89). Furthermore, a model trained on one cell line is applicable to another cell line suggesting generalizability of the model. Conclusions Here we present a first assessment of the predictability of genome-wide transcriptomic changes upon treatment with HDACi. Lack of appropriate omics data from clinical trials of epigenetic drugs currently hampers the assessment of applicability of our approach in clinical setting. " } { "@context": "http://schema.org", "@type": "BreadcrumbList", "itemListElement": [ { "@type": "ListItem", "position": "1", "item": { "@id": "https://f1000research.com/", "name": "Home" } }, { "@type": "ListItem", "position": "2", "item": { "@id": "https://f1000research.com/browse/articles", "name": "Browse" } }, { "@type": "ListItem", "position": "3", "item": { "@id": "https://f1000research.com/articles/12-1089/v3", "name": "Predicting gene expression changes upon epigenomic drug treatment" } } ] } Home Browse Predicting gene expression changes upon epigenomic drug treatment ALL Metrics - Views Downloads Get PDF Get XML Cite How to cite this article Agrawal P, Gopalan V, Rumi MA and Hannenhalli S. Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.12688/f1000research.140273.3 ) NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article. Close Copy Citation Details Export Export Citation Sciwheel EndNote Ref. Manager Bibtex ProCite Sente EXPORT Select a format first Track Share ▬ ✚ Research Article Revised Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] Piyush Agrawal https://orcid.org/0000-0003-2075-1111 1,2 , Vishaka Gopalan 2 , Monjura Afrin Rumi https://orcid.org/0009-0002-7353-8058 3 , Sridhar Hannenhalli 2 Piyush Agrawal https://orcid.org/0000-0003-2075-1111 1,2 , Vishaka Gopalan 2 , Monjura Afrin Rumi https://orcid.org/0009-0002-7353-8058 3 , Sridhar Hannenhalli 2 PUBLISHED 02 May 2025 Author details Author details 1 Division of Medical Research, SRM Medical College Hospital & Research Centre, SRMIST, Kattankulathur, Chennai, Tamil Nadu, India 2 Cancer Data Science Laboratory, National Cancer Institute, Bethesda, Maryland, 20814, USA 3 Faculty research assistant, Center for Bioinformatics and Computational Biology, University of Maryland, College Park, MD, 20742, USA Piyush Agrawal Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Supervision, Validation, Writing – Original Draft Preparation, Writing – Review & Editing Vishaka Gopalan Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Writing – Original Draft Preparation Monjura Afrin Rumi Roles: Data Curation, Formal Analysis, Methodology Sridhar Hannenhalli Roles: Conceptualization, Formal Analysis, Methodology, Supervision, Writing – Original Draft Preparation, Writing – Review & Editing OPEN PEER REVIEW DETAILS REVIEWER STATUS This article is included in the Genomics and Genetics gateway. This article is included in the Bioinformatics gateway. This article is included in the Bioinformatics in Cancer Research collection. Abstract Background Tumors are characterized by global changes in epigenetic modifications such as DNA methylation and histone modifications that are functionally linked to tumor progression. Accordingly, several drugs targeting the epigenome have been proposed for cancer therapy, notably, histone deacetylase inhibitors (HDACi) such as vorinostat and DNA methyltransferase inhibitors (DNMTi) such as zebularine. However, a fundamental challenge with such approaches is the lack of genomic specificity, i.e. , the transcriptional changes at different genomic loci can be highly variable, thus making it difficult to predict the consequences on the global transcriptome and drug response. For instance, treatment with DNMTi may upregulate the expression of not only a tumor suppressor but also an oncogene, leading to unintended adverse effect. Methods Given the pre-treatment transcriptome and epigenomic profile of a sample, we assessed the extent of predictability of locus-specific changes in gene expression upon treatment with HDACi using machine learning. Results We found that in two cell lines (HCT116 treated with Largazole at eight doses and RH4 treated with Entinostat at 1µM) where the appropriate data (pre-treatment transcriptome and epigenome as well as post-treatment transcriptome) is available, our model distinguished the post-treatment up versus downregulated genes with high accuracy (up to ROC of 0.89). Furthermore, a model trained on one cell line is applicable to another cell line suggesting generalizability of the model. Conclusions Here we present a first assessment of the predictability of genome-wide transcriptomic changes upon treatment with HDACi. Lack of appropriate omics data from clinical trials of epigenetic drugs currently hampers the assessment of applicability of our approach in clinical setting. READ ALL READ LESS Keywords Epigenetics, HDACi, DNMTi, Cancer therapy, Machine Learning, Transcriptomics Corresponding Author(s) Piyush Agrawal ( [email protected] ) Sridhar Hannenhalli ( [email protected] ) Close Corresponding authors: Piyush Agrawal, Sridhar Hannenhalli Competing interests: No competing interests were disclosed. Grant information: This work was supported by the National Cancer Institute (Intramural Research funding). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Copyright: © 2025 Agrawal P et al . This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The author(s) is/are employees of the US Government and therefore domestic copyright protection in USA does not apply to this work. The work may be protected under the copyright laws of other jurisdictions when used in those jurisdictions. How to cite: Agrawal P, Gopalan V, Rumi MA and Hannenhalli S. Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.12688/f1000research.140273.3 ) First published: 01 Sep 2023, 12 :1089 ( https://doi.org/10.12688/f1000research.140273.1 ) Latest published: 02 May 2025, 12 :1089 ( https://doi.org/10.12688/f1000research.140273.3 ) Revised Amendments from Version 2 In the revised manuscript, we have updated the gene ontology section providing an explanation about why genes are associated with different types of GO terms for different cell lines. Also, we have rewritten the results associated with the Figure 6. In the revised manuscript, we have updated the gene ontology section providing an explanation about why genes are associated with different types of GO terms for different cell lines. Also, we have rewritten the results associated with the Figure 6. See the authors' detailed response to the review by Sukhen Das Mandal See the authors' detailed response to the review by Angelika Merkel READ REVIEWER RESPONSES Introduction Phenotypic state of a cell or tissue, either in normal homeostasis or in disease such as cancer, is intimately linked to its transcriptional state, which in turn is profoundly determined by its global epigenome. 1 Cancer genomes display a substantially altered epigenome relative to their non-malignant counterparts. For instance, global DNA hypomethylation and focal hypermethylation, notably, at tumor suppressor gene promoters, have been noted as a general feature of many cancers. 2 Accordingly, drugs that alter the epigenome have emerged as potential candidates for cancer therapy. 3 Two of the most common classes of epigenome-altering drugs are DNA methyltransferase inhibitors (DNMTi), such as zebularine, and histone deacetylase inhibitors (HADCi), such as vorinostat. While DNMTi’s are standard of care in some hematological malignancies, across most cancers, the efficacy of epigenomic drugs has been mixed. One of the many reasons adversely affecting the success of epigenomic drugs is its lack of locus specificity. Note that the intent of DNMTi, for instance, is partly to reactivate aberrantly silenced genes by demethylating their (aberrantly methylated) promoters. 3 , 4 However, the drug is locus-agnostic and, a priori , can activate many other loci in the genome, some of which may have toxic side effects 5 and in the worst case, have pro-tumor effects; indeed, DNMTi’s are known to activate cancer testis antigens, which are known to be pro-tumor. Currently, we lack the sufficient knowledge to predict locus-specific effects of an epigenomic drug on the gene expression, to be able to develop more rational therapies. One of the facts complicating this understanding is incompletely understood interactions between different epigenomic marks. For instance, there is broad antagonism between two key mechanisms of transcriptional suppression, namely, histone modification H3K27me3 and DNA methylation 6 which may result in redistribution of one upon perturbation of the other. Because of the complex interactions between epigenomic marks as well as feed-forward and feed-back loops between the epigenome and the transcriptome, the ultimate effect of epigenomic perturbation on the global transcriptome may not be easily predictable, especially based only on the local genomic context. In previous studies, researchers tried to predict gene expression from the epigenome in a specific context. 7 – 9 However, we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. Hence, we believe, our work could substantially help assess the clinical efficacy of an epigenomic drug. Methods Processing of RH4 and HCT116 sequencing data FASTQ files were downloaded from Sequence Read Archive (SRA) (HCT116 dataset accession: SRP113250, RH4 accession: SRP151465). HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. We uniformly re-processed both RNA-seq and H3K27ac ChIP-seq data to minimize biases. We ran the fastqc toolkit (v0.11.9) to ensure quality. Trimgalore (v0.6.7) was run with default options to trim any adaptor sequence contamination in reads. For ChIP-seq data, bwa-mem2 (v2.2.1) 10 was used to align trimmed reads while salmon (v1.7.0) 11 was run to align trimmed reads with the –validateMappings option enabled. For ChIP-seq data, the read counts in each genomic bin (defined below) were normalized to TPM (transcripts per million) scale with genomic bin counts quantified using the Rsubread package. 12 Since the RH4 ChIP-seq data contained spike-in reads from the Drosophila melanogaster genome, bwa-mem2 was used to align reads using a joint BWA index of the hg38 and dm6 genome. Thus, for the RH4 data, in addition to the library size normalization that is applied to each sample, we additionally divided the TPM values by the total number of reads that aligned to the dm6 genome assembly following the recommendation in the source publication describing the Drosophila melanogaster spike-in protocol for ChIP-seq data. 13 All pseudo-aligned RNA-seq data from salmon was normalized to a TPM scale using the tximport function (v1.28.0). Distribution of histone marks in the genic region We identified 1000 most upregulated and downregulated genes post-drug treatment for HCT116 and RH4 cell line. The genes were selected based on Log2 Fold change i.e., Log2 Treated – Log2 Untreated. For every gene, we created 21 genomic bins to analyze the pattern of histone marks. The genomic bins include promoter region, transcription Start Site (TSS), and Gene Body (GB) region. TSS coordinates were obtained from the ENSEMBL Genes v101 database. 14 We defined the promoter as the 2kb region upstream to the TSS which was further divided into 10 equal-sized bins where the TSS was the single nucleotide position. Finally, the gene body was defined as the entire transcribed region and was also divided into 10 equal-sized bins. Overall, this resulted in a total of 21 bins for every gene. H3K27Ac read density was calculated in each of these 21 bins and was used to compare the up and down genes and as features for the prediction of up and down-regulated genes. Machine learning model to predict post-treatment transcriptional effect We used the histone mark distribution in 21 genic bins as features to develop machine learning models to distinguish up versus downregulated genes after HDACi-treatment separately for both HCT116 and RH4 cell lines. Using the conventional five-fold cross-validation we computed the area under curve (AUC) as performance measure. We used a python-based library known as Scikit-learn 15 and implemented three different machine learning techniques which include Support Vector Machine (SVM), Random Forest (RF), and Gradient Boosting. Models were developed in four different categories (i) using 10 Promoter features; (ii) using single TSS feature; (iii) using 10 GB features, and lastly (iv) using all 21 features. We further performed cross cell line prediction where a model trained on one cell line data was used to predict other cell line data. Gene Ontology analysis We used clusterProfiler 4.0 16 to identify biological processes associated with the identified up and downregulated genes. We used the following command to get the enriched significant processes: “ego <- enrichGO (de$Entrezid, OrgDb = “org. Hs.eg.db”, ont=“BP”, readable = TRUE, minGSSize = 10, maxGSSize = 500, keyType=“SYMBOL”)” As there are many redundant processes, we further obtained the parent processes using the following command: “ego2 <- simplify (ego, cutoff=0.8, by=“p.adjust”, select_fun=min, measure = “Wang”)” Dotplot of the above obtained processes were created using ggplot2 library in R. 17 Results Distinct patterns of epigenomic marks in gene locus between up- and down-regulated genes upon HDACi treatment For each of the two cell lines (HCT116 and RH4), for the respective dosage of HDACi drugs (eight doses of Largazole for HCT116, one dose of 1 μM Entinostat for RH4), we first identified the top 1000 up-regulated and 1000 down-regulated genes (Methods). Genes classified as up, down, and unchanged post-treatment for various doses in HCT116 and RH4 cell lines are provided in Tables S1 and S2 respectively of Extended data. 18 TPM value of each gene, untreated as well as treated for various concentrations in HCT116 cell line and single concentration for RH4 is also provided in Tables S3 and S4 of Extended data 18 respectively. We first identified enriched GO terms in each set of up and down-regulated genes (three pairs of gene sets for three representative doses in HCT116 [4.68 nM, 75 nM and 300 nM] and one pair for RH4). In general, the upregulated genes in both the cell lines were broadly enriched for the developmental and signaling processes ( Figure 1 ). The developmental process is in the direction of epithelial to mesenchymal transition (EMT). In the case of HCT116 cell line, we additionally observed response to hypoxia. Likewise, processes associated with downregulated genes are broadly associated with the cell cycle and cell division, whereas for RH4 cell line, additional processes such as histone modification and RNA splicing were also seen ( Figure 2 ). Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways, suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes, consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi. Some variability may also stem from the use of different HDAC inhibitors, which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling these two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design. Complete lists of processes associated with the upregulated and downregulated genes in HCT116 [4.68 nM, 75 nM and 300 nM] and RH4 cell lines are provided in Tables S5-S7 and S8 respectively of Extended data. 18 Figure 1. Upregulated genes are broadly enriched for the developmental and signaling processes. Top 20 enriched biological processes associated with upregulated genes in HCT116 cell line after treating with epigenetic drug largazole at 4.68 nM (A); 75 nM (B); 300 nM (C); and Top 20 enriched biological processes associated with upregulated genes in RH4 cell line after treating with epigenetic drug entinostat at 1 μM (D). Figure 2. Downregulated genes are broadly enriched for the cell cycle processes. Top 20 enriched biological processes associated with downregulated genes in HCT116 cell line after treating with epigenetic drug largazole at 4.68 nM (A); 75 nM (B); 300 nM (C); and Top 20 enriched biological processes associated with upregulated genes in RH4 cell line after treating with epigenetic drug entinostat at 1 μM (D). Next, we compared the pre-treatment epigenomic profiles of the up- and downregulated genes by plotting the H3K27Ac mark intensity (normalized read counts) in the pre-treatment sample along 21 genic bins (Methods). Distributions for three representative doses for HCT116 (4.68 nM (lowest), 75 nM, and 300 nM (highest)) and 1 μM dose for RH4 are included in Figure 3 ; all other distributions for the HCT116 cell line are provided in Figures S1-S5 of Extended data. 18 Overall, the following general trends emerged: (1) There was substantial variability across the bins around the genic locus in the upregulated versus downregulated H3K27Ac mark density, (2) in the upstream regions downregulated genes had a higher H3K27Ac pre-treatment; (3) this trend was also true in gene body but only at mid and higher dosage, while (4) at low dosage the trend was opposite in gene body where the downregulated genes had lower H3K27Ac; (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116. Overall, while there was a variable pre-treatment epigenomic pattern within the gene body across cell lines, drug, and dosages, there were nevertheless sufficient differences between up- and downregulated genes, motivating us to develop machine learning models to predict transcription effects given the H3K27Ac pattern at a gene locus. Figure 3. H3K27Ac mark distribution across genomic bins. Boxplot distribution of H3K27Ac marks across 21 genomic bins (10 equal sized bins of Promoter, Gene body and 1 bin of TSS) associated with upregulated (blue bars) and downregulated (brown bars) genes when HCT116 cell line is treated with largazole at concentrations 4.68 nM (A); 75 nM (B); 300 nM (C); and when RH4 cell line is treated with entinostat at 1 μM (D). Predicting HDACi treatment impact on gene expression from the epigenome Here, we assess whether the pre-treatment epigenetic profile at a gene locus can predict whether the gene will be upregulated or downregulated upon treatment with HDACi. The top 1000 upregulated and 1000 downregulated genes were compiled. For every gene, pre-treatment H3K27Ac read count in 21 regions relative to the gene (Methods) were used as features and three machine learning models – Support Vector Machine (SVM), Random Forest (RF), and Gradient Boosting (GB), were benchmarked based on five-fold cross-validation and accuracy was quantified as area under the ROC curve (AUC). A separate model was benchmarked for each of the eight drug dosages in HCT116 data. As shown in Table 1 , overall, various machine learning approaches performed comparably and using all features was preferable; specifically, the best performance was achieved by SVM for 75 nM dosage with AUC of 0.89 ( Figure 4 ). Analogous benchmarking for RH4 cell line data available at the single dosage using all features yielded comparable AUC ranging from 0.74-0.76 for the three machine learning methods. Overall, H3K27Ac signal near the gene is informative of the gene expression changes upon treatment with HDACi. Table 1. Performance of various machine learning models on HCT116 and RH4 cell line testing dataset. Here the concentration of drugs is in nM and μM. P stands for Promoter; TSS stands for Transcription Start Site; GB stands for Gene Body; and All is the combination of all three features. Support vector machine Random forest Gradient boosting Concentration P TSS GB All P TSS GB All P TSS GB All HCT116 (4.68 nM) 0.57 0.61 0.68 0.71 0.61 0.68 0.69 0.74 0.62 0.69 0.65 0.73 HCT116 (9.37 nM) 0.74 0.73 0.77 0.78 0.74 0.76 0.82 0.80 0.72 0.77 0.75 0.79 HCT116 (18.75 nM) 0.60 0.74 0.74 0.83 0.68 0.74 0.74 0.85 0.66 0.76 0.74 0.80 HCT116 (30 nM) 0.74 0.77 0.82 0.86 0.75 0.76 0.82 0.84 0.74 0.78 0.73 0.80 HCT116 (37.5 nM) 0.82 0.80 0.85 0.87 0.83 0.74 0.85 0.85 0.82 0.79 0.78 0.80 HCT116 (75 nM) 0.74 0.78 0.84 0.89 0.72 0.75 0.83 0.87 0.72 0.77 0.78 0.81 HCT116 (150 nM) 0.73 0.79 0.86 0.87 0.72 0.76 0.85 0.87 0.73 0.79 0.77 0.82 HCT116 (300 nM) 0.75 0.75 0.86 0.88 0.76 0.74 0.86 0.86 0.76 0.75 0.76 0.80 RH4 (1μM) 0.74 0.67 0.71 0.75 0.71 0.63 0.73 0.76 0.70 0.66 0.70 0.74 Figure 4. Performance of Support Vector Machine (SVM) model. Performance of various SVM based models in terms of Area Under Curve (AUC) at different concentrations when HCT116 cell line was treated with 8 different largazole concentration and RH4 cell line was treated with entinostat. Prediction model is generalizable across cell lines Next, we assessed whether a model trained on one cell line to predict the transcriptional effect of a certain epigenomic drug can predict the effect in a different cell line treated with a different drug, albeit also HDACi. Toward this, first, a SVM (75 nM) model trained on HCT116 cell line data was able to achieve an AUC value of 0.71 when applied to RH4 cell line data ( Figure 5A ). Likewise, the model trained on RH4 cell line data when applied to HCT116 data achieved an AUC of 0.81 ( Figure 5B ), supporting the cross-context generalizability of the model, consistent with similarity of epigenomic profile trends between the two cell lines as shown above ( Figure 3B and 3D ). Figure 5. Cross cell line prediction: (A) Performance of HCT116 data trained Support Vector Machine (SVM) model on RH4 cell line used as testing dataset. (B) Performance of RH4 data trained SVM model on HCT116 cell line used as testing dataset. Cell line-specific expression changes are reflected in their epigenome Next, we specifically assessed whether the context-specific differences across the two cell lines in their HDACi-induced gene expression were reflected in their context-specific pre-treatment epigenomic profile in the gene locus. Toward this we compared the data for HCT116 treated with 75 nM largazole with RH4 treated with 1 μM of entinostat. For each cell line we applied stringent criteria to identify genes which were upregulated in one cell line and downregulated in another cell line. We selected those genes whose fold change >3 in one cell line and <1/3 in another. This resulted in two gene sets: (1) 73 genes upregulated in HCT116 and downregulated in RH4, and (2) 184 genes upregulated in RH4 and downregulated in HCT116. To normalize for cell line-specific differences in H3K27Ac, we z-scored the cross-bin H3K27Ac signal for each gene. Then for these two gene sets, we plotted the normalized H3K27Ac intensities along the 21 genic bins, comparing two cell lines features. As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes ( Figure 6A & B ). Figure 6. H3K27Ac mark distribution across genomic bins in two cell lines. Boxplot distribution of H3K27Ac marks across 21 genomic bins (10 equal sized bins of Promoter, Gene body and 1 bin of TSS) associated with genes with positive log fold change in HCT116 but negative log fold change in RH4 cell line (A) and Genes with positive log fold change in RH4 but negative log fold change in HCT116 cell line (B). Discussion Epigenetic dysregulation is a key characteristic of cancers. A number of mutations have been observed in the genes encoding epigenetic modifiers such as DNA methylation and histone modification enzymes. 19 Accordingly, efforts have been made in targeting epigenetic regulators. 20 At present, seven epigenetics-targeting drugs have been approved by the FDA. 21 However, there are certain challenges associated with this class of drugs, limiting their success. Some of the key challenges include (i) Different epigenetic mutations are associated with different cancer types; (ii) The same gene may have opposite function in tumorigenesis of different cancers. For example, EZH2 deficiency causes myeloid malignancies 22 whereas gain-of-function causes B cell lymphomas 23 ; (iii) Another major issue is the selectivity of these drugs. For example, 30 enzymes of the KDM family with similar JMJC domain belong to five subfamilies. These enzymes demethylate different histone residues. Hence, drugs targeting these are broad-spectrum, affecting multiple KDM subfamilies and histone marks with potentially unintended consequences 24 ; (iv) Yet another issue with epigenetics-targeting drugs, focused on in this work, is the selectivity of genomic loci. For instance, a HDACi can both increase as well as decrease histone acetylation in different genomic loci and can thus upregulate certain genes while downregulating others, again with unintended consequences. Here, we tried to address the selectivity issue by developing a machine learning model based on pretreatment histone mark. In two cell lines, we established that the locus-specific effect of HDACi treatment on gene expression can be predicted to a reasonable accuracy from the pre-treatment histone acetylation pattern at a gene locus, and the model appears to be generalizable across cell lines. While the current study is promising and may potentially be applied to personalized therapy by predicting the transcriptomic consequence of HDACi treatment, there are a few limitations which need to be addressed. Our predictive model is based only on the H3K27ac mark. Several other marks such as H3K9ac, H3K4me3, H3K27me3, among others, should be incorporated in such modeling approaches in the future as and when such data become available. Our model was assessed only in cell lines and its efficacy in bulk tumor data representing the tumor microenvironment remains to be assessed. Last but not the least, pre- and post-treatment tumor epigenetic and transcriptomic data in clinical and pre-clinical models are still lacking, required for assessing the clinical applicability of our approach. Author contribution VG download and processed the data. PA, and SH perform the analysis. PA and SH perform the statistical analysis. PA, VG and SH wrote the manuscript. PA and SH supervised the study. All authors read the article and approved the submitted version. Declarations We declare that no third-party material was used in this study and also, we have not used AI tools at any point in the preparation of the manuscript. Data availability Underlying data Sequence Read Archive: Genome-wide Dose-dependent Inhibition of Histone Deacetylases Reveals Their Roles in Enhancer Remodeling and Suppression of Oncogenic Super-enhancers, https://identifiers.org/insdc.sra:SRP113250 . 25 Sequence Read Archive: Genome-wide Dose-dependent Inhibition of Histone Deacetylases Reveals Their Roles in Enhancer Remodeling and Suppression of Oncogenic Super-enhancers, https://identifiers.org/insdc.sra:SRP151465 . 26 Extended data Figshare: Supplementary_Figures, https://doi.org/10.6084/m9.figshare.23736882.v1 . 18 This project contains the following extended data: - Supplementary.xlsx - Supplementary_Figures.docx Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0). Analysis code Analysis code available from: https://github.com/hannenhalli-lab/Epigenetic_Project/ Archived analysis code at time of publication: https://zenodo.org/record/8212782 . 27 License: MIT Acknowledgements This work utilized the computational resources of the NIH HPC Biowulf cluster. Authors are also thankful to the other lab members for their valuable suggestions. This work was supported by the Intramural Research Program of the National Cancer Institute, Center for Cancer Research. References 1. Henikoff S, Greally JM: Epigenetics, cellular memory and gene regulation. Curr. Biol. 2016 Jul 25 [cited 2023 Jul 17]; 26 (14): R644–R648. Publisher Full Text Reference Source 2. Dawson MA: The cancer epigenome: Concepts, challenges, and therapeutic opportunities. Science. 2017 Mar 17 [cited 2023 Jul 17]; 355 (6330): 1147–1152. PubMed Abstract | Publisher Full Text 3. Jones PA, Issa JPJ, Baylin S: Targeting the cancer epigenome for therapy. Nat. Rev. Genet. 2016 Oct 1 [cited 2023 Jul 17]; 17 (10): 630–641. Publisher Full Text Reference Source 4. Amatori S, Bagaloni I, Donati B, et al. : DNA demethylating antineoplastic strategies: a comparative point of view. Genes Cancer. 2010 [cited 2023 Jul 17]; 1 (3): 197–209. PubMed Abstract | Publisher Full Text | Free Full Text 5. Feehley T, O’Donnell CW, Mendlein J, et al. : Drugging the epigenome in the age of precision medicine. Clin. Epigenetics. 2023 Dec 1 [cited 2023 Jul 17]; 15 (1): 6. PubMed Abstract | Publisher Full Text | Free Full Text 6. Reddington JP, Perricone SM, Nestor CE, et al. : Redistribution of H3K27me3 upon DNA hypomethylation results in de-repression of Polycomb target genes. Genome Biol. 2013 Mar 25 [cited 2023 Jul 17]; 14 (3): R25. PubMed Abstract | Publisher Full Text | Free Full Text 7. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning.Front. Genet [Internet].2022 Dec 14 [cited 2023 Dec 7]; 13 : 1081842. PubMed Abstract | Publisher Full Text 8. Singh R, Lanchantin J, Robins G, et al. : DeepChrome: deep-learning for predicting gene expression from histone modifications.Bioinformatics [Internet].2016 Sep 1 [cited 2023 Dec 7]; 32 (17): i639–i648. PubMed Abstract | Publisher Full Text 9. Karlić R, Chung HR, Lasserre J, et al. : Histone modification levels are predictive for gene expression.Proc. Natl. Acad. Sci. U S A. [Internet].2010 Feb 16 [cited 2023 Dec 7]; 107 (7): 2926–2931. PubMed Abstract | Publisher Full Text | Free Full Text 10. Md V, Misra S, Li H, et al. : Efficient architecture-aware acceleration of BWA-MEM for multicore systems. Proceedings - 2019 IEEE 33rd International Parallel and Distributed Processing Symposium, IPDPS 2019. 2019 May 1; pp. 314–324. 11. Patro R, Duggal G, Love MI, et al. : Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods. 2017 [cited 2023 Jul 17]; 14 (4): 417–419. PubMed Abstract | Publisher Full Text | Free Full Text 12. Liao Y, Smyth GK, Shi W: The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019 May 1 [cited 2023 Jul 17]; 47 (8): e47. PubMed Abstract | Publisher Full Text | Free Full Text 13. Orlando DA, Chen MW, Brown VE, et al. : Quantitative ChIP-Seq normalization reveals global modulation of the epigenome. Cell Rep. 2014 [cited 2023 Jul 17]; 9 (3): 1163–1170. PubMed Abstract | Publisher Full Text 14. Cunningham F, Allen JE, Allen J, et al. : Ensembl 2022. Nucleic Acids Res. 2022 Jan 7 [cited 2023 Jul 19]; 50 (D1): D988–D995. PubMed Abstract | Publisher Full Text | Free Full Text 15. Pedregosa FABIANPEDREGOSAF, Michel V, Grisel OLIVIERGRISELO, et al. : Scikit-learn: Machine Learning in Python Gaël Varoquaux Bertrand Thirion Vincent Dubourg Alexandre Passos PEDREGOSA, VAROQUAUX, GRAMFORT ET AL. Matthieu Perrot. J. Mach. Learn. Res. 2011 [cited 2023 Jul 17]; 12 : 2825–2830. Reference Source 16. Wu T, Hu E, Xu S, et al. : clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Cambridge (Mass)). 2021 Aug 28 [cited 2023 Jul 17]; 2 (3): 100141. Publisher Full Text Reference Source 17. Wickham H: ggpolt2 Elegant Graphics for Data Analysis. Use R! Series. 2016; 211. Publisher Full Text 18. Agrawal P: Supplementary_Figures.docx. [Data set]. figshare. 2023. Publisher Full Text 19. Morel D, Jeffery D, Aspeslagh S, et al. : Combining epigenetic drugs with other therapies for solid tumours - past lessons and future promise. Nat. Rev. Clin. Oncol. 2020 Feb 1 [cited 2023 Jul 17]; 17 (2): 91–107. PubMed Abstract | Publisher Full Text 20. Wang N, Ma T, Yu B: Targeting epigenetic regulators to overcome drug resistance in cancers. Signal Transduct. Target. Ther. 2023 Dec 1 [cited 2023 Jul 17]; 8 (1). Publisher Full Text Reference Source 21. Nepali K, Liou JP: Recent developments in epigenetic cancer therapeutics: clinical advancement and emerging trends. J. Biomed. Sci. 2021 Dec 1 [cited 2023 Jul 17]; 28 (1): 27. PubMed Abstract | Publisher Full Text | Free Full Text 22. Rinke J, Chase A, Cross NCP, et al. : EZH2 in Myeloid Malignancies. Cells. 2020 Jul 8 [cited 2023 Jul 17]; 9 (7). PubMed Abstract | Publisher Full Text | Free Full Text 23. Chu L, Tan D, Zhu M, et al. : EZH2 W113C is a gain-of-function mutation in B-cell lymphoma enabling both PRC2 methyltransferase activation and tazemetostat resistance. J. Biol. Chem. 2023 Apr 1 [cited 2023 Jul 17]; 299 (4): 103073. PubMed Abstract | Publisher Full Text | Free Full Text 24. Pfister SX, Ashworth A: Marked for death: targeting epigenetic changes in cancer. Nat. Rev. Drug Discov. 2017 Apr 1 [cited 2023 Jul 17]; 16 (4): 241–263. PubMed Abstract | Publisher Full Text 25. Chemistry and Biochemistry, University of Colorado Boulder: Genome-wide Dose-dependent Inhibition of Histone Deacetylases Reveals Their Roles in Enhancer Remodeling and Suppression of Oncogenic Super-enhancers, Sequence Read Archive. [Dataset]. 2017. Reference Source 26. Khan J; Genetics Branch, NCI, NIH: Genome-wide Dose-dependent Inhibition of Histone Deacetylases Reveals Their Roles in Enhancer Remodeling and Suppression of Oncogenic Super-enhancers, Sequence Read Archive. [Dataset]. 2019. Reference Source 27. Agrawal P, et al. : Github. Predicting gene expression changes upon epigenomic drug treatment.2023 [cited 2023 Aug 2]. Reference Source Comments on this article Comments (0) Version 3 VERSION 3 PUBLISHED 01 Sep 2023 ADD YOUR COMMENT Comment Author details Author details 1 Division of Medical Research, SRM Medical College Hospital & Research Centre, SRMIST, Kattankulathur, Chennai, Tamil Nadu, India 2 Cancer Data Science Laboratory, National Cancer Institute, Bethesda, Maryland, 20814, USA 3 Faculty research assistant, Center for Bioinformatics and Computational Biology, University of Maryland, College Park, MD, 20742, USA Piyush Agrawal Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Supervision, Validation, Writing – Original Draft Preparation, Writing – Review & Editing Vishaka Gopalan Roles: Conceptualization, Data Curation, Formal Analysis, Methodology, Writing – Original Draft Preparation Monjura Afrin Rumi Roles: Data Curation, Formal Analysis, Methodology Sridhar Hannenhalli Roles: Conceptualization, Formal Analysis, Methodology, Supervision, Writing – Original Draft Preparation, Writing – Review & Editing Competing interests No competing interests were disclosed. Grant information This work was supported by the National Cancer Institute (Intramural Research funding). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Article Versions (3) version 3 Revised Published: 02 May 2025, 12:1089 https://doi.org/10.12688/f1000research.140273.3 version 2 Revised Published: 19 Dec 2023, 12:1089 https://doi.org/10.12688/f1000research.140273.2 version 1 Published: 01 Sep 2023, 12:1089 https://doi.org/10.12688/f1000research.140273.1 Copyright © 2025 Agrawal P et al . This is an open access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The author(s) is/are employees of the US Government and therefore domestic copyright protection in USA does not apply to this work. The work may be protected under the copyright laws of other jurisdictions when used in those jurisdictions. Download Export To Sciwheel Bibtex EndNote ProCite Ref. Manager (RIS) Sente metrics Views Downloads F1000Research - - PubMed Central info_outline Data from PMC are received and updated monthly. - - Citations open_in_new 0 open_in_new 0 open_in_new SEE MORE DETAILS CITE how to cite this article Agrawal P, Gopalan V, Rumi MA and Hannenhalli S. Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.12688/f1000research.140273.3 ) NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS track receive updates on this article Track an article to receive email alerts on any updates to this article. TRACK THIS ARTICLE Share Open Peer Review Current Reviewer Status: ? Key to Reviewer Statuses VIEW HIDE Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Version 3 VERSION 3 PUBLISHED 02 May 2025 Revised Views 0 Cite How to cite this report: Kumar A. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383150 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383150 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 24 May 2025 Ashwini Kumar , Sharda University, Greater Noida, Uttar Pradesh, India Approved VIEWS 0 https://doi.org/10.5256/f1000research.181146.r383150 The article titled "Predicting gene expression changes upon epigenomic drug treatment" is an interesting approach to show the effect of drug treatment on the genomic and epigenomic changes in cancer cells and prepare a model that can predict the post-treatment ... Continue reading READ ALL The article titled "Predicting gene expression changes upon epigenomic drug treatment" is an interesting approach to show the effect of drug treatment on the genomic and epigenomic changes in cancer cells and prepare a model that can predict the post-treatment changes across the cancer types. Though drugs that target the epigenomic alterations in cancer have been clinically used, however they lack locus-specific action that may result in unpredictable adverse effects. Thus. the approach presented in this manuscript tackles pre- and post-treatment gene expression and pre-treatment epigenomic profile, and paves a way to predict the locus-specific effect(s) of the drugs that may lead us to personalize the therapy, and also lead us to develop more rational drugs. Although such models can be made on in vitro studies involving various cancer cells, however this manuscript presents a potential new approach to tackle the issue of epigenomic changes that have substantial effect pre- and post-treatment. These models can help the researchers develop better drug development strategies while clinicians to have better therapeutic strategies minimizing the adverse effects. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests: No competing interests were disclosed. Reviewer Expertise: Drug delivery and drug discovery I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Kumar A. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383150 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383150 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Dwivedi VD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383152 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383152 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 24 May 2025 Vivek Dhar Dwivedi , Quanta Calculus, Greater Noida, India Approved VIEWS 0 https://doi.org/10.5256/f1000research.181146.r383152 The current study proposes a novel and interesting strategy to predict gene expression changes using histone modification patterns, with a focus on H3K27 acetylation, in response to drug exposure. The findings may play a key role in enhancing personalized treatment ... Continue reading READ ALL The current study proposes a novel and interesting strategy to predict gene expression changes using histone modification patterns, with a focus on H3K27 acetylation, in response to drug exposure. The findings may play a key role in enhancing personalized treatment strategies for cancer patients undergoing histone deacetylase inhibitor therapy. Still, a few minor but important concerns need to be addressed prior to its acceptance. 1. Authors should explain why they use only the mentioned 2 cell lines in the current study. 2. Please provide an explanation if this kind of approach will also work for other histone marks? 3. Reason behind selecting the SVM method in the case of Figure 5. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests: No competing interests were disclosed. Reviewer Expertise: Bioinformatics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Dwivedi VD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383152 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383152 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Mandal SD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r382129 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-382129 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 24 May 2025 Sukhen Das Mandal , Ghani Khan Choudhary Institute of Engineering and Technology, Malda, West Bengal, India Approved VIEWS 0 https://doi.org/10.5256/f1000research.181146.r382129 All my concerns have been addressed ... Continue reading READ ALL All my concerns have been addressed satisfactorily. Therefore, this article may be approved. Competing Interests: No competing interests were disclosed. Reviewer Expertise: RNA biology, Machine learning I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Mandal SD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r382129 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-382129 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Version 2 VERSION 2 PUBLISHED 19 Dec 2023 Revised Views 0 Cite How to cite this report: Singh IK. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r267047 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-267047 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 07 May 2025 Indrakant K. Singh , Deshbandhu College, University of Delhi, Delhi, India Approved VIEWS 0 https://doi.org/10.5256/f1000research.159982.r267047 The study presents a novel machine learning approach to predict gene expression changes from histone modification data, specifically H3K27 acetylation, after drug treatment. This could have significant implications for personalized cancer treatment with histone deacetylase inhibitors. However, ... Continue reading READ ALL The study presents a novel machine learning approach to predict gene expression changes from histone modification data, specifically H3K27 acetylation, after drug treatment. This could have significant implications for personalized cancer treatment with histone deacetylase inhibitors. However, there are some concerns that need to be addressed before accepting. 1. The rationale behind choosing different parameters (e.g., z-scores of normalized read counts in Figure 6 versus normalized read counts in Figure 3) is not adequately explained. Consistency in data presentation or a clear justification for differences is needed. 2. The manuscript needs more detailed descriptions of certain methods. For example, the process of peak calling in ChIP-seq analysis and the significance of the log10 fold change in differential expression analysis are not clearly described. Detailed methodological transparency is essential for replication and understanding 3. Some results and their descriptions appear contradictory or unclear. For instance, the authors state that HCT116 genes showed higher concentrations of H3K27Ac marks, but this is not consistent with what is shown in Figure 6. Such inconsistencies need to be addressed and clarified. I recommend minor revision. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests: No competing interests were disclosed. Reviewer Expertise: Bioinformatics and Cancer Biology I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Singh IK. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r267047 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-267047 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Mandal SD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r243027 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-243027 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 01 Aug 2024 Sukhen Das Mandal , Ghani Khan Choudhary Institute of Engineering and Technology, Malda, West Bengal, India Approved with Reservations VIEWS 0 https://doi.org/10.5256/f1000research.159982.r243027 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. ... Continue reading READ ALL In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and down regulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. 6. The reference of the figure 5A is not there in the text. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? No Are the conclusions drawn adequately supported by the results? No Competing Interests: No competing interests were disclosed. Reviewer Expertise: RNA biology, Machine learning I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Mandal SD. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r243027 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-243027 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Author Response 07 May 2025 Piyush Agrawal , Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA 07 May 2025 Author Response Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change ... Continue reading Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and downregulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. Response: We thank the reviewer for this insightful comment. In response, we have expanded the discussion in the Results section to clarify the biological relevance of the distinct GO term enrichments observed across different cell lines following HDACi treatment. Specifically, we observed that upregulated genes were consistently enriched for developmental and cell signaling processes , while downregulated genes were predominantly associated with cell cycle-related functions . This pattern suggests that HDAC inhibition may generally suppress proliferation while activating differentiation or lineage-specifying pathways. However, the specific GO terms varied somewhat across different cell lines, which may either reflect context-specific effects of HDACi, or the differences in the specific HDAC inhibitors used ; however, since different HDACi compounds were used in different cell lines, it is not possible to distinguish between these two possibilities. We have now added this interpretation to the Results section (page 5, lines 139–149). The following section has been added in the revised manuscript: “Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways , suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes , consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi . Some variability may also stem from the use of different HDAC inhibitors , which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling thse two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design.” 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. Response: We completely agree and had considered this. Unfortunately, however the post-treatment H3K27Ac data is currently not available. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. Response: We regret the lack of clarity. For point 3, we are talking about higher distribution of H3K27Ac epigenomic marks in gene body regions of the downregulated genes for HCT cell line with concentration of 75nM and 300nM. In the figure, in case of 75nM, gene body bin 1-3, 6 and 8-10 shows higher H3K27Ac read density concentration, where in remaining bins its nearly equal among up and downregulated genes. Likewise, in case of 300nM, majority of the gene body shows higher H3K27Ac distribution then promoter region except bin 1, and 10. In revised manuscript, we have now revised the statement for clarity. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. Response: We appreciate this. However, we would like to clarify that our analysis already incorporates normalization that accounts for these differences . Specifically, as mentioned in the Methods section, we normalized the read counts within each genomic bin to TPM (Transcripts Per Million) using the Rsubread package. TPM normalization adjusts for both sequencing depth and bin length , thereby mitigating the potential bias introduced by longer bins containing more nucleotides. Thus, the enrichment patterns we observed are not confounded by gene length variation and reflect true biological signal distribution rather than technical artifacts. We have now made this point more explicit in the manuscript for clarity. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again, these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. Response: We thank the reviewer for pointing out this mistake and we also regret the misinterpretation of the results. In the revised manuscript, we have clarified this for the better explanation. Page 7; Lines 217:220. “As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes (Figure 6A&B)” . 6. The reference of the figure 5A is not there in the text. Response: Figure 5A is referenced on Page No. 7, Line No 201. Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and downregulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. Response: We thank the reviewer for this insightful comment. In response, we have expanded the discussion in the Results section to clarify the biological relevance of the distinct GO term enrichments observed across different cell lines following HDACi treatment. Specifically, we observed that upregulated genes were consistently enriched for developmental and cell signaling processes , while downregulated genes were predominantly associated with cell cycle-related functions . This pattern suggests that HDAC inhibition may generally suppress proliferation while activating differentiation or lineage-specifying pathways. However, the specific GO terms varied somewhat across different cell lines, which may either reflect context-specific effects of HDACi, or the differences in the specific HDAC inhibitors used ; however, since different HDACi compounds were used in different cell lines, it is not possible to distinguish between these two possibilities. We have now added this interpretation to the Results section (page 5, lines 139–149). The following section has been added in the revised manuscript: “Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways , suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes , consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi . Some variability may also stem from the use of different HDAC inhibitors , which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling thse two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design.” 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. Response: We completely agree and had considered this. Unfortunately, however the post-treatment H3K27Ac data is currently not available. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. Response: We regret the lack of clarity. For point 3, we are talking about higher distribution of H3K27Ac epigenomic marks in gene body regions of the downregulated genes for HCT cell line with concentration of 75nM and 300nM. In the figure, in case of 75nM, gene body bin 1-3, 6 and 8-10 shows higher H3K27Ac read density concentration, where in remaining bins its nearly equal among up and downregulated genes. Likewise, in case of 300nM, majority of the gene body shows higher H3K27Ac distribution then promoter region except bin 1, and 10. In revised manuscript, we have now revised the statement for clarity. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. Response: We appreciate this. However, we would like to clarify that our analysis already incorporates normalization that accounts for these differences . Specifically, as mentioned in the Methods section, we normalized the read counts within each genomic bin to TPM (Transcripts Per Million) using the Rsubread package. TPM normalization adjusts for both sequencing depth and bin length , thereby mitigating the potential bias introduced by longer bins containing more nucleotides. Thus, the enrichment patterns we observed are not confounded by gene length variation and reflect true biological signal distribution rather than technical artifacts. We have now made this point more explicit in the manuscript for clarity. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again, these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. Response: We thank the reviewer for pointing out this mistake and we also regret the misinterpretation of the results. In the revised manuscript, we have clarified this for the better explanation. Page 7; Lines 217:220. “As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes (Figure 6A&B)” . 6. The reference of the figure 5A is not there in the text. Response: Figure 5A is referenced on Page No. 7, Line No 201. Competing Interests: No competing interests were disclosed. Close Report a concern Respond or Comment COMMENTS ON THIS REPORT Author Response 07 May 2025 Piyush Agrawal , Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA 07 May 2025 Author Response Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change ... Continue reading Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and downregulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. Response: We thank the reviewer for this insightful comment. In response, we have expanded the discussion in the Results section to clarify the biological relevance of the distinct GO term enrichments observed across different cell lines following HDACi treatment. Specifically, we observed that upregulated genes were consistently enriched for developmental and cell signaling processes , while downregulated genes were predominantly associated with cell cycle-related functions . This pattern suggests that HDAC inhibition may generally suppress proliferation while activating differentiation or lineage-specifying pathways. However, the specific GO terms varied somewhat across different cell lines, which may either reflect context-specific effects of HDACi, or the differences in the specific HDAC inhibitors used ; however, since different HDACi compounds were used in different cell lines, it is not possible to distinguish between these two possibilities. We have now added this interpretation to the Results section (page 5, lines 139–149). The following section has been added in the revised manuscript: “Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways , suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes , consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi . Some variability may also stem from the use of different HDAC inhibitors , which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling thse two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design.” 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. Response: We completely agree and had considered this. Unfortunately, however the post-treatment H3K27Ac data is currently not available. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. Response: We regret the lack of clarity. For point 3, we are talking about higher distribution of H3K27Ac epigenomic marks in gene body regions of the downregulated genes for HCT cell line with concentration of 75nM and 300nM. In the figure, in case of 75nM, gene body bin 1-3, 6 and 8-10 shows higher H3K27Ac read density concentration, where in remaining bins its nearly equal among up and downregulated genes. Likewise, in case of 300nM, majority of the gene body shows higher H3K27Ac distribution then promoter region except bin 1, and 10. In revised manuscript, we have now revised the statement for clarity. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. Response: We appreciate this. However, we would like to clarify that our analysis already incorporates normalization that accounts for these differences . Specifically, as mentioned in the Methods section, we normalized the read counts within each genomic bin to TPM (Transcripts Per Million) using the Rsubread package. TPM normalization adjusts for both sequencing depth and bin length , thereby mitigating the potential bias introduced by longer bins containing more nucleotides. Thus, the enrichment patterns we observed are not confounded by gene length variation and reflect true biological signal distribution rather than technical artifacts. We have now made this point more explicit in the manuscript for clarity. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again, these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. Response: We thank the reviewer for pointing out this mistake and we also regret the misinterpretation of the results. In the revised manuscript, we have clarified this for the better explanation. Page 7; Lines 217:220. “As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes (Figure 6A&B)” . 6. The reference of the figure 5A is not there in the text. Response: Figure 5A is referenced on Page No. 7, Line No 201. Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and downregulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. Response: We thank the reviewer for this insightful comment. In response, we have expanded the discussion in the Results section to clarify the biological relevance of the distinct GO term enrichments observed across different cell lines following HDACi treatment. Specifically, we observed that upregulated genes were consistently enriched for developmental and cell signaling processes , while downregulated genes were predominantly associated with cell cycle-related functions . This pattern suggests that HDAC inhibition may generally suppress proliferation while activating differentiation or lineage-specifying pathways. However, the specific GO terms varied somewhat across different cell lines, which may either reflect context-specific effects of HDACi, or the differences in the specific HDAC inhibitors used ; however, since different HDACi compounds were used in different cell lines, it is not possible to distinguish between these two possibilities. We have now added this interpretation to the Results section (page 5, lines 139–149). The following section has been added in the revised manuscript: “Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways , suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes , consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi . Some variability may also stem from the use of different HDAC inhibitors , which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling thse two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design.” 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. Response: We completely agree and had considered this. Unfortunately, however the post-treatment H3K27Ac data is currently not available. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. Response: We regret the lack of clarity. For point 3, we are talking about higher distribution of H3K27Ac epigenomic marks in gene body regions of the downregulated genes for HCT cell line with concentration of 75nM and 300nM. In the figure, in case of 75nM, gene body bin 1-3, 6 and 8-10 shows higher H3K27Ac read density concentration, where in remaining bins its nearly equal among up and downregulated genes. Likewise, in case of 300nM, majority of the gene body shows higher H3K27Ac distribution then promoter region except bin 1, and 10. In revised manuscript, we have now revised the statement for clarity. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. Response: We appreciate this. However, we would like to clarify that our analysis already incorporates normalization that accounts for these differences . Specifically, as mentioned in the Methods section, we normalized the read counts within each genomic bin to TPM (Transcripts Per Million) using the Rsubread package. TPM normalization adjusts for both sequencing depth and bin length , thereby mitigating the potential bias introduced by longer bins containing more nucleotides. Thus, the enrichment patterns we observed are not confounded by gene length variation and reflect true biological signal distribution rather than technical artifacts. We have now made this point more explicit in the manuscript for clarity. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again, these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. Response: We thank the reviewer for pointing out this mistake and we also regret the misinterpretation of the results. In the revised manuscript, we have clarified this for the better explanation. Page 7; Lines 217:220. “As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes (Figure 6A&B)” . 6. The reference of the figure 5A is not there in the text. Response: Figure 5A is referenced on Page No. 7, Line No 201. Competing Interests: No competing interests were disclosed. Close Report a concern COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Liao D and Huo Z. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r239600 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-239600 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 22 Feb 2024 Daiqing Liao , Department of Anatomy and Cell Biology, College of Medicine & UF Health Cancer Center, University of Florida, Gainesville, USA Zhiguang Huo , Department of Biostatistics, University of Florida, Gainesville, Florida, USA Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.159982.r239600 Cancer epigenome underpins oncogenic gene expression. Histone acetylation is a major epigenomic event that impacts the transcriptome. The dynamic action of histone acetyltransferases (HATs) and deacetylases (HDACs) shape the acetylome. Small molecule inhibitors of both HATs and HDACs are promising ... Continue reading READ ALL Cancer epigenome underpins oncogenic gene expression. Histone acetylation is a major epigenomic event that impacts the transcriptome. The dynamic action of histone acetyltransferases (HATs) and deacetylases (HDACs) shape the acetylome. Small molecule inhibitors of both HATs and HDACs are promising agents for cancer therapy. Several HDAC inhibitors (HDACi) have been approved for treating cutaneous T-cell lymphoma and peripheral T-cell lymphoma. HDACi profoundly perturbs the transcriptome 1 . Thousands of genes can be downregulated and upregulated after treatment with an HDACi 2,3 . Pretreatment epigenome could influence transcriptomic perturbation by HDACi, although this has not been thoroughly examined. The manuscript titled “Predicting gene expression changes upon epigenomic drug treatment” by Agrawal, Gopalan, and Hannenhalli assessed whether the pretreatment epigenomic profile of cells could predict gene expression alterations induced by HDACi. The authors used the levels of a single histone mark, namely H3K27ac, in divided bins (21 genomic bins in total) from 2 kb upstream of the transcription start site (TSS) to the end of gene body and data of differentially expressed genes (DEGs) as inputs to study the relationship between the pretreatment H3K27ac levels and the transcription outputs in several machine learning models. The authors found that the pretreatment H3K27ac patterns could predict the post-treatment up versus downregulated genes with reasonable accuracy. Essentially, gene expression levels due to HDACi treatment seem to be inversely correlated with the levels of H3K27ac before treatment, which is observed in two cell lines treated with two HDACi of different chemical classes (entinostat and largazole). While the results are interesting, several significant issues should be addressed. The authors should also generate the feature importance score, e.g., the importance score for the random forest model (all gene models). This would be informative for readers to understand which genomic bin is most important to predict gene expression regulation. The authors should show the p-value for each comparison pair in Figure 3. The authors’ description “with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks, and for the second gene set, the opposite was true (Figure 6B), consistent with the patterns in Figure 3” (page 9) seems contradictory to what is shown in Figure 6. This must be clarified. In Figure 6, the z-scores of normalized read counts are used, while normalized read counts are shown in Figure 3. What is the rationale for using different parameters? Also, the p-value for each pair comparison should be shown in Figure 6. Although the authors acknowledged that different epigenomic marks should be included in future studies, they are encouraged to analyze at least whether a different histone mark associated with active transcription, such as H3K4me3, would exhibit a similar predictive power as H3K27ac. Another major limitation of the manuscript is that data based on only two cell lines and two HDACi are analyzed. An expanded study of data from a diverse set of samples (cells and tumors) treated with additional HDACi should improve the validity and rigor of the study. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Gryder BE, Wu L, Woldemichael GM, Pomella S, et al.: Chemical genomics reveals histone deacetylases are required for core regulatory transcription. Nat Commun . 2019; 10 (1): 3004 PubMed Abstract | Publisher Full Text 2. Lauffer BE, Mintzer R, Fong R, Mukund S, et al.: Histone deacetylase (HDAC) inhibitor kinetic rate constants correlate with cellular histone acetylation but not transcription and cell viability. J Biol Chem . 2013; 288 (37): 26926-43 PubMed Abstract | Publisher Full Text 3. Xiao Y, Hale S, Awasthee N, Meng C, et al.: HDAC3 and HDAC8 PROTAC dual degrader reveals roles of histone acetylation in gene regulation. Cell Chem Biol . 2023; 30 (11): 1421-1435.e12 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests were disclosed. Reviewer Expertise: bioinformatics (Zhiguang Huo) and epigenetics (Daiqing Liao) We confirm that we have read this submission and believe that we have an appropriate level of expertise to state that we do not consider it to be of an acceptable scientific standard, for reasons outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Liao D and Huo Z. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r239600 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-239600 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Respond or Comment COMMENT ON THIS REPORT Version 1 VERSION 1 PUBLISHED 01 Sep 2023 Views 0 Cite How to cite this report: Merkel A. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.153610.r203232 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v1#referee-response-203232 NOTE: it is important to ensure the information in square brackets after the title is included in this citation. Close Copy Citation Details Reviewer Report 20 Sep 2023 Angelika Merkel , Bioinformatics Unit, Josep Carreras Leukemia Research Institute (IJC), Barcelona, Spain Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.153610.r203232 Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone ... Continue reading READ ALL Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Was peak calling performed to eliminate background noise in the Chip-seq analysis? Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? The GO analysis is repeated in the methods and results section Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests were disclosed. Reviewer Expertise: transcriptomics, epigenetics, bioinformatics I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Merkel A. Reviewer Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.153610.r203232 ) The direct URL for this report is: https://f1000research.com/articles/12-1089/v1#referee-response-203232 NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article. COPY CITATION DETAILS Report a concern Author Response 17 Jan 2024 Piyush Agrawal , Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA 17 Jan 2024 Author Response Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful ... Continue reading Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. Response: We thank the reviewers for finding our work and interesting and potentially useful. Although we agree with reviewer regarding the importance of addressing variations in tissues, individual or disease type, unfortunately, lack of appropriate data currently prohibits such analysis. Besides the 2 cell lines used in this study, we didn’t come across any dataset (clinical/cell line) where for a given patient, both epigenomic marks and gene expression data is provided for both pre and post drug treatment. We would be happy to do the additional analysis if the reviewer can kindly point us to such data. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Response: We thank the reviewer for pointing out these studies, however, we want to clarify that the above-mentioned studies address a different question than ours, which is “to predict gene expression from the epigenome in a specific context”. Whereas we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. In the revised manuscript, we have highlighted this difference in the Introduction. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Response: Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Response: We regret lack of clarity. HCT116 and RH4 cell lines were chosen because the data we required for our study was available only for these 2 cell lines, namely, pre-treatment epigenomic profile and pre- and post-treatment gene expression data. HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. In the revised manuscript, we have added the description of these cell lines in the Methods section. Was peak calling performed to eliminate background noise in the Chip-seq analysis? Response: No, peak calling wasn’t performed in the current study. The goal of our analysis was to compare changes in H3K27ac read density upon HDAC treatment. Peak calling represents a 0/1 binarization of a continuous variable (in this case, normalized read counts in promoters and gene bodies). Our model uses the epigenomic signal intensity in specific windows around the gene magnitude, making peak-calling unnecessary. Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? Response: The mean Log-TPM values of each gene was computed for both treated and untreated conditions, with the Log-FC computed as the difference between treated and untreated conditions. The top 1000 most up-regulated and 1000 most down-regulated genes were considered to be differentially expressed. We mistakenly mentioned Log10, which is now corrected to Log2. We have corrected this in the revised version . The GO analysis is repeated in the methods and results section Response: We want to clarify that in Methods Section, we mentioned the software and command only used for GO analysis, however, in Result section, we discussed the enriched terms associated with the up and downregulated genes. Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Response: We agree with the reviewer that scale among plots differ and gene bins are of variable size. The only reason for plotting them together was to analyze all the 21 features together in single image i.e. change in distribution of the pre-treatment epigenomic marks in up and down genes post treatment. Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. Response: We thank the reviewers for finding our work and interesting and potentially useful. Although we agree with reviewer regarding the importance of addressing variations in tissues, individual or disease type, unfortunately, lack of appropriate data currently prohibits such analysis. Besides the 2 cell lines used in this study, we didn’t come across any dataset (clinical/cell line) where for a given patient, both epigenomic marks and gene expression data is provided for both pre and post drug treatment. We would be happy to do the additional analysis if the reviewer can kindly point us to such data. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Response: We thank the reviewer for pointing out these studies, however, we want to clarify that the above-mentioned studies address a different question than ours, which is “to predict gene expression from the epigenome in a specific context”. Whereas we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. In the revised manuscript, we have highlighted this difference in the Introduction. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Response: Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Response: We regret lack of clarity. HCT116 and RH4 cell lines were chosen because the data we required for our study was available only for these 2 cell lines, namely, pre-treatment epigenomic profile and pre- and post-treatment gene expression data. HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. In the revised manuscript, we have added the description of these cell lines in the Methods section. Was peak calling performed to eliminate background noise in the Chip-seq analysis? Response: No, peak calling wasn’t performed in the current study. The goal of our analysis was to compare changes in H3K27ac read density upon HDAC treatment. Peak calling represents a 0/1 binarization of a continuous variable (in this case, normalized read counts in promoters and gene bodies). Our model uses the epigenomic signal intensity in specific windows around the gene magnitude, making peak-calling unnecessary. Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? Response: The mean Log-TPM values of each gene was computed for both treated and untreated conditions, with the Log-FC computed as the difference between treated and untreated conditions. The top 1000 most up-regulated and 1000 most down-regulated genes were considered to be differentially expressed. We mistakenly mentioned Log10, which is now corrected to Log2. We have corrected this in the revised version . The GO analysis is repeated in the methods and results section Response: We want to clarify that in Methods Section, we mentioned the software and command only used for GO analysis, however, in Result section, we discussed the enriched terms associated with the up and downregulated genes. Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Response: We agree with the reviewer that scale among plots differ and gene bins are of variable size. The only reason for plotting them together was to analyze all the 21 features together in single image i.e. change in distribution of the pre-treatment epigenomic marks in up and down genes post treatment. Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests Close Report a concern Respond or Comment COMMENTS ON THIS REPORT Author Response 17 Jan 2024 Piyush Agrawal , Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA 17 Jan 2024 Author Response Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful ... Continue reading Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. Response: We thank the reviewers for finding our work and interesting and potentially useful. Although we agree with reviewer regarding the importance of addressing variations in tissues, individual or disease type, unfortunately, lack of appropriate data currently prohibits such analysis. Besides the 2 cell lines used in this study, we didn’t come across any dataset (clinical/cell line) where for a given patient, both epigenomic marks and gene expression data is provided for both pre and post drug treatment. We would be happy to do the additional analysis if the reviewer can kindly point us to such data. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Response: We thank the reviewer for pointing out these studies, however, we want to clarify that the above-mentioned studies address a different question than ours, which is “to predict gene expression from the epigenome in a specific context”. Whereas we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. In the revised manuscript, we have highlighted this difference in the Introduction. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Response: Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Response: We regret lack of clarity. HCT116 and RH4 cell lines were chosen because the data we required for our study was available only for these 2 cell lines, namely, pre-treatment epigenomic profile and pre- and post-treatment gene expression data. HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. In the revised manuscript, we have added the description of these cell lines in the Methods section. Was peak calling performed to eliminate background noise in the Chip-seq analysis? Response: No, peak calling wasn’t performed in the current study. The goal of our analysis was to compare changes in H3K27ac read density upon HDAC treatment. Peak calling represents a 0/1 binarization of a continuous variable (in this case, normalized read counts in promoters and gene bodies). Our model uses the epigenomic signal intensity in specific windows around the gene magnitude, making peak-calling unnecessary. Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? Response: The mean Log-TPM values of each gene was computed for both treated and untreated conditions, with the Log-FC computed as the difference between treated and untreated conditions. The top 1000 most up-regulated and 1000 most down-regulated genes were considered to be differentially expressed. We mistakenly mentioned Log10, which is now corrected to Log2. We have corrected this in the revised version . The GO analysis is repeated in the methods and results section Response: We want to clarify that in Methods Section, we mentioned the software and command only used for GO analysis, however, in Result section, we discussed the enriched terms associated with the up and downregulated genes. Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Response: We agree with the reviewer that scale among plots differ and gene bins are of variable size. The only reason for plotting them together was to analyze all the 21 features together in single image i.e. change in distribution of the pre-treatment epigenomic marks in up and down genes post treatment. Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. Response: We thank the reviewers for finding our work and interesting and potentially useful. Although we agree with reviewer regarding the importance of addressing variations in tissues, individual or disease type, unfortunately, lack of appropriate data currently prohibits such analysis. Besides the 2 cell lines used in this study, we didn’t come across any dataset (clinical/cell line) where for a given patient, both epigenomic marks and gene expression data is provided for both pre and post drug treatment. We would be happy to do the additional analysis if the reviewer can kindly point us to such data. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Response: We thank the reviewer for pointing out these studies, however, we want to clarify that the above-mentioned studies address a different question than ours, which is “to predict gene expression from the epigenome in a specific context”. Whereas we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. In the revised manuscript, we have highlighted this difference in the Introduction. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Response: Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Response: We regret lack of clarity. HCT116 and RH4 cell lines were chosen because the data we required for our study was available only for these 2 cell lines, namely, pre-treatment epigenomic profile and pre- and post-treatment gene expression data. HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. In the revised manuscript, we have added the description of these cell lines in the Methods section. Was peak calling performed to eliminate background noise in the Chip-seq analysis? Response: No, peak calling wasn’t performed in the current study. The goal of our analysis was to compare changes in H3K27ac read density upon HDAC treatment. Peak calling represents a 0/1 binarization of a continuous variable (in this case, normalized read counts in promoters and gene bodies). Our model uses the epigenomic signal intensity in specific windows around the gene magnitude, making peak-calling unnecessary. Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? Response: The mean Log-TPM values of each gene was computed for both treated and untreated conditions, with the Log-FC computed as the difference between treated and untreated conditions. The top 1000 most up-regulated and 1000 most down-regulated genes were considered to be differentially expressed. We mistakenly mentioned Log10, which is now corrected to Log2. We have corrected this in the revised version . The GO analysis is repeated in the methods and results section Response: We want to clarify that in Methods Section, we mentioned the software and command only used for GO analysis, however, in Result section, we discussed the enriched terms associated with the up and downregulated genes. Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Response: We agree with the reviewer that scale among plots differ and gene bins are of variable size. The only reason for plotting them together was to analyze all the 21 features together in single image i.e. change in distribution of the pre-treatment epigenomic marks in up and down genes post treatment. Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests Close Report a concern COMMENT ON THIS REPORT Comments on this article Comments (0) Version 3 VERSION 3 PUBLISHED 01 Sep 2023 ADD YOUR COMMENT Comment keyboard_arrow_left keyboard_arrow_right Open Peer Review Reviewer Status info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Reviewer Reports Invited Reviewers 1 2 3 4 5 6 Version 3 (revision) 02 May 25 read read read Version 2 (revision) 19 Dec 23 read read read Version 1 01 Sep 23 read Angelika Merkel , Josep Carreras Leukemia Research Institute (IJC), Barcelona, Spain Daiqing Liao , University of Florida, Gainesville, USA Zhiguang Huo , University of Florida, Gainesville, USA Sukhen Das Mandal , Ghani Khan Choudhary Institute of Engineering and Technology, Malda, West Bengal, India Indrakant K. Singh , Deshbandhu College, University of Delhi, Delhi, India Vivek Dhar Dwivedi , Quanta Calculus, Greater Noida, India Ashwini Kumar , Sharda University, Greater Noida, India Comments on this article All Comments (0) Add a comment Sign up for content alerts Sign Up You are now signed up to receive this alert Browse by related subjects keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Kumar A. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 24 May 2025 | for Version 3 Ashwini Kumar , Sharda University, Greater Noida, Uttar Pradesh, India 0 Views copyright © 2025 Kumar A. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The article titled "Predicting gene expression changes upon epigenomic drug treatment" is an interesting approach to show the effect of drug treatment on the genomic and epigenomic changes in cancer cells and prepare a model that can predict the post-treatment changes across the cancer types. Though drugs that target the epigenomic alterations in cancer have been clinically used, however they lack locus-specific action that may result in unpredictable adverse effects. Thus. the approach presented in this manuscript tackles pre- and post-treatment gene expression and pre-treatment epigenomic profile, and paves a way to predict the locus-specific effect(s) of the drugs that may lead us to personalize the therapy, and also lead us to develop more rational drugs. Although such models can be made on in vitro studies involving various cancer cells, however this manuscript presents a potential new approach to tackle the issue of epigenomic changes that have substantial effect pre- and post-treatment. These models can help the researchers develop better drug development strategies while clinicians to have better therapeutic strategies minimizing the adverse effects. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests No competing interests were disclosed. Reviewer Expertise Drug delivery and drug discovery I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. reply Respond to this report Responses (0) Kumar A. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383150) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383150 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Dwivedi V. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 24 May 2025 | for Version 3 Vivek Dhar Dwivedi , Quanta Calculus, Greater Noida, India 0 Views copyright © 2025 Dwivedi V. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The current study proposes a novel and interesting strategy to predict gene expression changes using histone modification patterns, with a focus on H3K27 acetylation, in response to drug exposure. The findings may play a key role in enhancing personalized treatment strategies for cancer patients undergoing histone deacetylase inhibitor therapy. Still, a few minor but important concerns need to be addressed prior to its acceptance. 1. Authors should explain why they use only the mentioned 2 cell lines in the current study. 2. Please provide an explanation if this kind of approach will also work for other histone marks? 3. Reason behind selecting the SVM method in the case of Figure 5. Is the work clearly and accurately presented and does it cite the current literature? Yes Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests No competing interests were disclosed. Reviewer Expertise Bioinformatics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. reply Respond to this report Responses (0) Dwivedi VD. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r383152) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-383152 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Mandal S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 24 May 2025 | for Version 3 Sukhen Das Mandal , Ghani Khan Choudhary Institute of Engineering and Technology, Malda, West Bengal, India 0 Views copyright © 2025 Mandal S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions All my concerns have been addressed satisfactorily. Therefore, this article may be approved. Competing Interests No competing interests were disclosed. Reviewer Expertise RNA biology, Machine learning I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. reply Respond to this report Responses (0) Mandal SD. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.181146.r382129) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v3#referee-response-382129 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2025 Singh I. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 07 May 2025 | for Version 2 Indrakant K. Singh , Deshbandhu College, University of Delhi, Delhi, India 0 Views copyright © 2025 Singh I. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions The study presents a novel machine learning approach to predict gene expression changes from histone modification data, specifically H3K27 acetylation, after drug treatment. This could have significant implications for personalized cancer treatment with histone deacetylase inhibitors. However, there are some concerns that need to be addressed before accepting. 1. The rationale behind choosing different parameters (e.g., z-scores of normalized read counts in Figure 6 versus normalized read counts in Figure 3) is not adequately explained. Consistency in data presentation or a clear justification for differences is needed. 2. The manuscript needs more detailed descriptions of certain methods. For example, the process of peak calling in ChIP-seq analysis and the significance of the log10 fold change in differential expression analysis are not clearly described. Detailed methodological transparency is essential for replication and understanding 3. Some results and their descriptions appear contradictory or unclear. For instance, the authors state that HCT116 genes showed higher concentrations of H3K27Ac marks, but this is not consistent with what is shown in Figure 6. Such inconsistencies need to be addressed and clarified. I recommend minor revision. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Yes Competing Interests No competing interests were disclosed. Reviewer Expertise Bioinformatics and Cancer Biology I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard. reply Respond to this report Responses (0) Singh IK. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r267047) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-267047 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2024 Mandal S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 01 Aug 2024 | for Version 2 Sukhen Das Mandal , Ghani Khan Choudhary Institute of Engineering and Technology, Malda, West Bengal, India 0 Views copyright © 2024 Mandal S. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (1) Approved With Reservations info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and down regulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. 6. The reference of the figure 5A is not there in the text. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Not applicable Are all the source data underlying the results available to ensure full reproducibility? No Are the conclusions drawn adequately supported by the results? No Competing Interests No competing interests were disclosed. Reviewer Expertise RNA biology, Machine learning I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above. reply Respond to this report Responses (1) Author Response 07 May 2025 Piyush Agrawal, Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA Reviewer 3 In the paper authors developed machine learning models to predict gene expression changes for HDACi treatment. Though it will enrich current understanding about how epigenomic drugs can change expression of genes, there are multiple/some issues those need to be addressed. Suggestions: 1. GO analysis of upregulated and downregulated genes for HDACi treatment showed that genes are associated with different types of GO terms for different cell lines. Whether it is because different cell lines respond differently to HDACi or different types of HDACi act through different mechanisms, you should discuss it more elaborately in the result section. Response: We thank the reviewer for this insightful comment. In response, we have expanded the discussion in the Results section to clarify the biological relevance of the distinct GO term enrichments observed across different cell lines following HDACi treatment. Specifically, we observed that upregulated genes were consistently enriched for developmental and cell signaling processes , while downregulated genes were predominantly associated with cell cycle-related functions . This pattern suggests that HDAC inhibition may generally suppress proliferation while activating differentiation or lineage-specifying pathways. However, the specific GO terms varied somewhat across different cell lines, which may either reflect context-specific effects of HDACi, or the differences in the specific HDAC inhibitors used ; however, since different HDACi compounds were used in different cell lines, it is not possible to distinguish between these two possibilities. We have now added this interpretation to the Results section (page 5, lines 139–149). The following section has been added in the revised manuscript: “Gene Ontology (GO) enrichment analysis revealed that genes upregulated by HDACi treatment were broadly involved in developmental processes and cell signaling pathways , suggesting a potential shift toward differentiation or activation of lineage-specific programs. In contrast, downregulated genes were strongly enriched for cell cycle-related processes , consistent with the known anti-proliferative effects of HDAC inhibitors. Notably, while this overall pattern was observed across all cell lines, the specific GO terms and pathways varied between them. This indicates that cell line-specific factors—such as differences in chromatin accessibility, baseline gene expression, or mutational background—may influence the transcriptional response to HDACi . Some variability may also stem from the use of different HDAC inhibitors , which may target different HDAC isoforms and thus modulate distinct regulatory networks. Disentangling thse two possibilities however will necessitate a broader set of experiments involving multiple cell lines and multiple HDACi with proper experimental design.” 2. H3K27Ac profile for treated vs untreated should be checked for up and down regulated genes. This can verify whether change of H3K27Ac causes transcriptional change which leads to expression change or post transcriptional regulation like RNA stability causes the change of the expression. Response: We completely agree and had considered this. Unfortunately, however the post-treatment H3K27Ac data is currently not available. 3. "(3) this trend was also true in gene body but only at mid and higher dosage" or " (5) RH4 trends at 1 μm dose of Entinostat most resembled the patterns at 75 nM dose of largazole in HCT116" - These statements do not match with the figures. Kindly explain these more clearly. Response: We regret the lack of clarity. For point 3, we are talking about higher distribution of H3K27Ac epigenomic marks in gene body regions of the downregulated genes for HCT cell line with concentration of 75nM and 300nM. In the figure, in case of 75nM, gene body bin 1-3, 6 and 8-10 shows higher H3K27Ac read density concentration, where in remaining bins its nearly equal among up and downregulated genes. Likewise, in case of 300nM, majority of the gene body shows higher H3K27Ac distribution then promoter region except bin 1, and 10. In revised manuscript, we have now revised the statement for clarity. 4. The length of GB is different for every gene. The bin size of the genes having larger GB size is higher compared to genes having shorter GB. As a longer bin having more numbers of nucleotides gives the possibility of higher numbers of H3K27Ac count. So, authors should consider the length of GB into their analysis and should check whether the pattern persists after normalisation by the length of the GB of a gene. Response: We appreciate this. However, we would like to clarify that our analysis already incorporates normalization that accounts for these differences . Specifically, as mentioned in the Methods section, we normalized the read counts within each genomic bin to TPM (Transcripts Per Million) using the Rsubread package. TPM normalization adjusts for both sequencing depth and bin length , thereby mitigating the potential bias introduced by longer bins containing more nucleotides. Thus, the enrichment patterns we observed are not confounded by gene length variation and reflect true biological signal distribution rather than technical artifacts. We have now made this point more explicit in the manuscript for clarity. 5. "As shown below, with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks and for the second gene set, the opposite was true (Figure 6B)". Again, these statements are not matching with the figures rather the opposite may be true or there is no clear pattern to mention. Kindly recheck your results and explain it properly. Response: We thank the reviewer for pointing out this mistake and we also regret the misinterpretation of the results. In the revised manuscript, we have clarified this for the better explanation. Page 7; Lines 217:220. “As shown below, we do not observe any clear pattern for the H3K27Ac read density among the two cell lines, except at TSS and gene body 1 region, where we observed lower H3K27Ac read distribution for the upregulated genes (Figure 6A&B)” . 6. The reference of the figure 5A is not there in the text. Response: Figure 5A is referenced on Page No. 7, Line No 201. View more View less Competing Interests No competing interests were disclosed. reply Respond Report a concern Mandal SD. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r243027) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-243027 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2024 Liao D et al. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 22 Feb 2024 | for Version 2 Daiqing Liao , Department of Anatomy and Cell Biology, College of Medicine & UF Health Cancer Center, University of Florida, Gainesville, USA Zhiguang Huo , Department of Biostatistics, University of Florida, Gainesville, Florida, USA 0 Views copyright © 2024 Liao D et al. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (0) Not Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Cancer epigenome underpins oncogenic gene expression. Histone acetylation is a major epigenomic event that impacts the transcriptome. The dynamic action of histone acetyltransferases (HATs) and deacetylases (HDACs) shape the acetylome. Small molecule inhibitors of both HATs and HDACs are promising agents for cancer therapy. Several HDAC inhibitors (HDACi) have been approved for treating cutaneous T-cell lymphoma and peripheral T-cell lymphoma. HDACi profoundly perturbs the transcriptome 1 . Thousands of genes can be downregulated and upregulated after treatment with an HDACi 2,3 . Pretreatment epigenome could influence transcriptomic perturbation by HDACi, although this has not been thoroughly examined. The manuscript titled “Predicting gene expression changes upon epigenomic drug treatment” by Agrawal, Gopalan, and Hannenhalli assessed whether the pretreatment epigenomic profile of cells could predict gene expression alterations induced by HDACi. The authors used the levels of a single histone mark, namely H3K27ac, in divided bins (21 genomic bins in total) from 2 kb upstream of the transcription start site (TSS) to the end of gene body and data of differentially expressed genes (DEGs) as inputs to study the relationship between the pretreatment H3K27ac levels and the transcription outputs in several machine learning models. The authors found that the pretreatment H3K27ac patterns could predict the post-treatment up versus downregulated genes with reasonable accuracy. Essentially, gene expression levels due to HDACi treatment seem to be inversely correlated with the levels of H3K27ac before treatment, which is observed in two cell lines treated with two HDACi of different chemical classes (entinostat and largazole). While the results are interesting, several significant issues should be addressed. The authors should also generate the feature importance score, e.g., the importance score for the random forest model (all gene models). This would be informative for readers to understand which genomic bin is most important to predict gene expression regulation. The authors should show the p-value for each comparison pair in Figure 3. The authors’ description “with few exceptions, for the first gene set, HCT116 genes showed higher concentration of H3K27Ac marks, and for the second gene set, the opposite was true (Figure 6B), consistent with the patterns in Figure 3” (page 9) seems contradictory to what is shown in Figure 6. This must be clarified. In Figure 6, the z-scores of normalized read counts are used, while normalized read counts are shown in Figure 3. What is the rationale for using different parameters? Also, the p-value for each pair comparison should be shown in Figure 6. Although the authors acknowledged that different epigenomic marks should be included in future studies, they are encouraged to analyze at least whether a different histone mark associated with active transcription, such as H3K4me3, would exhibit a similar predictive power as H3K27ac. Another major limitation of the manuscript is that data based on only two cell lines and two HDACi are analyzed. An expanded study of data from a diverse set of samples (cells and tumors) treated with additional HDACi should improve the validity and rigor of the study. Is the work clearly and accurately presented and does it cite the current literature? Partly Is the study design appropriate and is the work technically sound? Yes Are sufficient details of methods and analysis provided to allow replication by others? Yes If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Gryder BE, Wu L, Woldemichael GM, Pomella S, et al.: Chemical genomics reveals histone deacetylases are required for core regulatory transcription. Nat Commun . 2019; 10 (1): 3004 PubMed Abstract | Publisher Full Text 2. Lauffer BE, Mintzer R, Fong R, Mukund S, et al.: Histone deacetylase (HDAC) inhibitor kinetic rate constants correlate with cellular histone acetylation but not transcription and cell viability. J Biol Chem . 2013; 288 (37): 26926-43 PubMed Abstract | Publisher Full Text 3. Xiao Y, Hale S, Awasthee N, Meng C, et al.: HDAC3 and HDAC8 PROTAC dual degrader reveals roles of histone acetylation in gene regulation. Cell Chem Biol . 2023; 30 (11): 1421-1435.e12 PubMed Abstract | Publisher Full Text Competing Interests No competing interests were disclosed. Reviewer Expertise bioinformatics (Zhiguang Huo) and epigenetics (Daiqing Liao) We confirm that we have read this submission and believe that we have an appropriate level of expertise to state that we do not consider it to be of an acceptable scientific standard, for reasons outlined above. reply Respond to this report Responses (0) Liao D and Huo Z. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.159982.r239600) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v2#referee-response-239600 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2023 Merkel A. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 20 Sep 2023 | for Version 1 Angelika Merkel , Bioinformatics Unit, Josep Carreras Leukemia Research Institute (IJC), Barcelona, Spain 0 Views copyright © 2023 Merkel A. This is an open access peer review report distributed under the terms of the Creative Commons Attribution License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. format_quote Cite this report speaker_notes Responses (1) Not Approved info_outline Alongside their report, reviewers assign a status to the article: Approved The paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved Fundamental flaws in the paper seriously undermine the findings and conclusions Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Was peak calling performed to eliminate background noise in the Chip-seq analysis? Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? The GO analysis is repeated in the methods and results section Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text Competing Interests No competing interests were disclosed. Reviewer Expertise transcriptomics, epigenetics, bioinformatics I confirm that I have read this submission and believe that I have an appropriate level of expertise to state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above. reply Respond to this report Responses (1) Author Response 17 Jan 2024 Piyush Agrawal, Cancer Data Science Laboratory, National Cancer Institute, Bethesda, 20814, USA Agrawal et al have presented a comparison of machine learning approaches to predict gene expression from the distribution of histone marks, specifically H3K27 acetylation. Albeit interesting and potentially very useful for cancer specific and personalized treatment with histone deacetylase inhibitors as the authors suggest, the latter has actually not been addressed. The authors present a proof-of-concept study for predicting gene expression from histone marks but not do not address variations in tissues, individual or disease type. Response: We thank the reviewers for finding our work and interesting and potentially useful. Although we agree with reviewer regarding the importance of addressing variations in tissues, individual or disease type, unfortunately, lack of appropriate data currently prohibits such analysis. Besides the 2 cell lines used in this study, we didn’t come across any dataset (clinical/cell line) where for a given patient, both epigenomic marks and gene expression data is provided for both pre and post drug treatment. We would be happy to do the additional analysis if the reviewer can kindly point us to such data. As such, this is not a novel approach as there have been numerous studies before addressing the prediction of gene expression (e.g. Karlic et 2010 1 , Singh et al 2016 2 , Chen et al 2022 3 ). The authors have failed to cite any of the current literature on the topic. Response: We thank the reviewer for pointing out these studies, however, we want to clarify that the above-mentioned studies address a different question than ours, which is “to predict gene expression from the epigenome in a specific context”. Whereas we are interested in predicting the effect on the expression upon drug treatment, given the epigenomic profile in the pre-treatment sample. Our question necessitates availability of pre- and post-treatment gene expression and pre-treatment epigenomic profile, while these previous approaches are not concerned with the changes upon drug treatment. In the revised manuscript, we have highlighted this difference in the Introduction. Methods and results are presented in a somewhat rough manner and need improvement: Abstract and introduction contain poor grammar and sentence structure ('genomic specificity" is not a new term to be defined, better to to use target specificity) Response: Why were the HCT166 and RH4 cell lines chosen and what are these exactly (need description)? Response: We regret lack of clarity. HCT116 and RH4 cell lines were chosen because the data we required for our study was available only for these 2 cell lines, namely, pre-treatment epigenomic profile and pre- and post-treatment gene expression data. HCT116 is human colorectal carcinoma cell line initiated from an adult male whereas RH4 cell line is for studying alveolar rhabdomyosarcoma and is belongs to soft tissue lineage. In the revised manuscript, we have added the description of these cell lines in the Methods section. Was peak calling performed to eliminate background noise in the Chip-seq analysis? Response: No, peak calling wasn’t performed in the current study. The goal of our analysis was to compare changes in H3K27ac read density upon HDAC treatment. Peak calling represents a 0/1 binarization of a continuous variable (in this case, normalized read counts in promoters and gene bodies). Our model uses the epigenomic signal intensity in specific windows around the gene magnitude, making peak-calling unnecessary. Was differential expression analysis performed (normalization, statistical analysis, etc.)? What is the significance of the log10fold change? Response: The mean Log-TPM values of each gene was computed for both treated and untreated conditions, with the Log-FC computed as the difference between treated and untreated conditions. The top 1000 most up-regulated and 1000 most down-regulated genes were considered to be differentially expressed. We mistakenly mentioned Log10, which is now corrected to Log2. We have corrected this in the revised version . The GO analysis is repeated in the methods and results section Response: We want to clarify that in Methods Section, we mentioned the software and command only used for GO analysis, however, in Result section, we discussed the enriched terms associated with the up and downregulated genes. Figure 3: The scale among plots differ; gene bins are of variable size and larger promoter bins. Why have they been plotted them together? Response: We agree with the reviewer that scale among plots differ and gene bins are of variable size. The only reason for plotting them together was to analyze all the 21 features together in single image i.e. change in distribution of the pre-treatment epigenomic marks in up and down genes post treatment. Is the work clearly and accurately presented and does it cite the current literature? No Is the study design appropriate and is the work technically sound? Partly Are sufficient details of methods and analysis provided to allow replication by others? Partly If applicable, is the statistical analysis and its interpretation appropriate? Partly Are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions drawn adequately supported by the results? Partly References 1. Karlić R, Chung HR, Lasserre J, Vlahovicek K, et al.: Histone modification levels are predictive for gene expression. Proc Natl Acad Sci U S A . 2010; 107 (7): 2926-31 PubMed Abstract | Publisher Full Text 2. Singh R, Lanchantin J, Robins G, Qi Y: DeepChrome: deep-learning for predicting gene expression from histone modifications. Bioinformatics . 2016; 32 (17): i639-i648 PubMed Abstract | Publisher Full Text 3. Chen Y, Xie M, Wen J: Predicting gene expression from histone modifications with self-attention based neural networks and transfer learning. Front Genet . 2022; 13 : 1081842 PubMed Abstract | Publisher Full Text View more View less Competing Interests No competing interests reply Respond Report a concern Merkel A. Peer Review Report For: Predicting gene expression changes upon epigenomic drug treatment [version 3; peer review: 4 approved, 2 not approved] . F1000Research 2025, 12 :1089 ( https://doi.org/10.5256/f1000research.153610.r203232) NOTE: it is important to ensure the information in square brackets after the title is included in this citation. The direct URL for this report is: https://f1000research.com/articles/12-1089/v1#referee-response-203232 Alongside their report, reviewers assign a status to the article: Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit. Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions Adjust parameters to alter display View on desktop for interactive features Includes Interactive Elements View on desktop for interactive features Competing Interests Policy Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list: Examples of 'Non-Financial Competing Interests' Within the past 4 years, you have held joint grants, published or collaborated with any of the authors of the selected paper. You have a close personal relationship (e.g. parent, spouse, sibling, or domestic partner) with any of the authors. You are a close professional associate of any of the authors (e.g. scientific mentor, recent student). You work at the same institute as any of the authors. You hope/expect to benefit (e.g. favour or employment) as a result of your submission. You are an Editor for the journal in which the article is published. Examples of 'Financial Competing Interests' You expect to receive, or in the past 4 years have received, any of the following from any commercial organisation that may gain financially from your submission: a salary, fees, funding, reimbursements. You expect to receive, or in the past 4 years have received, shared grant support or other funding with any of the authors. You hold, or are currently applying for, any patents or significant stocks/shares relating to the subject matter of the paper you are commenting on. Stay Updated Sign up for content alerts and receive a weekly or monthly email with all newly published articles Register with F1000Research Already registered? Sign in Not now, thanks close PLEASE NOTE If you are an AUTHOR of this article, please check that you signed in with the account associated with this article otherwise we cannot automatically identify your role as an author and your comment will be labelled as a “User Comment”. If you are a REVIEWER of this article, please check that you have signed in with the account associated with this article and then go to your account to submit your report, please do not post your review here. If you do not have access to your original account, please contact us . All commenters must hold a formal affiliation as per our Policies . The information that you give us will be displayed next to your comment. User comments must be in English, comprehensible and relevant to the article under discussion. We reserve the right to remove any comments that we consider to be inappropriate, offensive or otherwise in breach of the User Comment Terms and Conditions . Commenters must not use a comment for personal attacks. When criticisms of the article are based on unpublished data, the data should be made available. I accept the User Comment Terms and Conditions Please confirm that you accept the User Comment Terms and Conditions. Affiliation ✕ refresh Please enter your institution. Note: To add your institution or organisation, start typing the name and then select the correct name from the list. Where applicable, the name will appear in both the original language and in English. Do not paste in the name. If the name does not appear in the drop-down list, we will display the information you have entered. ✕ refresh Country/Region * USA UK Canada China France Germany Afghanistan Aland Islands Albania Algeria American Samoa Andorra Angola Anguilla Antarctica Antigua and Barbuda Argentina Armenia Aruba Australia Austria Azerbaijan Bahamas Bahrain Bangladesh Barbados Belarus Belgium Belize Benin Bermuda Bhutan Bolivia Bosnia and Herzegovina Botswana Bouvet Island Brazil British Indian Ocean Territory British Virgin Islands Brunei Bulgaria Burkina Faso Burundi Cambodia Cameroon Canada Cape Verde Cayman Islands Central African Republic Chad Chile China Christmas Island Cocos (Keeling) Islands Colombia Comoros Congo Cook Islands Costa Rica Cote d'Ivoire Croatia Cuba Cyprus Czech Republic Democratic Republic of the Congo Denmark Djibouti Dominica Dominican Republic Ecuador Egypt El Salvador Equatorial Guinea Eritrea Estonia Ethiopia Falkland Islands Faroe Islands Federated States of Micronesia Fiji Finland France French Guiana French Polynesia French Southern Territories Gabon Georgia Germany Ghana Gibraltar Greece Greenland Grenada Guadeloupe Guam Guatemala Guernsey Guinea Guinea-Bissau Guyana Haiti Heard Island and Mcdonald Islands Holy See (Vatican City State) Honduras Hong Kong Hungary Iceland India Indonesia Iran Iraq Ireland Israel Italy Jamaica Japan Jersey Jordan Kazakhstan Kenya Kiribati Kosovo (Serbia and Montenegro) Kuwait Kyrgyzstan Lao People's Democratic Republic Latvia Lebanon Lesotho Liberia Libya Liechtenstein Lithuania Luxembourg Macao Madagascar Malawi Malaysia Maldives Mali Malta Marshall Islands Martinique Mauritania Mauritius Mayotte Mexico Minor Outlying Islands of the United States Moldova Monaco Mongolia Montenegro Montserrat Morocco Mozambique Myanmar Namibia Nauru Nepal Netherlands Antilles New Caledonia New Zealand Nicaragua Niger Nigeria Niue Norfolk Island North Korea North Macedonia Northern Mariana Islands Norway Oman Pakistan Palau Palestinian Territory Panama Papua New Guinea Paraguay Peru Philippines Pitcairn Poland Portugal Puerto Rico Qatar Reunion Romania Russian Federation Rwanda Saint Helena Saint Kitts and Nevis Saint Lucia Saint Pierre and Miquelon Saint Vincent and the Grenadines Samoa San Marino Sao Tome and Principe Saudi Arabia Senegal Serbia Seychelles Sierra Leone Singapore Slovakia Slovenia Solomon Islands Somalia South Africa South Georgia and the South Sandwich Is South Korea South Sudan Spain Sri Lanka Sudan Suriname Svalbard and Jan Mayen Swaziland Sweden Switzerland Syria Taiwan Tajikistan Tanzania Thailand The Gambia The Netherlands Timor-Leste Togo Tokelau Tonga Trinidad and Tobago Tunisia Turkey Turkmenistan Turks and Caicos Islands Tuvalu UK USA Uganda Ukraine United Arab Emirates United States Virgin Islands Uruguay Uzbekistan Vanuatu Venezuela Vietnam Wallis and Futuna West Bank and Gaza Strip Western Sahara Yemen Zambia Zimbabwe Please select your country/region. You must enter a comment. Competing Interests Please disclose any competing interests that might be construed to influence your judgment of the article's or peer review report's validity or importance. Competing Interests Policy Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list: Examples of 'Non-Financial Competing Interests' Within the past 4 years, you have held joint grants, published or collaborated with any of the authors of the selected paper. You have a close personal relationship (e.g. parent, spouse, sibling, or domestic partner) with any of the authors. You are a close professional associate of any of the authors (e.g. scientific mentor, recent student). You work at the same institute as any of the authors. You hope/expect to benefit (e.g. favour or employment) as a result of your submission. You are an Editor for the journal in which the article is published. Examples of 'Financial Competing Interests' You expect to receive, or in the past 4 years have received, any of the following from any commercial organisation that may gain financially from your submission: a salary, fees, funding, reimbursements. You expect to receive, or in the past 4 years have received, shared grant support or other funding with any of the authors. You hold, or are currently applying for, any patents or significant stocks/shares relating to the subject matter of the paper you are commenting on. Please state your competing interests The comment has been saved. An error has occurred. Please try again. Cancel Post var lTitle = "Predicting gene expression changes upon epigenomic...".replace("'", ''); var linkedInUrl = "http://www.linkedin.com/shareArticle?url=https://f1000research.com/articles/12-1089/v3" + "&title=" + encodeURIComponent(lTitle) + "&summary=" + encodeURIComponent('Read the article by '); var deliciousUrl = "https://del.icio.us/post?url=https://f1000research.com/articles/12-1089/v3&title=" + encodeURIComponent(lTitle); var redditUrl = "http://reddit.com/submit?url=https://f1000research.com/articles/12-1089/v3" + "&title=" + encodeURIComponent(lTitle); linkedInUrl += encodeURIComponent('Agrawal P et al.'); var offsetTop = /chrome/i.test( navigator.userAgent ) ? 4 : -10; var addthis_config = { ui_offset_top: offsetTop, services_compact : "facebook,twitter,www.linkedin.com,www.mendeley.com,reddit.com", services_expanded : "facebook,twitter,www.linkedin.com,www.mendeley.com,reddit.com", services_custom : [ { name: "LinkedIn", url: linkedInUrl, icon:"/img/icon/at_linkedin.svg" }, { name: "Mendeley", url: "http://www.mendeley.com/import/?url=https://f1000research.com/articles/12-1089/v3/mendeley", icon:"/img/icon/at_mendeley.svg" }, { name: "Reddit", url: redditUrl, icon:"/img/icon/at_reddit.svg" }, ] }; var addthis_share = { url: "https://f1000research.com/articles/12-1089", templates : { twitter : "Predicting gene expression changes upon epigenomic drug treatment. Agrawal P et al., published by " + "@F1000Research" + ", https://f1000research.com/articles/12-1089/v3" } }; if (typeof(addthis) != "undefined"){ addthis.addEventListener('addthis.ready', checkCount); addthis.addEventListener('addthis.menu.share', checkCount); } $(".f1r-shares-twitter").attr("href", "https://twitter.com/intent/tweet?text=" + addthis_share.templates.twitter); $(".f1r-shares-facebook").attr("href", "https://www.facebook.com/sharer/sharer.php?u=" + addthis_share.url); $(".f1r-shares-linkedin").attr("href", addthis_config.services_custom[0].url); $(".f1r-shares-reddit").attr("href", addthis_config.services_custom[2].url); $(".f1r-shares-mendelay").attr("href", addthis_config.services_custom[1].url); function checkCount(){ setTimeout(function(){ $(".addthis_button_expanded").each(function(){ var count = $(this).text(); if (count !== "" && count != "0") $(this).removeClass("is-hidden"); else $(this).addClass("is-hidden"); }); }, 1000); } close How to cite this report {{reportCitation}} Cancel Copy Citation Details $(function(){R.ui.buttonDropdowns('.dropdown-for-downloads');}); $(function(){R.ui.toolbarDropdowns('.toolbar-dropdown-for-downloads');}); $.get("/articles/acj/140273/181146") new F1000.Clipboard(); new F1000.ThesaurusTermsDisplay("articles", "article", "181146"); $(document).ready(function() { $( "#frame1" ).on('load', function() { var mydiv = $(this).contents().find("div"); var h = mydiv.height(); console.log(h) }); var tooltipLivingFigure = jQuery(".interactive-living-figure-label .icon-more-info"), titleLivingFigure = tooltipLivingFigure.attr("title"); tooltipLivingFigure.simpletip({ fixed: true, position: ["-115", "30"], baseClass: 'small-tooltip', content:titleLivingFigure + " " }); tooltipLivingFigure.removeAttr("title"); $("body").on("click", ".cite-living-figure", function(e) { e.preventDefault(); var ref = $(this).attr("data-ref"); $(this).closest(".living-figure-list-container").find("#" + ref).fadeIn(200); }); $("body").on("click", ".close-cite-living-figure", function(e) { e.preventDefault(); $(this).closest(".popup-window-wrapper").fadeOut(200); }); $(document).on("mouseup", function(e) { var metricsContainer = $(".article-metrics-popover-wrapper"); if (!metricsContainer.is(e.target) && metricsContainer.has(e.target).length === 0) { $(".article-metrics-close-button").click(); } }); var articleId = $('#articleId').val(); if($("#main-article-count-box").attachArticleMetrics) { $("#main-article-count-box").attachArticleMetrics(articleId, { articleMetricsView: true }); } }); var figshareWidget = $(".new_figshare_widget"); if (figshareWidget.length > 0) { window.figshare.load("f1000", function(Widget) { // Select a tag/tags defined in your page. In this tag we will place the widget. _.map(figshareWidget, function(el){ var widget = new Widget({ articleId: $(el).attr("figshare_articleId") //height:300 // this is the height of the viewer part. [Default: 550] }); widget.initialize(); // initialize the widget widget.mount(el); // mount it in a tag that's on your page // this will save the widget on the global scope for later use from // your JS scripts. This line is optional. //window.widget = widget; }); }); } close Error Close Add Reset F1000.MICROSERVICES.AFFILIATION = ''; $(document).ready(function () { $('.js-affiliations-form').each((index, form) => { new AffiliationForm({ formId: form.id, institutionErrorSelector: '.comment-enter-institution', departmentErrorSelector: '.comment-enter-department', placeSelector: '.js-add-comment-place', stateSelector: '.js-add-comment-state', zipCodeSelector: '.js-add-comment-zipcode', countrySelector: '.js-add-comment-country', countryErrorSelector: '.comment-enter-country', }); }); }); $(document).ready(function () { var reportIds = { "254720": 0, "254727": 0, "254726": 0, "254725": 0, "254730": 0, "231438": 0, "208151": 0, "267046": 0, "267047": 9, "383149": 0, "267052": 0, "383148": 0, "383151": 0, "383150": 8, "267055": 0, "267048": 0, "267049": 0, "267051": 0, "267060": 0, "382129": 6, "383152": 5, "267057": 0, "382131": 0, "267058": 0, "382130": 0, "243027": 29, "243031": 0, "243030": 0, "243029": 0, "243035": 0, "243034": 0, "243033": 0, "243032": 0, "203231": 0, "203230": 0, "203229": 0, "243037": 0, "243036": 0, "279396": 0, "279397": 0, "279398": 0, "203232": 70, "279399": 0, "258919": 0, "258918": 0, "279394": 0, "258917": 0, "279395": 0, "258916": 0, "258923": 0, "239595": 0, "258922": 0, "258921": 0, "258920": 0, "279400": 0, "239599": 0, "239598": 0, "258925": 0, "239597": 0, "258924": 0, "239596": 0, "239603": 0, "239602": 0, "239601": 0, "239600": 39, "254711": 0, "239604": 0, "254708": 0, "254714": 0, "254718": 0, "254716": 0, }; $(".referee-response-container,.js-referee-report").each(function(index, el) { var reportId = $(el).attr("data-reportid"), reportCount = reportIds[reportId] || 0; $(el).find(".comments-count-container,.js-referee-report-views").html(reportCount); }); var uuidInput = $("#article_uuid"), oldUUId = uuidInput.val(), newUUId = "5e5fadf2-8d6c-4b7b-bff9-92a9bcd4d31b"; uuidInput.val(newUUId); $("a[href*='article_uuid=']").each(function(index, el) { var newHref = $(el).attr("href").replace(oldUUId, newUUId); $(el).attr("href", newHref); }); }); An innovative open access publishing platform offering rapid publication and open peer review, whilst supporting data deposition and sharing. Browse Gateways Collections How it Works Contact For Developers Cookie Notice Privacy Notice RSS Submit Your Research Follow us © 2012-2026 F1000 Research Ltd. ISSN 2046-1402 | Legal | Partner of Research4Life • CrossRef • ORCID • FAIRSharing R.templateTests.simpleTemplate = R.template(' $text $text $text $text $text '); R.templateTests.runTests(); var F1000platform = new F1000.Platform({ name: "f1000research", displayName: "F1000Research", hostName: "f1000research.com", id: "1", editorialEmail: "[email protected]", infoEmail: "[email protected]", usePmcStats: true }); $(function(){R.ui.dropdowns('.dropdown-for-authors, .dropdown-for-about, .dropdown-for-myresearch');}); // $(function(){R.ui.dropdowns('.dropdown-for-referees');}); $(document).ready(function () { if ($(".cookie-warning").is(":visible")) { $(".sticky").css("margin-bottom", "35px"); $(".devices").addClass("devices-and-cookie-warning"); } $(".cookie-warning .close-button").click(function (e) { $(".devices").removeClass("devices-and-cookie-warning"); $(".sticky").css("margin-bottom", "0"); }); $("#tweeter-feed .tweet-message").each(function (i, message) { var self = $(message); self.html(linkify(self.html())); }); $(".partner").on("mouseenter mouseleave", function() { $(this).find(".gray-scale, .colour").toggleClass("is-hidden"); }); }); Sign In Remember me Forgotten your password? Sign In Cancel Email or password not correct. Please try again Please wait... $(function(){ // Note: All the setup needs to run against a name attribute and *not* the id due the clonish // nature of facebox... $("a[id=googleSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("GOOGLE"); $("form[id=oAuthForm]").submit(); }); $("a[id=facebookSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("FACEBOOK"); $("form[id=oAuthForm]").submit(); }); $("a[id=orcidSignInButton]").click(function(event){ event.preventDefault(); $("input[id=oAuthSystem]").val("ORCID"); $("form[id=oAuthForm]").submit(); }); }); If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password. The email address should be the one you originally registered with F1000. Email address not valid, please try again You registered with F1000 via Google, so we cannot reset your password. To sign in, please click here . If you still need help with your Google account password, please click here . You registered with F1000 via Facebook, so we cannot reset your password. To sign in, please click here . If you still need help with your Facebook account password, please click here . Code not correct, please try again Reset password Cancel Email us for further assistance. Server error, please try again. If your email address is registered with us, we will email you instructions to reset your password. If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance. Please wait... Register $(document).ready(function () { signIn.createSignInAsRow($("#sign-in-form-gfb-popup")); $(".target-field").each(function () { var uris = $(this).val().split("/"); if (uris.pop() === "login") { $(this).val(uris.toString().replace(",","/")); } }); });