Analyze genome collinearity and conserved gene order between species using MCScanX, SyRI, and JCVI for evolutionary and comparative genomics studies.
# MCScanX (compile from source)
git clone https://github.com/wyp1125/MCScanX
cd MCScanX && make
export PATH=$PATH:$(pwd)
# SyRI for structural variants
pip install syri
# JCVI for visualization
pip install jcvi
# BLAST for homology search
conda install -c bioconda blastTell your AI agent what you want to do:
- "Find syntenic blocks between human and mouse genomes"
- "Detect whole-genome duplications in my plant genome"
- "Identify chromosomal rearrangements between two assemblies"
"Find collinear gene blocks between Arabidopsis and rice"
"Create a synteny dot plot comparing my two genomes"
"Look for whole-genome duplication signatures in this genome"
"Calculate Ks distribution for syntenic gene pairs"
"Identify inversions and translocations between reference and query"
"Find structural rearrangements between chromosome assemblies"
- Prepare genome and annotation files for analysis
- Run all-vs-all BLASTP for homology detection
- Identify collinear gene blocks with MCScanX
- Classify syntenic relationships (1:1, 1:many)
- Detect structural variants with SyRI if requested
- Calculate Ks for dating duplications
- Generate visualization plots
- Block size - Minimum 5 genes per block reduces noise; use 10+ for stringent analysis
- E-value - Use 1e-10 for close species, 1e-5 for distant comparisons
- Ks saturation - Values >2 are unreliable; filter before interpretation
- WGD peaks - Multiple Ks peaks suggest multiple WGD events
- Visualization - JCVI produces publication-quality synteny plots