@sjcrh/proteinpaint-server 2.145.1-0 → 2.145.1

package/dataset/termdb.test.js

@@ -303,7 +303,6 @@ function termdb_test_default() {
       geneExpression: {
         src: "native",
         file: "files/hg38/TermdbTest/rnaseq/TermdbTest.fpkm.matrix.new.h5",
-        newformat: true,
         unit: "FPKM"
       },
       ssGSEA: {

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@sjcrh/proteinpaint-server",
-  "version": "2.145.1-0",
+  "version": "2.145.1",
   "type": "module",
   "description": "a genomics visualization tool for exploring a cohort's genotype and phenotype data",
   "main": "src/app.js",
@@ -63,9 +63,9 @@
     "@sjcrh/augen": "2.143.0",
     "@sjcrh/proteinpaint-python": "2.144.0",
     "@sjcrh/proteinpaint-r": "2.145.1-0",
-    "@sjcrh/proteinpaint-rust": "2.145.1-0",
+    "@sjcrh/proteinpaint-rust": "2.145.1",
     "@sjcrh/proteinpaint-shared": "2.145.1-0",
-    "@sjcrh/proteinpaint-types": "2.145.1-0",
+    "@sjcrh/proteinpaint-types": "2.145.1-1",
     "@types/express": "^5.0.0",
     "@types/express-session": "^1.18.1",
     "better-sqlite3": "^9.4.1",

package/routes/termdb.cluster.js

@@ -213,22 +213,6 @@ async function validate_query_geneExpression(ds, genome) {
   }
   throw "unknown queries.geneExpression.src";
 }
-async function queryGeneExpression(hdf5_file, geneNames) {
-  const jsonInput = JSON.stringify({
-    hdf5_file,
-    gene: geneNames
-  });
-  try {
-    const result = await run_rust("readHDF5", jsonInput);
-    if (!result || result.length === 0) {
-      throw new Error("Failed to retrieve expression data: Empty or missing response");
-    }
-    return result;
-  } catch (error) {
-    console.error(`Error querying gene expression for ${geneNames}`);
-    throw error;
-  }
-}
 async function queryHDF5(hdf5_file, query) {
   const jsonInput = JSON.stringify({
     hdf5_file,
@@ -250,34 +234,16 @@ async function validateNative(q, ds) {
   q.samples = [];
   try {
     await utils.file_is_readable(q.file);
-    let tmp;
-    if (q.newformat) {
-      tmp = await run_rust("readH5", JSON.stringify({ hdf5_file: q.file, validate: true }));
-    } else {
-      tmp = await run_rust("validateHDF5", JSON.stringify({ hdf5_file: q.file }));
-    }
+    const tmp = await run_rust("readH5", JSON.stringify({ hdf5_file: q.file, validate: true }));
     const vr = JSON.parse(tmp);
     if (vr.status !== "success") throw vr.message;
-    if (q.newformat) {
-      if (!vr.samples?.length) throw "HDF5 file has no samples, please check file.";
-      for (const sn of vr.samples) {
-        const si = ds.cohort.termdb.q.sampleName2id(sn);
-        if (si == void 0) throw `unknown sample ${sn} from HDF5 ${q.file}`;
-        q.samples.push(si);
-      }
-      console.log(`${ds.label}: geneExpression HDF5 file validated. Format: ${vr.format}, Samples:`, vr.samples.length);
-    } else {
-      if (!vr.sampleNames?.length) throw "HDF5 file has no samples, please check file.";
-      for (const sn of vr.sampleNames) {
-        const si = ds.cohort.termdb.q.sampleName2id(sn);
-        if (si == void 0) throw `unknown sample ${sn} from HDF5 ${q.file}`;
-        q.samples.push(si);
-      }
-      console.log(
-        `${ds.label}: geneExpression HDF5 file validated. Format: ${vr.format}, Samples:`,
-        vr.sampleNames.length
-      );
+    if (!vr.samples?.length) throw "HDF5 file has no samples, please check file.";
+    for (const sn of vr.samples) {
+      const si = ds.cohort.termdb.q.sampleName2id(sn);
+      if (si == void 0) throw `unknown sample ${sn} from HDF5 ${q.file}`;
+      q.samples.push(si);
     }
+    console.log(`${ds.label}: geneExpression HDF5 file validated. Format: ${vr.format}, Samples:`, vr.samples.length);
   } catch (error) {
     throw `${ds.label}: Failed to validate geneExpression HDF5 file: ${error}`;
   }
@@ -310,63 +276,32 @@ async function validateNative(q, ds) {
       return { term2sample2value, byTermId };
     }
     const time1 = Date.now();
-    let geneData;
-    if (q.newformat) {
-      geneData = JSON.parse(await queryHDF5(q.file, geneNames));
-    } else {
-      geneData = JSON.parse(await queryGeneExpression(q.file, geneNames));
-    }
+    const geneData = JSON.parse(await queryHDF5(q.file, geneNames));
     mayLog("Time taken to run gene query:", formatElapsedTime(Date.now() - time1));
-    if (q.newformat) {
-      const genesData = geneData.query_output || {};
-      if (!genesData) throw "No expression data returned from HDF5 query";
-      for (const tw of param.terms) {
-        if (!tw.term.gene) continue;
-        const geneResult = genesData[tw.term.gene];
-        if (!geneResult) {
-          console.warn(`No data found for gene ${tw.term.gene} in the response`);
-          continue;
-        }
-        const samplesData = geneResult.samples || {};
-        const s2v = {};
-        for (const sampleName in samplesData) {
-          const sampleId = ds.cohort.termdb.q.sampleName2id(sampleName);
-          if (!sampleId) continue;
-          if (limitSamples && !limitSamples.has(sampleId)) continue;
-          s2v[sampleId] = samplesData[sampleName];
-        }
-        if (Object.keys(s2v).length) {
-          term2sample2value.set(tw.$id, s2v);
-        }
+    const genesData = geneData.query_output || {};
+    if (!genesData) throw "No expression data returned from HDF5 query";
+    for (const tw of param.terms) {
+      if (!tw.term.gene) continue;
+      const geneResult = genesData[tw.term.gene];
+      if (!geneResult) {
+        console.warn(`No data found for gene ${tw.term.gene} in the response`);
+        continue;
       }
-      if (term2sample2value.size == 0) {
-        throw "No data available for the input " + param.terms?.map((tw) => tw.term.gene).join(", ");
+      const samplesData = geneResult.samples || {};
+      const s2v = {};
+      for (const sampleName in samplesData) {
+        const sampleId = ds.cohort.termdb.q.sampleName2id(sampleName);
+        if (!sampleId) continue;
+        if (limitSamples && !limitSamples.has(sampleId)) continue;
+        s2v[sampleId] = samplesData[sampleName];
       }
-    } else {
-      const genesData = geneData.genes || { [geneNames[0]]: geneData };
-      for (const tw of param.terms) {
-        if (!tw.term.gene) continue;
-        const geneResult = genesData[tw.term.gene];
-        if (!geneResult) {
-          console.warn(`No data found for gene ${tw.term.gene} in the response`);
-          continue;
-        }
-        const samplesData = geneResult.samples || {};
-        const s2v = {};
-        for (const [sampleName, value] of Object.entries(samplesData)) {
-          const sampleId = ds.cohort.termdb.q.sampleName2id(sampleName);
-          if (!sampleId) continue;
-          if (limitSamples && !limitSamples.has(sampleId)) continue;
-          s2v[sampleId] = value;
-        }
-        if (Object.keys(s2v).length) {
-          term2sample2value.set(tw.$id, s2v);
-        }
-      }
-      if (term2sample2value.size == 0) {
-        throw "No data available for the input " + param.terms?.map((tw) => tw.term.gene).join(", ");
+      if (Object.keys(s2v).length) {
+        term2sample2value.set(tw.$id, s2v);
       }
     }
+    if (term2sample2value.size == 0) {
+      throw "No data available for the input " + param.terms?.map((tw) => tw.term.gene).join(", ");
+    }
     return { term2sample2value, byTermId, bySampleId };
   };
 }