@platforma-open/milaboratories.mixcr-amplicon-alignment.workflow 1.19.9 → 1.20.0

package/.turbo/turbo-build.log

@@ -1,6 +1,6 @@
  WARN  Issue while reading "/home/runner/work/mixcr-amplicon-alignment/mixcr-amplicon-alignment/.npmrc". Failed to replace env in config: ${NPMJS_TOKEN}
 
-> @platforma-open/milaboratories.mixcr-amplicon-alignment.workflow@1.19.9 build /home/runner/work/mixcr-amplicon-alignment/mixcr-amplicon-alignment/workflow
+> @platforma-open/milaboratories.mixcr-amplicon-alignment.workflow@1.20.0 build /home/runner/work/mixcr-amplicon-alignment/mixcr-amplicon-alignment/workflow
 > rm -rf dist && pl-tengo check && pl-tengo build
 
 Processing "src/aggregate-by-clonotype-key.tpl.tengo"...

package/CHANGELOG.md

@@ -1,5 +1,11 @@
 # @platforma-open/milaboratories.mixcr-amplicon-alignment.workflow
 
+## 1.20.0
+
+### Minor Changes
+
+- e2b65c7: Support custom reference library file
+
 ## 1.19.9
 
 ### Patch Changes

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@platforma-open/milaboratories.mixcr-amplicon-alignment.workflow",
-  "version": "1.19.9",
+  "version": "1.20.0",
   "description": "MiXCR Amplicon Alignment Workflow",
   "type": "module",
   "dependencies": {

package/src/main.tpl.tengo

@@ -41,16 +41,25 @@ wf.body(func(args) {
 		"TCRDelta": { mixcrFilter: "TRG", name: "TCR Delta", shortName: "Delta" }
 }
 
-    // Generate reference library using repseqio if both V and J genes are available
-  
-	repseqioResults := render.create(repseqioLibraryTpl, {
-		vGenes: args.vGenes,
-		jGenes: args.jGenes,
-		chains: chainInfos[chains].mixcrFilter
-	})
-	
-	referenceLibrary := repseqioResults.output("referenceLibrary")
-	debugOutput := repseqioResults.output("debugOutput")
+	// Generate reference library using repseqio, or use directly imported library file
+	isLibraryFileGzipped := false
+	referenceLibrary := false
+	debugOutput := undefined
+	libraryImportHandle := undefined
+
+	if !is_undefined(args.libraryFile) {
+		fImport := file.importFile(args.libraryFile)
+    	libraryImportHandle = fImport.handle
+    	referenceLibrary = fImport.file
+	} else {
+		repseqioResults := render.create(repseqioLibraryTpl, {
+			vGenes: args.vGenes,
+			jGenes: args.jGenes,
+			chains: chainInfos[chains].mixcrFilter
+		})
+		referenceLibrary = repseqioResults.output("referenceLibrary")
+		debugOutput = repseqioResults.output("debugOutput")
+	}
     
 
 	runMixcr := render.createEphemeral(processTpl, {
@@ -71,6 +80,7 @@ wf.body(func(args) {
 			assemblingFeature: args.assemblingFeature,
 			imputeGermline: args.imputeGermline,
 			badQualityThreshold: args.badQualityThreshold,
+			isLibraryFileGzipped: isLibraryFileGzipped,
 			stopCodonTypes: args.stopCodonTypes,
 			stopCodonReplacements: args.stopCodonReplacements
 		}, { removeUndefs: true }))
@@ -92,6 +102,23 @@ wf.body(func(args) {
 		}
 	}
 
+	if !is_undefined(args.libraryFile) {
+		exports.library = {
+			data: referenceLibrary,
+			spec: {
+				kind: "File",
+				name: "pl7.app/vdj/library",
+				domain: {
+					"pl7.app/vdj/libraryId": blockId
+				},
+				annotations: {
+					"pl7.app/vdj/isLibrary": "true",
+					"pl7.app/vdj/libraryFormat": (isLibraryFileGzipped ? "repseqio.json.gz" : "repseqio.json")
+				}
+			}
+		}
+	}
+
 	outputs := {
 		qc: pframes.exportColumnData(runMixcr.output("qc.data")),
 		reports: pframes.exportColumnData(runMixcr.output("reports.data")),
@@ -99,9 +126,15 @@ wf.body(func(args) {
 		clones: pframes.exportFrame(runMixcr.output("clones")),
 		clns: runMixcr.output("clns.data"),
 		progress: runMixcr.output("progress.data"),
-		clonotypeTables: pframes.exportFrame(runMixcr.output("clonotypeTables")),
-		debugOutput: debugOutput
+		clonotypeTables: pframes.exportFrame(runMixcr.output("clonotypeTables"))
+	}
+	outputs.referenceLibrary = file.exportFile(referenceLibrary)
+	if !is_undefined(debugOutput) {
+		outputs.debugOutput = debugOutput
 	}
+	if !is_undefined(libraryImportHandle) {
+    	outputs.libraryImportHandle = libraryImportHandle
+    }
 
 	qcReportTable := runMixcr.output("qcReportTable")
 	outputs.qcReportTable = pframes.exportFrame(qcReportTable)

package/src/mixcr-analyze.tpl.tengo

@@ -81,11 +81,16 @@ self.body(func(inputs) {
 		return "{" + parts[0] + "Begin:" + parts[1] + "End}"
 	}
 
+	libraryFileName := "library.json"
+	if !is_undefined(params.isLibraryFileGzipped) && params.isLibraryFileGzipped {
+		libraryFileName = "library.json.gz"
+	}
+
     mixcrCmdBuilder.
         arg("--assemble-clonotypes-by").arg(formatAssemblingFeature(params.assemblingFeature)).
         arg("--species").arg("custom").
-		arg("--library").arg("library.json").
-		addFile("library.json", params.referenceLibrary).
+		arg("--library").arg(libraryFileName).
+		addFile(libraryFileName, params.referenceLibrary).
         arg("--rna").
 		arg("--force-overwrite").
 		arg("--rigid-left-alignment-boundary").

package/src/mixcr-export.tpl.tengo

@@ -216,12 +216,17 @@ self.body(func(inputs) {
 
 		additionalAction(mixcrCmdBuilder)
 
+		libraryFileName := "library.json"
+		if !is_undefined(params.isLibraryFileGzipped) && params.isLibraryFileGzipped {
+			libraryFileName = "library.json.gz"
+		}
+
 		return mixcrCmdBuilder.
 			arg("clones.clns").
 			addFile("clones.clns", clnsFile).
 			arg("clones.tsv").
 			saveFile("clones.tsv").
-			addFile("library.json", params.referenceLibrary).
+			addFile(libraryFileName, params.referenceLibrary).
 			cacheHours(3).
 			run()
 	}

package/src/process.tpl.tengo

@@ -22,8 +22,6 @@ calculateExportSpecs := import(":calculate-export-specs")
 self.awaitState("InputsLocked")
 self.awaitState("params", "ResourceReady")
 self.awaitState("inputSpec", "ResourceReady")
-self.awaitState("referenceLibrary", "ResourceReady")
-self.awaitState("cdr3Sequences", "ResourceReady")
 
 self.body(func(inputs) {
 
@@ -223,7 +221,8 @@ self.body(func(inputs) {
 					tagPattern: tagPattern,
 					assemblingFeature: params.assemblingFeature,
 					imputeGermline: params.imputeGermline,
-					badQualityThreshold: params.badQualityThreshold
+					badQualityThreshold: params.badQualityThreshold,
+					isLibraryFileGzipped: params.isLibraryFileGzipped
 				}, { removeUndefs: true }),
 				limitInput: limitInput
 			},
@@ -287,7 +286,8 @@ self.body(func(inputs) {
 					cdr3SeqColumns: cdr3SeqColumns,
 					stopCodonTypes: params.stopCodonTypes,
 					stopCodonReplacements: params.stopCodonReplacements,
-					perProcessMemGB: perProcessMemGB
+					perProcessMemGB: perProcessMemGB,
+					isLibraryFileGzipped: params.isLibraryFileGzipped
 				}, { removeUndefs: true })
 			}
 		}
@@ -369,7 +369,7 @@ self.body(func(inputs) {
 		sampleIdAxisSpec: sampleIdAxisSpec,
 		chains: [chains],
 		library: referenceLibrary,
-		isLibraryFileGzipped: false,
+		isLibraryFileGzipped: params.isLibraryFileGzipped,
 		clonotypeTablesData: clonotypeTablesData,
 		hasUmi: hasUMI,
 		umiTags: umiTags,