Skip to content

Commit

Permalink
Remove RNA-SeQC from RNA-seq workflows
Browse files Browse the repository at this point in the history
  • Loading branch information
marcellevstek committed Mar 14, 2024
1 parent 24b5840 commit 7b31ad8
Show file tree
Hide file tree
Showing 4 changed files with 4 additions and 61 deletions.
1 change: 1 addition & 0 deletions docs/CHANGELOG.rst
Original file line number Diff line number Diff line change
Expand Up @@ -19,6 +19,7 @@ Added
Changed
-------
- Remove ``rnaseqc-qc`` from MultiQC report
- Remove ``rnaseqc-qc`` from RNA-seq workflows

Fixed
-----
Expand Down
20 changes: 1 addition & 19 deletions resolwe_bio/processes/workflows/bbduk_salmon_qc.py
Original file line number Diff line number Diff line change
Expand Up @@ -42,7 +42,7 @@ class WorkflowBbdukSalmonQc(Process):
entity = {
"type": "sample",
}
version = "4.3.1"
version = "4.4.0"
process_type = "data:workflow:rnaseq:salmon"
category = "Pipeline"

Expand Down Expand Up @@ -449,22 +449,4 @@ def run(self, inputs, outputs):
]
}

# RNA-SeQC tool is initiated only if annotation source is ENSEMBL
if inputs.annotation.output.source == "ENSEMBL":
input_rnaseqc = {
"alignment": alignment_qc,
"annotation": inputs.annotation,
"strand_detection_options": {
"stranded": "auto",
"cdna_index": inputs.salmon_index,
"n_reads": 5000000,
},
}
rnaseqc = Data.create(
process=BioProcess.get_latest(slug="rnaseqc-qc"),
input=input_rnaseqc,
name=f"RNA-SeQC QC report ({inputs.reads.name})",
)
input_multiqc["data"].append(rnaseqc)

Data.create(process=BioProcess.get_latest(slug="multiqc"), input=input_multiqc)
22 changes: 1 addition & 21 deletions resolwe_bio/processes/workflows/bbduk_star.py
Original file line number Diff line number Diff line change
Expand Up @@ -41,7 +41,7 @@ class WorkflowSTAR(Process):
"expression-engine": "jinja",
}
data_name = "{{ reads|name|default('?') }}"
version = "1.4.0"
version = "1.5.0"
entity = {
"type": "sample",
}
Expand Down Expand Up @@ -731,24 +731,4 @@ def run(self, inputs, outputs):
"advanced": {"dirs": True, "config": True},
}

# RNA-SeQC tool is initiated only if annotation source is ENSEMBL
if inputs.annotation.output.source == "ENSEMBL":
input_rnaseqc = {
"alignment": alignment_downsampled,
"annotation": inputs.annotation,
"strand_detection_options": {"stranded": inputs.assay_type},
}

if inputs.cdna_index:
input_rnaseqc["strand_detection_options"][
"cdna_index"
] = inputs.cdna_index

rnaseqc = Data.create(
process=BioProcess.get_latest(slug="rnaseqc-qc"),
input=input_rnaseqc,
name=f"RNA-SeQC QC report ({inputs.reads.name})",
)
input_multiqc["data"].append(rnaseqc)

Data.create(process=BioProcess.get_latest(slug="multiqc"), input=input_multiqc)
22 changes: 1 addition & 21 deletions resolwe_bio/processes/workflows/bbduk_star_featurecounts_qc.py
Original file line number Diff line number Diff line change
Expand Up @@ -46,7 +46,7 @@ class WorkflowBBDukStarFcQC(Process):
entity = {
"type": "sample",
}
version = "6.2.0"
version = "6.3.0"
process_type = "data:workflow:rnaseq:featurecounts:qc"
category = "Pipeline"

Expand Down Expand Up @@ -764,24 +764,4 @@ def run(self, inputs, outputs):
"advanced": {"dirs": True, "config": True},
}

# RNA-SeQC tool is initiated only if annotation source is ENSEMBL
if inputs.annotation.output.source == "ENSEMBL":
input_rnaseqc = {
"alignment": alignment_downsampled,
"annotation": inputs.annotation,
"strand_detection_options": {"stranded": inputs.assay_type},
}

if inputs.cdna_index:
input_rnaseqc["strand_detection_options"][
"cdna_index"
] = inputs.cdna_index

rnaseqc = Data.create(
process=BioProcess.get_latest(slug="rnaseqc-qc"),
input=input_rnaseqc,
name=f"RNA-SeQC QC report ({inputs.reads.name})",
)
input_multiqc["data"].append(rnaseqc)

Data.create(process=BioProcess.get_latest(slug="multiqc"), input=input_multiqc)

0 comments on commit 7b31ad8

Please sign in to comment.