@platforma-open/milaboratories.top-antibodies.workflow 1.2.1 → 1.3.0

package/.turbo/turbo-build.log

@@ -1,6 +1,6 @@
  WARN  Issue while reading "/home/runner/work/top-antibodies/top-antibodies/.npmrc". Failed to replace env in config: ${NPMJS_TOKEN}
 
-> @platforma-open/milaboratories.top-antibodies.workflow@1.2.1 build /home/runner/work/top-antibodies/top-antibodies/workflow
+> @platforma-open/milaboratories.top-antibodies.workflow@1.3.0 build /home/runner/work/top-antibodies/top-antibodies/workflow
 > rm -rf dist && pl-tengo check && pl-tengo build
 
 Processing "src/main.tpl.tengo"...

package/CHANGELOG.md

@@ -1,5 +1,16 @@
 # @platforma-open/milaboratories.top-antibodies.workflow
 
+## 1.3.0
+
+### Minor Changes
+
+- 2e24f7a: Disable default normalization in VJ usage plot and change spectratype/VJ usage script to run on top clonotypes if provided
+
+### Patch Changes
+
+- Updated dependencies [2e24f7a]
+  - @platforma-open/milaboratories.top-antibodies.spectratype@1.3.0
+
 ## 1.2.1
 
 ### Patch Changes

package/package.json

@@ -1,17 +1,17 @@
 {
   "name": "@platforma-open/milaboratories.top-antibodies.workflow",
-  "version": "1.2.1",
+  "version": "1.3.0",
   "type": "module",
   "description": "Block Workflow",
   "dependencies": {
-    "@platforma-sdk/workflow-tengo": "^4.5.1",
+    "@platforma-sdk/workflow-tengo": "^4.6.1",
+    "@platforma-open/milaboratories.top-antibodies.spectratype": "1.3.0",
     "@platforma-open/milaboratories.top-antibodies.sample-clonotypes": "1.0.1",
-    "@platforma-open/milaboratories.top-antibodies.umap": "1.0.1",
-    "@platforma-open/milaboratories.top-antibodies.spectratype": "1.2.0"
+    "@platforma-open/milaboratories.top-antibodies.umap": "1.0.1"
   },
   "devDependencies": {
-    "@platforma-sdk/tengo-builder": "^2.1.3",
-    "@platforma-sdk/test": "^1.30.39",
+    "@platforma-sdk/tengo-builder": "^2.1.5",
+    "@platforma-sdk/test": "^1.31.16",
     "vitest": "^2.1.8"
   },
   "scripts": {

package/src/main.tpl.tengo

@@ -19,7 +19,8 @@ wf.prepare(func(args){
 	bundleBuilder.addAnchor("main", args.inputAnchor) 
 	
 	for col in args.rankingOrder {
-		bundleBuilder.addSingle(col)
+		bundleBuilder.addAnchor(col.value.anchorName, col.value.anchorRef)
+		bundleBuilder.addSingle(col.value.column)
 	}
 
 	// Add linker column
@@ -91,14 +92,14 @@ wf.body(func(args) {
 	////////// Clonotype Filtering //////////
 	
 	// Build clonotype table
-	cloneTable := columns.xsvTableBuilder()
+	cloneTable := pframes.csvFileBuilder() 
 	cloneTable.setAxisHeader(datasetSpec.axesSpec[1].name, "clonotypeKey")
 	
 	for i, col in args.rankingOrder {
-		cloneTable.add(col, {header: "Col" + string(i)})
+		cloneTable.add(columns.getColumn(col.value.column), {header: "Col" + string(i)})
 	
 		// If column does not have main anchor axis we have to include theirs
-		colsSpec := columns.getSpec(col)
+		colsSpec := columns.getSpec(col.value.column)
 		axesNames := slices.map(colsSpec.axesSpec, func (a) { return a.name})
 		if !slices.hasElement(axesNames, datasetSpec.axesSpec[1].name) {
 			for na, ax in colsSpec.axesSpec {
@@ -113,16 +114,18 @@ wf.body(func(args) {
 	linkerAxisSpec := {}
 	for i, col in columns.getColumns("linkers") {
 		if datasetSpec.axesSpec[1].name == col.spec.axesSpec[1].name {
-			cloneTable.add(col.key, {header: "linker." + string(i)})
+			cloneTable.add(col, {header: "linker." + string(i)})
 			cloneTable.setAxisHeader(col.spec.axesSpec[0].name, "cluster_" + string(i))
 			linkerAxisSpec["cluster_" + string(i)] = col.spec.axesSpec[0]
 		} else if datasetSpec.axesSpec[1].name == col.spec.axesSpec[0].name {
-			cloneTable.add(col.key, {header: "linker." + string(i)})
+			cloneTable.add(col, {header: "linker." + string(i)})
 			cloneTable.setAxisHeader(col.spec.axesSpec[1].name, "cluster_" + string(i))
 			linkerAxisSpec["cluster_" + string(i)] = col.spec.axesSpec[1]
 		}
 	}
-	cloneTable = cloneTable.build("csv")
+	cloneTable = cloneTable.build()
+
+	topClonotypesCsv := undefined
 
 	if topClonotypes != undefined {
 		// Run sampling script
@@ -137,33 +140,36 @@ wf.body(func(args) {
 			saveStdoutContent().
 			cache(24 * 60 * 60 * 1000).
 			run()
-		
+
+		// Save top clonotypes CSV file
+		topClonotypesCsv = sampleClones.getFile("sampledClonotypes_top.csv")
+
 		// Store outputs
 		sampledColsParams := sampledColsConv.getColumns(datasetSpec, linkerAxisSpec)
-		sampledColumnsPf := xsv.importFile(sampleClones.getFile("sampledClonotypes_top.csv"), "csv", sampledColsParams)
+		sampledColumnsPf := xsv.importFile(topClonotypesCsv, "csv", sampledColsParams)
 		outputs["sampledRows"] = pframes.exportFrame(sampledColumnsPf)
 
     	// Prepare filter col subset for UMAP
     	// Avoid taking cluster axis using other params
     	sampledColsUmapParams := sampledColsUmapConv.getColumns(datasetSpec) 
-    	sampledColsUmapPf := xsv.importFile(sampleClones.getFile("sampledClonotypes_top.csv"), "csv", sampledColsUmapParams)
+    	sampledColsUmapPf := xsv.importFile(topClonotypesCsv, "csv", sampledColsUmapParams)
 		outputs["sampledRowsUmap"] = pframes.exportFrame(sampledColsUmapPf)
 	}
 	
 	////////// UMAP //////////
 	// Generate input TSV with Clonotype ID and aa sequence
-	umapTable := columns.xsvTableBuilder()
+	umapTable := pframes.tsvFileBuilder()
 	umapTable.setAxisHeader(datasetSpec.axesSpec[1].name, "clonotypeKey")
 	for col in columns.getColumns("aaSequence") {
 		if isSingleCell {
 			chainLabel := col.spec.domain["pl7.app/vdj/scClonotypeChain"]
-			umapTable.add(col.key, {header: "aaSequence." + chainLabel})
+			umapTable.add(col, {header: "aaSequence." + chainLabel})
 		} else {
 			chainLabel := col.spec.axesSpec[0].domain["pl7.app/vdj/chain"]
-			umapTable.add(col.key, {header: "aaSequence." + chainLabel})
+			umapTable.add(col, {header: "aaSequence." + chainLabel})
 		}
 	}
-	umapTable = umapTable.build("tsv")
+	umapTable = umapTable.build()
 
 	// UMAP script should go here
 	umapClones := exec.builder().
@@ -183,7 +189,7 @@ wf.body(func(args) {
 
 	////////// CDR3 Length Calculation //////////
 	
-	cdr3SeqTable := columns.xsvTableBuilder()
+	cdr3SeqTable := pframes.tsvFileBuilder()
 	cdr3SeqTable.setAxisHeader(datasetSpec.axesSpec[1].name, "clonotypeKey")
 
     // Must deal with multiple CDR3 sequences (two for each cell in single cell data)
@@ -218,7 +224,7 @@ wf.body(func(args) {
 
 	for col in cdr3Sequences {
 		headerName := makeHeaderName(col, "cdr3Sequence", isSingleCell)
-		cdr3SeqTable.add(col.key, {header: headerName})
+		cdr3SeqTable.add(col, {header: headerName})
 	}
 
 	// Process V genes
@@ -226,7 +232,7 @@ wf.body(func(args) {
 
 	for col in vGenes {
 		headerName := makeHeaderName(col, "vGene", isSingleCell)
-		cdr3SeqTable.add(col.key, {header: headerName})
+		cdr3SeqTable.add(col, {header: headerName})
 	}
 
 	// Process J genes
@@ -234,17 +240,26 @@ wf.body(func(args) {
 
 	for col in jGenes {
 		headerName := makeHeaderName(col, "jGene", isSingleCell)
-		cdr3SeqTable.add(col.key, {header: headerName})
+		cdr3SeqTable.add(col, {header: headerName})
 	}
 
-	cdr3SeqTableBuilt := cdr3SeqTable.build("tsv")
+	cdr3SeqTableBuilt := cdr3SeqTable.build()
 
 	cdr3VspectratypeCmd := exec.builder().
 		software(assets.importSoftware("@platforma-open/milaboratories.top-antibodies.spectratype:main")).
 		addFile("cdr3_sequences_input.tsv", cdr3SeqTableBuilt).
 		arg("--input_tsv").arg("cdr3_sequences_input.tsv").
 		arg("--spectratype_tsv").arg("spectratype.tsv").
-		arg("--vj_usage_tsv").arg("vj_usage.tsv").
+		arg("--vj_usage_tsv").arg("vj_usage.tsv") // no dot here
+
+	// Add top clonotypes argument and file to the builder if provided
+	if topClonotypes != undefined {
+		cdr3VspectratypeCmd = cdr3VspectratypeCmd.
+			arg("--top_clonotypes_csv").arg("topClonotypes.csv").
+			addFile("topClonotypes.csv", topClonotypesCsv)
+	}
+
+	cdr3VspectratypeCmd = cdr3VspectratypeCmd. // continue building the command
 		saveFile("spectratype.tsv").
 		saveFile("vj_usage.tsv").
 		printErrStreamToStdout().