Explain comparative genomics — Biology notes

What IS Comparative Genomics?

The Foundation: How Evolution Creates Comparable Genomes

Derivation from First Principles

Starting point: All life shares a common ancestor (LUCA - Last Universal Common Ancestor, ~3.5 billion years ago).

Step 1 - Descent with Modification:

When species diverge, they inherit a copy of the ancestral genome
Over time: Point mutations, insertions, deletions, duplications accumulate
Rate of change depends on: mutation rate (μ), generation time, population size, selection pressure

Step 2 - The Molecular Clock: For neutral mutations (no selection): $K =2\mu t$

Where:

$K$ = number of substitutions per site
$\mu$ = mutation rate per site per generation
$t$ = time since divergence (in generations)

WHY the factor of 2? Each lineage accumulates mutations independently after the split, so total divergence is the sum of changes in BOTH lineages.

Step 3 - Selection Creates Conservation: For functionally important sequences, the substitution rate is reduced: $K_a = K_0 \cdot (1 - s)$

Where:

$K_a$ = actual substitution rate
$K_0$ = neutral rate
$s$ = selection coefficient (0= neutral, 1 = completely conserved)

This explains why some regions evolve slowly: strong purifying selection eliminates harmful mutations.

Key Concepts in Comparative Genomics

1. Homology: The Bridge Between Genomes

WHY this matters: Orthologues let us transfer functional knowledge across species. If we know a mouse gene causes a disease, its human orthologue is a candidate for the human version of that disease.

2. Sequence Alignment: The Core Tool

HOW alignment works:

Given two sequences, find the arrangement that maximizes similarity while accounting for:

Matches: Same nucleotide/amino acid at corresponding positions (+score)
Mismatches: Different nucleotide/amino acid (−score)
Gaps: Insertions or deletions in one sequence (−gap penalty)

The scoring function: $S = \sum_{i=1}^{L} s(a_i, b_i) - \sum_{j=1}^{G} g(l_j)$

Where:

$s(a_i, b_i)$ = match/mismatch score at position $i$
$g(l_j)$ = gap penalty for gap $j$ of length $l_j$
$L$ = alignment length
$G$ = number of gaps

WHY gap penalties? Indels (insertions/deletions) are rarer than point mutations, so we need to penalize opening gaps to avoid artificial gaps that inflate similarity.

Worked example Alignment Example

Sequences:

Species A: ATGCATGC
Species B: ATGC--GC Alignment:

A: ATGCATGC
B: ATGC--GC
   ||||  ||

Analysis:

6 matches (score: +6if match = +1)
2 gap positions (score: −2 if gap = −1)
Total score: 6 − 2 = 4 WHY this step? The alignment reveals that two nucleotides were deleted in species B (or inserted in A), but the flanking sequences are conserved, suggesting functional importance of the ATGC and GC regions.

3. Synteny: Chromosomal Organization

Quantifying synteny: $C = \frac{N_{syntenic}}{N_{total}} \times 100\%$

Where $C$ is the percentage of genes in conserved syntenic blocks.

Example: Human chr. 17 vs. mouse chr. 11 show ~80% synteny, while human vs. chicken show ~40% (more distant evolutionary relationship).

4. Ka/Ks Ratio: Measuring Selection Pressure

Derivation of interpretation:

Under neutral evolution, synonymous and non-synonymous sites should accumulate substitutions at equal rates (both proportional to mutation rate). Thus $K_a/K_s = 1$ is the null expectation.

If $K_a < K_s$ ( $\omega < 1$ ), non-synonymous mutations are being removed by selection faster than they arise → function is important.

If $K_a > K_s$ ( $\omega > 1$ ), non-synonymous changes are being FIXED by selection faster than expected → advantageous amino acid changes → adaptation.

Applications of Comparative Genomics

1. Gene Discovery and Functional Annotation

Method: Phylogenetic footprinting

HOW: Align genomes from multiple species → identify regions with unexpectedly high conservation → these are likely functional elements (genes, regulatory elements).

Example: The ENCODE project used29 mamalian genomes. Regions conserved across all mammals are enriched for:

Protein-coding exons (as expected)
Promoters and enhancers (regulatory DNA)
Structural RNAs (tRNAs, miRNAs)

Discovery: ~8% of the human genome shows conservation, but only ~1.5% codes for proteins → 6.5% is functional non-coding DNA (regulatory elements, structural elements).

2. Understanding Disease

Logic: If a human gene orthologue in mice causes a disease when mutated, humans with mutations in that gene likely have similar disease.

3. Evolutionary Insights

Comparative genomics reveals:

Speciation timing: More divergence → older split
Adaptive radiations: Rapid accumulation of differences in specific gene families
Gene loss: Genes lost in certain lineages (e.g., humans lost ability to synthesize vitamin C)

Example: Comparing primate genomes shows FOXP2 gene (speech-related) has accelerated evolution in humans after split from chimps → adaptive evolution related to language.

Common Mistakes and Fixes

Practical Workflow

Worked example Complete Comparative Genomics Analysis

Goal: Characterize an unknown human gene Step 1 - Find orthologues:

BLAST search against multiple species (mouse, zebrafish, fly, yeast)
Result: Orthologues found in all vertebrates, absent in invertebrates
Interpretation: Vertebrate-specific innovation or extreme divergence invertebrates

Step 2 - Multiple sequence alignment:

Align human, mouse, chicken, frog, zebrafish sequences
Identify conserved domains using tools like Pfam
Result: N-terminal domain highly conserved (80%+ identity), C-terminal variable
Interpretation: N-terminal domain is functionally critical, C-terminal may be regulatory/species-specific

Step 3 - Synteny analysis:

Check gene order around the target gene
Result: Syntenic block conserved in mammals, rearranged in fish
Interpretation: Mamalian chromosomal organization constraint or functional gene cluster

Step 4 - Ka/Ks calculation:

Calculate pairwise Ka/Ks for all species pairs
Result: Average $\omega = 0.15$ for N-terminal domain, $\omega = 0.8$ for C-terminal
Interpretation: Strong purifying selection on N-terminal (essential function), near-neutral evolution of C-terminal (less constrained)

Step 5 - Regulatory element prediction:

Align5 kb upstream regions
Find conserved non-coding sequences (CNS)
Result: Three CNS blocks, each ~50-100 bp, with 70%+ conservation
Interpretation: Likely transcription factor binding sites (enhancers/promoters)

Step 6 - Functional prediction:

Conserved N-terminal domain matches known kinase domain
Mouse orthologue knockout causes embryonic lethality (literature search)
Conclusion: Essential kinase involved in early development, with conserved catalytic activity but species-specific regulation

Deep Connections

Connections:

Molecular Evolution - The theoretical foundation for how sequences diverge
Phylogenetics - Building evolutionary trees from comparative genomics data
Functional Genomics - Using comparative data to predict gene function
BLAST and Sequence Similarity - The computational tool for finding homologues
Gene Duplication and Evolution - How genomes expand and paralogues arise
Regulatory Evolution - How non-coding sequence changes drive phenotypic diversity
Selection Types - Purifying, positive, and neutral selection evidenced by Ka/Ks
Whole Genome Duplication - Large-scale events creating opportunities for comparative analysis
Metagenomics - Comparative genomics applied to entire microbial communities
Human Genome Project - Foundation for modern comparative genomics

Recall Explain to a 12-year-old

Imagine you and your cousin both got the same LEGO instruction booklet from your grandma. Over years, you both built your sets, but sometimes you lost pieces, sometimes you added custom pieces, and sometimes you made mistakes and fixed them differently.

Now, if I put your two LEGO creations side by side, I can figure out:

What was in the original instruction booklet (the pieces you BOTH have are probably from grandma)
What's really important (if you BOTH kept the same pieces despite losing others, those must be essential)
How you're related (the more similar your LEGOs, the more recently you both got the instructions from grandma)
What makes you unique (the pieces only YOU have are your special additions)

Comparative genomics is exactly this, but with DNA instead of LEGO! We compare the "instruction manuals" (genomes) of different animals to figure out how life works, how species are related, and what makes each species special. When we find DNA that's EXACTLY the same in humans, mice, and fish, we know "this must be super important—evolution kept it the same for millions of years!"

Flashcards

#flashcards/biology

What is comparative genomics? :: The analysis and comparison of genome sequences across different organisms to identify conserved regions, understand evolutionary relationships, discover functional elements, and predict gene function.

Why does sequence conservation indicate functional importance?

Evolution removes harmful mutations through purifying selection. If a sequence is conserved across species, it means mutations that change it are harmful and eliminated, indicating the sequence performs a critical function.

What is the difference between orthologues and paralogues?

Orthologues are genes in different species that evolved from a common ancestral gene through speciation (usually same function). Paralogues are genes related by duplication within a genome (often different/modified functions).

What does the molecular clock equation K = 2μt represent?

K is the number of substitutions per site, μ is the mutation rate per site per generation, and t is time since divergence. The factor of 2 accounts for mutations accumulating independently in both lineages after they split.

What is synteny?

The preserved order of genes along chromosomes across different species, suggesting functional linkage, structural constraints, or recent divergence.

What does Ka/Ks ratio measure and how is it interpreted?

Ka/Ks (ω) measures the ratio of non-synonymous to synonymous substitution rates. ω < 1 indicates purifying selection (conserved function), ω = 1 indicates neutral evolution, and ω > 1 indicates positive selection (adaptive evolution).

Why are synonymous sites used as neutral baseline in Ka/Ks analysis?

Synonymous mutations don't change the amino acid, so they're usually not affected by selection. They accumulate at the neutral mutation rate, providing a baseline to compare against non-synonymous changes that DO affect protein function.

What is phylogenetic footprinting?

A method that aligns genomes from multiple species to identify regions with unexpectedly high conservation, which likely represent functional elements like genes, promoters, enhancers, or structural RNAs.

Why would a gene have high Ka/Ks ratio (ω close to or above 1)?

High Ka/Ks suggests either relaxed purifying selection (function is less critical) or positive selection (amino acid changes are advantageous). Common in immune system genes and genes involved in species-specific adaptations.

How does comparative genomics help with disease gene discovery?

If synteny and conservation identify a gene's location and importance, and the orthologue in model organisms (like mice) causes disease when mutated, then mutations in the human orthologue are likely candidates for causing similar human diseases.

What does it mean if a gene is conserved in all vertebrates but absent in invertebrates?

The gene likely represents a vertebrate-specific innovation, arose after the vertebrate-invertebrate split, or has diverged so much invertebrates that it's unrecognizable by sequence similarity alone.

Why use multiple species at different phylogenetic distances?

Close relatives reveal recent adaptations and lineage-specific features. Distant relatives identify deeply conserved essential functions. Multiple distances help resolve the evolutionary timing of when features arose or were lost.

Concept Map

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, comparative genomics ek bahut powerful technique hai jisme hum different species ke genomes ko compare karte hain taki yeh samajh sakein ki evolution kaise kaam karta hai aur genes ka function kya hai. Soch lo tumhare pas ek purani family recipe book hai jo tumhari dadi ne likhi thi. Ab tum aur tumhare cousins sab ke pas usi book ki copy hai, lekin saalon mein kuch changes ho gaye—kisi ne ingredients change kiye, kisi ne steps modify kiye. Agar

Test yourself — Genomics

Connections