This is a test version of Biostars. For the public version, visit https://www.biostars.org.
scrna pipeline

I have 2 fastq files and did their quality check through fastqc tool. Which other steps do I need to do to complete the pipeline?

scrna-seq

What have you tried? There is documentation on complete pipelines available online. Did you search show you any useful link?

Here's a list of steps you might want to try assuming library prep was done with 10x Genomics:

  1. Run cellranger to get to counts files (matrix, barcodes, genes).
  2. Load the data using scanpy.
  3. Compute QC metrics and filter the cells and genes.
  4. Run normalization, log transformation (log1p), scaling, PCA, neighborhood graph, and UMAP (in this order). I assumed no batch effect correction is needed.
  5. Run a clustering algorithm such as leiden.
  6. Based on your expected cell composition, you can plot marker genes to identify cell types. You can also color your UMAP by gene list scores to identify cell types.
  7. Perform additional analyses such as cell composition with scCODA, pseudotime ordering, cell-cell communication, perturbation modelling, etc.

I hope this gives you an idea of where you need to look.

0 answers

No answers yet.

Log in to answer this question.