Merge pull request #38 from platforma-open/julenmendieta/MILAB-6617_adaptToSingleCellVhh

Julenmendieta/MILAB-6617_adaptToSingleCellVhh

Author: julenmendieta

.changeset/vhh-presence-driven-chain-menu.md

@@ -0,0 +1,12 @@
+---
+"@platforma-open/milaboratories.vj-usage.model": patch
+"@platforma-open/milaboratories.vj-usage.ui": patch
+"@platforma-open/milaboratories.vj-usage": patch
+---
+
+Only offer single-cell chains that actually have columns. The chain selector was
+built from the receptor type alone, so on heavy-chain-only (VHH) single-cell input
+it still offered "Light" — picking it failed the workflow with `expected exactly 1
+V gene column for chain, got 0`. A new `availableScChains` model output reports the
+chains present for the selected dataset; the selector now lists only those (and is
+hidden when a single chain leaves nothing to choose).

model/src/index.ts

@@ -3,8 +3,8 @@ import { BlockModelV3, createPFrameForGraphs } from "@platforma-sdk/model";
 import { blockDataModel } from "./dataModel";
 import type { BlockArgs } from "./types";
 
-export { getDefaultBlockLabel } from "./label";
 export { blockDataModel } from "./dataModel";
+export { getDefaultBlockLabel } from "./label";
 export type { BlockArgs, BlockData } from "./types";
 
 export const platforma = BlockModelV3.create(blockDataModel)
@@ -45,6 +45,43 @@ export const platforma = BlockModelV3.create(blockDataModel)
     return ctx.resultPool.getPColumnSpecByRef(ctx.data.datasetRef);
   })
 
+  // Single-cell IG chain letters ("A" = heavy, "B" = light) that actually have columns
+  // for the selected dataset. Returns undefined — meaning "don't filter" — for bulk
+  // data, non-IG receptors, and while the pool is resolving.
+  .output("availableScChains", (ctx) => {
+    const ref = ctx.data.datasetRef;
+    if (ref === undefined) return undefined;
+
+    const spec = ctx.resultPool.getPColumnSpecByRef(ref);
+    // Only single-cell IG can be single-chain (heavy-only VHH): bulk has no chain axis,
+    // and single-cell TCR is always paired — nothing to filter in those cases.
+    if (spec?.axesSpec[1]?.name !== "pl7.app/vdj/scClonotypeKey") return undefined;
+    if (spec.axesSpec[1]?.domain?.["pl7.app/vdj/receptor"] !== "IG") return undefined;
+
+    // Ask for the V-gene-hit column per chain
+    const vGeneCols = ctx.resultPool.getAnchoredPColumns(
+      { main: ref },
+      [
+        {
+          axes: [{ anchor: "main", idx: 1 }],
+          name: "pl7.app/vdj/geneHit",
+          domain: { "pl7.app/vdj/reference": "VGene" },
+        },
+      ],
+      { ignoreMissingDomains: true },
+    );
+    if (vGeneCols === undefined) return undefined; // pool still resolving
+
+    const chains = new Set<string>();
+    for (const col of vGeneCols) {
+      const domain = col.spec.domain;
+      if (domain?.["pl7.app/vdj/scClonotypeChain/index"] !== "primary") continue;
+      const letter = domain?.["pl7.app/vdj/scClonotypeChain"];
+      if (letter) chains.add(letter);
+    }
+    return [...chains].sort();
+  })
+
   .outputWithStatus("pf", (ctx) => {
     const pCols = ctx.outputs?.resolve("pf")?.getPColumns();
     if (pCols === undefined) {