biopipen 0.24.1__py3-none-any.whl → 0.24.2__py3-none-any.whl

biopipen/__init__.py

@@ -1 +1 @@
-__version__ = "0.24.1"
+__version__ = "0.24.2"

biopipen/ns/cellranger.py

@@ -50,6 +50,7 @@ class CellRangerCount(Proc):
     script = "file://../scripts/cellranger/CellRangerCount.py"
     plugin_opts = {
         "report": "file://../reports/cellranger/CellRangerCount.svelte",
+        "report_paging": 5,
     }
 
 
@@ -98,4 +99,34 @@ class CellRangerVdj(Proc):
     script = "file://../scripts/cellranger/CellRangerVdj.py"
     plugin_opts = {
         "report": "file://../reports/cellranger/CellRangerVdj.svelte",
+        "report_paging": 5,
+    }
+
+
+class CellRangerSummary(Proc):
+    """Summarize cellranger metrics
+
+    Input:
+        indirs: The directories containing cellranger results
+            from `CellRangerCount`/`CellRangerVdj`.
+
+    Output:
+        outdir: The output directory
+
+    Envs:
+        group (type=auto): The group of the samples for boxplots.
+            If `None`, don't do boxplots.
+            It can be a dict of group names and sample names, e.g.
+            `{"group1": ["sample1", "sample2"], "group2": ["sample3"]}`
+            or a file containing the group information, with the first column
+            being the sample names and the second column being the group names.
+            The file should be tab-delimited with no header.
+    """
+    input = "indirs:dirs"
+    output = "outdir:dir:{{in.indirs | first | stem | append: '-etc.summary'}}"
+    lang = config.lang.rscript
+    script = "file://../scripts/cellranger/CellRangerSummary.R"
+    envs = {"group": None}
+    plugin_opts = {
+        "report": "file://../reports/cellranger/CellRangerSummary.svelte",
     }

biopipen/ns/cellranger_pipeline.py

@@ -0,0 +1,95 @@
+"""The cellranger pipelines
+
+Primarily cellranger process plus summary for summarizing the metrics for
+multiple samples.
+"""
+from __future__ import annotations
+from typing import TYPE_CHECKING
+
+from diot import Diot
+from pipen_args.procgroup import ProcGroup
+
+if TYPE_CHECKING:
+    from pipen import Proc
+
+
+class CellRangerCountPipeline(ProcGroup):
+    """The cellranger count pipeline
+
+    Run cellranger count for multiple samples and summarize the metrics.
+    """
+    DEFAULTS = Diot(input=None)
+
+    def post_init(self):
+        """Check if the input is a list of fastq files"""
+        if (
+            not isinstance(self.opts.input, (list, tuple))
+            or len(self.opts.input) == 0
+            or not isinstance(self.opts.input[0], (list, tuple))
+        ):
+            raise TypeError(
+                "The input of `CellRangerCountPipeline` should be a list of lists of "
+                "fastq files."
+            )
+
+    @ProcGroup.add_proc
+    def p_cellranger_count(self) -> Proc:
+        """Build CellRangerCount process"""
+        from .cellranger import CellRangerCount as _CellRangerCount
+
+        class CellRangerCount(_CellRangerCount):
+            input_data = self.opts.input
+
+        return CellRangerCount
+
+    @ProcGroup.add_proc
+    def p_cellranger_count_summary(self) -> Proc:
+        """Build CellRangerCountSummary process"""
+        from .cellranger import CellRangerSummary
+
+        class CellRangerCountSummary(CellRangerSummary):
+            requires = self.p_cellranger_count
+            input_data = lambda ch: [list(ch.iloc[:, 0])]
+
+        return CellRangerCountSummary
+
+
+class CellRangerVdjPipeline(ProcGroup):
+    """The cellranger vdj pipeline
+
+    Run cellranger vdj for multiple samples and summarize the metrics.
+    """
+    DEFAULTS = Diot(input=None)
+
+    def post_init(self):
+        """Check if the input is a list of fastq files"""
+        if (
+            not isinstance(self.opts.input, (list, tuple))
+            or len(self.opts.input) == 0
+            or not isinstance(self.opts.input[0], (list, tuple))
+        ):
+            raise TypeError(
+                "The input of `CellRangerVdjPipeline` should be a list of lists of "
+                "fastq files."
+            )
+
+    @ProcGroup.add_proc
+    def p_cellranger_vdj(self) -> Proc:
+        """Build CellRangerVdj process"""
+        from .cellranger import CellRangerVdj as _CellRangerVdj
+
+        class CellRangerVdj(_CellRangerVdj):
+            input_data = self.opts.input
+
+        return CellRangerVdj
+
+    @ProcGroup.add_proc
+    def p_cellranger_vdj_summary(self) -> Proc:
+        """Build CellRangerVdjSummary process"""
+        from .cellranger import CellRangerSummary
+
+        class CellRangerVdjSummary(CellRangerSummary):
+            requires = self.p_cellranger_vdj
+            input_data = lambda ch: [list(ch.iloc[:, 0])]
+
+        return CellRangerVdjSummary

biopipen/reports/cellranger/CellRangerCount.svelte

@@ -1,11 +1,14 @@
 {% from "utils/misc.liq" import report_jobs, table_of_images -%}
+<script>
+    import { Iframe } from "$libs";
+</script>
 
 {%- macro report_job(job, h=1) -%}
-    <iframe
+    <Iframe
         src="{{job.out.outdir}}/outs/web_summary.html"
         width="100%"
         frameborder="0"
-        style="min-height: 60vh"></iframe>
+        style="min-height: 60vh" />
 {%- endmacro -%}
 
 {%- macro head_job(job) -%}

biopipen/reports/cellranger/CellRangerSummary.svelte

@@ -0,0 +1,16 @@
+{% from "utils/misc.liq" import report_jobs -%}
+
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+
+{%- macro head_job(job) -%}
+    <h1>{{job.out.outdir | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen/reports/cellranger/CellRangerVdj.svelte

@@ -1,11 +1,14 @@
 {% from "utils/misc.liq" import report_jobs, table_of_images -%}
+<script>
+    import { Iframe } from "$libs";
+</script>
 
 {%- macro report_job(job, h=1) -%}
-    <iframe
+    <Iframe
         src="{{job.out.outdir}}/outs/web_summary.html"
         width="100%"
         frameborder="0"
-        style="min-height: 60vh"></iframe>
+        style="min-height: 60vh" />
 {%- endmacro -%}
 
 {%- macro head_job(job) -%}

biopipen/scripts/cellranger/CellRangerSummary.R

@@ -0,0 +1,189 @@
+source("{{biopipen_dir}}/utils/misc.R")
+library(rlang)
+library(dplyr)
+library(ggplot2)
+library(ggprism)
+
+theme_set(theme_prism())
+
+indirs <- {{in.indirs | r}}
+outdir <- {{out.outdir | r}}
+joboutdir <- {{job.outdir | r}}
+group <- {{envs.group | r}}
+
+if (is.character(group)) {
+    group <- read.csv(group, header = FALSE, row.names = NULL)
+    colnames(group) <- c("Sample", "Group")
+} else if (is.list(group)) {
+    group <- do_call(
+        rbind,
+        lapply(names(group), function(n) data.frame(Sample = group[[n]], Group = n))
+    )
+} else if (!is.null(group)) {
+    stop(paste0("Invalid group: ", paste0(group, collapse = ", ")))
+}
+
+cellranger_type <- NULL
+log_info("Reading and merging metrics for each sample ...")
+metrics <- NULL
+for (indir in indirs) {
+    sample <- basename(indir)
+    log_debug("- Reading metrics for sample: ", sample)
+    metric <- read.csv(
+        file.path(indir, "outs", "metrics_summary.csv"),
+        header = TRUE, row.names = NULL, check.names = FALSE)
+    metric$Sample <- sample
+    sample_cellranger_type <- ifelse(
+        file.exists(file.path(indir, "outs", "clonotypes.csv")),
+        "vdj",
+        "count"  # support more types in the future
+    )
+    cellranger_type <- cellranger_type %||% sample_cellranger_type
+    if (cellranger_type != sample_cellranger_type) {
+        stop("Multiple types of CellRanger output detected. Should be either count or vdj.")
+    }
+    if (!is.null(metrics)) {
+        missing_cols <- setdiff(colnames(metrics), colnames(metric))
+        if (length(missing_cols) > 0) {
+            log_warn('- Missing columns: {paste0(missing_cols, collapse = ", ")} in sample: {sample}')
+            metric[missing_cols] <- NA
+        }
+        missing_cols <- setdiff(colnames(metric), colnames(metrics))
+        if (length(missing_cols) > 0) {
+            log_warn('- Missing columns: {paste0(missing_cols, collapse = ", ")} in samples before {sample}')
+            metrics[missing_cols] <- NA
+        }
+    }
+    metrics <- rbind(metrics, metric)
+}
+
+if (is.null(metrics)) {
+    stop("No samples found, check the input directories.")
+}
+
+percent_columns <- sapply(colnames(metrics), function(x) {
+    is.character(metrics[[x]]) && grepl("%", metrics[[x]][1]) && x != "Sample"
+})
+percent_columns <- colnames(metrics)[percent_columns]
+# Remove %
+metrics <- metrics %>%
+    mutate(across(all_of(percent_columns), ~as.numeric(gsub("%", "", .x)))) %>%
+    rename_with(.fn = function(x) { paste0(x, " (%)") }, .cols = percent_columns) %>%
+    mutate(across(-Sample, ~as.numeric(gsub(",", "", .x))))
+
+# Save metrics
+write.table(
+    metrics,
+    file.path(outdir, "metrics.txt"),
+    sep = "\t",
+    quote = FALSE,
+    row.names = FALSE
+)
+
+add_report(
+    list(kind = "descr", content = "Metrics for all samples"),
+    list(kind = "table", src = file.path(outdir, "metrics.txt")),
+    h1 = "Metrics of all samples"
+)
+
+if (cellranger_type == "vdj") {
+    METRIC_DESCR = list(
+        `Estimated Number of Cells` = "The number of barcodes estimated to correspond to GEMs containing cells. See VDJ Cell Calling Algorithm.",
+        `Mean Read Pairs per Cell` = "Number of input read pairs divided by the estimated number of cells.",
+        `Number of Cells With Productive V-J Spanning Pair` = "Number of cell barcodes for which at least one productive sequence was found for each of TRA and TRB (or heavy and light chains, for Ig).",
+        `Number of Read Pairs` = "Total number of read pairs that were assigned to this library in demultiplexing.",
+        `Valid Barcodes` = "Fraction of reads with barcodes that match the whitelist after barcode correction.",
+        `Q30 Bases in Barcode` = "Fraction of cell barcode bases with Q-score greater than or equal to 30.",
+        `Q30 Bases in RNA Read 1` = "Fraction of read 1 bases with Q-score greater than or equal to 30. (Likewise for read 2.)",
+        `Q30 Bases in Sample Index` = "Fraction of sample index bases with Q-score greater than or equal to 30.",
+        `Q30 Bases in UMI` = "Fraction of UMI bases with Q-score ≥ 30.",
+        `Reads Mapped to Any V(D)J Gene` = "Fraction of reads that partially or wholly map to a V(D)J gene segment.",
+        `Reads Mapped to TRA` = "Fraction of reads that map partially or wholly to a TRA gene segment.",
+        `Mean Used Read Pairs per Cell` = "Mean number of read pairs used in assembly per cell barcode. These reads must have a cell barcode, map to a V(D)J gene, and have a UMI with sufficient read support, counted after subsampling.",
+        `Fraction Reads in Cells` = "Number of reads with cell barcodes divided by the number of reads with valid barcodes.",
+        `Median TRA UMIs per Cell` = "Median number of UMIs assigned to a TRA contig per cell.",
+        `Paired Clonotype Diversity` = "Effective diversity of the paired clonotypes, computed as the Inverse Simpson Index of the clonotype frequencies. A value of 1 indicates a minimally diverse sample - only one distinct clonotype was detected. A value equal to the estimated number of cells indicates a maximally diverse sample.",
+        `Cells With TRA Contig` = "Fraction of cell barcodes with at least one TRA contig annotated as a full or partial V(D)J gene.",
+        `Cells With CDR3-annotated TRA Contig` = "Fraction of cell barcodes with at least one TRA contig where a CDR3 was detected.",
+        `Cells With V-J Spanning Contig` = "Fraction of cell barcodes with at least one full-length contig.",
+        `Cells With V-J Spanning TRA Contig` = "Fraction of cell barcodes with at least one full-length TRA contig.",
+        `Cells With Productive TRA Contig` = "Fraction of cell barcodes with at least one full-length TRA contig that is productive.",
+        `Cells With Productive V-J Spanning Pair` = "Fraction of cell barcodes with at least one contig for each chain of the receptor pair that is productive."
+    )
+} else {
+    METRIC_DESCR = list(
+        `Estimated Number of Cells` = "The number of barcodes associated with cell-containing partitions, estimated from the barcode UMI count distribution.",
+        `Mean Reads per Cell` = "The total number of reads divided by the estimated number of cells.",
+        `Median Genes per Cell` = "Median number of read pairs sequenced from the cells assigned to this sample. In case of multiplexing, only cell-associated barcodes assigned exactly one CMO can be assigned to a sample.",
+        `Number of Reads` = "Total number of sequencing reads.",
+        `Valid Barcodes` = "Fraction of reads with cell-barcodes that match the whitelist.",
+        `Sequencing Saturation` = 'Fraction of reads originating from an already-observed UMI. This is a function of library complexity and sequencing depth. More specifically, this is a ratio where: the denominator is the number of confidently-mapped reads with a valid cell-barcode and valid UMI, and the numerator is the subset of those reads that had a non-unique combination of (cell-barcode, UMI, gene). This metric was called "cDNA PCR Duplication" in versions of Cell Ranger prior to 1.2.',
+        `Q30 Bases in Barcode` = "Fraction of bases with Q-score at least 30 in the cell barcode sequences. This is the i7 index (I1) read for the Single Cell 3' v1 chemistry and the R1 read for the Single Cell 3' v2 chemistry.",
+        `Q30 Bases in RNA` = "Fraction of bases with Q-score at least 30 in the RNA read sequences. This is Illumina R1 for the Single Cell 3' v1 chemistry and Illumina R2 for the Single Cell 3' v2 chemistry.",
+        `Q30 Bases in UMI` = "Fraction of bases with Q-score at least 30 in the UMI sequences. This is the R2 read for the Single Cell 3' v1 chemistry and the R1 read for the Single Cell 3' v2 chemistry.",
+        `Reads Mapped to Genome` = "Fraction of reads that mapped to the genome.",
+        `Reads Mapped Confidently to Genome` = "Fraction of reads that mapped uniquely to the genome. If a read mapped to exonic loci from a single gene and also to non-exonic loci, it is considered uniquely mapped to one of the exonic loci.",
+        `Reads Mapped Confidently to Intergenic Regions` = "Fraction of reads that mapped to the intergenic regions of the genome with a high mapping quality score as reported by the aligner.",
+        `Reads Mapped Confidently to Intronic Regions` = "Fraction of reads that mapped to the intronic regions of the genome with a high mapping quality score as reported by the aligner.",
+        `Reads Mapped Confidently to Exonic Regions` = "Fraction of reads that mapped to the exonic regions of the genome with a high mapping quality score as reported by the aligner.",
+        `Reads Mapped Confidently to Transcriptome` = "Fraction of reads that mapped to a unique gene in the transcriptome with a high mapping quality score as reported by the aligner. The read must be consistent with annotated splice junctions when include-introns=false. These reads are considered for UMI counting.",
+        `Reads Confidently Mapped Antisense` = "Fraction of reads confidently mapped to the transcriptome, but on the opposite strand of their annotated gene. A read is counted as antisense if it has any alignments that are consistent with an exon of a transcript but antisense to it, and has no sense alignments.",
+        `Total Genes Detected Median UMI Counts per Cell` = "The number of genes with at least one UMI count in any cell."
+    )
+}
+log_info("Plotting metrics ...")
+for (metric in colnames(metrics)) {
+    if (metric == "Sample") { next }
+    metric_name <- sub(" \\(%\\)$", "", metric)
+    log_info("- {metric_name}")
+
+    add_report(
+        list(
+            kind = "descr",
+            content = METRIC_DESCR[[metric_name]] %||% paste0("Metric: ", metric)
+        ),
+        h1 = metric
+    )
+
+    # barplot
+    p <- ggplot(metrics, aes(x = Sample, y = !!sym(metric))) +
+        geom_bar(stat = "identity", fill = "steelblue") +
+        labs(x = "Sample", y = metric) +
+        theme(axis.text.x = element_text(angle = 90, hjust = 1))
+
+    figfile <- file.path(outdir, paste0(slugify(metric), ".barplot.png"))
+    png(figfile, height = 600, res = 100, width = nrow(metrics) * 30 + 200)
+    print(p)
+    dev.off()
+
+    add_report(
+        list(src = figfile, name = "By Sample"),
+        ui = "table_of_images",
+        h1 = metric
+    )
+
+    if (is.null(group)) { next }
+    # boxplot, if group is provided
+    # group: Sample, Group
+    pdata <- group %>%
+        left_join(metrics, by = "Sample") %>%
+        mutate(Group = factor(Group, levels = unique(Group)))
+
+    p <- ggplot(pdata, aes(x = Group, y = !!sym(metric))) +
+        geom_boxplot(fill = "steelblue") +
+        labs(x = "Group", y = metric) +
+        theme(axis.text.x = element_text(angle = 90, hjust = 1))
+
+    figfile <- file.path(outdir, paste0(slugify(metric), ".boxplot.png"))
+    png(figfile, height = 600, res = 100, width = length(unique(pdata$Group)) * 30 + 200)
+    print(p)
+    dev.off()
+
+    add_report(
+        list(src = figfile, name = "By Group"),
+        ui = "table_of_images",
+        h1 = metric
+    )
+}
+
+save_report(joboutdir)

biopipen/scripts/scrna/CellsDistribution.R

@@ -8,7 +8,6 @@ library(dplyr)
 library(ggplot2)
 library(ggVennDiagram)
 library(UpSetR)
-library(slugify)
 library(circlize)
 library(ComplexHeatmap)
 
@@ -127,8 +126,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$sec_dir <- file.path(outdir, out$section_slug)
     if (create) {
@@ -139,7 +138,7 @@ casename_info <- function(casename, create = FALSE) {
 
 plot_heatmap <- function(group, meta, case, info, cluster_order_val) {
     log_info(paste("- Running heatmap for case:", info$casename, "group:", group))
-    hmfile <- file.path(info$sec_dir, paste0(info$case_slug, ".", slugify(group, tolower = FALSE), ".heatmap.png"))
+    hmfile <- file.path(info$sec_dir, paste0(info$case_slug, ".", slugify(group), ".heatmap.png"))
 
     # A matrix: 10 × 8 of type int
     #                   g3	g6	g0	g1	g7	g5	g4	g8
@@ -435,7 +434,7 @@ do_overlap <- function(section) {
         stop(paste0("Not enough cases for overlap for section: ", section))
     }
 
-    sec_dir <- file.path(outdir, slugify(section, tolower = FALSE))
+    sec_dir <- file.path(outdir, slugify(section))
     venn_plot <- file.path(sec_dir, "venn.png")
     venn_p <- ggVennDiagram(overlap_cases, label_percent_digit = 1) +
         scale_fill_distiller(palette = "Reds", direction = 1) +

biopipen/scripts/scrna/MarkersFinder.R

@@ -14,7 +14,6 @@ library(future)
 library(tidyseurat)
 library(ggVennDiagram)
 library(UpSetR)
-library(slugify)
 
 log_info("Setting up EnrichR ...")
 setEnrichrSite("Enrichr")
@@ -205,8 +204,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$casedir <- file.path(outdir, out$section_slug, out$case_slug)
     if (create) {

biopipen/scripts/scrna/MetaMarkers.R

@@ -12,7 +12,6 @@ library(ggplot2)
 library(ggprism)
 library(parallel)
 library(tidyseurat)
-library(slugify)
 
 setEnrichrSite("Enrichr")
 
@@ -133,8 +132,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$casedir <- file.path(outdir, out$section_slug, out$case_slug)
     if (create) {
@@ -283,7 +282,7 @@ do_case <- function(casename) {
         # Plot the top 10 genes in each group with violin plots
         geneplots <- list()
         for (gene in markers$gene) {
-            outfile <- file.path(plotdir, paste0(slugify(gene, tolower = FALSE), ".png"))
+            outfile <- file.path(plotdir, paste0(slugify(gene), ".png"))
             p <- ggplot(df, aes_string(x="GROUP", y=bQuote(gene), fill="GROUP")) +
                 geom_violin(alpha = .8) +
                 geom_boxplot(width=0.1, fill="white") +

biopipen/scripts/scrna/RadarPlots.R

@@ -7,7 +7,6 @@ library(tidyr)
 library(tibble)
 library(ggplot2)
 library(ggradar)
-library(slugify)
 library(ggprism)
 
 # input/output
@@ -143,8 +142,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$casedir <- file.path(outdir, out$section_slug, out$case_slug)
     if (create) {

biopipen/scripts/scrna/ScFGSEA.R

@@ -4,7 +4,6 @@ source("{{biopipen_dir}}/utils/mutate_helpers.R")
 library(rlang)
 library(Seurat)
 library(tidyseurat)
-library(slugify)
 
 srtfile <- {{in.srtobj | r}}  # nolint
 outdir <- {{out.outdir | r}}  # nolint
@@ -122,8 +121,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$casedir <- file.path(outdir, out$section_slug, out$case_slug)
     if (create) {

biopipen/scripts/scrna/SeuratClusterStats-hists.R

@@ -64,7 +64,7 @@ if (is.null(hists) || length(hists) == 0) {
         hists[[name]] <- list_update(hists_defaults, hists[[name]])
         case <- hists[[name]]
 
-        odir <- file.path(outdir, "hists", slugify(name, tolower = FALSE))
+        odir <- file.path(outdir, "hists", slugify(name))
         dir.create(odir, recursive=TRUE, showWarnings=FALSE)
 
         h1 <- name

biopipen/scripts/scrna/SeuratClusterStats.R

@@ -1,7 +1,6 @@
 source("{{biopipen_dir}}/utils/misc.R")
 source("{{biopipen_dir}}/utils/mutate_helpers.R")
 source("{{biopipen_dir}}/utils/plot.R")
-library(slugify)
 library(Seurat)
 library(rlang)
 library(dplyr)

biopipen/scripts/scrna/SeuratPreparing.R

@@ -5,7 +5,6 @@ library(future)
 library(bracer)
 library(ggplot2)
 library(tidyseurat)
-library(slugify)
 
 metafile = {{in.metafile | quote}}
 rdsfile = {{out.rdsfile | quote}}
@@ -17,13 +16,17 @@ options(future.globals.maxSize = 80000 * 1024^2)
 options(Seurat.object.assay.version = "v5")
 plan(strategy = "multicore", workers = envs$ncores)
 
+.stringify_list <- function(x) {
+    paste(sapply(names(x), function(n) paste(n, x[[n]], sep = " = ") ), collapse = "; ")
+}
+
 add_report(
     list(
         kind = "descr",
         name = "Filters applied",
         content = paste0(
             "<p>Cell filters: ", html_escape(envs$cell_qc), "</p>",
-            "<p>Gene filters: ", html_escape(envs$gene_qc), "</p>"
+            "<p>Gene filters: ", html_escape(.stringify_list(envs$gene_qc)), "</p>"
         )
     ),
     h1 = "Filters and QC"
@@ -177,7 +180,7 @@ for (feat in feats) {
             position = position_jitterdodge(jitter.width = 0.4, dodge.width = 0.9)
         ) + scale_color_manual(values = c("#181818", pal_biopipen()(1)), breaks = c(TRUE, FALSE))
 
-    vlnplot = file.path(plotsdir, paste0(slugify(feat, tolower = FALSE), ".vln.png"))
+    vlnplot = file.path(plotsdir, paste0(slugify(feat), ".vln.png"))
     png(
         vlnplot,
         width = 800 + length(samples) * 15, height = 600, res = 100
@@ -221,7 +224,7 @@ for (feat in setdiff(feats, "nCount_RNA")) {
     NoLegend() +
     scale_color_manual(values = c("#181818", pal_biopipen()(1)), breaks = c(TRUE, FALSE))
 
-    scatfile = file.path(plotsdir, paste0(slugify(feat, tolower = FALSE), "-nCount_RNA.scatter.png"))
+    scatfile = file.path(plotsdir, paste0(slugify(feat), "-nCount_RNA.scatter.png"))
     png(scatfile, width = 800, height = 600, res = 100)
     print(scat_p)
     dev.off()
@@ -302,7 +305,7 @@ if (envs$use_sct) {
     log_info("- Running SCTransform ...")
     SCTransformArgs <- envs$SCTransform
     # log to stdout but don't populate it to running log
-    print("  SCTransform: {.formatArgs(SCTransformArgs)}")
+    print(paste0("  SCTransform: ", .formatArgs(SCTransformArgs)))
     log_debug("  SCTransform: {.formatArgs(SCTransformArgs)}")
     SCTransformArgs$object <- sobj
     sobj <- do_call(SCTransform, SCTransformArgs)
@@ -310,21 +313,21 @@ if (envs$use_sct) {
 } else {
     log_info("- Running NormalizeData ...")
     NormalizeDataArgs <- envs$NormalizeData
-    print("  NormalizeData: {.formatArgs(NormalizeDataArgs)}")
+    print(paste0("  NormalizeData: ", .formatArgs(NormalizeDataArgs)))
     log_debug("  NormalizeData: {.formatArgs(NormalizeDataArgs)}")
     NormalizeDataArgs$object <- sobj
     sobj <- do_call(NormalizeData, NormalizeDataArgs)
 
     log_info("- Running FindVariableFeatures ...")
     FindVariableFeaturesArgs <- envs$FindVariableFeatures
-    print("  FindVariableFeatures: {.formatArgs(FindVariableFeaturesArgs)}")
+    print(paste0("  FindVariableFeatures: ", .formatArgs(FindVariableFeaturesArgs)))
     log_debug("  FindVariableFeatures: {.formatArgs(FindVariableFeaturesArgs)}")
     FindVariableFeaturesArgs$object <- sobj
     sobj <- do_call(FindVariableFeatures, FindVariableFeaturesArgs)
 
     log_info("- Running ScaleData ...")
     ScaleDataArgs <- envs$ScaleData
-    print("  ScaleData: {.formatArgs(ScaleDataArgs)}")
+    print(paste0("  ScaleData: ", .formatArgs(ScaleDataArgs)))
     log_debug("  ScaleData: {.formatArgs(ScaleDataArgs)}")
     ScaleDataArgs$object <- sobj
     sobj <- do_call(ScaleData, ScaleDataArgs)
@@ -333,7 +336,7 @@ if (envs$use_sct) {
 log_info("- Running RunPCA ...")
 RunPCAArgs <- envs$RunPCA
 RunPCAArgs$npcs <- if (is.null(RunPCAArgs$npcs)) { 50 } else { min(RunPCAArgs$npcs, ncol(sobj) - 1) }
-print("  RunPCA: {.formatArgs(RunPCAArgs)}")
+print(paste0("  RunPCA: ", .formatArgs(RunPCAArgs)))
 log_debug("  RunPCA: {.formatArgs(RunPCAArgs)}")
 RunPCAArgs$object <- sobj
 sobj <- do_call(RunPCA, RunPCAArgs)
@@ -367,7 +370,7 @@ if (!envs$no_integration) {
     if (is.null(IntegrateLayersArgs$new.reduction)) {
         IntegrateLayersArgs$new.reduction <- new_reductions[[method]]
     }
-    print("  IntegrateLayers: {.formatArgs(IntegrateLayersArgs)}")
+    print(paste0("  IntegrateLayers: ", .formatArgs(IntegrateLayersArgs)))
     log_debug("  IntegrateLayers: {.formatArgs(IntegrateLayersArgs)}")
     IntegrateLayersArgs$object <- sobj
     sobj <- do_call(IntegrateLayers, IntegrateLayersArgs)

biopipen/scripts/scrna/TopExpressingGenes.R

@@ -5,7 +5,6 @@ library(tibble)
 library(enrichR)
 library(rlang)
 library(dplyr)
-library(slugify)
 library(ggprism)
 
 setEnrichrSite("Enrichr")
@@ -134,8 +133,8 @@ casename_info <- function(casename, create = FALSE) {
         casename = casename,
         section = sec_case_names[1],
         case = cname,
-        section_slug = slugify(sec_case_names[1], tolower = FALSE),
-        case_slug = slugify(cname, tolower = FALSE)
+        section_slug = slugify(sec_case_names[1]),
+        case_slug = slugify(cname)
     )
     out$casedir <- file.path(outdir, out$section_slug, out$case_slug)
     if (create) {

biopipen/scripts/scrna_metabolic_landscape/MetabolicFeatures.R

@@ -3,7 +3,6 @@ source("{{biopipen_dir}}/utils/gsea.R")
 
 library(parallel)
 library(Seurat)
-library(slugify)
 
 sobjfile <- {{ in.sobjfile | r }}
 outdir <- {{ out.outdir | r }}
@@ -44,7 +43,7 @@ sobj <- readRDS(sobjfile)
 do_one_group <- function(obj, features, group, outputdir, h1) {
     log_info(paste("- Processing group", grouping, ":", group))
     groupname = paste0(grouping_prefix, group)
-    odir = file.path(outputdir, slugify(groupname, tolower = FALSE))
+    odir = file.path(outputdir, slugify(groupname))
     dir.create(odir, showWarnings = FALSE)
 
     classes = as.character(obj@meta.data[[grouping]])
@@ -109,7 +108,7 @@ do_one_subset <- function(s, subset_col, subset_prefix) {
         outputdir <- file.path(outdir, "ALL")
         subset_obj <- sobj
     } else {
-        outputdir <- file.path(outdir, slugify(paste0(subset_prefix, s), tolower = FALSE))
+        outputdir <- file.path(outdir, slugify(paste0(subset_prefix, s)))
         subset_code <- paste0("subset(sobj, subset = ", subset_col, "=='", s, "')")
         subset_obj <- eval(parse(text = subset_code))
     }

biopipen/scripts/scrna_metabolic_landscape/MetabolicFeaturesIntraSubset.R

@@ -4,7 +4,6 @@ source("{{biopipen_dir}}/utils/gsea.R")
 library(parallel)
 library(scater)
 library(Seurat)
-library(slugify)
 
 sobjfile <- {{ in.sobjfile | r }}
 outdir <- {{ out.outdir | r }}
@@ -139,7 +138,7 @@ do_one_group <- function(group) {
     )
     obj = eval(parse(text = group_code))
     groupname = paste0(grouping_prefix, group)
-    groupdir = file.path(outdir, slugify(groupname, tolower = FALSE))
+    groupdir = file.path(outdir, slugify(groupname))
     dir.create(groupdir, showWarnings = FALSE)
 
     report = list()

biopipen/scripts/scrna_metabolic_landscape/MetabolicPathwayHeterogeneity.R

@@ -8,7 +8,6 @@ library(ggprism)
 library(Matrix)
 library(sparseMatrixStats)
 library(Seurat)
-library(slugify)
 
 sobjfile <- {{ in.sobjfile | r }}
 outdir <- {{ out.outdir | r }}
@@ -53,7 +52,7 @@ do_one_subset <- function(s, subset_col, subset_prefix) {
         subset_dir = file.path(outdir, "ALL")
         subset_obj = sobj
     } else {
-        subset_dir = file.path(outdir, slugify(paste0(subset_prefix, s), tolower = FALSE))
+        subset_dir = file.path(outdir, slugify(paste0(subset_prefix, s)))
         subset_code = paste0("subset(sobj, subset = ", subset_col, " == '", s, "')")
         subset_obj = eval(parse(text = subset_code))
     }

biopipen/scripts/tcr/CDR3AAPhyschem.R

@@ -9,7 +9,6 @@ library(hash)
 library(glmnet)
 library(broom.mixed)
 library(stringr)
-library(slugify)
 
 immdatafile = {{in.immdata | quote}}
 srtobjfile = {{in.srtobj | r}}

biopipen/scripts/tcr/CloneResidency.R

@@ -10,7 +10,6 @@ library(ggplot2)
 library(ggprism)
 library(ggVennDiagram)
 library(UpSetR)
-library(slugify)
 
 theme_set(theme_prism())
 
@@ -163,19 +162,19 @@ perpare_case <- function(casename, case) {
 
     # Create case-specific directories
     # Scatter plots
-    scatter_dir <- file.path(outdir, slugify(casename, tolower = FALSE), "scatter")
+    scatter_dir <- file.path(outdir, slugify(casename), "scatter")
     dir.create(scatter_dir, recursive = TRUE, showWarnings = FALSE)
 
     # Venn diagrams
-    venn_dir <- file.path(outdir, slugify(casename, tolower = FALSE), "venn")
+    venn_dir <- file.path(outdir, slugify(casename), "venn")
     dir.create(venn_dir, recursive = TRUE, showWarnings = FALSE)
 
     # Upset plots
-    upset_dir <- file.path(outdir, slugify(casename, tolower = FALSE), "upset")
+    upset_dir <- file.path(outdir, slugify(casename), "upset")
     dir.create(upset_dir, recursive = TRUE, showWarnings = FALSE)
 
     # Counts
-    counts_dir <- file.path(outdir, slugify(casename, tolower = FALSE), "counts")
+    counts_dir <- file.path(outdir, slugify(casename), "counts")
     dir.create(counts_dir, recursive = TRUE, showWarnings = FALSE)
 
     case
@@ -396,7 +395,7 @@ headings <- function(section, casename, subject) {
 }
 
 handle_subject <- function(i, subjects, casename, case) {
-    casedir = file.path(outdir, slugify(casename, tolower = FALSE))
+    casedir = file.path(outdir, slugify(casename))
     # Generate a residency table
     # |    CDR3.aa    | Tumor | Normal |
     # | SEABESRWEFAEF | 0     | 10     |

biopipen/scripts/tcr/Immunarch-diversity.R

@@ -657,9 +657,9 @@ run_div_case = function(casename) {
     log_info("Processing case: {casename} ...")
     case = div_cases[[casename]]
     if (case$method == "raref") {
-        ddir = file.path(outdir, "rarefraction", slugify(casename, tolower = FALSE))
+        ddir = file.path(outdir, "rarefraction", slugify(casename))
     } else {
-        ddir = file.path(div_dir, slugify(casename, tolower = FALSE))
+        ddir = file.path(div_dir, slugify(casename))
     }
     dir.create(ddir, recursive = TRUE, showWarnings = FALSE)
 

biopipen/scripts/tcr/Immunarch-kmer.R

@@ -87,7 +87,7 @@ for (name in names(cases)) {
 do_one_case_kmer = function(name, case, kmer_dir) {
     # print(paste0("  Case: ", name))
     log_info("Processing case: {name} ...")
-    odir = file.path(kmer_dir, slugify(name, tolower = FALSE))
+    odir = file.path(kmer_dir, slugify(name))
     dir.create(odir, showWarnings = FALSE)
 
     if (!is.null(case$subset)) {

biopipen/scripts/tcr/Immunarch-spectratyping.R

@@ -44,7 +44,7 @@ if (is.null(spects$cases) || length(spects$cases) == 0) {
 do_one_case_spectratyping = function(name, case, spect_dir) {
     # print(paste0("  Case: ", name))
     log_info("- Processing case: {name} ...")
-    odir = file.path(spect_dir, slugify(name, tolower = FALSE))
+    odir = file.path(spect_dir, slugify(name))
     dir.create(odir, showWarnings = FALSE)
 
     if (!is.null(case$subset)) {
@@ -61,7 +61,7 @@ do_one_case_spectratyping = function(name, case, spect_dir) {
             .quant = case$quant,
             .col = case$col
         )
-        spectfile = file.path(odir, paste0(slugify(sample, tolower = FALSE), ".spect.png"))
+        spectfile = file.path(odir, paste0(slugify(sample), ".spect.png"))
         png(
             spectfile,
             res = case$devpars$res,

biopipen/scripts/tcr/Immunarch-vjjunc.R

@@ -47,7 +47,7 @@ dir.create(vjjunc_dir, showWarnings = FALSE)
 
 do_one_case_vjjunc <- function(name, case) {
     log_info("Processing case: {name} ...")
-    odir = file.path(vjjunc_dir, slugify(name, tolower = FALSE))
+    odir = file.path(vjjunc_dir, slugify(name))
     dir.create(odir, showWarnings = FALSE)
 
     if (!is.null(case$subset)) {

biopipen/scripts/tcr/Immunarch.R

@@ -12,7 +12,6 @@ library(glue)
 library(tidyr)
 library(tibble)
 library(logger)
-library(slugify)
 
 log_info("Loading arguments ...")
 theme_set(theme_prism())

biopipen/scripts/tcr/TCRClusterStats.R

@@ -6,7 +6,6 @@ library(dplyr)
 library(rlang)
 library(immunarch)
 library(ggprism)
-library(slugify)
 
 immfile = {{in.immfile | quote}}
 outdir = {{out.outdir | quote}}
@@ -56,7 +55,7 @@ sample_diversity_cases = expand_cases(sample_diversity_envs)
 cluster_size_distribution = function(name) {
     log_info("- Working on cluster size distribution: {name}")
 
-    odir = file.path(outdir, "ClusterSizeDistribution", slugify(name, tolower = FALSE))
+    odir = file.path(outdir, "ClusterSizeDistribution", slugify(name))
     dir.create(odir, showWarnings = FALSE, recursive = TRUE)
     case = cluster_size_cases[[name]]
 
@@ -100,7 +99,7 @@ cluster_size_distribution = function(name) {
 shared_clusters = function(name) {
     log_info("- Working on shared clusters: {name}")
 
-    odir = file.path(outdir, "SharedClusters", slugify(name, tolower = FALSE))
+    odir = file.path(outdir, "SharedClusters", slugify(name))
     dir.create(odir, showWarnings = FALSE, recursive = TRUE)
     case = shared_clusters_cases[[name]]
     if (!is.null(case$grouping)) {
@@ -169,7 +168,7 @@ shared_clusters = function(name) {
 }
 
 shared_clusters_by_grouping = function(name) {
-    odir = file.path(outdir, "SharedClusters", slugify(name, tolower = FALSE))
+    odir = file.path(outdir, "SharedClusters", slugify(name))
     case = shared_clusters_cases[[name]]
 
     data = list()
@@ -241,7 +240,7 @@ div_methods = list(
 sample_diversity = function(name) {
     log_info("- Working on sample diversity: {name}")
 
-    odir = file.path(outdir, "SampleDiversity", slugify(name, tolower = FALSE))
+    odir = file.path(outdir, "SampleDiversity", slugify(name))
     dir.create(odir, showWarnings = FALSE, recursive = TRUE)
     case = sample_diversity_cases[[name]]
 

biopipen/utils/gsea.R

@@ -1,8 +1,13 @@
 library(ggplot2)
 library(dplyr)
 library(tibble)
-library(slugify)
 
+.slugify <- function(x, non_alphanum_replace="-", collapse_replace=TRUE, tolower=FALSE) {
+    x <- gsub("[^[:alnum:]_]", non_alphanum_replace, x)
+    if(collapse_replace) x <- gsub(paste0(non_alphanum_replace, "+"), non_alphanum_replace, x)
+    if(tolower) x <- tolower(x)
+    x
+}
 
 localizeGmtfile <- function(gmturl, cachedir = tempdir()) {
     # Download the GMT file and save it to cachedir
@@ -176,7 +181,7 @@ runFGSEA = function(
     dev.off()
 
     for (pathway in topPathways) {
-        enrfig = file.path(outdir, paste0("fgsea_", slugify(pathway), ".png"))
+        enrfig = file.path(outdir, paste0("fgsea_", .slugify(pathway), ".png"))
         png(enrfig, res=100, width=1000, height=800)
         print(plotEnrichment(
             envs$pathways[[pathway]],

biopipen/utils/misc.R

@@ -22,14 +22,17 @@ bQuote <- function(x) {
     }
 }
 
-slugify <- function(x, non_alphanum_replace="", space_replace="_", tolower=TRUE) {
-  x <- gsub("[^[:alnum:] ]", non_alphanum_replace, x)
-  x <- trimws(x)
-  x <- gsub("[[:space:]]", space_replace, x)
-
-  if(tolower) { x <- tolower(x) }
-
-  return(x)
+#' Slugify a string
+#' @param x A string
+#' @param non_alphanum_replace Replace non-alphanumeric characters
+#' @param collapse_replace Collapse consecutive non-alphanumeric character replacements
+#' @param tolower Convert to lowercase
+#' @return A slugified string
+slugify <- function(x, non_alphanum_replace="-", collapse_replace=TRUE, tolower=FALSE) {
+    x <- gsub("[^[:alnum:]_]", non_alphanum_replace, x)
+    if(collapse_replace) x <- gsub(paste0(non_alphanum_replace, "+"), non_alphanum_replace, x)
+    if(tolower) x <- tolower(x)
+    x
 }
 
 do_call <- function (what, args, quote = FALSE, envir = parent.frame())  {

{biopipen-0.24.1.dist-info → biopipen-0.24.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: biopipen
-Version: 0.24.1
+Version: 0.24.2
 Summary: Bioinformatics processes/pipelines that can be run from `pipen run`
 License: MIT
 Author: pwwang

{biopipen-0.24.1.dist-info → biopipen-0.24.2.dist-info}/RECORD

@@ -1,4 +1,4 @@
-biopipen/__init__.py,sha256=23Q5wqrB_jTfEuw8bupXFhwqZO0_-oiISvjn1gCVUsA,23
+biopipen/__init__.py,sha256=UBF3OYTcBAovta3ux5ybxb0MZYAGpGO79WH_ax2NGeI,23
 biopipen/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/core/config.py,sha256=edK5xnDhM8j27srDzsxubi934NMrglLoKrdcC8qsEPk,1069
 biopipen/core/config.toml,sha256=Rn7Cta7WsMtmQkKGC4h9d5dU_STaIVBgR8UliiGgL6o,1757
@@ -10,7 +10,8 @@ biopipen/ns/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/ns/bam.py,sha256=5AsYrB0mtr_mH6mCL6gjJ5rC4NywpjFkpFjUrBGp7Fk,9301
 biopipen/ns/bcftools.py,sha256=puCDfIL-1z6cz2y1Rlz-ESNIr8xJgeIjEQ440qicCvM,3467
 biopipen/ns/bed.py,sha256=UN38qUChDeE-ipuSBY8RVLwvJqM2wxSRmlhOiDo4JG0,5395
-biopipen/ns/cellranger.py,sha256=vwyPNaiRIFMp1e8kSJ1UNsE9ZcUykNniakKE-BqZR60,3681
+biopipen/ns/cellranger.py,sha256=2FwE8YNsGEXn-vS0AmO7m5xmlTgNQjL1lhWUeJ6XNJA,4796
+biopipen/ns/cellranger_pipeline.py,sha256=D6gvIeasHjDCdro7f4wjomxRYTtsJT77Ld47XzlU3Fw,2891
 biopipen/ns/cnv.py,sha256=vq6dZfEOyuVuqg3nP6FQtNmQ-JocpBJMX9IYlZ0OPD0,6803
 biopipen/ns/cnvkit.py,sha256=5mA2Q8-YDs4g1HoxtpB_NWnyZYwEThNr3s3wlubLQrQ,31130
 biopipen/ns/cnvkit_pipeline.py,sha256=2fJLn70L2jJ81ZMNdnU84Sf3HoKA2CSnHuDzLGR8jmw,36854
@@ -29,8 +30,9 @@ biopipen/ns/web.py,sha256=3zucrDo-IVsSnIvlw-deoScuxqWa6OMTm8Vo-R4E44Q,2224
 biopipen/reports/bam/CNAClinic.svelte,sha256=D4IxQcgDCPQZMbXog-aZP5iJEQTK2N4i0C60e_iXyfs,213
 biopipen/reports/bam/CNVpytor.svelte,sha256=s03SlhbEPd8-_44Dy_cqE8FSErhUdqStLK39te5o7ZE,1364
 biopipen/reports/bam/ControlFREEC.svelte,sha256=OwN96RW0dN-gtQ1zWKbXYZCYkkrOC0RQmP3UG4x7zqU,837
-biopipen/reports/cellranger/CellRangerCount.svelte,sha256=l153cN-4F1aTN6TP9t-BkxMFiIyTsD2uOnRQOcWefqo,418
-biopipen/reports/cellranger/CellRangerVdj.svelte,sha256=l153cN-4F1aTN6TP9t-BkxMFiIyTsD2uOnRQOcWefqo,418
+biopipen/reports/cellranger/CellRangerCount.svelte,sha256=InkkGold9Vta8vSPeE7WgFuFWuMVBkQvwc2yxWFGVM4,466
+biopipen/reports/cellranger/CellRangerSummary.svelte,sha256=WEYPwmcmgtSqpFpTvNPV40yZR-7ERa5LgZni2RXxUZQ,348
+biopipen/reports/cellranger/CellRangerVdj.svelte,sha256=InkkGold9Vta8vSPeE7WgFuFWuMVBkQvwc2yxWFGVM4,466
 biopipen/reports/cnv/AneuploidyScore.svelte,sha256=x0LbhqjauZpqMzmzDWmYgx-rEh5Tzo8qBrXcLcM0h78,1020
 biopipen/reports/cnv/AneuploidyScoreSummary.svelte,sha256=AWlns70mChJGhH3z8r5uXuI4tc7VbVN_cOUdqBr3ZKg,4414
 biopipen/reports/cnv/TMADScoreSummary.svelte,sha256=tJutaMOqeXxKroAosOIqOJVyhTTFet-soMwuOYVHTYU,2060
@@ -78,6 +80,7 @@ biopipen/scripts/bed/BedConsensus.py,sha256=gfAxuIalvCEpS0tiOyAJGPYGgHN0L-hm0K37
 biopipen/scripts/bed/BedLiftOver.sh,sha256=Y4gBsz9w4zhE29UmWojO6F4PXMMMWC1uCzjrxa19eOs,256
 biopipen/scripts/bed/BedtoolsMerge.py,sha256=TjKO5MpUzDj931bQAWku2660MVSiZzdMHt_v2Xbt0IE,355
 biopipen/scripts/cellranger/CellRangerCount.py,sha256=ZDcry8suLhulXiTsl01LGKmSJkewJ-TgHazLtfsBr6U,2516
+biopipen/scripts/cellranger/CellRangerSummary.R,sha256=8EqP4j074RLi5Adw47xyqQwkDEcvSS6Yfk0UgKKYOuo,11045
 biopipen/scripts/cellranger/CellRangerVdj.py,sha256=-QbhPKqFBZ15Es6NJaU7Lwf1KQW_3Lyv0aISh-Urk2M,2504
 biopipen/scripts/cnv/AneuploidyScore.R,sha256=liAN8u8_lj8voJ01oBW9Dw09yi388KF5f_gwPOv0wdE,8437
 biopipen/scripts/cnv/AneuploidyScoreSummary.R,sha256=9Zni5zqYfzevs5XSAt3fqD9WZ_RWr_ByUnXReKLLWoY,12337
@@ -113,47 +116,47 @@ biopipen/scripts/scrna/CellTypeAnnotation-hitype.R,sha256=6_DBAlLKcHqaMyWGZWvTd4
 biopipen/scripts/scrna/CellTypeAnnotation-sccatch.R,sha256=1ejye0hs-EOwzzdP9gFWSLPcF6dOAA6VmNKXEjmS11E,1654
 biopipen/scripts/scrna/CellTypeAnnotation-sctype.R,sha256=1nKRtzhVoJ9y0yMg1sgI6u7czsrk2cN0FvNUCZo8l-o,3878
 biopipen/scripts/scrna/CellTypeAnnotation.R,sha256=6Le1SvZcKI8D0SLkFZ5SibGsW9ZWqirnBl3Q1BNZOuU,513
-biopipen/scripts/scrna/CellsDistribution.R,sha256=UDsiZp94oxyENpDS6H649SfeaTFMX5Ip1gG0mBJukkA,16647
+biopipen/scripts/scrna/CellsDistribution.R,sha256=r9XJwA-VZ2JP8TZEe4OS8Ah1b06yd58HrfAyPi4Be1M,16562
 biopipen/scripts/scrna/DimPlots.R,sha256=-mXOTMnpPxvR30XLjwcohFfFx7xTqWKKiICwJiD6yEo,1554
 biopipen/scripts/scrna/ExprImpution-alra.R,sha256=w3W1txJcdWg52-SETY2Z0lO7maDNfiMjBYIGN587YW0,843
 biopipen/scripts/scrna/ExprImpution-rmagic.R,sha256=jYIfqZpnvjKJkvItLnemPVtUApHBYQi1_L8rHVbEe1M,735
 biopipen/scripts/scrna/ExprImpution-scimpute.R,sha256=mg40qCUW7-nP5oHPvARq7dmtoahM0GRFWXQpum0BXVk,1082
 biopipen/scripts/scrna/ExprImpution.R,sha256=7768ezrr59xUZDXq8lO9jj2XhnkSsx-xxBmOD9_DO7c,313
-biopipen/scripts/scrna/MarkersFinder.R,sha256=60a4IVT2wmNa9ofuFRmJpw23XOh7v6uLyDy_VZkkoJY,26345
-biopipen/scripts/scrna/MetaMarkers.R,sha256=kgjk65EmewZ1uh8AoiENDLmNFe4lzmYvssJeINZ1JDo,12479
+biopipen/scripts/scrna/MarkersFinder.R,sha256=82XTp_8o4nBQ7yklqI5NoYW_4StLIdBj60M-AqhZnw4,26294
+biopipen/scripts/scrna/MetaMarkers.R,sha256=CyKRW81WFMv00W5mHbloaYrsq7B1niAFFtypf9Meezo,12411
 biopipen/scripts/scrna/ModuleScoreCalculator.R,sha256=JSHd-_-KiFqW8avCGxgU4T-C5BtDr2u0kwIvEu2lFIg,4188
-biopipen/scripts/scrna/RadarPlots.R,sha256=iR4JKtO2b3hGfqv_KAI7BR9tq02EAYfKeqp7tzAicKs,14808
+biopipen/scripts/scrna/RadarPlots.R,sha256=0A4rX0ngx3S0TDx-HMgLNIa6w5551vPvpTQZuZSb7J0,14757
 biopipen/scripts/scrna/SCImpute.R,sha256=dSJOHhmJ3x_72LBRXT72dbCti5oiB85CJ-OjWtqONbk,2958
-biopipen/scripts/scrna/ScFGSEA.R,sha256=E6Rx-0TjplP_nptDCE-LqMyipsOFaMU3hM7FjjFFFDY,7538
+biopipen/scripts/scrna/ScFGSEA.R,sha256=POyFFDmtDteVw0wN3ADkh85JhQvMPw3diW_Dg9NXMVM,7487
 biopipen/scripts/scrna/SeuratClusterStats-dimplots.R,sha256=pZKv1SnSNEGXDeE0_2VYp0GAikYitohW2FR5YGKjs8Q,2351
 biopipen/scripts/scrna/SeuratClusterStats-features.R,sha256=SaKTJloP1fttRXZQeb2ApX0ej7al13wOoEYkthSk13k,15489
-biopipen/scripts/scrna/SeuratClusterStats-hists.R,sha256=GTZfs1yOkuoMUM1Zb19i_My8B8b1Qtve8je55pU_w-g,5054
+biopipen/scripts/scrna/SeuratClusterStats-hists.R,sha256=YhuD-GePjJPSkR0iLRgV_hiGHD_bnOIKp-LB6GCwquo,5037
 biopipen/scripts/scrna/SeuratClusterStats-ngenes.R,sha256=GVKIXFNS_syCuSN8oxoBkjxxAeI5LdSxh-qLVkUsbDA,2146
 biopipen/scripts/scrna/SeuratClusterStats-stats.R,sha256=CE9989SaO75_KYEEVqivEbUoTcUOtiTkRGWLNtWzxI8,6450
-biopipen/scripts/scrna/SeuratClusterStats.R,sha256=PeXa3r2VTo0Q1rdXpSfAOIbYSJAcA8WUNM-tJkNUcPg,1284
+biopipen/scripts/scrna/SeuratClusterStats.R,sha256=ouWoj7Q644uG3MUlT23AES8f74g38-jPtPhINSvoUas,1267
 biopipen/scripts/scrna/SeuratClustering.R,sha256=kAvQq3RV86_KSv9NlUtUeQrPKkbhSsnv6Q4DoiTu8M0,6403
 biopipen/scripts/scrna/SeuratFilter.R,sha256=BrYK0MLdaTtQvInMaQsmOt7oH_hlks0M1zykkJtg2lM,509
 biopipen/scripts/scrna/SeuratLoading.R,sha256=ekWKnHIqtQb3kHVQiVymAHXXqiUxs6KKefjZKjaykmk,900
 biopipen/scripts/scrna/SeuratMap2Ref.R,sha256=tisYmoSaCX8Kl8y6euuuUroWdDsJ2NGI27J5AWr9Niw,4392
 biopipen/scripts/scrna/SeuratMetadataMutater.R,sha256=Pp4GsF3hZ6ZC2vroC3LSBmVa4B1p2L3hbh981yaAIeQ,1093
-biopipen/scripts/scrna/SeuratPreparing.R,sha256=KNrBPibqANhCsCctOGCOMPLTjDDd-iDRpMSCXkPuUQU,12994
+biopipen/scripts/scrna/SeuratPreparing.R,sha256=cgXWon2it6g4y-yrYk_zhivViOX8ZVf36u3wb9lKtj0,13133
 biopipen/scripts/scrna/SeuratSplit.R,sha256=vdK11V39_Uo_NaOh76QWCtxObGaEr5Ynxqq0hTiSvsU,754
 biopipen/scripts/scrna/SeuratSubClustering.R,sha256=L1SwKhNNKvsQGrcj0ZjScW9BLuvdO2pg7U48Ospsot8,6096
 biopipen/scripts/scrna/SeuratSubset.R,sha256=yVA11NVE2FSSw-DhxQcJRapns0tNNHdyDYi5epO6SKM,1776
 biopipen/scripts/scrna/SeuratTo10X.R,sha256=T2nJBTwOe12AIKC2FZsMSv6xx3s-67CYZokpz5wshqY,2679
-biopipen/scripts/scrna/TopExpressingGenes.R,sha256=m1BEpDfrzD9rJ2wEwA355GMqBq8Wfc_PzCjtLavB68Q,7754
+biopipen/scripts/scrna/TopExpressingGenes.R,sha256=J6tM54HStWkguUwFoIHRo_EAtCUKOuCBEfInzRfDYMQ,7703
 biopipen/scripts/scrna/Write10X.R,sha256=OMhXvJwvaH-aWsMpijKrvXQVabc1qUu5ZEwiLAhkDeY,285
 biopipen/scripts/scrna/sctype.R,sha256=NaUJkABwF5G1UVm1CCtcMbwLSj94Mo24mbYCKFqo1Bw,6524
-biopipen/scripts/scrna_metabolic_landscape/MetabolicFeatures.R,sha256=UbNoyspFB8166VgbkkFyTS6tkA-7ybylfvO50eqNnBU,4841
-biopipen/scripts/scrna_metabolic_landscape/MetabolicFeaturesIntraSubset.R,sha256=Kv69CkVm67K3w5EEcPQr7SlaWFpPwqTI0DVSUn62Rj0,5230
+biopipen/scripts/scrna_metabolic_landscape/MetabolicFeatures.R,sha256=b77yG5FeRse3bNfFgLIEYGHNZzydAn1OeyyR_n5Ju60,4790
+biopipen/scripts/scrna_metabolic_landscape/MetabolicFeaturesIntraSubset.R,sha256=ic8Fy8QqYDGh_izmvZVJ3KL66podg_CSF5ITL3FZsvo,5196
 biopipen/scripts/scrna_metabolic_landscape/MetabolicPathwayActivity.R,sha256=95DLX1Rz0tobOuDZ8V9YdGgO0KiNthhccoeeOK21tno,16216
-biopipen/scripts/scrna_metabolic_landscape/MetabolicPathwayHeterogeneity.R,sha256=0MIRrUjCtTDkUiL6GfKbtjC7bLlt50AGWFv4fl4e72k,9856
+biopipen/scripts/scrna_metabolic_landscape/MetabolicPathwayHeterogeneity.R,sha256=EpV4Wqj1dLqtaJyIJo6dN6vxeRH6elLibbHpK_N37bE,9822
 biopipen/scripts/tcgamaf/Maf2Vcf.py,sha256=Cxh7fiSNCxWDTfIJqZDOOnaSrw-85S_fH2U-PWY03hc,704
 biopipen/scripts/tcgamaf/MafAddChr.py,sha256=V10HMisl12O3ZfXuRmFNdy5p-3mr43WCvy0GHxSpwfA,494
 biopipen/scripts/tcgamaf/maf2vcf.pl,sha256=hJKcH-NbgWK6fmK7f3qex7ozJJl-PqCNPXqpwfcHwJg,22707
 biopipen/scripts/tcr/Attach2Seurat.R,sha256=C91TAh1cLSxWkdFPf84pbxlpTYMuWq_rduG4eiIkXZI,1345
-biopipen/scripts/tcr/CDR3AAPhyschem.R,sha256=rvXa9z7EY745Su0bzz0WpBd6QxXMGrzV_Eq7Ur9xb_g,16645
-biopipen/scripts/tcr/CloneResidency.R,sha256=C0XlkKOvWL0SeSCfhJ4BvqmxUQjx7Se6_PdzZk2HvgQ,21749
+biopipen/scripts/tcr/CDR3AAPhyschem.R,sha256=-0BS6cdt5GfQJphA3HlDgGjWr4XFF-7INLJyMBHQNAc,16628
+biopipen/scripts/tcr/CloneResidency.R,sha256=yE1Why0spX557D9aIQummoxXIdXBC6jt3YB5kIL1jt0,21647
 biopipen/scripts/tcr/CloneSizeQQPlot.R,sha256=5FPfWQjxTsv59KSDQaDWj3C95zPQMngKG7qOf95NEzI,4527
 biopipen/scripts/tcr/GIANA/GIANA.py,sha256=0qLhgCWxT8K-4JvORA03CzBPTT5pd4Di5B_DgrHXbFA,47198
 biopipen/scripts/tcr/GIANA/GIANA4.py,sha256=Z7Q3cUr1Pvmy4CFADN0P7i9g1-HbzWROMqk5HvL_F1Q,45762
@@ -161,20 +164,20 @@ biopipen/scripts/tcr/GIANA/Imgt_Human_TRBV.fasta,sha256=XUwDPXJxVH5O9Q0gCL6EILKX
 biopipen/scripts/tcr/GIANA/query.py,sha256=5NWSEDNrJomMt48tzLGpRwJTZB0zQxvMVTilyG8osX8,7298
 biopipen/scripts/tcr/Immunarch-basic.R,sha256=g64RXmiPw73vbnwrKoRaNs1d3O6mRV9uhcAFkBV0g3U,3121
 biopipen/scripts/tcr/Immunarch-clonality.R,sha256=48rbPCWka4eNEy-fjM0BlKDkMYdG2zlB8Sly1B4xdUI,3858
-biopipen/scripts/tcr/Immunarch-diversity.R,sha256=BDIhh1lbiaAEm-AQO0nb0mBmtQBA9yG_AP6N6zW8QOc,28330
+biopipen/scripts/tcr/Immunarch-diversity.R,sha256=C1qx50K2eIZ2c590jumGe1BpMmgohlaP-wVkmwzHa2E,28296
 biopipen/scripts/tcr/Immunarch-geneusage.R,sha256=c0C8-KtKI2q6O9xZ9f5COefQbPlshT2hz1f36qpnW34,6817
-biopipen/scripts/tcr/Immunarch-kmer.R,sha256=4BUDiclarawetjUdBmrCTfHZOTjRpE5TDxDGI_CtYds,6229
+biopipen/scripts/tcr/Immunarch-kmer.R,sha256=ttDKwPTde_DN2xekRF2LBrc_o8wQjxXa49VAzmUzQl8,6212
 biopipen/scripts/tcr/Immunarch-overlap.R,sha256=GVt-qJPtd6NEe5njAqNStf2AP6pLkv7Ittw0YT_qdNY,7465
-biopipen/scripts/tcr/Immunarch-spectratyping.R,sha256=Ii1E3QM0ZyUlmGUosHFQX1HXm93sJ-wqWiMKSbUiTgo,2916
+biopipen/scripts/tcr/Immunarch-spectratyping.R,sha256=jm8qI_A6HBgr1RT6-kq13z46_3aXHwz_LbQk4Ggeq9Y,2882
 biopipen/scripts/tcr/Immunarch-tracking.R,sha256=0tiywpGhd3H0REp4xrhOlkWzJM4ntrQjrVesHbEWT40,4374
-biopipen/scripts/tcr/Immunarch-vjjunc.R,sha256=gnLYaUS9uhwTcTqqba2ZmsepuBWcAw3pr1q3FUFQJFs,4451
-biopipen/scripts/tcr/Immunarch.R,sha256=DjSLeEkB01inDS6iSCnQqQhRBJ20IbZn55nefL0gmQU,3054
+biopipen/scripts/tcr/Immunarch-vjjunc.R,sha256=vSVPdRPSlpAfuX4ynA2UjfUC0_61yMAaPieJM-5J4sw,4434
+biopipen/scripts/tcr/Immunarch.R,sha256=FVPK5qu6twSSJJNDCAZO4wNned6uIOBZHgBNeIgLmKw,3037
 biopipen/scripts/tcr/Immunarch2VDJtools.R,sha256=QB9ILGbnsfoWaRANK6ceb14wpSWy8F1V1EdEmfIqiks,706
 biopipen/scripts/tcr/ImmunarchFilter.R,sha256=o25O36FwH_0w6F8DFQ0SfpcwDzlzaGefXqr9ESrvb4k,3974
 biopipen/scripts/tcr/ImmunarchLoading.R,sha256=l_l-gojiCKI_MWgIUe2zG5boVtNipBv4rACRJEcrnFE,5734
 biopipen/scripts/tcr/ImmunarchSplitIdents.R,sha256=FGCeGV0uSmFU91lKkldUAeV4A2m3hHw5X4GNi8ffGzI,1873
 biopipen/scripts/tcr/SampleDiversity.R,sha256=jQ1OU3b8vswD8tZhLt3fkcqJKrl2bhQX0giHM2rXz3Y,2643
-biopipen/scripts/tcr/TCRClusterStats.R,sha256=3YxIfsTBbFFI6fBTU3gM60bGuVv52PmL7bs16_WciGw,12089
+biopipen/scripts/tcr/TCRClusterStats.R,sha256=D7q1svXQxl1uOya8bePvR9e6NJXjCjXbPsXnEPTWdlE,12004
 biopipen/scripts/tcr/TCRClustering.R,sha256=-BWbDqvDBEpfVaxrVvzVHK5bm6FCOFmGHydg1c3EgAM,8747
 biopipen/scripts/tcr/TESSA.R,sha256=bfOixWLZy8yi0MzXncP67KjtCukwXEzsK5fCdMzB5VM,6822
 biopipen/scripts/tcr/TESSA_source/Atchley_factors.csv,sha256=SumqDOqP67P54uM7Cuc5_O_rySTWcGo7eX3psMSPX9s,763
@@ -210,9 +213,9 @@ biopipen/utils/caching.R,sha256=qANQqH8p-VpvD8V4VSoqSfp0TFr4esujC7x3OFZsJMw,1687
 biopipen/utils/common_docstrs.py,sha256=ro9dUXeHMXBaqb-hTafwrG6xW5IOBEeiUM2_REjFoCo,5842
 biopipen/utils/gene.R,sha256=BzAwlLA8hO12vF-3t6IwEuTEeLa_jBll4zm_5qe3qoE,1243
 biopipen/utils/gene.py,sha256=qE_BqTayrJWxRdniffhcz6OhZcw9GUoOrj2EtFWH9Gw,2246
-biopipen/utils/gsea.R,sha256=rB_wHUBrn9SihMy_ICSRSgj-e48SGm_ZYs1oZtFgroA,7284
+biopipen/utils/gsea.R,sha256=UMQOlWGstQTOBScvy1wIzrB7I3CE28Xo2v1sy4lmJ-M,7549
 biopipen/utils/io.R,sha256=jIYdqdn0iRWfQYAZa5CjXi3fikqmYvPPLIXhobRe8sw,537
-biopipen/utils/misc.R,sha256=nkjiAsEsilq0AeiKRDNqrhTx-1Grqg-rFlkjOEOEDYg,5224
+biopipen/utils/misc.R,sha256=sB2cM-sJn9H90FAZD89uNUYPU5k_6BT8d_nG33fI7d4,5520
 biopipen/utils/misc.py,sha256=Pmh3CBiKJ3vC_RqorfOfRAvTVKXrGDJT8DMLfYbTivs,3055
 biopipen/utils/mutate_helpers.R,sha256=Bqy6Oi4rrPEPJw0Jq32bVAwwBfZv7JJL9jFcK5x-cek,17649
 biopipen/utils/plot.R,sha256=pzl37PomNeUZPxohHZ2w93j3Fc4T0Qrc62FF-9MTKdw,4417
@@ -220,7 +223,7 @@ biopipen/utils/reference.py,sha256=6bPSwQa-GiDfr7xLR9a5T64Ey40y24yn3QfQ5wDFZkU,4
 biopipen/utils/rnaseq.R,sha256=Ro2B2dG-Z2oVaT5tkwp9RHBz4dp_RF-JcizlM5GYXFs,1298
 biopipen/utils/single_cell.R,sha256=bKduqOQjSC8BtZJuwfUShR49omoEMbB57n3Gi6dYlqA,4147
 biopipen/utils/vcf.py,sha256=ajXs0M_QghEctlvUlSRjWQIABVF02wPdYd-0LP4mIsU,9377
-biopipen-0.24.1.dist-info/METADATA,sha256=ZzmqmOreZpC57jB-7zZEeu5vFfpKlR9LOhUCa5kauW4,932
-biopipen-0.24.1.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
-biopipen-0.24.1.dist-info/entry_points.txt,sha256=16Apdku3RFwghe1nb0JR7eVo4IzLae6hCWjU1VxYUn0,525
-biopipen-0.24.1.dist-info/RECORD,,
+biopipen-0.24.2.dist-info/METADATA,sha256=it_bxwfsavzwDtAkIWf2RxQIcmbr0810VrAEonBE9DQ,932
+biopipen-0.24.2.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
+biopipen-0.24.2.dist-info/entry_points.txt,sha256=O7K7PHoJi-4tERYpqmryuTvrFNEsptMpqliZzB2oBQQ,577
+biopipen-0.24.2.dist-info/RECORD,,

{biopipen-0.24.1.dist-info → biopipen-0.24.2.dist-info}/entry_points.txt

@@ -3,6 +3,7 @@ bam=biopipen.ns.bam
 bcftools=biopipen.ns.bcftools
 bed=biopipen.ns.bed
 cellranger=biopipen.ns.cellranger
+cellranger_pipeline=biopipen.ns.cellranger_pipeline
 cnv=biopipen.ns.cnv
 cnvkit=biopipen.ns.cnvkit
 cnvkit_pipeline=biopipen.ns.cnvkit_pipeline