Protein Phylogeny
by @billwanttobetop
Comprehensive protein family phylogenetic analysis workflow with quality control, conservation analysis, coevolution network analysis, and publication-ready...
clawhub install protein-phylogenyπ About This Skill
name: protein-phylogeny description: "Comprehensive protein family phylogenetic analysis workflow with quality control, conservation analysis, coevolution network analysis, and publication-ready visualization. Use when: (1) analyzing protein family evolution, (2) building phylogenetic trees from sequences, (3) identifying conserved/coevolved residues, (4) generating publication-quality figures and reports, (5) quality-controlling sequence datasets, or (6) performing systematic evolutionary analysis of enzyme families, protein superfamilies, or any homologous protein groups."
Protein Family Phylogenetic Analysis
Complete workflow for protein family evolutionary analysis: quality control β conservation β coevolution β phylogeny β publication report.
Quick Start
Input: FASTA file with protein sequences (any family, any size) Output: Publication-ready report with phylogenetic tree, conservation analysis, coevolution networks, and high-quality figures
Typical workflow:
# 1. Quality control (removes low-quality sequences)
bash scripts/01_quality_control.sh input.fasta output_dir/2. Conservation analysis
bash scripts/02_conservation.sh output_dir/qc/final.fasta output_dir/3. Coevolution analysis
bash scripts/03_coevolution.sh output_dir/qc/final.fasta output_dir/4. Phylogenetic tree
bash scripts/04_phylogeny.sh output_dir/qc/final.fasta output_dir/5. Generate figures
bash scripts/05_visualize.sh output_dir/6. Create report
bash scripts/06_report.sh output_dir/ "Family Name"
Workflow Overview
Stage 1: Quality Control (references/01-quality-control.md)
Purpose: Filter raw sequences to high-quality, non-redundant dataset
Steps: 1. Literature validation (remove predicted sequences) 2. Length filtering (remove fragments/fusions) 3. CD-HIT redundancy removal (90% identity) 4. Complexity check (remove low-complexity regions) 5. Motif validation (confirm family membership) 6. MAFFT alignment (high accuracy mode) 7. trimAl trimming (automatic strategy) 8. Final validation (gap ratio, coverage)
Key parameters:
Output: qc/final.fasta (high-quality aligned sequences)
Stage 2: Conservation Analysis (references/02-conservation.md)
Purpose: Identify functionally important conserved residues
Method: Shannon entropy
Output:
Stage 3: Coevolution Analysis (references/03-coevolution.md)
Purpose: Identify residue pairs that evolve together
Method: Normalized Mutual Information (NMI)
Output:
Stage 4: Phylogenetic Analysis (references/04-phylogeny.md)
Purpose: Reconstruct evolutionary relationships
Method: IQ-TREE maximum likelihood
Output:
Stage 5: Visualization (references/05-visualization.md)
Purpose: Generate publication-quality figures (300 DPI)
Figures: 1. Workflow diagram 2. Conservation heatmap 3. Coevolution network 4. Hub analysis 5. Quality metrics 6. Phylogenetic tree 7. Bootstrap distribution 8. Supplementary plots
Style: Clean, colorblind-friendly, Nature/Science standards
Stage 6: Report Generation (references/06-report.md)
Purpose: Create comprehensive analysis report
Sections: 1. Overview (dataset summary) 2. Quality control (methods + results) 3. Conservation analysis (algorithms + findings) 4. Coevolution analysis (networks + hubs) 5. Phylogenetic analysis (tree + support) 6. Quality assessment (standards comparison) 7. Conclusions (biological insights)
Format: Markdown β Feishu/Word/PDF
Key Features
AI-Friendly Design
Token Efficiency
Professional Quality
Dependencies
Required tools:
Installation:
bash scripts/install_dependencies.sh
Common Pitfalls
1. Low Sequence Similarity (< 25%)
Problem: Alignment unreliable, phylogeny uncertain Solution:
2. High Gap Ratio (> 30%)
Problem: Many unreliable positions Solution:
-gt 0.8)3. Bootstrap Convergence Failure (< 0.99)
Problem: Tree topology unstable Solution:
4. No Conserved Motifs
Problem: Family definition unclear Solution:
Advanced Usage
Custom Quality Control
Edit scripts/01_quality_control.sh parameters:
CDHIT_THRESHOLD=0.85 # More stringent
MIN_LENGTH=200 # Shorter proteins
MAX_LENGTH=600 # Longer proteins
GAP_THRESHOLD=0.25 # Stricter gap cutoff
Alternative Phylogeny Methods
See references/04-phylogeny.md for:
Custom Visualization
Edit scripts/05_visualize.sh for:
Troubleshooting
Issue: CD-HIT crashes with large datasets Fix: Split input, process in batches, merge results
Issue: IQ-TREE runs forever
Fix: Use -fast mode or reduce bootstrap replicates
Issue: Figures look pixelated
Fix: Increase DPI in scripts/05_visualize.sh (default 300)
Issue: Report generation fails Fix: Check all intermediate files exist, rerun failed stages
References
For detailed methodology, see:
Citation
If you use this workflow, please cite:
Example Usage
# Download your sequences
(from UniProt, NCBI, or your own database)
Run full workflow
bash scripts/run_full_workflow.sh sequences.fasta analysis_output/ "Your Family Name"Results in analysis_output/:
- qc/final.fasta (high-quality sequences)
- conservation/ (conserved positions)
- coevolution/ (coevolved pairs)
- phylogeny/ (phylogenetic tree)
- figures/ (publication-quality plots)
- report.md (complete analysis)
π‘ Examples
Input: FASTA file with protein sequences (any family, any size) Output: Publication-ready report with phylogenetic tree, conservation analysis, coevolution networks, and high-quality figures
Typical workflow:
# 1. Quality control (removes low-quality sequences)
bash scripts/01_quality_control.sh input.fasta output_dir/2. Conservation analysis
bash scripts/02_conservation.sh output_dir/qc/final.fasta output_dir/3. Coevolution analysis
bash scripts/03_coevolution.sh output_dir/qc/final.fasta output_dir/4. Phylogenetic tree
bash scripts/04_phylogeny.sh output_dir/qc/final.fasta output_dir/5. Generate figures
bash scripts/05_visualize.sh output_dir/6. Create report
bash scripts/06_report.sh output_dir/ "Family Name"
π Tips & Best Practices
Issue: CD-HIT crashes with large datasets Fix: Split input, process in batches, merge results
Issue: IQ-TREE runs forever
Fix: Use -fast mode or reduce bootstrap replicates
Issue: Figures look pixelated
Fix: Increase DPI in scripts/05_visualize.sh (default 300)
Issue: Report generation fails Fix: Check all intermediate files exist, rerun failed stages