@sjcrh/proteinpaint-server 2.63.0 → 2.63.3-3

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@sjcrh/proteinpaint-server",
-  "version": "2.63.0",
+  "version": "2.63.3-3",
   "type": "module",
   "description": "a genomics visualization tool for exploring a cohort's genotype and phenotype data",
   "main": "src/app.js",
@@ -47,7 +47,7 @@
     "@typescript-eslint/eslint-plugin": "^5.60.0",
     "babel-loader": "^8.2.2",
     "esbuild": "^0.19.12",
-    "glob": "^7.2.3",
+    "glob": "^10.4.1",
     "node-watch": "^0.7.1",
     "nodemon": "^2.0.19",
     "prettier": "^2.8.8",

package/routes/genesetEnrichment.js

@@ -0,0 +1,110 @@
+import fs from "fs";
+import { spawn } from "child_process";
+import { Readable } from "stream";
+import path from "path";
+import serverconfig from "../src/serverconfig.js";
+const api = {
+  endpoint: "genesetEnrichment",
+  methods: {
+    all: {
+      init,
+      request: {
+        typeId: "genesetEnrichmentRequest"
+      },
+      response: {
+        typeId: "genesetEnrichmentResponse"
+        // will combine this with type checker
+        //valid: (t) => {}
+      }
+    }
+  }
+};
+function init({ genomes }) {
+  return async (req, res) => {
+    try {
+      const results = await run_genesetEnrichment_analysis(req.query, genomes);
+      res.send(results);
+    } catch (e) {
+      res.send({ status: "error", error: e.message || e });
+    }
+  };
+}
+async function run_genesetEnrichment_analysis(q, genomes) {
+  if (!genomes[q.genome].termdbs)
+    throw "termdb database is not available for " + q.genome;
+  const genesetenrichment_input = {
+    genes: q.genes,
+    fold_change: q.fold_change,
+    db: genomes[q.genome].termdbs.msigdb.cohort.db.connection.name,
+    gene_set_group: q.geneSetGroup
+  };
+  const gsea_output = await run_gsea(
+    path.join(serverconfig.binpath, "../python/src", "gsea.py"),
+    "/" + JSON.stringify(genesetenrichment_input)
+    // "/" is needed for python to accept the bracket "{" as a bracket
+  );
+  let result;
+  for (const line of gsea_output.split("\n")) {
+    if (line.startsWith("result: ")) {
+      result = JSON.parse(line.replace("result: ", ""));
+    } else {
+      console.log(line);
+    }
+  }
+  return result;
+}
+async function run_gsea(path2, data) {
+  try {
+    await fs.promises.stat(path2);
+  } catch (e) {
+    throw `${path2} does not exist`;
+  }
+  return new Promise((resolve, reject) => {
+    const _stdout = [];
+    const _stderr = [];
+    let python_path = "python3";
+    if (serverconfig.py_path)
+      python_path = serverconfig.py_path;
+    const sp = spawn(python_path, [path2]);
+    if (data) {
+      try {
+        const input = data.endsWith("\n") ? data : data + "\n";
+        Readable.from(input).pipe(sp.stdin);
+      } catch (e) {
+        sp.kill();
+        let errmsg = e;
+        const stderr = _stderr.join("").trim();
+        if (stderr)
+          errmsg += `
+python stderr: ${stderr}`;
+        reject(errmsg);
+      }
+    }
+    sp.stdout.on("data", (data2) => _stdout.push(data2));
+    sp.stderr.on("data", (data2) => _stderr.push(data2));
+    sp.on("error", (err) => reject(err));
+    sp.on("close", (code) => {
+      const stdout = _stdout.join("").trim();
+      const stderr = _stderr.join("").trim();
+      if (code !== 0) {
+        let errmsg = `python process exited with non-zero status code=${code}`;
+        if (stdout)
+          errmsg += `
+python stdout: ${stdout}`;
+        if (stderr)
+          errmsg += `
+python stderr: ${stderr}`;
+        reject(errmsg);
+      }
+      if (stderr) {
+        const errmsg = `python process emitted standard error
+python stderr: ${stderr}`;
+        reject(errmsg);
+      }
+      resolve(stdout);
+    });
+  });
+}
+export {
+  api
+};