bioc_sc
WDL pipeline for standard Bioconductor/OSCA single-cell RNA-seq post-processing: load matrix, QC filter, normalize, select highly variable genes, run dimensionality reduction, cluster, identify marker genes, and annotate cell types
WDL Version
1.0
Nested Inputs Not Allowed
RUN WITH
sprocket run --target bioc_sc ww-bioc-sc.wdlww-bioc-sc.wdl [INPUTS]...
authorurlhttps://raw.githubusercontent.com/getwilds/wilds-wdl-library/refs/heads/main/pipelines/ww-bioc-sc/ww-bioc-sc.wdlInputs
Required Inputs
Name
Type
Description
samplesArray[SingleCellSample]List of single-cell sample objects, each with a name and either a pre-filtered 1...
List of single-cell sample objects, each with a name and either a pre-filtered 10x feature-barcode matrix (h5_matrix) or a raw, unfiltered one (raw_h5_matrix). Exactly one of the two must be provided per sample; if raw_h5_matrix is given, real cells are called with emptyDrops before QC
Other Inputs
Name
Type
Default
Description
mito_patternString"^MT-"Regex pattern identifying mitochondrial gene symbols, passed to ww-scran
Regex pattern identifying mitochondrial gene symbols, passed to ww-scran
nmadsFloat3.0Number of median absolute deviations from the median used to flag low-quality cells, passed to ww-scran
Number of median absolute deviations from the median used to flag low-quality cells, passed to ww-scran
n_hvgsInt2000Number of top highly variable genes to select in ww-scran and use for PCA in ww-scater
Number of top highly variable genes to select in ww-scran and use for PCA in ww-scater
n_pcsInt50Number of principal components to compute in ww-scater
Number of principal components to compute in ww-scater
reference_datasetString"HumanPrimaryCellAtlasData"Name of the celldex reference dataset function to fetch for cell-type annotation, passed to ww-singler
Name of the celldex reference dataset function to fetch for cell-type annotation, passed to ww-singler
reference_ensemblBooleantrueFetch the celldex reference with Ensembl gene IDs instead of gene symbols, passed to ww-singler
Fetch the celldex reference with Ensembl gene IDs instead of gene symbols, passed to ww-singler
label_columnString"label.main"Column name in the reference object's colData containing cell-type labels, passed to ww-singler
Column name in the reference object's colData containing cell-type labels, passed to ww-singler
cpu_coresInt2Number of CPU cores allocated per task
Number of CPU cores allocated per task
memory_gbInt8Memory allocated in GB for ww-dropletutils, ww-scran, and ww-scater tasks
Memory allocated in GB for ww-dropletutils, ww-scran, and ww-scater tasks
singler_memory_gbInt16Memory allocated in GB for ww-singler tasks (higher default since celldex reference loading and clustering need more headroom)
Memory allocated in GB for ww-singler tasks (higher default since celldex reference loading and clustering need more headroom)
Outputs
Name
Type
Expression
Description
sce_rdsArray[File]sample_sce_rdsArray of loaded SingleCellExperiment RDS objects for each sample, prior to QC/normalization
Array of loaded SingleCellExperiment RDS objects for each sample, prior to QC/normalization
empty_drops_csvArray[File?]empty_drops_filter.empty_drops_csvArray of per-barcode emptyDrops statistics for samples provided as a raw matrix (empty for samples provided pre-filtered)
Array of per-barcode emptyDrops statistics for samples provided as a raw matrix (empty for samples provided pre-filtered)
barcode_rank_plotArray[File?]empty_drops_filter.barcode_rank_pdfArray of barcode-rank QC plots for samples provided as a raw matrix (empty for samples provided pre-filtered)
Array of barcode-rank QC plots for samples provided as a raw matrix (empty for samples provided pre-filtered)
scran_sce_objectArray[File]run_scran.sce_objectArray of normalized SingleCellExperiment RDS objects with size factors and PCA for each sample
Array of normalized SingleCellExperiment RDS objects with size factors and PCA for each sample
scran_qc_plotArray[File]run_scran.qc_plotArray of per-cell QC scatter plots from scran for each sample
Array of per-cell QC scatter plots from scran for each sample
scran_size_factor_plotArray[File]run_scran.size_factor_plotArray of size factor histograms for each sample
Array of size factor histograms for each sample
scran_mean_variance_plotArray[File]run_scran.mean_variance_plotArray of mean-variance trend plots for each sample
Array of mean-variance trend plots for each sample
scran_hvg_tableArray[File]run_scran.hvg_tableArray of CSVs of genes ranked by biological variance for each sample
Array of CSVs of genes ranked by biological variance for each sample
scater_sce_objectArray[File]run_scater.sce_objectArray of SingleCellExperiment RDS objects with PCA and UMAP embeddings for each sample
Array of SingleCellExperiment RDS objects with PCA and UMAP embeddings for each sample
scater_pca_plotArray[File]run_scater.pca_plotArray of PCA scatter plots for each sample
Array of PCA scatter plots for each sample
scater_umap_plotArray[File]run_scater.umap_plotArray of UMAP scatter plots for each sample
Array of UMAP scatter plots for each sample
scater_qc_plotArray[File]run_scater.qc_plotArray of per-cell QC scatter plots from scater for each sample
Array of per-cell QC scatter plots from scater for each sample
singler_sce_objectArray[File]run_singler.sce_objectArray of final SingleCellExperiment RDS objects with cluster assignments and predicted cell types for each sample
Array of final SingleCellExperiment RDS objects with cluster assignments and predicted cell types for each sample
singler_cluster_tableArray[File]run_singler.cluster_tableArray of CSVs mapping cell barcodes to clusters for each sample
Array of CSVs mapping cell barcodes to clusters for each sample
singler_marker_tableArray[File]run_singler.marker_tableArray of CSVs of top marker genes per cluster for each sample
Array of CSVs of top marker genes per cluster for each sample
singler_prediction_tableArray[File]run_singler.prediction_tableArray of CSVs of SingleR cell-type predictions for each sample
Array of CSVs of SingleR cell-type predictions for each sample