• Product manuals

Browse the manual

• Introduction
  • The concept of CLC Single Cell Analysis Module
  • Contact information
  • System requirements and installation
    • System requirements
    • Installation of modules
    • Uninstalling modules
  • Licensing modules
    • Server License
  • The QIAGEN Cell Ontology
  • Reference data management
    • The Reference Data Manager
  • Running tools and workflows
• Data import
  • Import Expression Matrix
  • Import Cell Annotations
  • Import Cell Clusters
  • Import Cell Clonotypes
• Data export
  • Export Cell Ranger HDF5 expression matrix
  • Export Loom expression matrix
  • Export MEX expression matrix
• Creating a count matrix from FASTQ
  • Annotate Reads with Cell and UMI
    • Setting the sample name
    • Interpreting the output of Annotate Reads with Cell and UMI
  • Single Cell RNA-Seq Analysis
    • The Single Cell RNA-Seq Analysis report
    • The Single Cell RNA-Seq Analysis algorithm
• Cell preparation
  • QC for Single Cell
    • Empty droplets filter
    • Count-based and extra-chromosomal filters
    • Doublets filter
    • Interpreting the output of QC for Single Cell
    • Choosing barcodes to retain
    • The cell calling algorithm
    • The doublet calling algorithm
  • Normalize Single Cell Data
    • When is batch correction appropriate?
    • Interpreting the output of Normalize Single Cell Data
    • The Normalize Single Cell Data algorithm
  • Combine Cell Annotations
  • Combine Cell Clusters
  • Convert Metadata to Cell Annotations
  • Update to Common Sample Name
• Noise reduction through feature selection and PCA
  • Feature selection and PCA
    • Calculation of estimated biological variation
• Cell Annotation
  • Predict Cell Types
  • Train Cell Type Classifier
    • Features used for training and prediction
    • Interpreting the output of Train Cell Type Classifier
    • SVMs for cell type classification
  • Cluster Single Cell Data
    • The Cluster Single Cell Data algorithm
  • Create Heat Map for Cell Abundance
    • Interpreting the output of Create Heat Map for Cell Abundance
• Dimensionality reduction
  • Create UMAP Plot
  • Create tSNE Plot
  • Add Information to Plot
• Expression Analysis
  • Differential Expression for Single Cell
    • Interpreting the output of Differential Expression for Single Cell
    • The differential expression algorithm
  • Create Expression Plot
    • The Heat Map output of Create Expression Plot
    • The Dot Plot output of Create Expression Plot
    • The Violin Plot output of Create Expression Plot
• Immune Repertoire
  • Single Cell Immune Repertoire Analysis
    • The report output from Single Cell Immune Repertoire Analysis
    • The Cell Clonotypes element
    • The clonotype identification algorithm
  • Filter Cell Clonotypes
  • Combine Cell Clonotypes
  • Compare Single Cell Immune Repertoires
    • Interpreting the output of Compare Single Cell Immune Repertoires
  • Convert Clonotypes to Cell Annotations
• Single cell workflows from reads
  • Expression Analysis from Reads
    • Configuring the batch units
    • Output from Expression Analysis from Reads
  • Immune Repertoire Analysis from Reads (10xVDJ)
    • Output from Immune Repertoire Analysis from Reads (10xVDJ)
  • Immune Repertoire and Expression Analysis from Reads (10xVDJ)
    • Configuring the batch units with multiple input types
    • Output from Immune Repertoire and Expression Analysis from Reads (10xVDJ)
• Single cell workflows from imported data
  • Expression Analysis from Matrix
    • Output from Expression Analysis from Matrix
  • Immune Repertoire and Expression Analysis from Clonotypes and Matrix
    • Output from Immune Repertoire and Expression Analysis from Clonotypes and Matrix
• UMAP and tSNE plot functionality
  • Manual Annotation
  • Create Subset
  • Extract to Table
  • Launching of Create Expression Plot
  • Launching of Create Heat Map for Cell Abundance
  • Launching of Differential Expression for Single Cell
• Licensing requirements for the CLC Single Cell Analysis Module
  • Licensing modules on a Workbench
    • Request an evaluation license
    • Download a license using a license order ID
    • Import a license from a file
    • Configure License Server connection
    • Download a static license on a non-networked computer
  • Licensing Server Extensions on a CLC Server
    • Download a static license on a non-networked machine
• Bibliography

The cell calling algorithm

The algorithm for cell calling is based on EmptyDrops [Lun et al., 2019]:

• The ambient RNA profile is estimated from the barcodes with low number of reads. The expressions from these barcodes are added together and a proportion vector for the ambient profile is obtained using the Good Turing algorithm [Gale and Sampson, 1995].
• If Specify minimum number of reads for barcodes to be retained is disabled (see Empty droplets filter), the knee is automatically estimated from the log-log rank plot (figure 5.5). The knee estimation does not follow that from EmptyDrops. Instead, the knee is identified on a smoothed data using an adaptation of the [Satopaa et al., 2011] algorithm implemented in https://github.com/mariolpantunes/ml.
• Barcodes with fewer number of reads than the automatically inferred knee or a manually specified threshold, are tested for significant deviations from the ambient profile. For each barcode, the probability of obtaining its expression profile from the ambient is calculated. A p-value is obtained from the probabilities of ambient simulated barcodes containing the same total number of reads.
• FDR correction is applied to the p-values for barcodes that are not part of the ambient profile.
• Barcodes with FDR-corrected p-values below the provided value in FDR threshold (see Empty droplets filter) are retained as non-empty droplets.

---