HOME  ›   pipelines
If your question is not answered here, please email us at:  ${email.software}

Cell Ranger


Loupe

10x Genomics
Chromium Single Cell Gene Expression

Comparing Targeted and Parent Samples with cellranger targeted-compare

Table of Contents

Overview

Targeted Gene Expression libraries are derived by performing hybrid capture on a Single Cell Gene Expression library. This starting input library, referred to as the parent library, therefore contains the molecules present in the targeted library. While it is not necessary to sequence this parent sample when analyzing Targeted Gene Expression data, if both the targeted and parent samples have been sequenced, the targeted-compare pipeline may be used to directly compare the two datasets. This pipeline can be used to assess targeting performance with greater accuracy than when only the targeted data are known. This comparative analysis provides quality control metrics to verify the extent to which targeted genes were enriched and the parent sample data were recovered.

Requirements

First, run cellranger count on the parent sample (see Single-Library Analysis with Cell Ranger). Then run cellranger count --target-panel target_panel.csv (see Targeted Gene Expression Analysis). The targeted and parent samples must be run using the same transcriptome reference. For example:

cd /home/jdoe/runs
cellranger count --id=targeted_sample --target-panel target_panel.csv ...
... wait for pipeline to finish ...

cellranger count --id=parent_sample ... ... wait for pipeline to finish ...

After the above two analyses are complete, run targeted-compare as outlined in the following section.

Command line interface

These are the command line arguments (also accessible by running cellranger targeted-compare --help):

ArgumentDescription
--id=IDA unique run ID string: e.g. TGTCMP123
--targeted=MOL_INFO_H5Path to a molecule_info.h5 file from a Targeted Gene Expression Analysis.
--parent=MOL_INFO_H5Path to a molecule_info.h5 file from a Single Cell Gene Expression Analysis.
--target-panel=CSVTarget panel CSV file from the targeting experiment. Must be the same target set used as the --target-panel flag used in the targeted gene expression analysis.

Resuming the above example, run cellranger targeted-compare as follows:

cellranger targeted-compare --id=TGTCMP123
    --targeted=targeted_sample/outs/molecule_info.h5
    --parent=parent_sample/outs/molecule_info.h5
    --target-panel target_panel.csv

The pipeline will begin to run, creating a new folder named with the ID you specified (e.g. /home/jdoe/runs/TGTCMP123) for its output. If this folder already exists, cellranger will assume it is from an incomplete targeted-compare run and attempt to resume running it.

Output files

A successful run should conclude with a message similar to this:

2020-06-29 12:25:58 [runtime] (join_complete)   ID.TGTCMP123.TARGETED_COMPARE_CS.TARGETED_COMPARE.SUMMARIZE_COMPARISON
 
Outputs:
- Sample comparison summary HTML: /home/jdoe/runs/TGTCMP123/outs/web_summary.html
- Sample comparison summary CSV:  /home/jdoe/runs/TGTCMP123/outs/metrics_summary.csv
- Per barcode metrics CSV:        /home/jdoe/runs/TGTCMP123/outs/barcode_summary.csv
- Per feature metrics CSV:        /home/jdoe/runs/TGTCMP123/outs/feature_summary.csv
 
Pipestance completed successfully!

The output folder (e.g. TGTCMP123/runs/outs) will contain the following files:

File NameDescription
web_summary.htmlRun summary metrics and charts in HTML format.
metrics_summary.csvRun summary metrics in CSV format.
barcode_summary.csvPer-barcode data in the Targeted and Parent samples.
feature_summary.csvPer-gene data in the Targeted and Parent samples.

See Targeted-Compare Outputs for more information on these output files.