Unsupervised gene selection for predicting cell... | 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/9-124" }, "headline": "Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo", "datePublished": "2020-02-19T11:54:34", "dateModified": "2021-06-07T17:55:21", "author": [ { "@type": "Person", "name": "Yang Chen" }, { "@type": "Person", "name": "Disheng Mao" }, { "@type": "Person", "name": "Yuping Zhang" }, { "@type": "Person", "name": "Zhengqing Ouyang" } ], "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": "Analyzing single cell RNA-seq data is important for deciphering the spatial relationships, expression patterns, and developmental processes of cells. Combining in situ hybridization-based gene expression atlas images, some works have successfully recovered spatial locations of cells in zebrafish and Drosophila embryos. In this article, we describe a highly ranked method in the DREAM Single Cell Transcriptomics Challenge for predicting cell positions in the Drosophila embryo. The method performs unsupervised feature extraction to select a small number of driver genes and then uses them to predict gene expression and spatial position of each individual cell. First, hierarchical clustering is used to select a subset of driver genes. Second, the similarity matrix of single cells in the bins of the reference atlas is computed. Based on the similarity matrix, the spatial positions of cells are then determined by hierarchical clustering. This method is evaluated on the cell positions and gene expressions in the DREAM Single Cell Transcriptomics Challenge. The comparison with the “silver standard” suggests that our method is effective in reconstructing the cell spatial positions and gene expression patterns in tissues." } { "@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/9-124", "name": "Unsupervised gene selection for predicting cell spatial positions..." } } ] } Home Browse Unsupervised gene selection for predicting cell spatial positions... ALL Metrics - Views Downloads Get PDF Get XML Cite How to cite this article Chen Y, Mao D, Zhang Y and Ouyang Z. Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.12688/f1000research.20446.2 ) 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 ▬ ✚ Method Article Revised Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] Previously titled: An unsupervised learning method for reconstructing cell spatial organization with application to the DREAM Single Cell Transcriptomics Challenge Yang Chen https://orcid.org/0000-0003-0147-3443 1 , Disheng Mao 2 , Yuping Zhang 2 , Zhengqing Ouyang https://orcid.org/0000-0003-2842-8503 1 Yang Chen https://orcid.org/0000-0003-0147-3443 1 , Disheng Mao 2 , Yuping Zhang 2 , Zhengqing Ouyang https://orcid.org/0000-0003-2842-8503 1 PUBLISHED 07 Jun 2021 Author details Author details 1 Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA, USA 2 Department of Statistics, University of Connecticut, Storrs, CT, USA Yang Chen Roles: Data Curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – Original Draft Preparation Disheng Mao Roles: Investigation, Methodology Yuping Zhang Roles: Funding Acquisition, Investigation, Methodology, Resources, Supervision, Writing – Review & Editing Zhengqing Ouyang Roles: Conceptualization, Formal Analysis, Funding Acquisition, Investigation, Methodology, Project Administration, Resources, Supervision, Writing – Review & Editing OPEN PEER REVIEW DETAILS REVIEWER STATUS Abstract Analyzing single cell RNA-seq data is important for deciphering the spatial relationships, expression patterns, and developmental processes of cells. Combining in situ hybridization-based gene expression atlas images, some works have successfully recovered spatial locations of cells in zebrafish and Drosophila embryos. In this article, we describe a highly ranked method in the DREAM Single Cell Transcriptomics Challenge for predicting cell positions in the Drosophila embryo. The method performs unsupervised feature extraction to select a small number of driver genes and then uses them to predict gene expression and spatial position of each individual cell. First, hierarchical clustering is used to select a subset of driver genes. Second, the similarity matrix of single cells in the bins of the reference atlas is computed. Based on the similarity matrix, the spatial positions of cells are then determined by hierarchical clustering. This method is evaluated on the cell positions and gene expressions in the DREAM Single Cell Transcriptomics Challenge. The comparison with the “silver standard” suggests that our method is effective in reconstructing the cell spatial positions and gene expression patterns in tissues. READ ALL READ LESS Keywords Spatial organization, single cell RNA-seq, Drosophila embryo, clustering, DREAM challenge Corresponding Author(s) Zhengqing Ouyang ( [email protected] ) Close Corresponding author: Zhengqing Ouyang Competing interests: No competing interests were disclosed. Grant information: This work was partially supported by the Faculty Research Excellence Program Award at UConn (to YZ). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Copyright: © 2021 Chen Y 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. How to cite: Chen Y, Mao D, Zhang Y and Ouyang Z. Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.12688/f1000research.20446.2 ) First published: 19 Feb 2020, 9 :124 ( https://doi.org/10.12688/f1000research.20446.1 ) Latest published: 07 Jun 2021, 9 :124 ( https://doi.org/10.12688/f1000research.20446.2 ) Revised Amendments from Version 1 The title is modified to avoid ambiguity. The abstract is modified to improve the clarity of method summary. More citations and description of existing methods are added. All figures and legends are modified to describe the results more clearly. More details are added on binarization of scRNA-seq data and prediction of spatial gene expression. More discussions are added on the interpretation of the results as well as the advantage and limitation of the method. The title is modified to avoid ambiguity. The abstract is modified to improve the clarity of method summary. More citations and description of existing methods are added. All figures and legends are modified to describe the results more clearly. More details are added on binarization of scRNA-seq data and prediction of spatial gene expression. More discussions are added on the interpretation of the results as well as the advantage and limitation of the method. See the authors' detailed response to the review by Xianwen Ren See the authors' detailed response to the review by Mark S. Cembrowski and Larissa Kraus READ REVIEWER RESPONSES Introduction The development of single cell RNA sequencing (scRNA-seq) has provided a powerful solution for building a global transcription landscape of all cell types in tissues, finding new cell types, cell lineage tracing, spatial reconstruction, and combining with other omics 1 – 9 . Single cells are originally made from dissociated tissues without spatial information, and spatial gene expression pattern is unknown. In situ hybridization (ISH) and its variants can detect the spatial information of mRNA transcripts and produce gene expression reference atlas. Using enough marker genes, users can reconstruct the spatial position of single cell RNA-seq data by combing the ISH reference atlas 10 – 15 . Some works have also combined sequential fluorescence in situ hybridization (seqFISH) and multiplexed error robust fluorescence in situ hybridization (MERFISH) with scRNA-seq data to map cell types to the reference atlas 16 – 18 . Recent methods have successfully mapped cells from scRNA-seq data to the spatial positions using dozens of landmark genes 10 – 15 . Nikos et al. developed a DistMap method for mapping the ~1,300 Drosophila embryo cells into ~3,000 bins in the spatial position using 84 marker genes by computing the Matthews Correlation Coefficient (MCC) score between in situ gene expression and binarized scRNA-seq data 10 . Satija et al. mapped 851 cells of zebrafish embryo into 64 bins of the spatial embryo with 47 genes. They used a bimodal mixture model to binarize the scRNA-seq data and predict a cell originated from the bin by posterior probability 11 . Kaia et al. computed the correspondence score and mapped 139 cells into a spatial Platynereis dumerilii brain using a set of 98 genes by whole-mount in situ hybridization (WMISH). They binarized the scRNA-seq data and defined a special score to characterize the relation between every cell-voxel combination. They suggested that 50–100 in situ genes were needed to map cells to spatial position with high confidence using simulated spatial reference atlases 12 . Andreas et al. also reconstructed spatial single enterocytes along the villus axis in 1-D space using 102 bottom and top landmark genes from single-molecule fluorescence in situ hybridization (smFISH 19 ) for 1,383 single cells. They used the summed expression of landmark genes to infer the locations along the 1-D coordinates and tried 50 and 20 landmark genes for prediction. The reconstruction errors were relatively low for 50 genes, but the prediction errors increased remarkably using 20 landmark genes 13 . Nitzan et al. also proposed a new method named novoSpaRc for spatially mapping the scRNA-seq cells into an existing reference atlas. Based on the optimal transport method, it reported the gene expression and spatial position with at most 2 landmark genes for virtual scRNA-seq data and virtual Drosophila embryo 14 . The performance of novoSpaRc is still not good for real Drosophila embryo when the marker gene number is smaller than 30 14 , 15 . In the methods mentioned above, the dimension and resolution of the spatial region, as well as the number of marker genes, are key factors that affect the recovery of the spatial position. Although these methods demonstrated good performance, they do not focus on markers genes selection in which biological information were used to determine the spatial gene expression patterns. To explore new algorithms that include biological information contained in the genes to determine embryonic patterning, the DREAM Single Cell Transcriptomics Challenge asked participating teams to predict the positions in the embryo of 1,297 cells using the expression patterns of 60 (sub-challenge 1), 40 (sub-challenge 2), and 20 (sub-challenge 3) driver genes from the in situ hybridization data 20 . The challenge endeavors to use less marker genes to infer the spatial locations of cells. In this article, we assume that cells with similar gene expression will have close spatial distances. By clustering the possible spatial bins for a cell, we select the most possible cluster as the predicted spatial position. The details of our method are the following: First, we use hierarchical clustering for gene selection. Then, we map the cells to spatial positions by computing the mapping similarity between the binarized scRNA-seq and in situ hybridization data and cluster the bins. Finally, we select the optimal thresholds by computing the maximum MCC score for spatial gene expression prediction. We detail our solution for the DREAM challenge, and validate the results using the “silver standard” derived from DistMap. The paper is organized as follows: in Methods, we briefly describe the solutions for all the three sub-challenges; in Results, we present the results of the three sub-challenges in the DREAM Single Cell Transcriptomics Challenge; finally, we discuss and summarize our work. Methods Dataset The scRNA-seq dataset is from ~1,000 handpicked stage 6 fly embryos using Drop-seq 21 . It contains both raw and normalized UMI counts with 1,297 cells and 8,924 genes per cell. In situ hybridization expression patterns of 84 driver genes are from the Berkeley Drosophila Transcription Network Project (BDNTP). The BDNTP reference atlas are binarized. The bin number of one half of the symmetric Drosophila melanogaster embryo at stage 6-pre-gastrulation (which is used in this study) is 3,039. The spatial coordinates of these bins are also specified. The dataset files can be downloaded from the DREAM Single Cell Transcriptomics Challenge after registration with Synapse free of charge ( https://www.synapse.org/#!Synapse:syn16782360 ). We directly use the normalized scRNA-seq data, the in situ matrix and the geometry of the embryo. The gene names “ E.spl.m5.HLH ” and “ Blimp.1 ” are replaced by “ E(spl)m5-HLH ” and “ Blimp-1 ”. The DistMap “silver standard” is the best available cell position reference, which is determined by calculating the maximum similarity for cells and spatial bins using all the 84 in situ driver genes in BDNTP. It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells 10 . Gene selection We use a hierarchical clustering method to select 60, 40, and 20 driver genes from the 84 genes based on the normalized scRNA-seq data ( Figure 1 ). Specifically, based on the belief that the scRNA-seq gene expression pattern is driven by the driver genes’ expression pattern, we propose to select the essential driver genes based on the information provided by scRNA-seq data. If two genes have high correlation in the scRNA-seq data, we choose only one of them without losing too much of the information. To find the correlated genes, we perform hierarchical clustering on the normalized scRNA-seq data to separate all 84 genes into 60 clusters (with the Euclidean distance and the Mcquitty linkage). The Mcquitty linkage gives more weights for objects in small clusters than those in large clusters in calculating the distance between two clusters. Thus, it is suitable for situations with many small clusters. Since the numbers of clusters are fairly large in the sub-challenge 1 and 2, we opt to use the Mcquitty linkage for distance calculation. In sub-challenge 3, since the total number of clusters is shrunk to 20, which is smaller than sub-challenge 1 and 2, we choose to use the ward linkage in the hierarchical clustering part to obtain larger-sized clusters from the data. After this step, the gene selection process remains the same as sub-challenge 1 and 2. Figure 1. Workflow of our spatial position prediction method for scRNA-seq with in situ reference gene panel. ( 1 ) Use hierarchical clustering (HC) to select 60/40/20 driver genes from the given 84 driver genes based on scRNA-seq data with the Euclidean distance and the Mcquitty (for 60 and 40 genes) and ward linkage (for 20 genes). ( 2 ) Apply the binarization method for the selected genes by minimizing root-mean-square error between of Pearson correlation matrix between the binarized scRNA-seq and the in situ matrix, and use the same criteria to binarize the remaining 24/44/64 predicted genes. ( 3 ) Calculate mapping similarity matrix from single cells to spatial bins based on the proportion of matched gene expression between each cell and each bin. ( 4 ) For each cell, select top bins based on mapping similarity values by quantile thresholding. Then do HC on the spatial distance of these top bins and select the number of clusters based on the Silhouette score. The top 10 predicted bins are the ones surrounding the largest cluster center. ( 5 ) For each gene, we compute the MCC score between predicted gene expression and binarize gene expression given different thresholds. The optimal threshold is corresponding to the maximum MCC. Finally, we apply the optimal threshold for the computeVISH () method in DistMap to predict gene expression. The black points in the scatter plot are the predicted MCC. The red points are Pearson correlation between predicted in situ gene expression with in situ gene expression. After getting the clusters, we pick the most representative gene of each cluster by calculating the distance between each member gene and the cluster center based on the Euclidean distance and selecting the closest one. Binarization of scRNA-seq data We perform binarization on the normalized scRNA-seq data for the selected genes based on the “binarizeSingleCellData()” function in DistMap ( https://github.com/rajewsky-lab/distmap ). The details of binarization is as follows: for each quantiles from 0.15 to 0.5 with 0.01 as interval, we perform binarization on the scRNA-seq data for 60, 40, and 20 driver genes in sub-challenge 1, 2, 3, respectively. If the gene expression value is larger than the quantile gene expression value, it will be set as 1, otherwise it is 0. Then we compute the difference between the correlation matrix of the binarized scRNA-seq data and that of the in situ matrix based on the root-mean square error. Finally, we select the quantile threshold which has the smallest root mean square error. Next, to predict the spatial gene expression for all the 84 in situ genes, we used a similar method to binarize the remaining 24, 44, and 64 in situ genes in sub-challenges with the binarized 60, 40, and 20 binarized gene expression. First, we compute the Pearson correlation matrix C 1 between the selected genes (60, 40, and 20) and the remaining genes (24, 44, and 64) using the normalized scRNA-seq data. Second, we binarize the remaining genes based on a series of quantiles from 0.15 to 0.5 with 0.01 as interval. Furthermore, we calculate the Pearson correlation matrix C 2 between the binarized selected genes and the remaining genes. Finally, we select the best quantile for binarizing the remaining genes by minimizing the root mean square error between the Pearson correlation matrices C 1 and C 2 . Compute the similarity matrix between cells and bins Given the binarized scRNA-seq data and the in situ matrix, we calculate the similarity matrix between the cells and bins based on the selected driver genes. Here, we denote the number of selected driver genes is n g . The similarity between a cell c i (i ∈ [1,1,297]) and a bin b j (j ∈ [1,3,039]) is calculated as follows: p i j = n s i j n g , ( 1 ) where n s i j is the number of the same gene expression values (0 or 1) in the two binarized vectors corresponding to cell c i and bin b j , and n g is the in situ matrix for cell c i and bin b j . Select candidate the cell positions based on the similarity matrix Select top bins. The likelihood of a cell mapped to a bin is determined by the gene expression in the bin and the cell. The bins with the higher similarity are possibly the potential cell position. To generate stable results, we select enough bins (see below) based on the similarity values. Then we use clustering to determine a more refined cell position. To select the potential bins for predicting cell position, we check the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell from sub-challenge 1, 2, and 3 ( Figure 2 ). The average number of bins with mapping similarity > 0.8 increases as the number of driver genes decreases. In sub-challenge 3, there are more bins with high similarity values produced. The small similarity difference among these bins makes it harder to determine real spatial bins where the cell originated. Thus, as the number of driver genes decreases, the accuracy of bins with high similarity values also decreases. Here, we use the third quartile of the similarity values to select the top bins in sub-challenge 1 with 60 driver genes. We use the first quartile of similarity values in sub-challenge 2 with 40 driver genes. And we use all bins in sub-challenge 3 with 20 driver genes. If the number of selected bins is 0, the bin which has the maximum similarity will be the predicted position. If the number of the selected bins is larger than 100, then only the top 100 bins will be kept based on the similarity values. Figure 2. Compare the number of high similarity values (> 0.8) in the similarity matrix in each bin for sub-challenge 1, 2, and 3. The numbers of high similarity values in sub-challenge 1 and 2 are much smaller than those in sub-challenge 3. Silhouette score to determine the number of clusters of bins For bins with high similarity values, we need to perform clustering to select the cluster which has the maximum sum of similarity as the predicted cell position. Here, we use hierarchical clustering on the selected bins based on their three-dimensional distance. The cluster number is determined by the silhouette score, which measures the average distance of a point to other points in its cluster compared to the smallest average distance to other clusters. The silhouette score ranges from -1 to +1. The higher the silhouette score, the closer the point is closer to its own cluster and the farther it is away from other clusters. We use the average silhouette score across all selected bins to select the optimal clustering number. Here, we use the NbClust package 22 to perform hierarchical clustering with the “centroid” method. Based on the silhouette score, we obtain the number of clusters. Predict cell positions based on clustering of bins We compute the sum of similarities for each cluster and select the cluster with the maximum value. Finally, we assign 10 nearest bins to the selected cluster center as the top 10 most possible cell positions. Predict spatial gene expression based on the similarity matrix To predict the final spatial in situ gene expression, we use the method computeVISH () in DistMap. It multiplies the normalized gene expression matrix (gene X cell) with the mapping similarity matrix (cell X bin) and sets a threshold to binarize final spatial gene expression. For sub-challenge 1, 2, 3, we select an optimal threshold. First, we generate a series of thresholds between 0 and 1. For each threshold, we compute the predicted spatial gene expression using computeVISH (), then compute the MCC score between the predicted gene expression and the binarized expression of the driver genes. Finally, we select the best threshold based on the largest MCC score to compute the final spatial gene expression. To estimate the performance of our method to predict spatial gene expression, we compute the Pearson correlation between the in situ gene expression and the predicted gene expression. Performance evaluation To evaluate the performance of our method, we use the three performance scores in the DREAM challenge ( https://www.synapse.org/#!Synapse:syn17091286 , https://github.com/dream-sctc/Scoring/blob/master/dream_scoring_clean.R ). The first score (score 1) is the primary score to evaluate the precision of the assignment for the single cells. The second score (score 2) is the average of the relative assignment precision over all the single cells which is used when the first scores are equal for two methods. The third score (score 3) is to compare predictions of gene patterns. Ambiguous cells: If the predicted top 1 and top 2 positions are the same in the DistMap results, the prediction position will be ambiguous, and the cell will not be computed in the score 1–3. In this challenge, the number of ambiguous cells derived from DistMap are 287. Results Select driver genes We calculated the sums of gene expression values in the in situ matrix for each bin of the 3,039 bins. Each bin has at least one gene expressed in the in situ matrix. It suggests that our selected driver genes cover all bins as at least one driver gene is expressed in each bin. As the number of expressed driver genes decreases when the number of selected driver genes decreases (from 60 to 40 to 20), the prediction accuracy of cell spatial positions also decreases ( Figure 3(a) ). We also compared the overlaps when varying the number of the driver genes in Figure 3(b) . Among the 40 driver genes of the sub-challenge 2, only one driver gene is not in the selected 60 driver genes of the sub-challenge 1. There are 14 overlapped driver genes between the sub-challenge 1, 2 and 3. It suggests our method is consistent on selecting different number of driver genes. Figure 3. ( a ) Compare number of expressed driver genes in each bin in the in situ matrix in sub-challenge 1, 2, and 3 under different selected driver genes scenarios. ( b ) The Venn diagram of the selected 60, 40, and 20 gene in the three sub-challenges 1, 2, 3. There are 14 common genes in three sub-challenges 1, 2 and 3 (labeled as gene60, gene40, and gene20). Evaluate the predicted position and spatial gene expression We used the score 1–3 to evaluate our method for different number of selected driver genes. Figure 4(a) shows the scores of our submitted results for the sub-challenge 1, 2, 3. The blue bar is the score for the silver standard method using 84 driver genes from DistMap. For score 1, the primary score for cell position prediction, the performance of our method is close to the silver standard when using 60 driver genes. The performance decreases when using 40 and 20 driver genes. For score 2, our method shows high scores when using 60 and 40 driver genes and lower score when using 20 driver genes. As Score 2 is the average relative precision for all cells, it suggests that our method is robust for predicting the right position using 60 and 40 driver genes, probably because our algorithm uses the top 10 bins around the largest cluster centers. The score 3 estimates the predicted gene expression patterns. It shows a small difference in our method when using 60 and 40 driver genes. To check the effect of threshold on the prediction results, we test our method under different thresholds. Figure 4(b)–(d) shows the consistency of the score 1–3 over a range of thresholds in the scenarios of different numbers of driver genes. Overall, the results show that our method performs well when using 60 and 40 genes. Figure 4. ( a ) Comparing the score 1, 2, 3 for sub-challenge 1 (60 driver genes), 2 (40 driver genes) and 3 (20 driver genes) with the silver standard (84 driver genes). ( b - d ) Comparing the score 1, 2, 3 using different thresholds (min, 1st quantile, median, mean, 3rd quantile, max value) for ( b ) sub-challenge 1 with 60 genes; ( c ) sub challenge 2 with 40 genes; ( d ) sub-challenge 3 with 20 genes. As shown in Figure 5 , we used boxplots to present the spatial gene expression prediction accuracy of the 84 driver genes by computing the Pearson correlation between the predicted and in situ gene expression. Consistent with Figure 4(a) , the average Pearson correlation of our prediction (60 or 40 driver genes) is higher than that of the silver standard (84 driver genes by DisMap). When the number of driver genes decreases to 20, the average Pearson correlation of our prediction becomes smaller but close to that of the silver standard, which used 84 driver genes. From the results, it suggests that our optimal threshold selection may improve the spatial gene expression prediction. It is feasible to use less genes for spatial gene expression prediction. Figure 5. Comparing Pearson correlation of predicted spatial gene expression and in situ spatial gene expression for 84 genes between DistMap (84 genes) with our method using ( a ) 60 genes in subchallenge 1; ( b ) 40 genes in sub-challenge 2; and ( c ) 20 genes in sub-challenge 3. For each gene, we calculate the Pearson correlation between the predicted spatial gene expression and the in situ spatial gene expression. The average Pearson correlation value of DistMap (84 genes), sub-challenge 1, sub-challenge 2, and sub-challenge 3 are 0.532, 0.622, 0.561, and 0.422, respectively. Discussion With high quality Drosophila dataset, Karaiskos et al. mapped single cells to the 3D embryo positions using the DistMap method 10 . It assumed that transcriptomes of neighboring cells should show similar patterns in spatial tissues. We further assumed that cells with similar gene expression will have close spatial distance. We applied hierarchical clustering to select a small number of driver genes and group similar cells to be mapped to the 3D embryo. After selecting the largest cluster, we determined the most possible predicted spatial cell position around the cluster center. We assessed the performance of our method on the selected 60, 40 or 20 driver genes using score 1–3 and compared it with the silver standard (DistMap results). For 60 driver genes, our results show close performance to the silver standard in score 1 and 3, and higher performance in score 2. It suggests that our method applied to 60 driver genes can achieve the same or higher level of performance as the silver standard to predict cell position. Our method obtains better performance for gene expression pattern prediction with 60 and 40 genes. When using 40 driver genes, the score 1 of our method decreases a little but the score 3 is very close to the silver standard. It suggests that our method can still predict gene expression pattern reasonably well using 40 genes. Fore 20 driver genes, the performance of our method is lower than but still close to the silver standard in both cell position and gene expression prediction. Overall, our method achieves good performance to predict spatial cell position and gene expression with 60 and 40 driver genes. The primary reason may be that our method selects the most representing genes in the gene clusters for prediction. Also, for potential spatial bins, our method utilizes bin clusters and selects the cluster center as predicted position. The good performance of our method is consistent to our assumption that cells with close spatial position show similar gene expression patterns. The performance of our method decreases when the number of driver genes reduces to 20. This is probably because the similarities between cells and spatial bins have larger variations when measured using smaller number of genes. In this case, it becomes harder to find the most possible bin position based on cell-bin similarities. In Nitzan’s work 14 , 15 , they successfully reconstructed the virtual Drosophila embryo with one marker gene and the structural information between gene expression space and physical space. It suggested that gene-gene correlation and spatial correlation contain useful information for gene expression and spatial cell position prediction. In addition, our method can be extended using an iterative approach by updating the mapping similarity matrix at each step. Data availability The dataset associated with the DREAM Single Cell Transcriptomics Challenge is available at https://www.synapse.org/#!Synapse:syn16782360 and https://www.synapse.org/#!Synapse:syn18632189 . Limited by the sharing settings of Synapse, users should register in Synapse (free of charge; https://www.synapse.org/ ) using their email address, and agree to the dataset conditions of use ( https://www.synapse.org/#!Synapse:syn15665609/wiki/583240 ). Software availability Source code implementation for the method presented in this article and used in the DREAM Single Cell Transcriptomics Challenge is available at: https://github.com/ouyang-lab/SCTC-Challenge-zho_team . The scoring scripts are available at: https://github.com/dream-sctc/Scoring/blob/master/dream_scoring_clean.R . Archived source code as at time of publication can be accessed at: https://doi.org/10.5281/zenodo.3592532 License: GLP 3.0 Author information All authors actively contributed to the results and discussed the procedures. They have been involved in the test, design and evaluation of various approaches; many of the explored techniques were not used in the final workflow due to strong overfitting. DM: gene selection and coding. YC: data pre-processing, clustering, and writeup. ZO: supervision, realization of workflows, and writeup. YZ: supervision, realization of workflows, and writeup. Acknowledgements Data used in this publication were generated by Prof. Dr. Nikolaus Rajewsky, Max Delbrück at the Center for Molecular Medicine, and these results were obtained as part of the DREAM Single Cell Transcriptomics Challenge project through Synapse ID (syn15665609). We thank Pablo Meyer and Jovan Tanevski for providing the documents, code and suggestions for score 1–3. Faculty Opinions recommended References 1. Farrell JA, Wang Y, Riesenfeld SJ, et al. : Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis. Science. 2018; 360 (6392): eaar3131. PubMed Abstract | Publisher Full Text | Free Full Text 2. Butler A, Hoffman P, Smibert P, et al. : Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018; 36 (5): 411–420. PubMed Abstract | Publisher Full Text | Free Full Text 3. Colomé-Tatché M, Theis FJ: Statistical single cell multi-omics integration. Curr Opin Syst Biol. 2018; 7 : 54–59. Publisher Full Text 4. Gu C, Liu S, Wu Q, et al. : Integrative single-cell analysis of transcriptome, DNA methylome and chromatin accessibility in mouse oocytes. Cell Res. 2019; 29 (2): 110–123. PubMed Abstract | Publisher Full Text | Free Full Text 5. Bian S, Hou Y, Zhou X, et al. : Single-cell multiomics sequencing and analyses of human colorectal cancer. Science. 2018; 362 (6418): 1060–1063. PubMed Abstract | Publisher Full Text 6. Kulkarni A, Anderson AG, Merullo Dp, et al. : Beyond bulk: a review of single cell transcriptomics methodologies and applications. Curr Opin Biotechnol. 2019; 58 : 129–136. PubMed Abstract | Publisher Full Text | Free Full Text 7. Savulescu AF, Jacobs C, Negishi Y, et al. : Pinpointing cell identity in time and space. Front Mol Biosci. 2020; 7 : 209. PubMed Abstract | Publisher Full Text | Free Full Text 8. Stuart T, Satija R: Integrative single-cell analysis. Nat Rev Genet. 2019; 20 (5): 257–272. PubMed Abstract | Publisher Full Text 9. Mayr U, Serra D, Liberali P: Exploring single cells in space and time during tissue development, homeostasis and regeneration. Development. 2019; 146 (12): dev176727. PubMed Abstract | Publisher Full Text 10. Karaiskos N, Wahle P, Alles J, et al. : The Drosophila embryo at single-cell transcriptome resolution. Science. 2017; 358 (6360): 194–199. PubMed Abstract | Publisher Full Text 11. Satija R, Farrell JA, Gennert D, et al. : Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015; 33 (5): 495–502. PubMed Abstract | Publisher Full Text | Free Full Text 12. Achim K, Pettit JB, Saraiva LR, et al. : High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin. Nat Biotechnol. 2015; 33 (5): 503–9. PubMed Abstract | Publisher Full Text 13. Moor AE, Harnik Y, Ben-Moshe S, et al. : Spatial reconstruction of single enterocytes uncovers broad zonation along the intestinal villus axis. Cell. 2018; 175 (4): 1156–1167.e15. PubMed Abstract | Publisher Full Text 14. Nitzan M, Karaiskos N, Friedman N, et al. : Charting a tissue from single-cell transcriptomes. bioRxiv. 2018; 456350. Publisher Full Text 15. Nitzan M, Karaiskos N, Friedman N, et al. : Gene expression cartography. Nature. 2019; 576 (7785): 132–137. PubMed Abstract | Publisher Full Text 16. Shah S, Lubeck E, Zhou W, et al. : In situ transcription profiling of single cells reveals spatial organization of cells in the mouse hippocampus. Neuron. 2016; 92 (2): 342–357. PubMed Abstract | Publisher Full Text | Free Full Text 17. Zhu Q, Shah S, Dries R, et al. : Identification of spatially associated subpopulations by combining scRNAseq and sequential fluorescence in situ hybridization data. Nat Biotechnol. 2018; 36 : 1183–1190. PubMed Abstract | Publisher Full Text | Free Full Text 18. Moffitt JR, Bambah-Mukku D, Eichhorn SW, et al. : Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region. Science. 2018; 362 (6416): eaau5324. PubMed Abstract | Publisher Full Text | Free Full Text 19. Halpern KB, Shenhav R, Matcovitch-Natan O, et al. : Single-cell spatial reconstruction reveals global division of labour in the mammalian liver. Nature. 2017; 542 (7641): 352–356. PubMed Abstract | Publisher Full Text | Free Full Text 20. Tanevski J, Nguyen T, Truong B, et al. : Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data. Life Sci Alliance. 2020; 3 (11): e202000867. PubMed Abstract | Publisher Full Text | Free Full Text 21. Macosko EZ, Basu A, Satija R, et al. : Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell. 2015; 161 (5): 1202–1214. PubMed Abstract | Publisher Full Text | Free Full Text 22. Charrad M, Ghazzali N, Boiteau V, et al. : NbClust: An R Package for Determining the Relevant Number of Clusters in a Data Set. J Stat Softw. 2014; 61 (6): 1–36. Publisher Full Text Comments on this article Comments (0) Version 2 VERSION 2 PUBLISHED 19 Feb 2020 ADD YOUR COMMENT Comment Author details Author details 1 Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA, USA 2 Department of Statistics, University of Connecticut, Storrs, CT, USA Yang Chen Roles: Data Curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – Original Draft Preparation Disheng Mao Roles: Investigation, Methodology Yuping Zhang Roles: Funding Acquisition, Investigation, Methodology, Resources, Supervision, Writing – Review & Editing Zhengqing Ouyang Roles: Conceptualization, Formal Analysis, Funding Acquisition, Investigation, Methodology, Project Administration, Resources, Supervision, Writing – Review & Editing Competing interests No competing interests were disclosed. Grant information This work was partially supported by the Faculty Research Excellence Program Award at UConn (to YZ). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Article Versions (2) version 2 Revised Published: 07 Jun 2021, 9:124 https://doi.org/10.12688/f1000research.20446.2 version 1 Published: 19 Feb 2020, 9:124 https://doi.org/10.12688/f1000research.20446.1 Copyright © 2021 Chen Y 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. 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 Chen Y, Mao D, Zhang Y and Ouyang Z. Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.12688/f1000research.20446.2 ) 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 2 VERSION 2 PUBLISHED 07 Jun 2021 Revised Views 0 Cite How to cite this report: Raja KKB. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290214 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290214 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 28 Jun 2024 Komal Kumar Bollepogu Raja , Baylor college of medicine, Houston, TX, USA Approved with Reservations VIEWS 0 https://doi.org/10.5256/f1000research.53961.r290214 In this manuscript Chen et al develop an unsupervised gene selection model that predicts cell positions. They have utilized published single cell RNA seq dataset from Drosophila embryo. Further, hierarchical clustering was used to select a set of driver genes ... Continue reading READ ALL In this manuscript Chen et al develop an unsupervised gene selection model that predicts cell positions. They have utilized published single cell RNA seq dataset from Drosophila embryo. Further, hierarchical clustering was used to select a set of driver genes to predict cell positions. The authors compare their model with Distmap and show that their model performs well with 60 genes. The manuscript is well written and the data is presented well. Spatial reconstruction of tissues using single cell and other omics approaches is an emerging field which allows identifying dynamic gene expression programs across tissues. Therefore, the current study can be a valuable resource. However, the authors have evaluated their model on only a single dataset. It will be worthwhile to evaluate this method on at least one more independent dataset to see its consistency. As of now, the method cannot be generalized to a wide range of tissues. Minor comment: 1. The last line in the Abstract should be changed: "The comparison with the âÂÂsilver standardâ suggests that our method is effective in reconstructing the cell spatial positions and gene expression patterns in tissues." Instead of tissues I suggest Drosophila embryo. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly Competing Interests: No competing interests were disclosed. Reviewer Expertise: Single cell omics, Genomics and 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, however I have significant reservations, as outlined above. Close READ LESS CITE CITE HOW TO CITE THIS REPORT Raja KKB. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290214 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290214 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: Choudhury SD. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290210 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290210 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 28 Jun 2024 Saumitra Dey Choudhury , All India Institute of Medical Sciences, New Delhi, India Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.53961.r290210 In this manuscript titled Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo, the authors have attempted to computationally deduce the spatial positions of individual cells in Drosophila embryos from their Drop-seq data in reference to the ... Continue reading READ ALL In this manuscript titled Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo, the authors have attempted to computationally deduce the spatial positions of individual cells in Drosophila embryos from their Drop-seq data in reference to the in situ hybridization expression patterns of 84 driver genes from the BDNTP. Though the method is well defined and could be useful in general to map spatial gene expression using transcriptomic data, some concerns exist: First, the Methods section does not have any explanation of the wet lab procedure of homogenization of embryos, making the single-cell suspension and cDNA library preparation. Second, the quality of the Drop-seq data cannot be established from the manuscript. Third, the word Unsupervised in the title is misleading. Fourth and most importantly, experimental validation of the performance of the method is required; at least for a few of the driver genes using antibodies/ immunofluorescence. In general, this manuscript brings forth a strong computational method to extrapolate spatial cell positions but lacks the experimental validation to support it. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? No Competing Interests: No competing interests were disclosed. Reviewer Expertise: Drosophila Neurogenetics, Single-cell RNA Transcriptomics, Confocal and Super-resolution Microscopy 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 Choudhury SD. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290210 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290210 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: Chen YC. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290211 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290211 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 27 Jun 2024 Yen-Chung Chen , New York University, New York, New York, USA Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.53961.r290211 Chen et al. described a novel method for feature selection for spatially mapping single-cell transcriptomic data. Originally a leading solution in DREAM Single-Cell Transcriptomics challenge [1], this method has demonstrated its relative edge in the challenge and superior performance in certain aspects against ... Continue reading READ ALL Chen et al. described a novel method for feature selection for spatially mapping single-cell transcriptomic data. Originally a leading solution in DREAM Single-Cell Transcriptomics challenge [1], this method has demonstrated its relative edge in the challenge and superior performance in certain aspects against distMap [2]. While being indeed an advance in methodology, the work left several aspect to be desired to generalize and address feature selection in spatial mapping. First, the original 84-gene pool was pre-selected based on prior knowledge of importance and spatial expression pattern in developing fly embryos [3] while most single-cell transcriptomic atlas are and will be generated from samples lacking prior knowledge with such breadth and resolution. It therefore remains open whether this proposed method will remain efficient when feature selection has to be performed on all genes or all variable genes. To address this issue, the method should be tested on more recent spatial and classic single-cell transcriptomic dataset to benchmark whether the method performs favorably when selecting features from single-cell RNA-seq without supervision from prior knowledge . Second, since assigning a cell to a bin can be considered as transferring the label of the bin to the cell of interest, potential users would be curious h ow the method holds up against the broad selection of existing label transfer methods , including Seurat [4], scVI [5], and classic classifiers like multimodal logistic regression, random forest classifier, and support vector machine for predicting which in situ bin a cell will match. It might be of interest that a previous iteration of Seurat used to demonstrate a similar spatial mapping task [6] and will be a great benchmarking dataset. Finally, there are minor methodological decisions that were reported where discussion of rationale could help a potential user and future method development. Specifically, in section Gene selection , the choice between Ward and Mcquitty linkage appeared arbitrary to me and made me wonder that for a user working on a different dataset, how the decision should be made (e.g., with regard to a fixed total number of clusters?). In section Compute the similarity matrix between cells and bins , it is unclear why p ij was chosen as the scoring scheme instead of correlation coefficient or mutual information. In sum, the work represents a novel and efficient method addressing the pressing need to systematically supplement existing single-cell transcriptomic atlases with spatial information, and the authors have demonstrated its power in answering a specific question. However, its general applicability remains to be tested on available spatial transcriptomic datasets collected in other contexts and tissues (e.g., [7-9]). Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Partly Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly References 1. Tanevski J, Nguyen T, Truong B, Karaiskos N, et al.: Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data. Life Sci Alliance . 2020; 3 (11). PubMed Abstract | Publisher Full Text 2. Karaiskos N, Wahle P, Alles J, Boltengagen A, et al.: The Drosophila embryo at single-cell transcriptome resolution. Science . 2017; 358 (6360): 194-199 PubMed Abstract | Publisher Full Text 3. Fowlkes CC, Hendriks CL, Keränen SV, Weber GH, et al.: A quantitative spatiotemporal atlas of gene expression in the Drosophila blastoderm. Cell . 2008; 133 (2): 364-74 PubMed Abstract | Publisher Full Text 4. Hao Y, Stuart T, Kowalski MH, Choudhary S, et al.: Dictionary learning for integrative, multimodal and scalable single-cell analysis. Nat Biotechnol . 2024; 42 (2): 293-304 PubMed Abstract | Publisher Full Text 5. Gayoso A, Lopez R, Xing G, Boyeau P, et al.: A Python library for probabilistic analysis of single-cell omics data. Nature Biotechnology . 2022; 40 (2): 163-166 Publisher Full Text 6. Mayer C, Hafemeister C, Bandler RC, Machold R, et al.: Developmental diversification of cortical inhibitory interneurons. Nature . 2018; 555 (7697): 457-462 PubMed Abstract | Publisher Full Text 7. Maniatis S, Äijö T, Vickovic S, Braine C, et al.: Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis. Science . 2019; 364 (6435): 89-93 PubMed Abstract | Publisher Full Text 8. Langlieb J, Sachdev NS, Balderrama KS, Nadaf NM, et al.: The molecular cytoarchitecture of the adult mouse brain. Nature . 2023; 624 (7991): 333-342 PubMed Abstract | Publisher Full Text 9. Wang M, Hu Q, Lv T, Wang Y, et al.: High-resolution 3D spatiotemporal transcriptomic maps of developing Drosophila embryos and larvae. Dev Cell . 2022; 57 (10): 1271-1283.e4 PubMed Abstract | Publisher Full Text Competing Interests: No competing interests were disclosed. Reviewer Expertise: Genomics, developmental biology, biostatistics 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 Chen YC. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290211 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-290211 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: Ren X. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r87000 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-87000 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 17 Jun 2021 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China Approved VIEWS 0 https://doi.org/10.5256/f1000research.53961.r87000 The authors have ... Continue reading READ ALL The authors have clarified my concerns. Competing Interests: No competing interests were disclosed. 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 Ren X. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r87000 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v2#referee-response-87000 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 19 Feb 2020 Views 0 Cite How to cite this report: Cembrowski MS and Kraus L. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r70705 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v1#referee-response-70705 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 23 Sep 2020 Mark S. Cembrowski , Dept. of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC, Canada Larissa Kraus , Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC, Canada Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.22477.r70705 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap ... Continue reading READ ALL In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Partly If any results are presented, are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions about the method and its performance adequately supported by the findings presented in the article? No Competing Interests: No competing interests were disclosed. Reviewer Expertise: Transcriptomics 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 Cembrowski MS and Kraus L. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r70705 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v1#referee-response-70705 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 Jun 2021 Zhengqing Ouyang , Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA 07 Jun 2021 Author Response Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors ... Continue reading Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. We thank the reviewer for the constructive comments. The “gold standard” is generated by DistMap, which is the best available single cell position reference used in the DREAM Challenge. DistMap determines single cell position by searching the maximum similarity for cells and spatial bins using the in situ expression of 84 driver genes in the Berkeley Drosophila Transcription Network Project (BDNTP). It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells (Karaiskos, et al. 2017). Thus, it is used in the DREAM Challenge to assess the performance of methods using smaller numbers (60, 40, and 20) of driver genes. In the revised manuscript, we changed “gold standard” to “silver standard”, as the cell positions generated by DistMap are not experimentally validated. We added the description of the DistMap “silver standard” in the Method section. We also added more introductions to the common methods for the cell position prediction from single cell transcriptome data in the Introduction section. We described the limitation of current methods and explain the motivation to study how the driver gene selection affecting the prediction of cell spatial position. The comparisons of our methods (Zho_team) and other top performing methods have been described in the DREAM Single Cell Transcriptomics Challenge consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867). Our manuscript is focused on describing the details of our method used in the DREAM Challenge. We also expanded the Discussion section of our manuscript. In addition to summarizing our method and results, we also pointed that gene-gene correlation and spatial correlation are useful information for spatial cell position and gene expression prediction. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. We thank the reviewer for the constructive comments. In the revised manuscript, we improved all figures and legends. We enlarged the font sizes for all figures. In Figure 1, we clarified the whole workflow to display the motivation of our method. In Figure 2, we changed it to display the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell when the number of driver genes is 60, 40, or 20. In Figure 3, we enlarged the font size and added more information in the figure labels. In Figure 4, we enlarged the font size. In figure 5, we changed it to the comparison of the Pearson correlation of predicted gene expression from DistMap (using 84 driver genes) and our method (using 60, 40, or 20 genes). Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. We thank the reviewer for the constructive comments. The “gold standard” is generated by DistMap, which is the best available single cell position reference used in the DREAM Challenge. DistMap determines single cell position by searching the maximum similarity for cells and spatial bins using the in situ expression of 84 driver genes in the Berkeley Drosophila Transcription Network Project (BDNTP). It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells (Karaiskos, et al. 2017). Thus, it is used in the DREAM Challenge to assess the performance of methods using smaller numbers (60, 40, and 20) of driver genes. In the revised manuscript, we changed “gold standard” to “silver standard”, as the cell positions generated by DistMap are not experimentally validated. We added the description of the DistMap “silver standard” in the Method section. We also added more introductions to the common methods for the cell position prediction from single cell transcriptome data in the Introduction section. We described the limitation of current methods and explain the motivation to study how the driver gene selection affecting the prediction of cell spatial position. The comparisons of our methods (Zho_team) and other top performing methods have been described in the DREAM Single Cell Transcriptomics Challenge consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867). Our manuscript is focused on describing the details of our method used in the DREAM Challenge. We also expanded the Discussion section of our manuscript. In addition to summarizing our method and results, we also pointed that gene-gene correlation and spatial correlation are useful information for spatial cell position and gene expression prediction. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. We thank the reviewer for the constructive comments. In the revised manuscript, we improved all figures and legends. We enlarged the font sizes for all figures. In Figure 1, we clarified the whole workflow to display the motivation of our method. In Figure 2, we changed it to display the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell when the number of driver genes is 60, 40, or 20. In Figure 3, we enlarged the font size and added more information in the figure labels. In Figure 4, we enlarged the font size. In figure 5, we changed it to the comparison of the Pearson correlation of predicted gene expression from DistMap (using 84 driver genes) and our method (using 60, 40, or 20 genes). Competing Interests: No competing interests were disclosed. Close Report a concern Respond or Comment COMMENTS ON THIS REPORT Author Response 07 Jun 2021 Zhengqing Ouyang , Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA 07 Jun 2021 Author Response Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors ... Continue reading Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. We thank the reviewer for the constructive comments. The “gold standard” is generated by DistMap, which is the best available single cell position reference used in the DREAM Challenge. DistMap determines single cell position by searching the maximum similarity for cells and spatial bins using the in situ expression of 84 driver genes in the Berkeley Drosophila Transcription Network Project (BDNTP). It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells (Karaiskos, et al. 2017). Thus, it is used in the DREAM Challenge to assess the performance of methods using smaller numbers (60, 40, and 20) of driver genes. In the revised manuscript, we changed “gold standard” to “silver standard”, as the cell positions generated by DistMap are not experimentally validated. We added the description of the DistMap “silver standard” in the Method section. We also added more introductions to the common methods for the cell position prediction from single cell transcriptome data in the Introduction section. We described the limitation of current methods and explain the motivation to study how the driver gene selection affecting the prediction of cell spatial position. The comparisons of our methods (Zho_team) and other top performing methods have been described in the DREAM Single Cell Transcriptomics Challenge consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867). Our manuscript is focused on describing the details of our method used in the DREAM Challenge. We also expanded the Discussion section of our manuscript. In addition to summarizing our method and results, we also pointed that gene-gene correlation and spatial correlation are useful information for spatial cell position and gene expression prediction. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. We thank the reviewer for the constructive comments. In the revised manuscript, we improved all figures and legends. We enlarged the font sizes for all figures. In Figure 1, we clarified the whole workflow to display the motivation of our method. In Figure 2, we changed it to display the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell when the number of driver genes is 60, 40, or 20. In Figure 3, we enlarged the font size and added more information in the figure labels. In Figure 4, we enlarged the font size. In figure 5, we changed it to the comparison of the Pearson correlation of predicted gene expression from DistMap (using 84 driver genes) and our method (using 60, 40, or 20 genes). Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. We thank the reviewer for the constructive comments. The “gold standard” is generated by DistMap, which is the best available single cell position reference used in the DREAM Challenge. DistMap determines single cell position by searching the maximum similarity for cells and spatial bins using the in situ expression of 84 driver genes in the Berkeley Drosophila Transcription Network Project (BDNTP). It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells (Karaiskos, et al. 2017). Thus, it is used in the DREAM Challenge to assess the performance of methods using smaller numbers (60, 40, and 20) of driver genes. In the revised manuscript, we changed “gold standard” to “silver standard”, as the cell positions generated by DistMap are not experimentally validated. We added the description of the DistMap “silver standard” in the Method section. We also added more introductions to the common methods for the cell position prediction from single cell transcriptome data in the Introduction section. We described the limitation of current methods and explain the motivation to study how the driver gene selection affecting the prediction of cell spatial position. The comparisons of our methods (Zho_team) and other top performing methods have been described in the DREAM Single Cell Transcriptomics Challenge consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867). Our manuscript is focused on describing the details of our method used in the DREAM Challenge. We also expanded the Discussion section of our manuscript. In addition to summarizing our method and results, we also pointed that gene-gene correlation and spatial correlation are useful information for spatial cell position and gene expression prediction. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. We thank the reviewer for the constructive comments. In the revised manuscript, we improved all figures and legends. We enlarged the font sizes for all figures. In Figure 1, we clarified the whole workflow to display the motivation of our method. In Figure 2, we changed it to display the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell when the number of driver genes is 60, 40, or 20. In Figure 3, we enlarged the font size and added more information in the figure labels. In Figure 4, we enlarged the font size. In figure 5, we changed it to the comparison of the Pearson correlation of predicted gene expression from DistMap (using 84 driver genes) and our method (using 60, 40, or 20 genes). Competing Interests: No competing interests were disclosed. Close Report a concern COMMENT ON THIS REPORT Views 0 Cite How to cite this report: Ren X. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r69520 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v1#referee-response-69520 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 02 Sep 2020 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China Not Approved VIEWS 0 https://doi.org/10.5256/f1000research.22477.r69520 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells ... Continue reading READ ALL Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly 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 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 Ren X. Reviewer Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r69520 ) The direct URL for this report is: https://f1000research.com/articles/9-124/v1#referee-response-69520 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 Jun 2021 Zhengqing Ouyang , Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA 07 Jun 2021 Author Response Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release ... Continue reading Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. 1. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". We thank the reviewer for the comments. To clarify, we used unsupervised hierarchical clustering for gene selection. Based on the selected genes, we predicted the spatial position of each individual cell. Unsupervised and supervised gene selection approaches are used by the methods in the DREAM Single Cell Transcriptomics Challenge for predicting cell spatial position. Table S2 of the consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867) listed the top performing methods (included ours) using either unsupervised or supervised gene selection. 2.Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. We thank the reviewer for the constructive comments. The results of our method (Zho_team) and other top performing methods have been systematically and independently evaluated on the Drosophila embryo dataset in the DREAM Single Cell Transcriptomics Challenge (https://www.synapse.org/#!Synapse:syn15665609/wiki/583234) and in the consortium paper aforementioned. Our manuscript, collected in the F1000Research DREAM Challenges Gateway (https://f1000research.com/gateways/dreamchallenges/singlecelltranscriptomics), is focused on describing the details of our method used in the DREAM Challenge. To avoid ambiguity, in the revised manuscript, we changed the title of our manuscript to “Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo”. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. We thank the reviewer for the comments. Our method is the only one that uses unsupervised hierarchical clustering for gene selection among the top performing methods as described in Table S2 of the DREAM Challenge consortium paper mentioned above. The power has been systematically demonstrated in our manuscript as well as the consortium paper using independent measurements and comprehensive comparisons. Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. 1. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". We thank the reviewer for the comments. To clarify, we used unsupervised hierarchical clustering for gene selection. Based on the selected genes, we predicted the spatial position of each individual cell. Unsupervised and supervised gene selection approaches are used by the methods in the DREAM Single Cell Transcriptomics Challenge for predicting cell spatial position. Table S2 of the consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867) listed the top performing methods (included ours) using either unsupervised or supervised gene selection. 2.Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. We thank the reviewer for the constructive comments. The results of our method (Zho_team) and other top performing methods have been systematically and independently evaluated on the Drosophila embryo dataset in the DREAM Single Cell Transcriptomics Challenge (https://www.synapse.org/#!Synapse:syn15665609/wiki/583234) and in the consortium paper aforementioned. Our manuscript, collected in the F1000Research DREAM Challenges Gateway (https://f1000research.com/gateways/dreamchallenges/singlecelltranscriptomics), is focused on describing the details of our method used in the DREAM Challenge. To avoid ambiguity, in the revised manuscript, we changed the title of our manuscript to “Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo”. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. We thank the reviewer for the comments. Our method is the only one that uses unsupervised hierarchical clustering for gene selection among the top performing methods as described in Table S2 of the DREAM Challenge consortium paper mentioned above. The power has been systematically demonstrated in our manuscript as well as the consortium paper using independent measurements and comprehensive comparisons. Competing Interests: No competing interests were disclosed. Close Report a concern Reviewer Response 17 Jun 2021 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China 17 Jun 2021 Reviewer Response The authors have clarified my concerns. Competing Interests: NA The authors have clarified my concerns. The authors have clarified my concerns. Competing Interests: NA Close Report a concern Respond or Comment COMMENTS ON THIS REPORT Author Response 07 Jun 2021 Zhengqing Ouyang , Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA 07 Jun 2021 Author Response Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release ... Continue reading Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. 1. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". We thank the reviewer for the comments. To clarify, we used unsupervised hierarchical clustering for gene selection. Based on the selected genes, we predicted the spatial position of each individual cell. Unsupervised and supervised gene selection approaches are used by the methods in the DREAM Single Cell Transcriptomics Challenge for predicting cell spatial position. Table S2 of the consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867) listed the top performing methods (included ours) using either unsupervised or supervised gene selection. 2.Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. We thank the reviewer for the constructive comments. The results of our method (Zho_team) and other top performing methods have been systematically and independently evaluated on the Drosophila embryo dataset in the DREAM Single Cell Transcriptomics Challenge (https://www.synapse.org/#!Synapse:syn15665609/wiki/583234) and in the consortium paper aforementioned. Our manuscript, collected in the F1000Research DREAM Challenges Gateway (https://f1000research.com/gateways/dreamchallenges/singlecelltranscriptomics), is focused on describing the details of our method used in the DREAM Challenge. To avoid ambiguity, in the revised manuscript, we changed the title of our manuscript to “Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo”. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. We thank the reviewer for the comments. Our method is the only one that uses unsupervised hierarchical clustering for gene selection among the top performing methods as described in Table S2 of the DREAM Challenge consortium paper mentioned above. The power has been systematically demonstrated in our manuscript as well as the consortium paper using independent measurements and comprehensive comparisons. Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. 1. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". We thank the reviewer for the comments. To clarify, we used unsupervised hierarchical clustering for gene selection. Based on the selected genes, we predicted the spatial position of each individual cell. Unsupervised and supervised gene selection approaches are used by the methods in the DREAM Single Cell Transcriptomics Challenge for predicting cell spatial position. Table S2 of the consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867) listed the top performing methods (included ours) using either unsupervised or supervised gene selection. 2.Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. We thank the reviewer for the constructive comments. The results of our method (Zho_team) and other top performing methods have been systematically and independently evaluated on the Drosophila embryo dataset in the DREAM Single Cell Transcriptomics Challenge (https://www.synapse.org/#!Synapse:syn15665609/wiki/583234) and in the consortium paper aforementioned. Our manuscript, collected in the F1000Research DREAM Challenges Gateway (https://f1000research.com/gateways/dreamchallenges/singlecelltranscriptomics), is focused on describing the details of our method used in the DREAM Challenge. To avoid ambiguity, in the revised manuscript, we changed the title of our manuscript to “Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo”. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. We thank the reviewer for the comments. Our method is the only one that uses unsupervised hierarchical clustering for gene selection among the top performing methods as described in Table S2 of the DREAM Challenge consortium paper mentioned above. The power has been systematically demonstrated in our manuscript as well as the consortium paper using independent measurements and comprehensive comparisons. Competing Interests: No competing interests were disclosed. Close Report a concern Reviewer Response 17 Jun 2021 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China 17 Jun 2021 Reviewer Response The authors have clarified my concerns. Competing Interests: NA The authors have clarified my concerns. The authors have clarified my concerns. Competing Interests: NA Close Report a concern COMMENT ON THIS REPORT Comments on this article Comments (0) Version 2 VERSION 2 PUBLISHED 19 Feb 2020 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 Version 2 (revision) 07 Jun 21 read read read read Version 1 19 Feb 20 read read Xianwen Ren , Peking University, Beijing, China Mark S. Cembrowski , University of British Columbia, Vancouver, Canada Larissa Kraus , University of British Columbia, Vancouver, Canada Yen-Chung Chen , New York University, New York, USA Saumitra Dey Choudhury , All India Institute of Medical Sciences, New Delhi, India Komal Kumar Bollepogu Raja , Baylor college of medicine, Houston, USA 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 © 2024 Raja K. 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. 28 Jun 2024 | for Version 2 Komal Kumar Bollepogu Raja , Baylor college of medicine, Houston, TX, USA 0 Views copyright © 2024 Raja K. 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 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 this manuscript Chen et al develop an unsupervised gene selection model that predicts cell positions. They have utilized published single cell RNA seq dataset from Drosophila embryo. Further, hierarchical clustering was used to select a set of driver genes to predict cell positions. The authors compare their model with Distmap and show that their model performs well with 60 genes. The manuscript is well written and the data is presented well. Spatial reconstruction of tissues using single cell and other omics approaches is an emerging field which allows identifying dynamic gene expression programs across tissues. Therefore, the current study can be a valuable resource. However, the authors have evaluated their model on only a single dataset. It will be worthwhile to evaluate this method on at least one more independent dataset to see its consistency. As of now, the method cannot be generalized to a wide range of tissues. Minor comment: 1. The last line in the Abstract should be changed: "The comparison with the âÂÂsilver standardâ suggests that our method is effective in reconstructing the cell spatial positions and gene expression patterns in tissues." Instead of tissues I suggest Drosophila embryo. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly Competing Interests No competing interests were disclosed. Reviewer Expertise Single cell omics, Genomics and 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, however I have significant reservations, as outlined above. reply Respond to this report Responses (0) Raja KKB. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290214) 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/9-124/v2#referee-response-290214 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2024 Choudhury 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. 28 Jun 2024 | for Version 2 Saumitra Dey Choudhury , All India Institute of Medical Sciences, New Delhi, India 0 Views copyright © 2024 Choudhury 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) 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 In this manuscript titled Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo, the authors have attempted to computationally deduce the spatial positions of individual cells in Drosophila embryos from their Drop-seq data in reference to the in situ hybridization expression patterns of 84 driver genes from the BDNTP. Though the method is well defined and could be useful in general to map spatial gene expression using transcriptomic data, some concerns exist: First, the Methods section does not have any explanation of the wet lab procedure of homogenization of embryos, making the single-cell suspension and cDNA library preparation. Second, the quality of the Drop-seq data cannot be established from the manuscript. Third, the word Unsupervised in the title is misleading. Fourth and most importantly, experimental validation of the performance of the method is required; at least for a few of the driver genes using antibodies/ immunofluorescence. In general, this manuscript brings forth a strong computational method to extrapolate spatial cell positions but lacks the experimental validation to support it. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? No Competing Interests No competing interests were disclosed. Reviewer Expertise Drosophila Neurogenetics, Single-cell RNA Transcriptomics, Confocal and Super-resolution Microscopy 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 (0) Choudhury SD. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290210) 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/9-124/v2#referee-response-290210 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2024 Chen Y. 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. 27 Jun 2024 | for Version 2 Yen-Chung Chen , New York University, New York, New York, USA 0 Views copyright © 2024 Chen Y. 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 Chen et al. described a novel method for feature selection for spatially mapping single-cell transcriptomic data. Originally a leading solution in DREAM Single-Cell Transcriptomics challenge [1], this method has demonstrated its relative edge in the challenge and superior performance in certain aspects against distMap [2]. While being indeed an advance in methodology, the work left several aspect to be desired to generalize and address feature selection in spatial mapping. First, the original 84-gene pool was pre-selected based on prior knowledge of importance and spatial expression pattern in developing fly embryos [3] while most single-cell transcriptomic atlas are and will be generated from samples lacking prior knowledge with such breadth and resolution. It therefore remains open whether this proposed method will remain efficient when feature selection has to be performed on all genes or all variable genes. To address this issue, the method should be tested on more recent spatial and classic single-cell transcriptomic dataset to benchmark whether the method performs favorably when selecting features from single-cell RNA-seq without supervision from prior knowledge . Second, since assigning a cell to a bin can be considered as transferring the label of the bin to the cell of interest, potential users would be curious h ow the method holds up against the broad selection of existing label transfer methods , including Seurat [4], scVI [5], and classic classifiers like multimodal logistic regression, random forest classifier, and support vector machine for predicting which in situ bin a cell will match. It might be of interest that a previous iteration of Seurat used to demonstrate a similar spatial mapping task [6] and will be a great benchmarking dataset. Finally, there are minor methodological decisions that were reported where discussion of rationale could help a potential user and future method development. Specifically, in section Gene selection , the choice between Ward and Mcquitty linkage appeared arbitrary to me and made me wonder that for a user working on a different dataset, how the decision should be made (e.g., with regard to a fixed total number of clusters?). In section Compute the similarity matrix between cells and bins , it is unclear why p ij was chosen as the scoring scheme instead of correlation coefficient or mutual information. In sum, the work represents a novel and efficient method addressing the pressing need to systematically supplement existing single-cell transcriptomic atlases with spatial information, and the authors have demonstrated its power in answering a specific question. However, its general applicability remains to be tested on available spatial transcriptomic datasets collected in other contexts and tissues (e.g., [7-9]). Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Partly Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly References 1. Tanevski J, Nguyen T, Truong B, Karaiskos N, et al.: Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data. Life Sci Alliance . 2020; 3 (11). PubMed Abstract | Publisher Full Text 2. Karaiskos N, Wahle P, Alles J, Boltengagen A, et al.: The Drosophila embryo at single-cell transcriptome resolution. Science . 2017; 358 (6360): 194-199 PubMed Abstract | Publisher Full Text 3. Fowlkes CC, Hendriks CL, Keränen SV, Weber GH, et al.: A quantitative spatiotemporal atlas of gene expression in the Drosophila blastoderm. Cell . 2008; 133 (2): 364-74 PubMed Abstract | Publisher Full Text 4. Hao Y, Stuart T, Kowalski MH, Choudhary S, et al.: Dictionary learning for integrative, multimodal and scalable single-cell analysis. Nat Biotechnol . 2024; 42 (2): 293-304 PubMed Abstract | Publisher Full Text 5. Gayoso A, Lopez R, Xing G, Boyeau P, et al.: A Python library for probabilistic analysis of single-cell omics data. Nature Biotechnology . 2022; 40 (2): 163-166 Publisher Full Text 6. Mayer C, Hafemeister C, Bandler RC, Machold R, et al.: Developmental diversification of cortical inhibitory interneurons. Nature . 2018; 555 (7697): 457-462 PubMed Abstract | Publisher Full Text 7. Maniatis S, Äijö T, Vickovic S, Braine C, et al.: Spatiotemporal dynamics of molecular pathology in amyotrophic lateral sclerosis. Science . 2019; 364 (6435): 89-93 PubMed Abstract | Publisher Full Text 8. Langlieb J, Sachdev NS, Balderrama KS, Nadaf NM, et al.: The molecular cytoarchitecture of the adult mouse brain. Nature . 2023; 624 (7991): 333-342 PubMed Abstract | Publisher Full Text 9. Wang M, Hu Q, Lv T, Wang Y, et al.: High-resolution 3D spatiotemporal transcriptomic maps of developing Drosophila embryos and larvae. Dev Cell . 2022; 57 (10): 1271-1283.e4 PubMed Abstract | Publisher Full Text Competing Interests No competing interests were disclosed. Reviewer Expertise Genomics, developmental biology, biostatistics 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 (0) Chen YC. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r290211) 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/9-124/v2#referee-response-290211 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2021 Ren X. 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. 17 Jun 2021 | for Version 2 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China 0 Views copyright © 2021 Ren X. 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 authors have clarified my concerns. Competing Interests No competing interests were disclosed. 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) Ren X. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.53961.r87000) 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/9-124/v2#referee-response-87000 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2020 Cembrowski M 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. 23 Sep 2020 | for Version 1 Mark S. Cembrowski , Dept. of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC, Canada Larissa Kraus , Department of Cellular and Physiological Sciences, University of British Columbia, Vancouver, BC, Canada 0 Views copyright © 2020 Cembrowski M 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 (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 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Partly If any results are presented, are all the source data underlying the results available to ensure full reproducibility? Yes Are the conclusions about the method and its performance adequately supported by the findings presented in the article? No Competing Interests No competing interests were disclosed. Reviewer Expertise Transcriptomics 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 (1) Author Response 07 Jun 2021 Zhengqing Ouyang, Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA Point-by-point response to the comments of Reviewer 2 Reviewer 2 In this manuscript, the authors describe an approach to infer spatial gene expression using scRNA-seq data. The authors aimed to identify spatial expression of cell clusters using 20, 40 or 60 marker genes, and compared to the “gold standard” Distmap method using 84 marker genes. The authors begin to explain the rationale of the manuscript by describing the concept of the DREAM challenge, which aims to introduce novel computational techniques to map spatial gene expression using transcriptomic data. Since the DREAM challenge aims to implement novel methods, it is imperative to introduce and carefully compare common methods with the novel proposed approach, which was underdeveloped by the authors. As an example, the “gold standard” method is mentioned several times, though it is not clear what exactly the “gold standard” is or how well-controlled it is. As another key example, a discussion of any kind is completely missing. We thank the reviewer for the constructive comments. The “gold standard” is generated by DistMap, which is the best available single cell position reference used in the DREAM Challenge. DistMap determines single cell position by searching the maximum similarity for cells and spatial bins using the in situ expression of 84 driver genes in the Berkeley Drosophila Transcription Network Project (BDNTP). It is shown that the 84 driver gens are sufficient to uniquely and individually label most of the 1,297 cells (Karaiskos, et al. 2017). Thus, it is used in the DREAM Challenge to assess the performance of methods using smaller numbers (60, 40, and 20) of driver genes. In the revised manuscript, we changed “gold standard” to “silver standard”, as the cell positions generated by DistMap are not experimentally validated. We added the description of the DistMap “silver standard” in the Method section. We also added more introductions to the common methods for the cell position prediction from single cell transcriptome data in the Introduction section. We described the limitation of current methods and explain the motivation to study how the driver gene selection affecting the prediction of cell spatial position. The comparisons of our methods (Zho_team) and other top performing methods have been described in the DREAM Single Cell Transcriptomics Challenge consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867). Our manuscript is focused on describing the details of our method used in the DREAM Challenge. We also expanded the Discussion section of our manuscript. In addition to summarizing our method and results, we also pointed that gene-gene correlation and spatial correlation are useful information for spatial cell position and gene expression prediction. The design of figures for a manuscript should be chosen to clearly emphasize the main points of the data. Here, the figures are difficult to interpret because key information and description is missing (e.g., legends are too brief, axis labels are missing, subplot headings are redundant, figure design is inconsistent). In addition, the design and presentation of the figures seem unmotivated and thereby difficult for the reader to see the importance of any particular data or easily draw conclusions from the results. In addition, consistent colour schemes and designs would help the reader to follow the manuscript and understand the results. We thank the reviewer for the constructive comments. In the revised manuscript, we improved all figures and legends. We enlarged the font sizes for all figures. In Figure 1, we clarified the whole workflow to display the motivation of our method. In Figure 2, we changed it to display the number of bins with high similarity values (> 0.8) in the similarity matrix for each cell when the number of driver genes is 60, 40, or 20. In Figure 3, we enlarged the font size and added more information in the figure labels. In Figure 4, we enlarged the font size. In figure 5, we changed it to the comparison of the Pearson correlation of predicted gene expression from DistMap (using 84 driver genes) and our method (using 60, 40, or 20 genes). View more View less Competing Interests No competing interests were disclosed. reply Respond Report a concern Cembrowski MS and Kraus L. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r70705) 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/9-124/v1#referee-response-70705 keyboard_arrow_left Back to all reports Reviewer Report 0 Views copyright © 2020 Ren X. 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. 02 Sep 2020 | for Version 1 Xianwen Ren , Biomedical Pioneering Innovation Center, Peking University, Beijing, China 0 Views copyright © 2020 Ren X. 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 (2) 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 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. Is the rationale for developing the new method (or application) clearly explained? Yes Is the description of the method technically sound? Yes Are sufficient details provided to allow replication of the method development and its use by others? Yes If any results are presented, are all the source data underlying the results available to ensure full reproducibility? No source data required Are the conclusions about the method and its performance adequately supported by the findings presented in the article? Partly 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 state that I do not consider it to be of an acceptable scientific standard, for reasons outlined above. reply Respond to this report Responses (2) Author Response 07 Jun 2021 Zhengqing Ouyang, Department of Biostatistics and Epidemiology, School of Public Health and Health Sciences, University of Massachusetts, Amherst, USA Point-by-point response to the comments of Reviewer 1 Reviewer 1 Reconstructing the spatial information of single cells from single-cell RNA-seq data is a pivotal question to further release the revolutionary power of the scRNA-seq technology. Here the authors propose a computational method to infer the spatial positions of single cells of Drosophila embryos based on scRNA-seq data and a reference spatial map based on in situ hybridization of 84 driver genes. While the method is finely tuned for this specific scientific question, some concerns exist. 1. First, the title is misleading. The authors claimed "an unsupervised learning method". This is not valid because this method uses the reference spatial map of Drosophila embryos. Unless the method only uses the scRNA-seq data, it cannot be claimed "an unsupervised learning method". We thank the reviewer for the comments. To clarify, we used unsupervised hierarchical clustering for gene selection. Based on the selected genes, we predicted the spatial position of each individual cell. Unsupervised and supervised gene selection approaches are used by the methods in the DREAM Single Cell Transcriptomics Challenge for predicting cell spatial position. Table S2 of the consortium paper (Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics data, Life Sci Alliance. 2020 Nov; 3(11): e202000867) listed the top performing methods (included ours) using either unsupervised or supervised gene selection. 2.Second, the current method has only been evaluated on one dataset. More independent evaluations should be added to demonstrate the generalizing capacity of this method. We thank the reviewer for the constructive comments. The results of our method (Zho_team) and other top performing methods have been systematically and independently evaluated on the Drosophila embryo dataset in the DREAM Single Cell Transcriptomics Challenge (https://www.synapse.org/#!Synapse:syn15665609/wiki/583234) and in the consortium paper aforementioned. Our manuscript, collected in the F1000Research DREAM Challenges Gateway (https://f1000research.com/gateways/dreamchallenges/singlecelltranscriptomics), is focused on describing the details of our method used in the DREAM Challenge. To avoid ambiguity, in the revised manuscript, we changed the title of our manuscript to “Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo”. Overall, this manuscript addresses an important scientific question. But the method needs more justifications and evidence to demonstrate its novelty and power. We thank the reviewer for the comments. Our method is the only one that uses unsupervised hierarchical clustering for gene selection among the top performing methods as described in Table S2 of the DREAM Challenge consortium paper mentioned above. The power has been systematically demonstrated in our manuscript as well as the consortium paper using independent measurements and comprehensive comparisons. View more View less Competing Interests No competing interests were disclosed. reply Respond Report a concern Reviewer Response 17 Jun 2021 Xianwen Ren, Biomedical Pioneering Innovation Center, Peking University, Beijing, China The authors have clarified my concerns. View more View less Competing Interests NA reply Respond Report a concern Ren X. Peer Review Report For: Unsupervised gene selection for predicting cell spatial positions in the Drosophila embryo [version 2; peer review: 1 approved, 1 approved with reservations, 3 not approved] . F1000Research 2021, 9 :124 ( https://doi.org/10.5256/f1000research.22477.r69520) 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/9-124/v1#referee-response-69520 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 = "Unsupervised gene selection for predicting...".replace("'", ''); var linkedInUrl = "http://www.linkedin.com/shareArticle?url=https://f1000research.com/articles/9-124/v2" + "&title=" + encodeURIComponent(lTitle) + "&summary=" + encodeURIComponent('Read the article by '); var deliciousUrl = "https://del.icio.us/post?url=https://f1000research.com/articles/9-124/v2&title=" + encodeURIComponent(lTitle); var redditUrl = "http://reddit.com/submit?url=https://f1000research.com/articles/9-124/v2" + "&title=" + encodeURIComponent(lTitle); linkedInUrl += encodeURIComponent('Chen Y 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/9-124/v2/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/9-124", templates : { twitter : "Unsupervised gene selection for predicting cell spatial positions.... Chen Y et al., published by " + "@F1000Research" + ", https://f1000research.com/articles/9-124/v2" } }; 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/20446/53961") new F1000.Clipboard(); new F1000.ThesaurusTermsDisplay("articles", "article", "53961"); $(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 = { "66315": 0, "66314": 0, "66317": 0, "66316": 0, "69519": 0, "66318": 0, "69518": 0, "69521": 0, "69520": 36, "69523": 0, "69522": 0, "62996": 0, "62997": 0, "86999": 0, "62998": 0, "62999": 0, "87000": 13, "290205": 0, "290206": 0, "290207": 0, "290212": 0, "290213": 0, "290214": 10, "290208": 0, "290209": 0, "290210": 16, "67495": 0, "290211": 14, "67494": 0, "67497": 0, "67496": 0, "67498": 0, "70701": 0, "70700": 0, "60269": 0, "70703": 0, "60270": 0, "70702": 0, "60271": 0, "60272": 0, "70705": 34, "70704": 0, "60273": 0, "64945": 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 = "f565fdf3-5534-4914-aa63-0035e91cfc75"; 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(",","/")); } }); });
Text is read by the "Ask this paper" AI Q&A widget below.
Extraction quality varies by source — PMC NXML preserves structure
cleanly, OA-HTML may include some navigation residue, and OA-PDF can
have broken hyphenation. The publisher copy
(via DOI)
is the canonical version.