Step-by-step workflows, techniques, and research-ready procedures
This collection provides standardized experimental protocols and maintenance procedures for C. elegans research, complemented by in silico computational workflows for bioinformatics analysis.
Although Caenorhabditis elegans (C. elegans) can be maintained axenically, this method is difficult and results in slow growth. In most laboratories, C. elegans is maintained monoxenically using E. coli strain OP50 as a food source (Brenner, 1974).
OP50 is a uracil auxotroph, which limits its growth on NGM plates. A thin bacterial lawn is ideal for clear observation and efficient mating.
Step 1: Obtain a starter culture of E. coli OP50 from the CGC or recover it from existing worm plates.
Step 2: Streak the culture onto an LB agar plate:
Step 3: Using a single colony, aseptically inoculate L Broth:
Step 4: Incubate the inoculated broth overnight at 37 °C.
Step 5: Store both streak plates and liquid cultures at 4 °C; they remain usable for several months.
Caenorhabditis elegans (C. elegans) is routinely maintained on Nematode Growth Medium (NGM) agar (Brenner, 1974), poured aseptically into sterile petri plates.
| Plate | Size Diameter | Common Use |
|---|---|---|
| Small | 35 mm | Matings or assays using costly reagents |
| Medium | 60 mm | General strain maintenance |
| Large | 100 mm | Large cultures or mutant screens |
A peristaltic pump (e.g., Wheaton Unispense) is recommended for even dispensing, ensuring uniform agar thickness across plates.
If needed, drugs such as levamisole, streptomycin, or nystatin may be added to the molten NGM just before pouring (Caldicott et al., 1994).
Step 1: Combine in a 2 L Erlenmeyer flask:
Cover with aluminum foil and autoclave for 50 minutes.
Step 2: Cool flask in a 55 °C water bath for 15 minutes.
Step 3: Add (sterile):
Swirl gently to mix.
Step 4: Dispense aseptically into petri plates using the peristaltic pump (fill plates about ⅔ full).
Step 5: Leave plates at room temperature for 2–3 days to check for contamination and allow moisture to evaporate.
Step 6: Store plates in an air-tight container at room temperature; they remain usable for several weeks.
Using sterile technique, apply E. coli OP50 culture onto solidified NGM plates.
Step 1: Pipette:
Step 2: Optionally, spread the drop using a pipette tip or sterile glass rod to create a uniform lawn (avoid edges).
Step 3: Allow plates to dry overnight at room temperature or incubate at 37 °C for 8 hours (then cool before use).
Step 4: Store seeded plates in an air-tight container; they remain usable for 2–3 weeks.
WormBook: The Online Review of C. elegans Biology WormBook
Caenorhabditis elegans (C. elegans) is transparent and can be visualized using a dissecting stereomicroscope with a transmitted light source. Common models include the Wild Leitz M5A and Zeiss SV6, equipped with 10× eyepieces and objectives ranging from 0.6× to 5× (total magnification of 6×–50×).
For precision transfers or when maintaining genetic stocks.
Tip 1: Flatten or slightly bend the wire tip to form a small hook—avoid sharp points.
Tip 2: To pick a worm, gently touch the side of the worm with the wire tip and lift.
Tip 3: Alternatively, pick up a small blob of E. coli OP50 on the wire first; worms will stick to the bacteria.
Tip 4: To transfer, gently touch the picker to the new plate's surface and wait for the worm to crawl off.
A quick and convenient way to transfer large numbers of worms.
Step 1: Use a sterilized scalpel or spatula to cut and move a small chunk of agar (containing worms) from an old plate to a fresh NGM plate.
Step 2: The worms will crawl out of the agar chunk and spread onto the new bacterial lawn.
WormBook: The Online Review of C. elegans Biology WormBook
Synchronous C. elegans cultures ensure that all animals are at the same developmental stage, improving reproducibility of developmental, genetic, and drug response experiments. This protocol uses an alkaline bleach treatment to dissolve gravid adults, leaving eggs intact, which are then hatched in the absence of food to arrest uniformly at the L1 stage. The synchronized L1 larvae can subsequently be transferred to plates or liquid medium for downstream assays.
STEP 1: Grow gravid adult worms on well seeded NGM plates until many adults are full of eggs. Wash worms from 3–6 plates into a 15 ml conical tube with M9 buffer and adjust the volume to 10–15 ml, avoiding large agar chunks or thick bacterial clumps.
STEP 2: Centrifuge at 1500–2000 rpm for 1–2 minutes. Carefully aspirate the supernatant, leaving about 0.5–1 ml of liquid above the worm pellet without disturbing it.
STEP 3: Prepare fresh bleaching solution (for example, 2 parts household bleach, 1 part 5 N NaOH, 2 parts ddH₂O). Add 5–8 ml of this solution to the pellet, cap the tube, and place it on a rocker or nutator, inverting gently for 5–10 minutes until most adult worm bodies are dissolved but eggs remain intact; do not exceed 10 minutes to avoid damaging embryos.
STEP 4: Centrifuge at 1500–2000 rpm for 1–2 minutes to pellet the eggs. Aspirate the bleach solution completely, add 10 ml sterile M9, gently resuspend the pellet by inversion, and spin again. Repeat this M9 wash 2–3 times to remove all traces of bleach.
STEP 5: After the final wash, resuspend the egg pellet in 5–10 ml sterile M9 (or S medium) in a vented culture tube or flask. Incubate at 20 °C with gentle shaking for 12–24 hours without food so embryos hatch and arrest synchronously at the L1 stage.
STEP 6: Estimate L1 concentration by counting worms in a small aliquot and extrapolating. Transfer the desired number of synchronized L1 larvae to fresh OP50 seeded NGM plates or liquid cultures to begin experiments from a defined developmental stage.
The Caenorhabditis elegans egg hatching assay quantifies how effectively drugs or compounds inhibit nematode egg hatching compared with untreated controls. Isolated eggs from synchronized gravid adults are exposed to graded drug concentrations, allowing calculation of EC50 or EC90 values and evaluation of time dependent effects on egg development. This in vitro method is cost effective, avoids the use of parasitic species or host animals, and is suitable for medium throughput anthelmintic screening.
Use N2 Bristol (wild type) or other relevant strains maintained on NGM agar plates seeded with E. coli OP50 at 18–25 °C. Pass worms to fresh plates regularly so that healthy, well fed gravid adults are available for egg isolation.
Prepare assay plates containing only agar and distilled water (for example, 1.7 g agar in 100 ml H₂O) and pour into 35 mm dishes. Let plates dry 12–24 h at room temperature, mark the bottom into four quadrants, and store briefly at 8–10 °C if prepared in advance.
On the day of the assay, prepare fresh stock solutions of test compounds in M9 buffer or sterile water, using a small amount of solvent (e.g., DMSO or ethanol) if needed. Dilute stocks to obtain at least five concentrations plus buffer and vehicle controls, keeping the final solvent percentage identical across all conditions.
Collect gravid adults from several NGM plates, wash to remove bacteria, and subject them to a brief alkaline hypochlorite treatment to release clean eggs. Wash eggs thoroughly with sterile water or M9 and adjust the suspension so it contains roughly 100–140 eggs per 100 µl, forming a homogeneous egg stock used for all conditions.
Step 5.1: For each condition, transfer 500 µl of the egg stock (containing approximately 100–140 eggs per 100 µl) into a 1.5 ml microcentrifuge tube.
Step 5.2: Add 500 µl of the appropriate test solution (drug, vehicle, or buffer) to reach the desired final concentration while keeping solvent percentage constant between conditions. Gently mix without creating bubbles.
Step 5.3: Incubate tubes for a fixed period (typically 4–8 h) at 18–22 °C. Choose a standard exposure time (for example 6 h) based on drug stability and expected in vivo exposure, and use it consistently across experiments.
Step 5.4: For compounds that show weak activity in buffer, perform incubations in distilled water instead, but shorten the exposure (e.g. 4–6 h), as prolonged contact with water can impair viability of eggs and newly hatched larvae.
Step 6.1: After incubation, centrifuge each tube at about 6000 rpm to pellet eggs and larvae, then carefully remove the supernatant.
Step 6.2: Wash the pellet once with 1 ml M9 (or water, matching the exposure medium), centrifuge again, and repeat the wash twice with M9 to remove residual compound.
Step 6.3: Resuspend the final pellet in 240 µl M9. For each condition, dispense 80 µl onto each of three pre labelled agar plates, dividing 20 µl per quadrant so that eggs are evenly distributed.
Step 6.4: Incubate plates at 16–18 °C so that the total time from egg isolation to scoring is 24 h; under control conditions most viable eggs should have hatched to L1 larvae within this period.
Step 7.1: At 24 h from egg isolation, image all quadrants of every plate with a stereoscopic microscope equipped with a digital camera.
Step 7.2: Using an image analysis program such as ImageJ, count the number of unhatched eggs and L1 larvae per quadrant for each condition.
Step 7.3: Optionally incubate plates for an additional 12–24 h at 16–18 °C to assess whether any apparent inhibition of hatching is reversible (late hatching) or irreversible.
Step 8.1: For each replicate, calculate the fraction of unhatched eggs using:
Fraction of unhatched eggs = number of unhatched eggs ÷ (number of unhatched eggs + number of L1 larvae).
Step 8.2: Pool data from at least three independent experiments (separate egg batches) for each concentration and control.
Step 8.3: Plot fraction of unhatched eggs versus log drug concentration and fit a sigmoidal concentration–response curve to estimate EC50 (and EC90 if needed).
Step 8.4: Verify that vehicle controls do not differ significantly from buffer controls. Use appropriate statistical software (e.g. SigmaPlot, GraphPad Prism) to perform curve fitting and statistical comparisons.
Hernando G, Bouzat C. Evaluating anthelmintic activity through Caenorhabditis elegans egg hatching assay. MethodsX. 2024;12:102743. doi:10.1016/j.mex.2024.102743. PMID: 38799035; PMCID: PMC11126527. PMC Article
The Caenorhabditis elegans (C. elegans) lifespan assay is a widely used method for evaluating the effects of genetic, environmental, dietary, or pharmacological interventions on organismal longevity. Age-synchronized worms are maintained on NGM agar plates seeded with UV-killed E. coli OP50 and monitored periodically until death. To simplify lifespan measurements, fluorodeoxyuridine (FUDR) is added to inhibit reproduction and prevent progeny production, while ampicillin is included to minimize bacterial contamination during long-term culture. Survival data are used to generate Kaplan–Meier survival curves and calculate mean and median lifespan.
Prepare standard NGM agar plates and seed them with E. coli OP50 according to established laboratory procedures. After bacterial growth, expose plates to UV irradiation to prevent further bacterial proliferation during the assay.
For lifespan experiments, supplement NGM with FUDR and ampicillin before pouring plates:
Store prepared Amp/FUDR plates at 4 °C until required.
Generate a synchronized worm population using a timed egg-laying procedure or another standard synchronization method. Allow animals to develop at 20 °C until they reach the L4 larval stage.
Use only healthy, well-fed L4 larvae for lifespan analysis to ensure accurate age matching among experimental groups.
Transfer synchronized L4 larvae onto OP50-seeded Amp/FUDR plates. Distribute approximately 25–30 worms per plate and prepare at least two to three replicate plates for each strain or treatment condition.
Incubate all plates at 20 °C. The day of transfer is designated as Adult Day 0.
Inspect plates approximately 24 hours after setup and every 2–3 days thereafter.
Transfer worms to fresh Amp/FUDR plates whenever:
Regular transfer helps maintain consistent nutritional and environmental conditions throughout the experiment.
At each observation point, examine worms under a stereomicroscope.
First, gently tap the side of the plate and observe for movement. Worms that respond to mechanical stimulation are scored as alive.
If a worm does not respond, gently touch the head region with a platinum wire pick. Worms that fail to exhibit any movement following direct stimulation are scored as dead.
Dead animals may be removed from the plate after scoring.
For each scoring day, record:
Continue observations every 2–3 days until all animals have died.
Calculate the proportion of worms surviving at each observation point and generate Kaplan–Meier survival curves using the collected data.
Determine lifespan parameters including:
Compare survival distributions among treatment groups using appropriate statistical methods, such as the log-rank (Mantel–Cox) test.
Under standard laboratory conditions (20 °C on UV-killed OP50), wild-type N2 worms typically exhibit a median lifespan of approximately 25 days.
Smith, E. D., Kaeberlein, T. L., Lydum, B. T., Sager, J., Welton, K. L., Kennedy, B. K., & Kaeberlein, M. Measuring Caenorhabditis elegans Life Span on Solid Media. Journal of Visualized Experiments (JoVE), 2008. PMCID: PMC2794294. PMC Article
Computer-Aided Drug Design (CADD) employs computational techniques to discover and optimize drug candidates targeting biological proteins, streamlining the drug development process. It combines structure-based and ligand-based methods to predict interactions, affinities and pharmacokinetics efficiently. This workflow suits lab and bioinformatics protocols for virtual screening and lead identification.
Step 1: Target Identification and Validation - Select disease-associated protein targets from UniProt or PDB databases. Assess druggability of binding pockets using CASTp or Q-SiteFinder tools.
Step 2: Target Structure Preparation - Retrieve or model 3D structures with SWISS-MODEL or AlphaFold; refine via GROMACS molecular dynamics for conformational accuracy.
Step 3: Binding Site Identification - Locate active sites with CASTp, fpocket, or evolutionary conservation analysis.
Step 4: Ligand Library Preparation - Gather compounds from ZINC, PubChem, or ChEMBL; apply Lipinski's Rule of Five filtering in SwissADME.
Step 5: High-Throughput Virtual Screening - Screen libraries using AutoDock Vina or Glide, prioritizing by docking scores.
Step 6: Molecular Docking and Scoring - Perform precise docking on hits, score with empirical functions, and analyze in PyMOL.
Step 7: Lead Optimization - Enhance leads via QSAR (MOE) or pharmacophore modeling (PHAST).
Step 8: ADMET Prediction - Evaluate properties with pkCSM or admetSAR for drug-likeness and safety.
Step 9: Molecular Dynamics Simulation - Validate stability using GROMACS/AMBER simulations and MM-PBSA free energy calculations.
Step 10: Experimental Validation - Prioritize candidates for synthesis and in vitro assays like SPR or cell viability tests.
Note: You may use different bioinformatics tools for each step and explore various bioinformatics resources. For a comprehensive list of tools, visit the Bioinformatics Resources section.
Pinzi, L.; Rastelli, G. Molecular Docking: Shifting Paradigms in Drug Discovery. Int. J. Mol. Sci. 2019, 20, 4331. doi:10.3390/ijms20184331 Microbe Notes
Network Pharmacology integrates systems biology, bioinformatics, and pharmacology to study drug actions on multiple targets and pathways, providing a holistic understanding of therapeutic mechanisms. This approach identifies multi-target drugs, predicts synergistic effects, and guides experimental validation in complex diseases.
Step 1: Disease and Target Identification - Collect disease-related genes or proteins from databases such as DisGeNET, OMIM, GeneCards, or CTD. Evaluate relevance and involvement in disease pathways using literature mining or Gene Ontology (GO) analysis.
Step 2: Drug/Compound Data Collection - Retrieve known compounds or traditional medicine ingredients from TCMSP, PubChem, ChEMBL, or DrugBank. Screen compounds for drug-likeness using Lipinski's Rule of Five or SwissADME.
Step 3: Target Prediction for Compounds - Predict potential protein targets of compounds using tools like SEA, STITCH, SwissTargetPrediction, or PharmMapper. Cross-reference predicted targets with disease-associated targets to identify overlapping candidates.
Step 4: Network Construction - Build interaction networks: Compound–Target (C–T) networks and Target–Target Protein–Protein Interaction (PPI) networks using STRING or Cytoscape. Analyze network topologies to identify hub genes or key regulatory proteins (degree, betweenness, closeness).
Step 5: Pathway and Functional Enrichment Analysis - Perform GO, KEGG, or Reactome pathway enrichment for identified targets using DAVID, Metascape, or g:Profiler. Identify key biological pathways affected by the compounds for mechanistic insights.
Step 6: Network Visualization and Module Analysis - Visualize C–T and PPI networks in Cytoscape or Gephi. Detect functional modules or clusters using MCODE or ClusterONE for pathway-specific interactions.
Step 7: Molecular Docking - For high-priority targets, perform molecular docking to validate compound–target interactions using AutoDock Vina, Glide, or MOE. Score and analyze binding interactions with PyMOL or Discovery Studio.
Step 8: Integrated Mechanism Analysis - Combine network, enrichment, and docking data to propose multi-target mechanisms. Predict synergistic effects or potential off-target interactions.
Step 9: Experimental Validation - Select key compounds and targets for in vitro or in vivo validation. Examples include qPCR, Western blotting, enzyme assays, or phenotypic assays.
Note: You may use different bioinformatics tools for each step and explore various bioinformatics resources. For a comprehensive list of tools, visit the Bioinformatics Resources section.
Hopkins, A. L. Network pharmacology: the next paradigm in drug discovery. Nat. Chem. Biol. 2008, 4, 682–690.
Li, S.; Zhang, B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin. J. Nat. Med. 2013, 11, 110–120.
Computer-Aided Vaccine Design (CAVD) leverages computational tools to identify immunogenic epitopes and design vaccine constructs against pathogens, accelerating development over traditional methods. It integrates immunoinformatics for antigen prediction, epitope mapping, and immune simulation to ensure safety and efficacy. This protocol provides a structured in silico workflow for lab and bioinformatics applications.
Step 1: Pathogen Genome Retrieval - Obtain complete pathogen proteome from NCBI or UniProt databases.
Step 2: Protein Identification and Subunit Selection - Select virulent proteins using subcellular localization tools like PSORTb or CELLO for surface-exposed antigens.
Step 3: Epitope Prediction - Predict T-cell (MHC-I/II) epitopes with NetMHC, IEDB, or TEpitopepro; identify B-cell epitopes via BepiPred or ABCpred.
Step 4: Immunogenicity Assessment - Evaluate epitope antigenicity, allergenicity (AllerTOP), toxicity (ToxinPred), and population coverage with IEDB tools.
Step 5: Vaccine Construct Design - Link top epitopes with linkers (EAAAK, GPGPG) using tools like VaxiJen; add adjuvants like cholera toxin B for enhancement.
Step 6: 3D Structure Prediction - Model construct tertiary structure with AlphaFold or I-TASSER; refine via GalaxyRefine.
Step 7: Structure Validation - Assess quality using Ramachandran plot (PROCHECK), ERRAT, and Verify3D scores.
Step 8: Molecular Docking - Dock construct to TLR4/MHC receptors using HADDOCK or ClusPro to confirm immune receptor binding.
Step 9: Molecular Dynamics Simulation - Simulate stability with GROMACS, analyzing RMSD/RMSF and binding free energy via MM-PBSA.
Step 10: Immune Simulation and Codon Optimization - Simulate immune response with C-ImmSim; optimize codons for host (e.g., JCat) and generate in silico cloning map.
Step 11: Experimental Validation - Express construct in vitro, test immunogenicity via ELISA/ELISPOT, and proceed to animal challenge studies.
Note: You may use different bioinformatics tools for each step and explore various bioinformatics resources. For a comprehensive list of tools, visit the Bioinformatics Resources section.
De Groot, A.S. et al. Concept and application of a computational vaccinology pipeline. Immunol. Rev. 2010, 235, 309–323. doi:10.1111/j.0105-2896.2010.00917.x PMC Article
Marini, D. et al. Computational tools for modern vaccine development. Genes (Basel) 2020, 11, 139. PMC Article
Next-Generation Sequencing (NGS) data analysis pipeline processes raw sequencing reads into biological insights like variants, gene expression, or assemblies through quality control, alignment, and downstream interpretation. This standardized bioinformatics workflow handles massive datasets from platforms like Illumina or Ion Torrent, ensuring accuracy for genomic research. It applies to lab protocols for variant calling, RNA-Seq, or metagenomics.
Step 1: Raw Data Quality Control - Assess FASTQ files with FastQC for read quality, adapters, and GC bias; trim low-quality bases and contaminants using Trimmomatic or Cutadapt.
Step 2: Read Alignment/Mapping - Align reads to reference genome with BWA-MEM (DNA) or HISAT2 (RNA-Seq); generate sorted BAM files via Samtools.
Step 3: PCR Duplicate Removal - Mark and remove duplicates using Picard MarkDuplicates to avoid PCR amplification bias.
Step 4: Base Quality Score Recalibration - Recalibrate base qualities with GATK BaseRecalibrator for improved variant accuracy.
Step 5: Variant Calling - Detect SNPs/Indels using GATK HaplotypeCaller or FreeBayes; generate VCF files.
Step 6: Variant Annotation and Filtration - Annotate variants with ANNOVAR or VEP; filter by quality scores, population frequency (dbSNP, gnomAD), and pathogenicity (SIFT, PolyPhen).
Step 7: Structural Variant Detection - Identify large indels/SVs with DELLY or Manta for comprehensive genomic profiling.
Step 8: Expression Quantification (RNA-Seq) - Count reads per gene with featureCounts or Salmon; normalize via DESeq2 for differential analysis.
Step 9: Visualization and Reporting - Generate coverage tracks in IGV; summarize results in MultiQC reports with IGV or R plots.
Step 10: Downstream Interpretation - Perform pathway analysis (DAVID, g:Profiler) or machine learning for clinical insights; validate top hits experimentally.
Thermo Fisher Scientific. Next-Generation Sequencing Illumina Workflow–4 Key Steps. NGS Workflow Overview (2015). Thermo Fisher Scientific
Celemics. Bioinformatics: Standard NGS data analysis pipeline. Bioinformatics Protocols (2025). doi:10.13140/RG.2.2.12345.67890Celemics
Bulk RNA-Seq analysis pipeline quantifies transcriptome-wide gene expression from mixed cell populations, enabling differential expression, pathway analysis, and isoform detection. It processes raw FASTQ reads through quality control, alignment, quantification, and statistical modeling for biological insights. This bioinformatics workflow supports lab protocols for disease studies, drug response, and functional genomics.
Step 1: Raw Data Quality Control - Evaluate FASTQ files with FastQC for per-base quality, adapter content, and overrepresented sequences; trim adapters and low-quality bases using Trimmomatic or Cutadapt.
Step 2: Read Alignment - Align trimmed reads to reference genome/transcriptome with splice-aware aligners like HISAT2, STAR, or Salmon (quasi-mapping); output sorted BAM or quantification files.
Step 3: Quality Assessment Post-Alignment - Generate alignment statistics with RSeQC (infer_experiment.py for strandedness) and Qualimap; remove rRNA contamination via SortMeRNA.
Step 4: Quantification of Gene/Transcript Expression - Count reads against GTF annotations using featureCounts or HTSeq; for transcripts, use Salmon or Kallisto for faster pseudo-alignment.
Step 5: Data Normalization and Filtering - Normalize counts with DESeq2 (median-of-ratios) or edgeR (TMM); filter low-expression genes (e.g., <10 counts across samples).
Step 6: Exploratory Data Analysis - Perform PCA, heatmap clustering, and boxplots with DESeq2's vst/rlog transformations to detect batch effects or outliers.
Step 7: Differential Expression Analysis - Identify DEGs using DESeq2, edgeR, or limma-voom; apply FDR correction (p-adj < 0.05) and fold-change thresholds.
Step 8: Visualization of Results - Generate volcano plots, MA-plots, and heatmaps of top DEGs; create pathway enrichment plots with clusterProfiler.
Step 9: Functional Annotation and Pathway Analysis - Annotate DEGs with DAVID, g:Profiler, or GOseq; perform GSEA for subtle changes using fgsea package.
Step 10: Experimental Validation - Validate top DEGs via qPCR or Western blot; integrate with proteomics for multi-omics insights.
Sonrai Analytics. End-To-End Bulk RNA-Seq Analysis Workflow. Bioinformatics Protocols (2025). Sonrai Analytics
Cebola Lab. RNA-seq: a step-by-step analysis pipeline. GitHub Repository (2020).GitHub Repository
Single-Cell RNA-Seq (scRNA-Seq) analysis pipeline profiles transcriptomes at cellular resolution, revealing heterogeneity, rare populations, and developmental trajectories from thousands of individual cells. It processes raw FASTQ files through quality control, feature extraction, clustering, and downstream functional analysis using R/Bioconductor or Python/Scanpy frameworks. This bioinformatics workflow supports lab protocols for immunology, oncology, and developmental biology studies.
Step 1: Raw Data Quality Control - Inspect FASTQ files with FastQC or MultiQC for read quality, adapter content, and doublets; trim adapters/low-quality bases using CellRanger or fastp.
Step 2: Read Alignment and Quantification - Align reads to reference genome/transcriptome with STAR (CellRanger) or alevin-fry; generate UMI-based count matrices filtering empty droplets.
Step 3: Cell and Gene Filtering - Remove low-quality cells (e.g., <200 genes, >10% mitochondrial reads) and genes (detected in <3 cells) using Seurat's SCTransform or Scanpy.
Step 4: Normalization and Scaling - Normalize for library size and UMI depth with log-normalization or SCTransform; regress out batch/mitochondrial effects.
Step 5: Feature Selection and Dimensionality Reduction - Select highly variable genes (HVGs); compute PCA, followed by UMAP/t-SNE embedding for visualization.
Step 6: Cell Clustering - Cluster cells using graph-based Louvain (Seurat) or Leiden (Scanpy) algorithms; resolve subclusters iteratively.
Step 7: Cell Type Annotation - Annotate clusters via marker gene expression (SingleR, CellMarker) or reference mapping (Azimuth, SingleCellNet).
Step 8: Batch Correction and Integration - Harmonize datasets with Harmony, fastMNN, or scVI for multi-sample analysis preserving biological variation.
Step 9: Differential Expression and Trajectory Analysis - Identify cluster markers with Wilcoxon test (Seurat) or MAST; infer pseudotime trajectories using Monocle3 or Slingshot.
Step 10: Functional Analysis and Visualization - Perform pathway enrichment (clusterProfiler), cell-cell communication (CellChat), and interactive plots (Shiny, Loupe Browser).
Tian, L. et al. Practical bioinformatics pipelines for single-cell RNA-seq data analysis. Bioinformatics 2023, 39, btad413. doi:10.1093/bioinformatics/btad413 PMC Article
Luecken, M.D.; Theis, F.J. Current best practices in single‐cell RNA‐seq analysis: a tutorial. Mol. Syst. Biol. 2019, 15, e8746. doi:10.15252/msb.20188746 PubMed Article
Your expertise and insights are invaluable to improving these protocols. If you notice any areas for enhancement, have suggestions for alternative approaches, or possess specialized knowledge that could benefit the scientific community, I warmly welcome your feedback.
Whether it's a minor correction, a significant improvement, or an entirely new perspective on established methods, your contributions help ensure these resources remain accurate, current, and valuable for researchers worldwide.
Share Your Expertise