Comparative Genomics
Evolutionary Track • Chapter 6
Estimated Study Time: 22 minutes
Learning Objectives
- Distinguish homologous gene classes: orthologs, paralogs (ohnologs, tandem duplicates), and xenologs.
- Identify syntenic chromosomal blocks and detect large-scale evolutionary inversions and translocations.
- Calculate and interpret Ka/Ks selection pressure ratios across homologous gene alignments.
- Construct population-scale pangenome reference variation graphs with Minigraph-Cactus.
1. Orthology, Paralogy & Evolutionary Homology
Comparative genomics leverages evolutionary divergence to understand genome structure, gene function, and adaptation. Central to this field is the classification of homologous genes (genes descended from a common ancestral sequence):
- Orthologs: Genes separated by a speciation event (e.g. human β-globin and mouse β-globin). Orthologs generally retain identical or closely related biological functions across species.
- Paralogs: Genes separated by a gene duplication event within a genome (e.g. α-globin and β-globin). Gene duplication frees one copy from purifying selection, allowing it to acquire novel functions (neofunctionalization) or divide ancestral subfunctions (subfunctionalization).
- Ohnologs: Paralogs that arose through ancestral Whole Genome Duplication (WGD) events (such as the two rounds of WGD, 2R, at the base of vertebrate evolution).
2. Measuring Selection Pressure: The Ka/Ks Ratio
By comparing the rate of non-synonymous amino acid-altering substitutions ($K_a$ or $dN$) to the rate of synonymous silent substitutions ($K_s$ or $dS$), we quantify natural selection acting upon protein-coding genes:
$$\omega = \frac{K_a}{K_s}$$
- $\omega < 1$ (Purifying / Negative Selection): Non-synonymous mutations are deleterious and eliminated by selection. Most essential structural and catalytic proteins exhibit $\omega pprox 0.05 - 0.2$.
- $\omega = 1$ (Neutral Evolution): Non-synonymous mutations accumulate at the same rate as synonymous mutations, indicating neutral genetic drift (common in pseudogenes).
- $\omega > 1$ (Positive / Diversifying Selection): Non-synonymous amino acid changes confer selective fitness advantages (e.g. host-pathogen arms races, immune receptors, viral capsid antigens).
3. Hands-on Orthology & Synteny Commands
# 1. High-throughput orthology inference with OrthoFinder
orthofinder -f fasta_proteomes/ -t 16 -M msa -T fasttree
# 2. Synteny block identification with MCScanX
# Combine BLASTP hits and GFF gene coordinates
MCScanX my_species_comparison -k 50 -s 5
# 3. Whole-genome progressive alignment with Cactus
cactus jobStore seqFile.txt alignment.hal --maxCores 16