Single Cell Analysis II
Applications of Single-Cell RNA Sequencing
Single-Cell RNA Sequencing Analysis Pipeline
Expression Matrix: Combines genes with cells based on barcodes and unique molecular identifiers.
Filtering Cells (Quality Control):
Removing dead cells, often indicated by high levels of mitochondrial genes.
Removing multiplets (cells clumped together).
Normalization: Adjusting for relative gene expression, especially when comparing different conditions or cells.
Variable Gene Identification: Determining which genes are expressed differently by individual cells.
Dimensionality Reduction: Using methods like TSNE plots or UMAP plots to visualize cell similarity.
Cells close to each other are more similar; distant cells are more different.
Clustering: Grouping cells based on their similarity in gene expression.
Exploring Known Gene Markers: Annotating clusters based on known marker genes (e.g., CD3E for T cells).
Differential Gene Expression Analysis: Identifying genes that are differentially expressed between clusters.
Cell Type Assignment: Assigning cell types to clusters based on gene expression profiles.
Functional Annotation: Identifying the functions associated with differentially expressed genes (e.g., pro-inflammatory state).
Pseudo Time Analysis: Tracking cellular differentiation or maturation by analyzing gene expression trajectories.
Following gene expression changes as cells differentiate (e.g., from stem cell to different cell types).
Superimposing disease conditions to observe changes in the transition.
Outputs and Strengths of Single-Cell RNA Sequencing
Cellular Resolution: Ability to identify expression profiles of individual cells and subclusters.
Heterogeneity: Revealing great heterogeneity even within one cell type.
Rare Cell Type Detection: Identifying rare cell types or subpopulations.
Dot Plots
Used to identify cells based on differentially expressed genes.
Show the top 5 most highly expressed genes for each cell type.
Enable identification of cell types based on unique gene expression patterns.
Example: Epicardial cells have a unique set of top 5 genes not expressed by other cells.
TSNE Plots
Visualize different clusters of cells.
Demonstrate cellular heterogeneity within populations like fibroblasts.
Example: In hypertensive heart disease, a specific fibroblast population appears that is not present in healthy hearts.
This population expresses new genes and may represent a new drug target for treating hypertensive heart failure.
Identification of New Cell Types
Single-cell RNA sequencing can uncover new cell states or types that are missed by bulk RNA sequencing.
Bulk RNA sequencing provides an average of gene expression, masking subtle changes in specific cell types.
*Example: A study in aorta of healthy and hypertensive mice revealed novel fibroblast cluster upregulated in hypertensive vessels
Quantifying Relative Abundance
Single-cell data allows quantification of cell type abundance.
Example: Stack graphs showing an increase in a novel fibroblast type (Fibro 1) in hypertensive vessels compared to healthy vessels.
Probing for Unique Genes
Dot plots can identify unique genes expressed by specific cell clusters.
Focusing on the top 10 genes expressed by a novel population (e.g., Fibro 1) helps identify its characteristics.
These genes may be associated with disease processes like fibrosis (upregulating fibrotic genes).
Validation
Using top 10 expressed genes to identify markers for specific cells.
Example: CTHRC1 is highly expressed in Fibro 1 but not in other fibroblasts, making it a potential marker.
Feature Plots
Highlight the expression of a specific gene.
Example: A feature plot showing CTHRC1 expression predominantly in the Fibro 1 population.
Challenging Current Paradigms
Single-cell RNA sequencing can challenge existing beliefs about cell function in disease.
Example: Myofibroblasts are traditionally thought to be the main cells driving collagen deposition in fibrotic diseases, but single-cell data suggests another cell type (Fibro CTHRC1) is also involved.
Functional Annotation
Gene ontology analysis and gene set enrichment analysis can identify functions or cellular pathways associated with specific cell types.
Pathways are associated with clusters of genes.
Example: Collagen or extracellular matrix molecules form a pathway associated with extracellular matrix deposition or fibrosis.
Generating gene ontology graphs:
The X-axis shows the number of genes associated with each pathway.
Graph visualizes genes that are upregulated in disease.
Single-cell RNA sequencing enables multiple comparisons, such as comparing a specific cell to an entire healthy organ.
Trajectory Analysis
Orders cells along a developmental pathway.
Reveals insights not obtainable with bulk RNA sequencing.
Identifies directionality of differentiating cells through genes that are switched on as they head toward a different cell type.
Limitations: Requires determination of a starting point (origin), which can greatly influence analysis.
Selecting genes is often done by reading the literature to look at genes turned on by certain cells during differentiation or maturation processes.
Cell-to-Cell Interaction
Bulk RNA sequencing only shows what ligands are up regulated in a tissue.
Single-cell RNA sequencing can describe what cells the ligands are coming from.
Ligand-receptor interactions:
Ligands (e.g., cytokines, cellular mediators) bind to receptors on other cells, initiating cell signaling.
Interactome:
Visualizes cell-to-cell communication networks.
Arrowheads show the direction of signals between cells, with colors indicating cell types.
Analyzing outgoing signals:
Allows for isolation of one cell type.
Allows you to see what cells the isolated cell is talking to.
Limitations: Reed depth can influence detection of ligands.
Outputs as Hypothesis-Generating Data
Major outputs from single-cell RNA sequencing are hypothesis generating data.
Identify novel cell types and molecular biomarkers.
Studying potential drug targets:
Genes and potentially target proteins can be studied as potential drug targets.
Ligands and receptors can be studied as potential drug targets.
Limitations of Single-Cell RNA Sequencing
Small Cell Populations: If the population of cells is small, or a unique cell type that's not at a high level, an imbalance in sequencing can occur.
Technical Noise: Noisy data can influence the data in the sense that amplification biases or even dropouts of signals can occur.
Cellular Heterogeneity: Can be problematic with imbalanced populations.
Lowly Expressed Genes: Unlikely to be detected with single-cell RNAC.
Computational Complexity: Requires specialized tools, software packages, computing, coding, and bioinformatic skills.
Importance of Validation
Artefacts In Silica: High risk of artefacts with analysis using a computer.
Relative Gene Expression: Lowly expressed genes may not be picked up, or really abundant cell information is picked up.
Validation Needed: Relative gene expression should be validated with validation experiments.
Techniques for validation:
QPCR.
Looking at protein after validation.
mRNA Expression: Not protein.
Location matters for validation:
Important to validate localization.
Think about anatomical barriers with communication.