Overview
For decades, the evolutionary history of squid and cuttlefish remained a puzzle. How did these intelligent, soft-bodied creatures survive the planet's most devastating mass extinctions? New genomic research has finally provided answers, revealing that these cephalopods originated in the deep ocean over 100 million years ago and survived catastrophic events by retreating to oxygen-rich deep-sea refuges. A period of stasis followed, until a post-extinction boom triggered rapid diversification into shallow waters.
Prerequisites
Knowledge Requirements
- Basic evolutionary biology: understanding of natural selection, speciation, and mass extinctions.
- Genomics fundamentals: familiar with DNA sequencing, genome assembly, and comparative genomics.
- Bioinformatics tools: ability to use command-line software like BLAST, MAFFT, RAxML, or BEAST.
Data Requirements
- Newly sequenced genomes from modern squid (e.g., Doryteuthis pealeii) and cuttlefish (e.g., Sepia officinalis).
- Global datasets of cephalopod fossils and oceanographic conditions.
Step-by-Step Instructions
Step 1: Collect High-Quality Genomes
Start by obtaining raw sequencing data from multiple cephalopod species. Use long-read platforms (e.g., PacBio, Oxford Nanopore) for continuity. Assemble genomes using tools like Flye or Canu. Ensure coverage depth >30x for accurate assembly.
# Example assembly command using Canu
canu -p cephalopod -d assembly genomeSize=2.5g useGrid=false -pacbio-raw reads.fastq
Step 2: Align Sequences Across Species
Extract orthologous genes (e.g., conserved single-copy genes) from each genome. Use OrthoFinder to identify orthogroups. Then align each gene with MAFFT using the L-INS-i algorithm for accuracy.
# Align a single gene
mafft --localpair --maxiterate 1000 gene.fasta > aligned.fasta
Step 3: Build a Phylogenetic Tree
Concatenate aligned genes into a supermatrix (partitioned by gene). Use RAxML-NG with a GTR+G model and 100 bootstrap replicates to infer a maximum-likelihood tree. This reveals the evolutionary relationships among modern squid and cuttlefish.
# Run RAxML-NG
raxml-ng --msa supermatrix.phy --model GTR+G --prefix tree --threads 4
Step 4: Estimate Divergence Times with Molecular Clocks
Use a relaxed molecular clock in BEAST 2. Calibrate the tree with fossil-based priors, e.g., the oldest known decapodiform fossil (~200 million years ago). Run for 100 million generations, then combine traces with LogCombiner and analyze with Tracer.
<!-- BEAST XML configuration (abbreviated) -->
<beast>
<run spec='MCMC' chainLength='100000000'>
<state>
<parameter id='tree' value='@startingTree'/>
</state>
</run>
</beast>
Step 5: Correlate Diversification with Mass Extinctions
Plot the timeline of cephalopod branching events against known extinction dates (e.g., Cretaceous-Paleogene 66 million years ago). Use the phytools R package to generate lineage-through-time (LTT) plots. Look for sudden rate shifts – the original study found a explosion after the end-Cretaceous.
# R example
library(phytools)
tree <- read.tree('tree.nex')
ltt.plot(tree, log='y')
add extinction event line at 66 Mya
Step 6: Reconstruct Ancestral Habitats
Characterize habitat preferences of each species (e.g., deep vs. shallow, oxygen levels). Use BioGeoBEARS to model ancestral ranges. The results confirm that early squid and cuttlefish lived in well-oxygenated deep waters, and only later moved to shallow, more diverse niches after extinctions cleared the way.
# BioGeoBEARS analysis (pseudocode)
run Biogeography using DEC model with stratigraphic time slices
Common Mistakes
Over-reliance on Single-Gene Trees
Using only one gene (e.g., COI) can produce misleading topologies. Always use multiple loci or whole-genome data to reduce stochastic error.
Ignoring Rate Heterogeneity
Not accounting for variable substitution rates across lineages can skew divergence time estimates. Use relaxed clock models (e.g., uncorrelated lognormal).
Misinterpreting Fossil Calibrations
Place fossil calibrations incorrectly or without proper justification may artificially push dates. Always use multiple well-constrained fossils and cross-validate with geological data.
Overlooking Ecological Niche Conservatism
The deep-sea refuge hypothesis assumes squid ancestors were already deep-dwelling. If your analysis doesn't include extinct lineages, you might miss that they could have been shallow before retreating. Use ancestral state reconstruction to test this.
Summary
By sequencing and analyzing modern squid and cuttlefish genomes, scientists have traced their roots to a deep-ocean origin over 100 million years ago. These cephalopods survived mass extinctions by retreating to oxygen-rich deep-sea refuges, underwent little evolutionary change for millions of years, and then rapidly diversified into shallow waters after the end-Cretaceous event. This genomic detective work demonstrates how combining high-throughput sequencing, phylogenetic inference, and ecological modeling can illuminate deep-time evolutionary history. Future studies should expand the taxon sampling and incorporate proteomic data to refine the scenario.