qc/protocols: implement new protocol for 10x Cytassist Visium PEX data.

fls-bioinformatics-core · Dec 14, 2023 · 7c25a1c · 7c25a1c
1 parent 9eb9c99
commit 7c25a1c
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 4 deletions.
diff --git a/auto_process_ngs/qc/protocols.py b/auto_process_ngs/qc/protocols.py
@@ -259,7 +259,7 @@ class in the ``check_outputs`` method; new modules must be added
     },
 
     "10x_Visium_FFPE": {
-        "description": "10xGenomics Visium FFPE spatial data",
+        "description": "10xGenomics Visium FFPE spatial RNA-seq/GEX",
         "reads": {
             "seq_data": ('r2:1-50',),
             "index": ('r1',)
@@ -274,6 +274,18 @@ class in the ``check_outputs`` method; new modules must be added
         ]
     },
 
+    "10x_Visium_FFPE_PEX": {
+        "description": "10xGenomics Visium FFPE spatial PEX",
+        "reads": {
+            "seq_data": ('r2:1-50',),
+            "index": ('r1',)
+        },
+        "qc_modules": [
+            'fastqc',
+            'sequence_lengths'
+        ]
+    },
+
     "ParseEvercode": {
         "description": "Parse Biosciences Evercode data",
         "reads": {
@@ -635,9 +647,10 @@ def determine_qc_protocol(project):
                                              "10xGenomics CytAssist Visium"):
         # 10xGenomics Visium spatial transcriptomics
         if project.info.library_type in ("FFPE Spatial RNA-seq",
-                                         "FFPE Spatial GEX",
-                                         "FFPE Spatial PEX"):
+                                         "FFPE Spatial GEX"):
             protocol = "10x_Visium_FFPE"
+        elif project.info.library_type == "FFPE Spatial PEX":
+            protocol = "10x_Visium_FFPE_PEX"
         else:
             protocol = "10x_Visium"
     # Multiome ATAC+GEX

diff --git a/auto_process_ngs/test/qc/test_protocols.py b/auto_process_ngs/test/qc/test_protocols.py
@@ -658,7 +658,7 @@ def test_determine_qc_protocol_10x_visium_cytassist_ffpe_pex(self):
         project = AnalysisProject("PJB",
                                   os.path.join(self.wd,"PJB"))
         self.assertEqual(determine_qc_protocol(project),
-                         "10x_Visium_FFPE")
+                         "10x_Visium_FFPE_PEX")
 
     def test_determine_qc_protocol_10x_multiome_atac(self):
         """determine_qc_protocol: single cell multiome ATAC run (10xGenomics Multiome ATAC)