@platforma-open/milaboratories.top-antibodies.workflow 4.0.1 → 4.0.2

package/.turbo/turbo-build.log

@@ -1,6 +1,6 @@
  WARN  Issue while reading "/home/runner/work/antibody-tcr-lead-selection/antibody-tcr-lead-selection/.npmrc". Failed to replace env in config: ${NPMJS_TOKEN}
 
-> @platforma-open/milaboratories.top-antibodies.workflow@4.0.1 build /home/runner/work/antibody-tcr-lead-selection/antibody-tcr-lead-selection/workflow
+> @platforma-open/milaboratories.top-antibodies.workflow@4.0.2 build /home/runner/work/antibody-tcr-lead-selection/antibody-tcr-lead-selection/workflow
 > shx rm -rf dist && pl-tengo check && pl-tengo build
 
 Processing "src/assembling-fasta.tpl.tengo"...

package/CHANGELOG.md

@@ -1,5 +1,13 @@
 # @platforma-open/milaboratories.top-antibodies.workflow
 
+## 4.0.2
+
+### Patch Changes
+
+- 2a2533d: Fix minor issues
+- 6042e4a: Minor fix
+- 461999c: Fix minor issues
+
 ## 4.0.1
 
 ### Patch Changes

package/dist/tengo/lib/utils.lib.tengo

@@ -1,5 +1,6 @@
 
 
+ll := import("@platforma-sdk/workflow-tengo:ll")
 slices := import("@platforma-sdk/workflow-tengo:slices")
 json := import("json")
 
@@ -466,11 +467,18 @@ initializeCloneTable := func(pframes, columns, args, datasetSpec) {
     }
 
 
+
     if !addedCols {
         cdr3Sequences := columns.getColumns("cdr3Sequences")
         if len(cdr3Sequences) > 0 {
-            cloneTable.add(cdr3Sequences[0], {header: "cdr3_fallback"})
+            cloneTable.add(cdr3Sequences[0], {header: "sequence_fallback"})
             addedCols = true
+        } else {
+            peptideMainSeqs := columns.getColumns("peptideMainSeqs")
+            if len(peptideMainSeqs) > 0 {
+                cloneTable.add(peptideMainSeqs[0], {header: "sequence_fallback"})
+                addedCols = true
+            }
         }
     }
 
@@ -601,11 +609,28 @@ detectBulkChain := func(seqCols) {
 
 
 
-initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell) {
+
+initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell, isScFv) {
     assemSeqTable := pframes.parquetFileBuilder()
     assemSeqTable.setAxisHeader(datasetSpec.axesSpec[1].name, "clonotypeKey")
 
-    seqCols := columns.getColumns("assemblingAaSeqs")
+    seqCols := undefined
+    if isScFv {
+        allowedFeatures := { "VDJRegion": true, "VDJRegionInFrame": true }
+        byChain := {}
+        for col in columns.getColumns("scFvPerChainSeqs") {
+            if is_undefined(allowedFeatures[col.spec.domain["pl7.app/vdj/feature"]]) { continue }
+            if col.spec.domain["pl7.app/vdj/scClonotypeChain/index"] != "primary" { continue }
+            sc := col.spec.domain["pl7.app/vdj/scClonotypeChain"]
+            if sc != "A" && sc != "B" { continue }
+            if is_undefined(byChain[sc]) { byChain[sc] = col }
+        }
+        seqCols = []
+        for _, c in byChain { seqCols = append(seqCols, c) }
+    } else {
+        seqCols = columns.getColumns("assemblingAaSeqs")
+    }
+
     for col in seqCols {
         headerName := makeHeaderName(col, "assemblingFeature", isSingleCell)
         assemSeqTable.add(col, {header: headerName})
@@ -613,7 +638,7 @@ initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell) {
 
     assemSeqTable.mem("16GiB")
     assemSeqTable.cpu(1)
-    
+
 
     bulkChain := undefined
     if !isSingleCell {

package/package.json

@@ -1,16 +1,16 @@
 {
   "name": "@platforma-open/milaboratories.top-antibodies.workflow",
-  "version": "4.0.1",
+  "version": "4.0.2",
   "type": "module",
   "description": "Block Workflow",
   "dependencies": {
     "@platforma-sdk/workflow-tengo": "5.20.0",
     "@platforma-open/milaboratories.software-anarci": "^0.0.3",
     "@platforma-open/milaboratories.top-antibodies.sample-clonotypes": "2.1.2",
-    "@platforma-open/milaboratories.top-antibodies.assembling-fasta": "1.3.2",
+    "@platforma-open/milaboratories.top-antibodies.spectratype": "1.8.3",
     "@platforma-open/milaboratories.top-antibodies.umap": "1.2.3",
     "@platforma-open/milaboratories.top-antibodies.anarci-kabat": "1.4.3",
-    "@platforma-open/milaboratories.top-antibodies.spectratype": "1.8.3"
+    "@platforma-open/milaboratories.top-antibodies.assembling-fasta": "1.3.2"
   },
   "devDependencies": {
     "@platforma-sdk/tengo-builder": "2.5.20"

package/src/main.tpl.tengo

@@ -111,12 +111,37 @@ wf.prepare(func(args){
 		}
 	}, "JGenes")
 
-	// Add assembling feature aminoacid sequences (bulk, sc, scFv)
+	// Add assembling feature aminoacid sequences (bulk, sc).
 	bundleBuilder.addMulti({
 		axes: [{ anchor: "main", idx: 1 }], // Clonotype axis
+		name: "pl7.app/vdj/sequence",
 		annotations: { "pl7.app/vdj/isAssemblingFeature": "true" },
 		domain: { "pl7.app/alphabet": "aminoacid" }
 	}, "assemblingAaSeqs")
+
+	// Detect scFv mode by presence of the combined construct column.
+	bundleBuilder.addMulti({
+		axes: [{ anchor: "main", idx: 1 }],
+		name: "pl7.app/vdj/scFv-sequence"
+	}, "scFvDetect")
+
+	// Per-chain VDJRegion AA sequences for ANARCI/Kabat in scFv mode.
+	bundleBuilder.addMulti({
+		axes: [{ anchor: "main", idx: 1 }],
+		name: "pl7.app/vdj/sequence",
+		domain: { "pl7.app/alphabet": "aminoacid" }
+	}, "scFvPerChainSeqs")
+
+	// Peptide main sequence — single-axis column on variantKey, used as
+	// modality-aware fallback when no filter/ranking columns load.
+	bundleBuilder.addMulti({
+		axes: [{ anchor: "main", idx: 1 }],
+		annotations: {
+			"pl7.app/isAssemblingFeature": "true",
+			"pl7.app/isMainSequence": "true"
+		},
+		domain: { "pl7.app/alphabet": "aminoacid" }
+	}, "peptideMainSeqs")
     
     return {
         columns: bundleBuilder.build()
@@ -226,7 +251,8 @@ wf.body(func(args) {
                 if args.kabatNumbering == true {
                     ////////// Assembling AA sequences //////////
                     // Initialize and build assembling sequence table
-                    assemInit := utils.initializeAssemSeqTable(pframes, columns, datasetSpec, isSingleCell)
+                    isScFv := len(columns.getColumns("scFvDetect")) > 0
+                    assemInit := utils.initializeAssemSeqTable(pframes, columns, datasetSpec, isSingleCell, isScFv)
                     assemSeqTableBuilt := assemInit.assemSeqTable
                     bulkChain := assemInit.bulkChain
                     seqCols := assemInit.seqCols

package/src/utils.lib.tengo

@@ -1,5 +1,6 @@
 // Utility functions for antibody-tcr-lead-selection workflow
 
+ll := import("@platforma-sdk/workflow-tengo:ll")
 slices := import("@platforma-sdk/workflow-tengo:slices")
 json := import("json")
 
@@ -465,12 +466,19 @@ initializeCloneTable := func(pframes, columns, args, datasetSpec) {
         }
     }
 
-    // Fallback: if no columns added, add at least one CDR3 sequence column
+    // Fallback: if no columns added, add a single-axis trunk-keyed sequence
+    // column. Try VDJ CDR3 first
     if !addedCols {
         cdr3Sequences := columns.getColumns("cdr3Sequences")
         if len(cdr3Sequences) > 0 {
-            cloneTable.add(cdr3Sequences[0], {header: "cdr3_fallback"})
+            cloneTable.add(cdr3Sequences[0], {header: "sequence_fallback"})
             addedCols = true
+        } else {
+            peptideMainSeqs := columns.getColumns("peptideMainSeqs")
+            if len(peptideMainSeqs) > 0 {
+                cloneTable.add(peptideMainSeqs[0], {header: "sequence_fallback"})
+                addedCols = true
+            }
         }
     }
 
@@ -594,18 +602,35 @@ detectBulkChain := func(seqCols) {
 
 /**
  * Initializes and builds assembling sequence table with assembling AA sequences.
- * 
+
  * @param pframes - PFrames import
  * @param columns - PBundle containing all columns
  * @param datasetSpec - Dataset specification with axes
  * @param isSingleCell - Whether the data is single cell
+ * @param isScFv - Whether the input is scFv data
  * @return Map with keys: assemSeqTable (built table), bulkChain, seqCols
  */
-initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell) {
+initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell, isScFv) {
     assemSeqTable := pframes.parquetFileBuilder()
     assemSeqTable.setAxisHeader(datasetSpec.axesSpec[1].name, "clonotypeKey")
 
-    seqCols := columns.getColumns("assemblingAaSeqs")
+    seqCols := undefined
+    if isScFv {
+        allowedFeatures := { "VDJRegion": true, "VDJRegionInFrame": true }
+        byChain := {}
+        for col in columns.getColumns("scFvPerChainSeqs") {
+            if is_undefined(allowedFeatures[col.spec.domain["pl7.app/vdj/feature"]]) { continue }
+            if col.spec.domain["pl7.app/vdj/scClonotypeChain/index"] != "primary" { continue }
+            sc := col.spec.domain["pl7.app/vdj/scClonotypeChain"]
+            if sc != "A" && sc != "B" { continue }
+            if is_undefined(byChain[sc]) { byChain[sc] = col }
+        }
+        seqCols = []
+        for _, c in byChain { seqCols = append(seqCols, c) }
+    } else {
+        seqCols = columns.getColumns("assemblingAaSeqs")
+    }
+
     for col in seqCols {
         headerName := makeHeaderName(col, "assemblingFeature", isSingleCell)
         assemSeqTable.add(col, {header: headerName})
@@ -613,7 +638,7 @@ initializeAssemSeqTable := func(pframes, columns, datasetSpec, isSingleCell) {
 
     assemSeqTable.mem("16GiB")
     assemSeqTable.cpu(1)
-    
+
     // Detect bulk chain if needed
     bulkChain := undefined
     if !isSingleCell {