\r

Volcano Plot Script Generator\r

\r A skill for generating publication-ready volcano plots from differential gene expression analysis results.\r \r

When to Use\r

\r

• Use this skill when the task is to Generate R/Python code for volcano plots from DEG (Differentially Expressed Genes) analysis results. Triggered when user needs visualization of gene expression data, p-value vs fold-change scatter plots, publication-ready figures for bioinformatics analysis.\r
• Use this skill for data analysis tasks that require explicit assumptions, bounded scope, and a reproducible output format.\r
• Use this skill when you need a documented fallback path for missing inputs, execution errors, or partial evidence.\r \r

Key Features\r

\r

• Scope-focused workflow aligned to: Generate R/Python code for volcano plots from DEG (Differentially Expressed Genes) analysis results. Triggered when user needs visualization of gene expression data, p-value vs fold-change scatter plots, publication-ready figures for bioinformatics analysis.\r
• Packaged executable path(s): scripts/main.py.\r
• Reference material available in references/ for task-specific guidance.\r
• Reusable packaged asset(s), including assets/example_volcano.R.\r
• Structured execution path designed to keep outputs consistent and reviewable.\r \r

Dependencies\r

\r

Example Usage\r

\r See ## Usage above for related details.\r \r

cd "20260318/scientific-skills/Data Analytics/volcano-plot-script"\r
python -m py_compile scripts/main.py\r
python scripts/main.py --help\r
```\r
\r
Example run plan:\r
1. Confirm the user input, output path, and any required config values.\r
2. Edit the in-file `CONFIG` block or documented parameters if the script uses fixed settings.\r
3. Run `python scripts/main.py` with the validated inputs.\r
4. Review the generated output and return the final artifact with any assumptions called out.\r
\r
## Implementation Details\r
\r
See `## Workflow` above for related details.\r
\r
- Execution model: validate the request, choose the packaged workflow, and produce a bounded deliverable.\r
- Input controls: confirm the source files, scope limits, output format, and acceptance criteria before running any script.\r
- Primary implementation surface: `scripts/main.py`.\r
- Reference guidance: `references/` contains supporting rules, prompts, or checklists.\r
- Packaged assets: reusable files are available under `assets/`.\r
- Parameters to clarify first: input path, output path, scope filters, thresholds, and any domain-specific constraints.\r
- Output discipline: keep results reproducible, identify assumptions explicitly, and avoid undocumented side effects.\r
\r
## Quick Check\r
\r
Use this command to verify that the packaged script entry point can be parsed before deeper execution.\r
\r
```bash\r
python -m py_compile scripts/main.py\r
```\r
\r
## Audit-Ready Commands\r
\r
Use these concrete commands for validation. They are intentionally self-contained and avoid placeholder paths.\r
\r
```bash\r
python -m py_compile scripts/main.py\r
python scripts/main.py --help\r
python scripts/main.py --input "Audit validation sample with explicit symptoms, history, assessment, and next-step plan."\r
```\r
\r
## Workflow\r
\r
1. Confirm the user objective, required inputs, and non-negotiable constraints before doing detailed work.\r
2. Validate that the request matches the documented scope and stop early if the task would require unsupported assumptions.\r
3. Use the packaged script path or the documented reasoning path with only the inputs that are actually available.\r
4. Return a structured result that separates assumptions, deliverables, risks, and unresolved items.\r
5. If execution fails or inputs are incomplete, switch to the fallback path and state exactly what blocked full completion.\r
\r
## Overview\r
\r
Volcano plots visualize the relationship between statistical significance (p-values) and magnitude of change (fold changes) in gene expression data. This skill generates customizable R or Python scripts for creating high-quality figures suitable for publications.\r
\r
## Use Cases\r
\r
- Visualize RNA-seq DEG analysis results\r
- Identify significantly upregulated and downregulated genes\r
- Highlight genes of interest (markers, pathways)\r
- Generate publication-quality figures for manuscripts\r
- Compare multiple experimental conditions\r
\r
## Input Requirements\r
\r
Required input data format:\r
- Gene identifier (gene symbol or ENSEMBL ID)\r
- Log2 fold change values\r
- Adjusted or raw p-values\r
- Optional: gene annotations, pathways\r
\r
## Output\r
\r
- Publication-ready volcano plot (PNG/PDF/SVG)\r
- Customizable R or Python script\r
- Optional: labeled significant gene lists\r
\r
## Usage\r
\r
```python\r
\r
# Example: Run the volcano plot generator\r
python scripts/main.py --input deg_results.csv --output volcano_plot.png\r
```\r
\r
## Parameters\r
\r
| Parameter | Description | Default |\r
|-----------|-------------|---------|\r
| `--input` | Path to DEG results CSV/TSV | required |\r
| `--output` | Output plot file path | volcano_plot.png |\r
| `--log2fc-col` | Column name for log2 fold change | log2FoldChange |\r
| `--pvalue-col` | Column name for p-value | padj |\r
| `--gene-col` | Column name for gene IDs | gene |\r
| `--log2fc-thresh` | Log2 FC threshold for significance | 1.0 |\r
| `--pvalue-thresh` | P-value threshold | 0.05 |\r
| `--label-genes` | File with genes to label | None |\r
| `--top-n` | Label top N significant genes | 10 |\r
| `--color-up` | Color for upregulated genes | #E74C3C |\r
| `--color-down` | Color for downregulated genes | #3498DB |\r
| `--color-ns` | Color for non-significant genes | #95A5A6 |\r
\r
## Technical Difficulty\r
\r
**Medium** - Requires understanding of:\r
- DEG analysis concepts (fold change, p-values, FDR)\r
- Data visualization principles\r
- Matplotlib/ggplot2 plotting libraries\r
\r
### Python\r
- pandas\r
- matplotlib\r
- seaborn\r
- numpy\r
\r
### R\r
- ggplot2\r
- dplyr\r
- ggrepel (for label positioning)\r
\r
## References\r
\r
- [Example datasets and templates](references/)\r
- Best practices for volcano plot visualization\r
- Color schemes for accessibility\r
\r
## Author\r
\r
Auto-generated skill for bioinformatics visualization.\r
\r
## Risk Assessment\r
\r
| Risk Indicator | Assessment | Level |\r
|----------------|------------|-------|\r
| Code Execution | Python/R scripts executed locally | Medium |\r
| Network Access | No external API calls | Low |\r
| File System Access | Read input files, write output plots | Medium |\r
| Instruction Tampering | Standard prompt guidelines | Low |\r
| Data Exposure | Output files saved to workspace | Low |\r
\r
## Security Checklist\r
\r
- [ ] No hardcoded credentials or API keys\r
- [ ] Input file paths validated (no ../ traversal)\r
- [ ] Output directory restricted to workspace\r
- [ ] Script execution in sandboxed environment\r
- [ ] Error messages sanitized (no stack traces exposed)\r
- [ ] Dependencies audited (pandas, matplotlib, seaborn, numpy)\r
\r
## Prerequisites\r
\r
```text\r
\r
# Python dependencies\r
pip install -r requirements.txt\r
\r
# R dependencies (if using R)\r
install.packages(c("ggplot2", "dplyr", "ggrepel"))\r
```\r
\r
## Evaluation Criteria\r
\r
### Success Metrics\r
- [ ] Successfully generates executable Python/R script\r
- [ ] Output plot is publication-ready quality\r
- [ ] Correctly identifies significant genes based on thresholds\r
- [ ] Handles missing or malformed data gracefully\r
- [ ] Color scheme is accessible (colorblind-friendly)\r
\r
### Test Cases\r
1. **Basic DEG Visualization**: Input standard DESeq2 results → Valid volcano plot\r
2. **Custom Thresholds**: Adjust log2FC and p-value thresholds → Correct gene classification\r
3. **Gene Labeling**: Specify genes to label → Labels appear correctly\r
4. **Large Dataset**: Input 20,000+ genes → Performance remains acceptable\r
5. **Malformed Data**: Input with missing values → Graceful error handling\r
\r
## Lifecycle Status\r
\r
- **Current Stage**: Draft\r
- **Next Review Date**: 2026-03-06\r
- **Known Issues**: None\r
- **Planned Improvements**:\r
  - Add interactive plot option (Plotly)\r
  - Support for multiple comparison groups\r
  - Integration with pathway enrichment tools\r
\r
## Output Requirements\r
\r
Every final response should make these items explicit when they are relevant:\r
\r
- Objective or requested deliverable\r
- Inputs used and assumptions introduced\r
- Workflow or decision path\r
- Core result, recommendation, or artifact\r
- Constraints, risks, caveats, or validation needs\r
- Unresolved items and next-step checks\r
\r
## Error Handling\r
\r
- If required inputs are missing, state exactly which fields are missing and request only the minimum additional information.\r
- If the task goes outside the documented scope, stop instead of guessing or silently widening the assignment.\r
- If `scripts/main.py` fails, report the failure point, summarize what still can be completed safely, and provide a manual fallback.\r
- Do not fabricate files, citations, data, search results, or execution outcomes.\r
\r
## Input Validation\r
\r
This skill accepts requests that match the documented purpose of `volcano-plot-script` and include enough context to complete the workflow safely.\r
\r
Do not continue the workflow when the request is out of scope, missing a critical input, or would require unsupported assumptions. Instead respond:\r
\r
> `volcano-plot-script` only handles its documented workflow. Please provide the missing required inputs or switch to a more suitable skill.\r
\r
## Response Template\r
\r
Use the following fixed structure for non-trivial requests:\r
\r
1. Objective\r
2. Inputs Received\r
3. Assumptions\r
4. Workflow\r
5. Deliverable\r
6. Risks and Limits\r
7. Next Checks\r
\r
If the request is simple, you may compress the structure, but still keep assumptions and limits explicit when they affect correctness.\r