Sticklebacks

Lake sticklebacks lose armour within one to two decades by reusing an Eda allele already present in the sea, and lose the pelvis through different, overlapping deletions of one fragile Pitx1 enhancer. No new adaptive mutation was watched arising. A measured mutation rate dates the reused allele to about two million generations, which young-earth models have not answered.

By AI (Claude Opus 5.5, thinking level unavailable) Draft:  September 2026

"Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment." — Lescak et al., PNAS, 20151

"The hypothesis does not address where or when an advantageous mutation arose" — Schluter and Conte on their own "transporter" model, PNAS, 20092

Summary

The threespine stickleback, Gasterosteus aculeatus, is a finger-length fish that lives in the sea and breeds in rivers.  Sea-run fish carry a full row of bony side plates and a spiny pelvis.  Thousands of lakes and streams hold resident freshwater populations with few plates and, in some lakes, no pelvis.  Sticklebacks are the best-documented case of adaptation happening fast in a wild vertebrate.  This article asks what was actually watched, what was inferred, and how each genetic change scores under the site's gain, loss, modification or transfer taxonomy.

The findings, in order of strength:

  1. Armour loss has been watched in real time, repeatedly. In Loberg Lake, Alaska, the low-plated form was absent in 1990 and made up 75% of the fish by 2001.3 In Cheney and Scout Lakes, stocked with about 3,000 sea-run fish each in 2009 and 2011, the freshwater Eda allele rose from under 1% to over 50% in eight years.4 4 Ponds lifted out of the sea by the 1964 Alaska earthquake hold freshwater-type fish that have diverged nearly as far as populations thought to be thousands of years old.1
  2. None of that change needed a new mutation. The freshwater Eda allele already sits in sea-run fish at 0.2% to 3.8%.5 In the three Alaskan lakes, 94% of the 344 regions that changed fastest were regions already known to differ between old freshwater and marine populations.4 No adaptive variant has been documented arising during any observed period. The study designs could not have detected one, so this is an absence of evidence, not a measured zero.
  3. The armour change is a shift in allele frequency. Its cause is a weakened switch, not a deleted one. The freshwater Eda allele carries one T to G change in an enhancer that lowers its response to a growth signal in armour and leaves it working elsewhere.6 That scores as a modification, down.
  4. The pelvis is lost by loss. Nine pelvic-reduced populations carry nine different deletions of one Pitx1 enhancer, overlapping in 488 base pairs.7 The region breaks easily because of a run of TG repeats.8 The deletions are independent and overlapping, not "the same deletions each time." Each scores as a loss.
  5. Stickleback adaptation is not all loss. Freshwater benthic fish grew more teeth by raising expression of Bmp6,9 and several freshwater populations carry extra copies of Fads2, which lets them make an essential fatty acid their diet lacks.10 Both are inferred from comparisons, and neither adds unique functional sequence: one strengthens an existing switch, the other copies an existing gene.
  6. The mainstream objection is strong, and the young-earth literature has not answered it. The standing variation had to come from somewhere. The Eda low and complete haplotypes differ at 1.9% of synonymous sites, more than two stickleback lineages dated at about 2 million years.5 A pedigree-measured stickleback mutation rate turns that difference into about 1.9 million generations.11 Producing it by mutation in 4,400 years would need a rate about 430 times the measured one.12 A young-earth model has to claim the difference was created, not accumulated. Answers in Genesis says exactly that, but no creationist source found tests it.

Verdict.  Sticklebacks are strong evidence for fast adaptation from variation that was already there, and that part fits the creation model of rapid speciation well.  They are weak evidence for the origin of anything new, because nothing new was seen arising.  But the young-earth reading has a hole of its own:  the reused variants look old by a measured clock, and the only published creationist answer, that they were created, has not been put to the test that could check it.

Part 1: What was watched

The membership rule of this series is that something must have been watched happening.  For sticklebacks, six things were.

Loberg Lake, 1982 to 2017

In 1982 the Alaska Department of Fish and Game poisoned Loberg Lake to clear it for sport fish.  Sea-run sticklebacks, probably from neighbouring Rabbit Slough, got in between 1983 and 1988.4  Michael Bell's group sampled the lake every year from 1990.  In 1990 the fish were "nearly monomorphic for the complete plate morph" and the low morph "was absent."  By 2001 complete morphs had fallen to 11% and low morphs had risen to 75%.3  By 2017 the armour-reduced form was about 95%, and the freshwater Eda allele was "approaching fixation in LB by 15 years."4 4

Bell and colleagues drew the conclusion that "threespine stickleback adapt to freshwater habitats within decades after invasion from the ocean."3

Cheney and Scout Lakes, 2009 to 2017

To watch the start of the process, Bell's group released about 3,000 Rabbit Slough fish into each of two poisoned lakes with no outlet:  Cheney Lake in 2009 and Scout Lake in 2011.  Low-plated fish appeared "in the second and third generation after founding" and reached 20 to 30% by 2017.4  Garrett Roberts Kingman and colleagues sequenced a pool of 100 fish from every year, plus 200 of the founders' source population.4  They found 344 regions whose allele frequencies moved fast, and fitted an average selection coefficient of 0.30 to them.4  The Eda allele that marks low armour rose "from an initial frequency of <1% to over 50% within 8 years."4

This is the one genome-wide time series of wild vertebrates adapting from the first generation, and it is the centre of this article.

Barrett's ponds, 2006 to 2007

Rowan Barrett, Sean Rogers and Dolph Schluter measured selection on Eda directly.  They trapped about 35,000 fish at Oyster Lagoon, British Columbia, kept the 354 with incomplete armour, and found 182 of them heterozygous at Eda.13  On 1 June 2006 they put 45 or 46 heterozygotes into each of four ponds and sampled 50 offspring from each pond ten times over a year.13

Two selection coefficients came out, and they point in opposite directions.  From October 2006 to May 2007, the complete allele fell from 67% to 49%, a coefficient of 0.52 against it.13  But from July to October 2006, while the fish were still small, selection ran the other way at 0.50 against the low allele, "which offset the gains occurring later in life."13  Across the whole generation the low allele went from 50% at birth to 51.6% in breeding adults, and the first sample of their offspring stood at 57.0%, a difference with P = 0.061.13

That tension matters, because the single number 0.52 circulates on its own.  Lescak and colleagues write that "The strength of selection (s) on the region containing Eda" in transplanted fish "was found to be 0.52," without the opposing 0.50.1  Barrett's own reading is more careful:  "countervailing selection pressures diminish the advantage of the low allele over the whole life span," and the fitness effects of Eda "do not seem to be determined solely by differences in lateral plate number."13  The honest figure for one generation is a small net gain for the low allele after two large and opposite pushes.  The same group's ponds also showed cold tolerance improving by 2.5 °C in three generations, a rate the authors call "among the most rapid to be observed in a natural population."14 14

Lake Washington, 1957 to 2005

Evolution in lakes does not only strip armour.  In Lake Washington, Seattle, only 6% of sticklebacks were completely plated in 1968 to 1969.  By 2005 the figure was 49%, with 35% partially plated and 16% low.15  Jun Kitano and colleagues tie the change to Eda and to trout that could see their prey once the polluted lake cleared in the early 1970s.  They estimate a selection coefficient of 0.58 to 0.72 for the plated form between 1969 and 1976.15  The plate counts were watched through museum and survey samples.  The coefficients are model fits, not measurements like Barrett's.

The earthquake islands, 1964 to 2012

On 27 March 1964 the Great Alaska Earthquake lifted Middleton Island 3.4 m and turned sea floor into freshwater ponds.1  Emily Lescak and colleagues genotyped more than 1,000 fish from three islands.  The freshwater fish came from the sea after 1964, more than once, and "have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago."1  The date of the ponds is known.  The colonisations themselves were not watched;  they are reconstructed from the genotypes.

What the watched cases share

Case Years Founders What changed Genetics
Loberg Lake 1983-88 to 2017 Unknown, sea-run Low morph 0% to 75% by 2001; armour-reduced about 95% by 2017 Eda allele near fixation
Cheney, Scout Lakes 2009, 2011 to 2017 About 3,000 each Low morph 20 to 30% Eda under 1% to over 50%; 344 regions moved
Barrett's ponds 2006 to 2007 45 or 46 per pond Opposing selection within one year Eda 50% to 51.6%
Barrett's ponds 2006 to 2008 45 or 46 per pond, three ponds Cold tolerance 2.5 °C better Not identified
Lake Washington 1957 to 2005 Resident Complete morph 6% to 49% Eda explains 75.2% of plate variance
Earthquake ponds 1964 to about 2012 Unknown, sea-run Freshwater form Genome-wide divergence

Every case is a change in how common an existing form is.  In every case where genes were checked, the variant that rose was already present when the watching began.  That is the fact both sides have to explain.

Part 2: Armour plates and Eda

What Colosimo and colleagues measured

Pamela Colosimo, David Kingsley and colleagues published the key paper in 2005.5  It did five separate things, and they should be kept apart.

  1. Mapped the armour gene. In a cross between a Japanese marine fish and a Paxton Lake freshwater fish, plate number mapped to Ectodysplasin (Eda), a signalling gene that in mammals is needed for hair, teeth and sweat glands. A mouse EDA transgene gave low-plated fish extra plates.
  2. Sequenced the difference. Against the marine sequence, the freshwater allele carried "numerous noncoding changes" and four amino-acid changes, none "at sites that are highly conserved between mammals and fish."5
  3. Showed the freshwater alleles are one family. Of 15 low-plated populations from across the Northern Hemisphere, every one but Nakagawa Creek in Japan carried alleles from one clade. Genes elsewhere in the genome did not group the low-plated populations together, so these are not one migrating lineage.
  4. Measured the allele in the sea. In completely plated marine fish, the low allele turned up at 3.8% at the Navarro River, California (8 of 218 chromosomes), and at 0.2% at the Little Campbell River, British Columbia (1 of 604).5 Two later samples fall in the same range. O'Brown found the key freshwater base in 12 of 263 Alaskan sea-run fish, 2.3%.6 Barrett's Oyster Lagoon trapping implies at least about 0.26%.16
  5. Estimated the allele's age. This step is an inference, and its method is set out below.

The conclusion follows from steps 3 and 4:  "worldwide evolution of the low-plated phenotype has occurred mainly by recurrent local selection for a family of alleles repeatedly brought into freshwater populations by marine founders."5  Lake populations did not invent low armour.  They inherited it from sea-run ancestors who carried it hidden, because heterozygotes are usually fully or partly plated.

How the age was estimated

Colosimo's team compared synonymous differences, the base changes that do not alter the protein and are assumed to build up at a steady rate.  The low and complete haplotypes differ at a fraction Ks = 0.01938 of synonymous sites.  That is more than the 0.01384 between Japanese Pacific and Sea of Japan sticklebacks, a pair the paper says was "thought to have diverged" about 2 million years ago "on the basis of genetic and geological estimates."  It is less than the 0.05875 between G. aculeatus and its sister species G. wheatlandi, dated at about ten million years "on the basis of the split between Pacific and Atlantic populations."  In the paper's words, "The common clade of low-plated alleles thus probably arose between 2 and 10 million years ago."5  The abstract rounds this to "more than two million years ago."5

The estimate is relative.  It depends on four things:

  1. The dates of the two calibration pairs, which come from geology and biogeography.
  2. A steady rate of synonymous change across the three comparisons.
  3. Synonymous sites being close to neutral.
  4. The difference having accumulated by mutation after the two haplotypes split, rather than being present from the start.

The third and fourth assumptions are where the two sides part, and Part 6 returns to them.  A second route to the age now exists that does not use either calibration pair:  a mutation rate measured in stickleback families.  Chaowei Zhang and colleagues sequenced 40 parent-offspring trios and counted 115 new mutations, a rate of 5.11 × 10⁻⁹ per base pair per generation.17  Dividing the Eda difference by twice that rate gives about 1.9 million generations, or 1.9 to 5.7 million years at one to three years a generation.11  The two routes agree.

What the freshwater allele does

Natasha O'Brown and colleagues found the change that matters in 2015.  In hybrids carrying one allele of each kind, the freshwater allele was expressed "almost fourfold lower" in the developing flanks and spines.  It was not lower in the fins or the lower jaw.6  Every low-plated population tested shares one base change, T to G, in an enhancer downstream of the gene.  Putting that one change into the marine enhancer "strongly reduces expression in posterior armor plates," because it "reduces the sensitivity of an EDA enhancer to Wnt signaling, and alters expression in developing armor plates while preserving expression in other tissues."6

The same study revised Colosimo's one exception.  The Nakagawa Creek fish carry the same T to G change on an otherwise marine-like background.  The authors offer two explanations and do not choose between them:  the same mutation arose twice, or a tiny piece of the usual freshwater haplotype was carried in by recombination.6

How the Eda change scores

The site's scoring procedure asks four questions in order.  Nothing arrived from another species.  No element stopped working:  the enhancer still drives expression in other tissues, and anterior plates still form.  No new element appeared.  So the T to G change is a modification, down.  The four amino-acid differences are modifications whose effect has not been measured.  Behe's scheme would reach the same verdict, since his modification-of-function covers a mutation that does "something else—either tweak an old FCT in a way that leaves it still working or affect some noncoded feature of a cell."18

Two claims on opposite sides fail this check.  Creation Ministries International wrote that "the mutated (i.e. broken) genetic switch simply prevents expression of a gene."19  It does not.  It weakens one input to one enhancer.  On the other side, Gregory Lang and Amber Rice cited Eda against Behe:  "none of the observed Eda mutations were predicted to damage the function of the protein."20  That is true, and it is not a gain.  Lower expression in armour is a reduction of an existing function.

What was watched in the lakes is simpler still.  An allele that already existed rose in frequency.  No sequence was gained or lost during the observation.

Part 3: The pelvis and Pitx1

Nine different deletions

Many freshwater sticklebacks have lost the pelvic girdle and its spines.  In crosses, most of that loss maps to the region of Pitx1, a gene that shapes the hind limb in vertebrates and is lethal when knocked out in mice.  Yingguang Frank Chan and colleagues found in 2010 that the gene itself is intact.  What is missing is a stretch of DNA upstream that switches Pitx1 on in the developing pelvis.7

Paxton Lake benthic fish carry a 1,868 base-pair deletion that removes this pelvic enhancer.  The Paxton version drove no pelvic expression in any of 8 transgenic fish.  Bear Paw Lake fish carry a different deletion of 757 base pairs, and Hump Lake fish a third of 973.  "The three different deletions in PAXB, BEPA and HUMP overlap in a 488bp region."7  A wider survey of 13 pelvic-reduced and 21 pelvic-complete populations found deletions in 9 of the reduced ones and none of the complete ones:  "nine different haplotypes with different staggered deletions, each consistently seen within a pelvic-reduced population."7  Putting the marine enhancer and Pitx1 back into pelvic-reduced fish restored pelvic structures.7  Four reduced populations had no detectable deletion, and three of those look as if they lost the pelvis another way.7

The catalogue once described this as "the same deletions each time."  That is wrong.  The deletions are different in size and endpoints, each is private to its own lake, and they overlap on one target.  Chan's team calls them "independent deletion mutations fixed in different pelvic-reduced stickleback populations."7  The Eda story is the opposite:  one old allele reused many times.

Why the same place breaks

Kathleen Xie and colleagues asked in 2019 why the same enhancer keeps being deleted.  In yeast artificial chromosomes, test chromosomes broke at a background rate of 3.37 per million divisions.  Chromosomes carrying the marine pelvic enhancer "broke ~25–50 times more frequently," a rate "even higher than previously analyzed human fragile sites."  Versions from pelvic-reduced lakes broke at the background rate, "suggesting that natural Pel mutations remove breakage-prone regions."8  The cause is about 170 base pairs of TG repeats.8

Their model is the key point for the ledger.  "Population genetic modeling indicates that new mutations occurring at the low rates of typical single nucleotide changes (~10⁻⁹) would rarely arise at a particular locus in postglacial stickleback populations, while mutations occurring at elevated rates (~10⁻⁵ for fragile sites) would arise often."8  So the Pitx1 deletions are the one class of stickleback adaptive change that the mainstream literature treats as new mutations arising inside young lake populations.  The breakage rates were measured in yeast and monkey cells.  The deletions in wild fish were not watched arising.

How the Pitx1 change scores

A coded element that worked, the pelvic enhancer, stopped working, and nothing new took its place.  That is a loss under the site's procedure and under Behe's.  It is also adaptive:  Chan found reduced diversity and an excess of derived alleles centred on the enhancer in pelvic-reduced lakes, the signature of selection.7  The sibling article Measuring Biological Information scores the same case the same way.  The fish gained fitness by losing a part.  Both facts are true, and neither cancels the other.

Part 4: The genome-wide picture, and the search for a gain

Old variants, mostly regulatory

Felicity Jones and colleagues sequenced a reference genome and 20 more fish from around the Northern Hemisphere in 2012.  About 0.5% of the genome groups fish by habitat rather than by geography:  242 regions by either of two methods, 147 by both.21  Three are chromosome inversions.  Of the 64 regions with the strongest evidence, 41% lie entirely outside genes, 17% change amino acids, and 43% contain genes but no habitat-specific amino-acid change.21  The authors conclude that "regulatory changes appear to predominate."21

"Regulatory" here means "not an amino-acid change."  It does not mean deleted.  Most of these regions are alternative versions of the same stretch of DNA, differing at many single bases.  The catalogue's phrase for this article, adaptive change "that works by removing regulatory DNA rather than adding any," fits Pitx1 and does not fit the genome as a whole.  Answers in Genesis summarised Jones fairly as "about 80% of the genetic differences arise from alterations in regulatory DNA and about 20% in actual coding DNA."22

Thomas Nelson and William Cresko dated these regions in 2018.  Using the ninespine stickleback as an outgroup and fixing the split between the two genera at 15 million years, they found the habitat-linked variation "averages six million years old, nearly twice the genome‐wide average."23 23  In some regions, marine and freshwater threespine sticklebacks "are as divergent as threespine and ninespine stickleback, which are classified into separate genera."23  Roberts Kingman's team, using three calibration species, put the reused regions at 1,600 thousand years against 700 thousand for the rest of the genome, "the majority (161 of 209)" at over a million years.4 4  The two papers' absolute ages differ severalfold.  Both find the adaptive regions about 1.6 to 2.3 times older than the genomic background.24

Candidates for a gain

The series rule is to look hard for a gain before saying there is none.  Four candidates turn up.

More teeth.  Freshwater benthic sticklebacks from Paxton Lake have more throat teeth than marine fish.  Phillip Cleves and colleagues mapped the largest effect to Bmp6, found no amino-acid differences, and found the freshwater allele expressed more strongly in developing tooth plates.9 9  They present it as a rare case of "how morphological features may be gained in natural species."  At the level of the organism it is a gain of structures.  At the level of the molecule, an existing gene is switched on more strongly:  modification, up.  The causal bases were not identified in that paper.

More copies of a fat gene.  Freshwater food is poor in DHA, a fatty acid fish need.  Asano Ishikawa and colleagues found that several freshwater populations carry extra copies of Fads2, a gene for making DHA, duplicated on the X chromosome, and that "fish with higher Fads2 copy number better synthesized DHA when fed with DHA-deficient diets."10  This is the clearest adaptive gain of function in the stickleback record.  Under the scoring procedure an element "present in more copies" is a modification, up.  It adds no unique functional sequence, because a duplicate is a copy.  Behe scores duplication the same way:  "no new functional coded element was gained or lost, just copied."18  The 2019 Science paper that reported the duplications was not obtained;  the summary here is the authors' own, from 2021.

Haplotypes that grew.  Roberts Kingman's team found that older reused regions are larger, and that the oldest bases sit at their centres.  Their reading:  "alleles originally based on an initial beneficial mutation accumulating additional linked favorable mutations, snowballing over time to form a finely tuned haplotype with multiple adaptive changes."4  If that is right, these haplotypes were built by mutation and selection, and they would be real gains of sequence-specific DNA.  It is inferred from age patterns over a million years.  Nothing like it was watched.

The Eda T to G itself.  On the mainstream account it arose by mutation once, or twice, long ago.  It scores as a modification that lowers activity, not a gain.

So nothing in the stickleback record is a newly arisen gain, if "newly arisen" means watched.  If it means inferred, the best candidates are the Fads2 copies, which are copies, and the snowballing haplotypes, which are a reconstruction.

The scored record

Change What changed Class Behe's scheme New sequence Watched?
Eda allele frequency in lakes Old allele from under 1% to near fixation No mutation; frequency shift Not a mutation 0 Observed
Eda T to G One base; enhancer answers Wnt less in armour Modification, down Modification 1 Inferred
Eda amino acids Four substitutions at unconserved sites Modification, unmeasured Modification 4 Inferred
Pitx1 enhancer Nine different deletions of one enhancer Loss Loss Negative Inferred
Bmp6 Stronger expression in tooth plates; more teeth Modification, up Modification Not identified Inferred
Fads2 Extra gene copies; more DHA made Modification, up, by copy Modification 0 Inferred
Lake Washington Eda Complete allele rose again No mutation; frequency shift Not a mutation 0 Observed
Cold tolerance 2.5 °C better in three generations Not scored; genes unknown Not scored Unknown Observed

"New sequence" counts bases changed by a mutation that are not a copy, following Measuring Biological Information.  The two observed genetic rows have zero because nothing new arose while anyone watched.

Part 5: The ledger

The series runs the same ledger for every organism:  individuals times generations times mutation rate gives the number of new mutations the population tried.  Then it lists every documented adaptive gain.  For sticklebacks the ledger can be run for the observed lakes, and it gives a clear answer about why standing variation matters.  It cannot be run for output.

Mutations tried in the watched lakes

Input Value Source
Founders, Cheney and Scout Lakes About 3,000 each Roberts Kingman et al. 20214
Years to 2017 8 and 6 Same
Genome size 463 Mb Jones et al. 201221
Mutation rate 5.11 × 10⁻⁹ per base pair per generation Zhang et al. 202517
New mutations per fish About 4.7 2 × 463 Mb × 5.11 × 10⁻⁹
Fish-generations, floor About 42,000 Founders × years, if the lakes never grew and bred every year
New mutations tried, floor About 2 × 10⁵ 42,000 × 4.7
Chance of one given base change arising About 1 in 7,000 84,000 gametes × 5.11 × 10⁻⁹ ÷ 3
Expected fragile-site deletions About 0.8 84,000 gametes × 10⁻⁵ (Xie et al. 2019)

The arithmetic is in the note.25  The populations surely grew past their founding size, so the true number of mutations tried is larger, perhaps by one or two orders of magnitude.  That does not change the conclusion.  A specific single-base change, such as the Eda T to G, had about one chance in 7,000 of arising anew in these lakes in the time observed.  It rose from under 1% to over 50% in eight years because it was already there.  Xie's team makes the same argument for postglacial lakes.8  Both sides agree on this point.  Fast adaptation in small lakes runs on standing variation, except at a fragile site, where a deletion could appear about once in a decade.

Adaptive variants that arose during observation

Observed system Adaptive variants tracked Already present at the start Shown to arise during observation
Cheney, Scout, Loberg 344 fast-moving regions 323 overlap known old regions (94%); the rest not shown to be new 0
Barrett's ponds Eda Founders chosen as carriers 0
Lake Washington Eda Present at 6% complete morphs in 1968 to 1969 0
Earthquake ponds Genome-wide Attributed to marine standing variation 0

The zeros need a caution.  A new mutation starts as one copy among thousands.  Pooled sequencing of 100 fish cannot see anything below a few percent, and a new allele could not climb that high in a decade unless its advantage were enormous.  These studies were designed to follow old variants, not to find new ones.  So the right entry is "none documented," not "none happened."

Gains that cannot be priced

The candidate gains in Part 4, the Bmp6 teeth, the Fads2 copies and the snowballing haplotypes, happened before anyone was watching and in populations of unknown size.  They cannot be put on this ledger at all.  The series treats such gains as a reason the ledger is a lower bound on what evolution produced, which cuts against the young-earth side.  Here the effect is small, because none of the three adds unique functional sequence in the observed period, but the principle stands.

Unique functional sequence per organism-generation

This cell is left empty, for two reasons.  The numerator from the watched lakes is zero, but it is not a measured zero, for the reason above.  The inferred numerator, the sequence built into old haplotypes over a million years or more, has never been counted in functional nucleotides.  The EcoPeaks are regions tens of thousands of bases long, and the bases that matter inside them are known for only one, the Eda T to G.  No denominator is offered either.  The long-term effective size of sea-run sticklebacks is about 160,000 by the standard formula,23 26 but the census population is far larger and unmeasured.  An invented denominator would sink the argument faster than a missing one.

Part 6: What each side makes of it

The creationist and design reading

Creation Ministries International's 2009 article is the fullest.  Its central claim is that "the so-called “rapid evolution” did not result in the addition of any new features but rather the reduction or even complete loss of existing features—bony plates and pelvic spines."19  Answers in Genesis goes further in 2012:  "The “ancient genes” the fish seems to switch on and off are the genetic information God provided in order to enable this fish to adapt to various environments."22  And in 2010 it used the pond experiments for the young-earth timescale:  "This three-year example of stickleback “evolution” is a soft echo of the adaptation and speciation that would have occurred in the centuries after the Ark landed."27  The Institute for Creation Research's current model, Continuous Environmental Tracking, holds that "sensors trigger an organism’s pre-programmed logic to select pre-existing information designed to respond to the changed condition," and cites sticklebacks as a case of traits "produced “on the fly” by changes to the “instruction code”" of DNA.28 28

On speed and on standing variation, this reading is right, and the mainstream papers say the same.  Lescak's team concluded that "most stickleback evolution in fresh water occurs within the first few decades."1  Roberts Kingman's team calls "the ancient SGV" (standing genetic variation) "the primary genomic feature enabling extremely fast evolution."4  A post-Flood model needs freshwater forms to appear within centuries.  Sticklebacks show them appearing within decades.

The creationist sources also make errors a skeptic will find.  CMI attributes armour loss to Pitx1:  "Specifically, it affects expression of the Pitx gene."19  Armour is Eda;  Pitx1 is the pelvis.  CMI and Answers in Genesis both give 200 founders for Barrett's ponds;  the paper had 182.19 29  Brian Thomas of ICR wrote that the sequences Jones compared "were very similar between all the fishes studied, showing that they have been resistant to mutational changes,"30 when Jones had found the opposite.  In 2010 Thomas described the cold-tolerant fish as descendants of "rare individuals" and then called "the reference to "selection" superfluous."31  Survival of rare existing variants is what selection means.

The mainstream objection, at full strength

Standing variation is not free.  Every allele that a lake population reuses had to arise by mutation at some earlier time, in some population.  Schluter and Conte's transporter model explains how freshwater alleles are kept in the sea and carried into new lakes, and they say plainly that it "does not address where or when an advantageous mutation arose."2  Their answer to the "when" is that "the original low-armor mutation occurred long ago."2  The measured numbers back them.  The Eda haplotypes differ by more than a pair of populations dated at about 2 million years.  The habitat-linked regions are older than the genome around them.  A measured family mutation rate gives about 1.9 million generations for Eda.  And measuring the rate directly made stickleback dates older, not younger:  the family rate is "3.1 times smaller than the commonly used substitution rate," which moved the threespine-ninespine split from 5.6 to 17.5 million years.32 32  This is the opposite of the human mitochondrial case, where family rates run faster than phylogenetic ones and creationists have used the gap to shorten dates.

On this view, "no new mutation in the lakes" is exactly what the mainstream predicts.  The new mutations happened over millions of years.  The lakes show only the harvest.

What the age depends on

The pedigree route to the Eda age rests on four inputs:  the difference between the haplotypes, 0.01938 per synonymous site;  the family mutation rate, 5.11 × 10⁻⁹;  the generation time, one to three years;  and the assumption that the difference built up by mutation after the haplotypes split.  The first two are measurements.  The third moves the answer threefold and no more.  The fourth is the whole question.

What a young-earth model has to claim

A model with about 4,400 years since the Flood has two options.

Option 1:  the difference arose by mutation, faster.  Producing a 1.9% synonymous difference between two lineages in 4,400 one-year generations needs a rate of about 2.2 × 10⁻⁶ per base per generation, 431 times the measured rate, and 1,293 times with three-year generations.12  Applied to the whole genome, that is about 2,000 new mutations per fish per generation instead of about 4.7.  No creationist source found proposes this.  The numbers are here for scale only.

Option 2:  the difference was created.  This is what Answers in Genesis says:  the ancient genes are "the genetic information God provided."22  On this option, the created stickleback kind carried both marine and freshwater haplotypes from the start, and the ages are artefacts of reading created differences as accumulated mutations.  This option is coherent.  It is also untested.

What would test each mechanism

Continuous Environmental Tracking has been tested at Eda, and failed there.  If fish changed their own DNA "on the fly," genotypes would shift toward the freshwater allele within individuals or at conception.  Barrett's F1 fish were born at the Mendelian 1:2:1 ratio in all four ponds.13  The frequency then moved because some genotypes died faster than others.  Roberts Kingman's team found the low-armour phenotype tracking the square of the allele frequency, the ordinary signature of a recessive allele.4  CET may still apply to traits that are plastic.  Armour plates are not one of them.

Created versus accumulated differences can be told apart.  Mutations leave a fingerprint.  Some kinds of base change are far commoner than others, and Zhang's team has now measured the stickleback pattern in families.17  O'Brown's team noted that the Eda T to G is a transversion, "one of the least common types of changes seen in large scale studies of spontaneous germ-line mutations."6  If the thousands of differences between marine and freshwater haplotypes carry the same fingerprint as new mutations, that supports accumulation.  A created set of differences has no reason to carry it.  The public genomes from Jones and Roberts Kingman would allow the test.  No study found has run it, on either side.

The age of the reused variants can be checked against the kind.  If the stickleback family was created as one kind, front-loaded haplotypes should be shared across its genera.  Nelson and Cresko note that genetic mapping of armour traits points "either to independent genetic architectures in the two species" or to repeated mutation at one locus.  That leans against a shared created pool, though it rests on a small number of mapping studies.23

Objections and responses

"Sticklebacks are evolution in action.  Why call it anything else?"  They are, in the population-genetic sense of a change in allele frequency.  Answers in Genesis concedes the definition:  "the stickleback population’s loss of armor was evolution, and now the reappearance of the armor is evolution."33  What was watched is the spread of existing alleles.  It is not evidence for how those alleles began.

"The standing variation had to come from somewhere, and the dates say millions of years."  This is the strongest objection, and Part 6 states it at full strength.  The young-earth answer is that the variation was created.  That answer is coherent but has not been tested, and a test is available.  Until someone compares the marine-freshwater differences with the measured mutation fingerprint, this objection stands unanswered.

"The Pitx1 deletions show evolution happening by new mutation, right now."  They show new mutation happening easily at a fragile site, and every such mutation is a loss.  They were inferred from lake populations, not watched.  They fit the claim that losses are cheap and gains are expensive.  They do not show new function arising.

"Creationists say it is all loss, but more teeth and extra Fads2 copies are gains."  At the level of the fish, yes, and this article says so.  At the level of the molecule, one is a stronger existing switch and the other is a copy.  Neither adds unique functional sequence.  CMI's line that there is "no new genetic information here"19 is right about the Pitx1 enhancer and wrong as a description of stickleback adaptation as a whole.

"Barrett measured a selection coefficient of 0.52.  That is enormous."  It was one half of a pair.  The low allele was selected against at 0.50 earlier in life.  Across one generation the net change was small.  The large figure describes part of a life cycle, not a generation.

"Fast change within decades proves the young-earth timescale."  It shows the timescale is not needed for this change.  It says nothing about the age of the variants that changed, and the one measured clock available says those are old.  The speed of the lake change and the age of the alleles are separate questions.

"The mutation rate may have been higher in the past."  It may have varied, and Zhang's team says so.  But the required rate is 300 to 1,300 times today's.  At that rate each fish would carry about 2,000 new mutations per generation.  A proposal of that size needs its own evidence, and none has been offered.

"Old dates assume common descent."  The pedigree route assumes only that the two Eda haplotypes share an ancestor and that their differences arose by mutation.  Within one interbreeding species, common ancestry of two versions of one gene is not in dispute.  The second assumption is the real point at issue.

Issues and unanswered questions

  • Loberg Lake population size and Bell 2004 full text. Only the abstract was read. The ledger uses Cheney and Scout founders instead.
  • Ishikawa et al. 2019 (Science 364:886) not obtained. The Fads2 account rests on the authors' 2021 summary. Whether the extra copies sit under new control, which would make them a gain by copy rather than a modification, is not known here.
  • The mutation-fingerprint test has not been run. It would decide between created and accumulated haplotype differences. It needs published genomes and the spectrum from Zhang 2025.
  • What the 6% of fast-moving regions outside known old regions are. Roberts Kingman does not say they are new. Nor does it say they are not.
  • Behe has not written on sticklebacks. Darwin Devolves does not mention them. The Behe scores in the table are this article's application of his published definitions.
  • The generation time of sea-run sticklebacks. Zhang assumed three years. Lake fish often breed at one. The Eda age moves threefold with it.

Verdict

Sticklebacks are a strong case for rapid adaptation from standing variation, and a weak case for the origin of new function.  Everything watched happening was the rise of alleles already present.  In the watched lakes, no adaptive mutation was seen arising, though the studies could not have seen one.  Of the inferred changes, the pelvis was lost by deletions, the armour by a weakened switch, and the two gains of function are a stronger switch and a copied gene.  That record supports the claim that fast adaptation mostly spends existing information.

The young-earth side should not overclaim it.  The variants being spent look old by a measured mutation rate, not only by an assumed calibration.  The one published creationist answer, that the differences were created, is coherent.  It has not been put to the test that could separate it from accumulation.  Until it is, this is the main open weakness of the creation model of rapid speciation in its best example.

Sources


  1. 1.

    Lescak, Emily A. et al.  “Evolution of Stickleback in 50 Years on Earthquake-Uplifted Islands.”  PNAS 112(52):E7204–E7212.  2015.  Read from the publisher PDF via the Wayback Machine.  Mirrors:  Archive.org | PMC

    1. 1a.

      Page E7204:  “We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.” ↩ ↩ ↩ ↩

    2. 1b.

      Page E7205:  “Middleton Island was uplifted 3.4 m, creating a new terrace with ponds from a previously submarine platform (1).” ↩

    3. 1c.

      Page E7209:  “The strength of selection (s) on the region containing Eda—the locus associated with lateral plate loss (35, 36, 70)—in oceanic stickleback transplanted to artificial freshwater ponds was found to be 0.52 (47, 69, 71).” ↩

  2. 2.

    Schluter, Dolph, and Gina L. Conte.  “Genetics and Ecological Speciation.”  PNAS 106(Suppl. 1):9955–9962.  2009.  Read from the PMC copy, which has no page numbers, so claims give section headings.  Mirror:  PMC

    1. 2a.

      Section Ecological Speciation from Standing Variation in Stickleback:  “The hypothesis does not address where or when an advantageous mutation arose, but proposes that once a mutation reaches appreciable frequency in a freshwater population it could participate in the transporter process.” ↩ ↩

    2. 2b.

      Section Ecological Speciation from Standing Variation in Stickleback:  “The origin of this clade predates the ages of all known freshwater populations, in most cases by orders of magnitude, suggesting that the original low-armor mutation occurred long ago.” ↩

  3. 3.

    Bell, Michael A. et al.  “Twelve Years of Contemporary Armor Evolution in a Threespine Stickleback Population.”  Evolution 58(4):814–824.  2004.  Abstract only, read from the Oxford Academic page and checked against PubMed 15154557.  The full text was not obtained.

    1. 3a.

      Section Abstract:  “The 1990 sample was nearly monomorphic for the complete plate morph, which is monomorphic in local sea-run populations; the low plate morph, which is usually monomorphic in local freshwater populations, was absent. By 2001, the frequency of completes had declined to 11%, and lows had increased to 75%.”  “Contemporary evolution in this population indicates that threespine stickleback adapt to freshwater habitats within decades after invasion from the ocean”. ↩ ↩ ↩

  4. 4.

    Roberts Kingman, Garrett A. et al.  “Predicting Future from Past:  The Genomic Basis of Recurrent and Rapid Stickleback Evolution.”  Science Advances 7(25).  Jun 18, 2021.  Read from the PMC copy, which has no page numbers, so claims give section headings.  CC BY-NC.  The supplement was not read.  Article number eabg5285.  Mirror:  PMC

    1. 4a.

      Section Rapid contemporary evolution and TempoPeak identification:  “In 1982, stickleback in Loberg Lake (LB) were exterminated to improve recreational fishing (17). Sometime between 1983 and 1988, LB was invaded by completely plated (~33 plates per side) anadromous stickleback [most likely from neighboring Rabbit Slough (RS)]. The characteristic freshwater, armor-reduced phenotype increased rapidly from ~16% in 1991 to ~50% by 1995 and to ~95% by 2017”. ↩ ↩

    2. 4b.

      Section Rapid contemporary evolution and TempoPeak identification:  “Bell et al. (47) introduced ~3000 anadromous RS fish into each of two other Cook Inlet lakes without outlets that had been similarly treated to exterminate fish: Cheney Lake (CH) in 2009 and Scout Lake (SC) in 2011. Low–armor-plated (~5 to 7 plates per side) stickleback began to appear in the second and third generation after founding in CH and SC respectively, and, by 2017, they had increased to 20 to 30% (Fig. 2B).” ↩ ↩ ↩ ↩

    3. 4c.

      Section Rapid contemporary evolution and TempoPeak identification:  “Each freshwater pool-seq experiment consisted of 100 individuals (with three exceptions), with mean coverage of 223× per pool. In addition, we resequenced a pool of 200 anadromous RS individuals used to found the CH population in 2009 (RS2009) to 585×.” ↩

    4. 4d.

      Section Rapid contemporary evolution and TempoPeak identification:  “we identified 524 sensitive and 344 specific TempoPeaks. Despite operating over very different time spans, the visual correspondence between the Pacific EcoPeaks in long-established populations and the TempoPeaks in recently established populations is notable, particularly for the specific TempoPeaks, of which 323 of 344 (94%) overlap with the sensitive Pacific EcoPeaks”. ↩

    5. 4e.

      Section Rapid contemporary evolution and TempoPeak identification:  “This marked congruity strongly suggests that the ancient SGV represented by Pacific EcoPeaks is the primary genomic feature enabling extremely fast evolution of freshwater phenotypes in stickleback from the northeast Pacific basin.” ↩

    6. 4f.

      Section Rapid contemporary evolution and TempoPeak identification:  “In both CH and SC, the G allele increases rapidly from an initial frequency of <1% to over 50% within 8 years, while approaching fixation in LB by 15 years.” ↩ ↩ ↩

    7. 4g.

      Section Rapid contemporary evolution and TempoPeak identification:  “Notably, the square of G-allele frequencies (i.e., the expected number of GG homozygotes) tracks closely with frequencies of the low–armor plate phenotype, consistent with almost complete recessiveness (h = 0.0) for the G allele for this phenotype (Fig. 2B).”  To fit the speed of the rise the authors nonetheless had to assume h = 1.0 and s = 0.55, and say so. ↩

    8. 4h.

      Section Rapid contemporary evolution and TempoPeak identification:  “When modeling each peak SNP as independent, we find an extremely high mean s of 0.30 (5th, 95th percentile 0.08 to 0.53) and h of 0.98 (5th, 95th percentile 0.95 to 1.0) for the 344 specific TempoPeaks found in CH + SC + LB.” ↩

    9. 4i.

      Section Features associated with EcoPeak evolution:  “Using Gasterosteus nipponicus, Gasterosteus wheatlandi, and Pungitius pungitius for calibration, we estimated molecular divergence time between a pair of freshwater (Little Campbell upstream) and marine (Little Campbell downstream) stickleback in windows tiled across the genome (section S11). We find that EcoPeaks as a whole are significantly older than the rest of the genome [1600 thousand years (ka) versus 700 ka”. ↩ ↩

    10. 4j.

      Section Features associated with EcoPeak evolution:  “We estimate that the majority (161 of 209) are over a million years old and have cycled between freshwater and marine environments many times during this long history”. ↩

    11. 4k.

      Section Features associated with EcoPeak evolution:  “This suggests that individual regions may grow over time, with alleles originally based on an initial beneficial mutation accumulating additional linked favorable mutations, snowballing over time to form a finely tuned haplotype with multiple adaptive changes.” ↩

  5. 5.

    Colosimo, Pamela F. et al.  “Widespread Parallel Evolution in Sticklebacks by Repeated Fixation of Ectodysplasin Alleles.”  Science 307(5717):1928–1933.  Mar 25, 2005.  Read from the Schluter lab reprint.  Pages anchored by extracting each PDF page and reading its folio.  Mirror:  Zoology.ubc.ca ↩

    1. 5a.

      Page 1928:  “parallel evolution of stickleback low-plated phenotypes at most freshwater locations around the world has occurred by repeated selection of Eda alleles derived from an ancestral low-plated haplotype that first appeared more than two million years ago.” ↩

    2. 5b.

      Page 1929:  “Comparison to the sequence of the completely plated SRMA marine population identified numerous noncoding changes (table S2), and four mutations that lead to amino acid changes in EDA protein (Fig. 2). None of the amino acid changes occur at sites that are highly conserved between mammals and fish or in residues that have previously been associated with defects in humans (25, 26).” ↩

    3. 5c.

      Page 1930:  “Eight of 109 completely plated fish from NAV were heterozygous for complete and low-morph alleles at both markers, which gave an estimated frequency of low-morph alleles of 3.8% (8/218 chromosomes). One of 302 completely plated marine sticklebacks collected from LITC was also heterozygous at both markers and gave an estimated allele frequency of 0.2% (1/604 chromosomes).” ↩ ↩

    4. 5d.

      Page 1930:  “suggest that worldwide evolution of the low-plated phenotype has occurred mainly by recurrent local selection for a family of alleles repeatedly brought into freshwater populations by marine founders.” ↩

    5. 5e.

      Page 1930:  “To estimate the age of the shared haplotype, we examined the rate of neutral sequence evolution at third–base pair positions from all genes (EDA to Garp) in the shared region.”  Pages 1930–1931:  “thought to have diverged”.  “on the basis of genetic and geological estimates”.  “on the basis of the split between Pacific and Atlantic populations”.  Page 1931:  “(Ks values of 0.01938, 0.01384, and 0.05875, respectively). The common clade of low-plated alleles thus probably arose between 2 and 10 million years ago.”  The printed approximately-equal sign before each calibration date does not survive text extraction, so the calibration sentences are quoted in pieces. ↩ ↩ ↩

  6. 6.

    O'Brown, Natasha M. et al.  “A Recurrent Regulatory Change Underlying Altered Expression and Wnt Response of the Stickleback Armor Plates Gene EDA.”  eLife 4:e05290.  Jan 28, 2015.  Read from the PMC copy.  CC BY.  Mirror:  PMC

    1. 6a.

      Section Abstract:  “Recreation of the T → G change in a marine enhancer strongly reduces expression in posterior armor plates.”  “Thus parallel evolution of low-plated sticklebacks has occurred through a shared DNA regulatory change, which reduces the sensitivity of an EDA enhancer to Wnt signaling, and alters expression in developing armor plates while preserving expression in other tissues.” ↩ ↩

    2. 6b.

      Section Results:  “the freshwater EDA allele was expressed almost fourfold lower than the marine allele in the developing anterior and posterior flanks”. ↩

    3. 6c.

      Section Discussion:  “Characteristic flanking SNPs are not shared between NAKA and other low-plated populations, suggesting that the same T → G mutation has likely occurred independently on two very different haplotypes.”  “It is possible that the shared T → G change arose not by independent mutation, but by extensive recombination or gene conversion from the typical freshwater EDA haplotype.” ↩

    4. 6d.

      Section Discussion:  “the recurrent regulatory mutation we have identified at the stickleback EDA locus is a T → G transversion substitution, one of the least common types of changes seen in large scale studies of spontaneous germ-line mutations in humans”. ↩

    5. 6e.

      Section Discussion:  “we have recently surveyed 263 migratory marine sticklebacks from Alaska and identified 12 completely plated individuals that are heterozygous for the T → G change in the EDA enhancer (minor allele frequency 2.3%).” ↩

  7. 7.

    Chan, Yingguang Frank et al.  “Adaptive Evolution of Pelvic Reduction in Sticklebacks by Recurrent Deletion of a Pitx1 Enhancer.”  Science 327(5963):302–305.  2010.  Read from the PMC author manuscript, which has no page numbers, so claims give its section headings.  Mirror:  PMC

    1. 7a.

      Section Abstract:  “Pelvic loss in different natural populations of threespine stickleback fish has occurred by regulatory mutations deleting a tissue-specific enhancer of the Pituitary homeobox transcription factor 1 (Pitx1) gene.” ↩

    2. 7b.

      Section Nature of mutations in pelvic-reduced fish:  “Bacterial Artificial Chromosome sequencing from the PAXB population identified a 1,868bp deletion present in the Pel-2.5kb region (Fig. S7). We cloned the PAXB deleted variant and found that it no longer drove expression in the developing pelvis (0 out of 8 transgenic animals; Fig. 2D, G)”.  “The three different deletions in PAXB, BEPA and HUMP overlap in a 488bp region”. ↩

    3. 7c.

      Section Nature of mutations in pelvic-reduced fish:  “Nine of the 13 pelvic-reduced stickleback populations—but zero out of 21 pelvic-complete populations—showed consistent missing genotypes for multiple consecutive SNP markers located in and around the Pel enhancer”.  “The larger genotyping survey identified a total of nine different haplotypes with different staggered deletions, each consistently seen within a pelvic-reduced population, and each overlapping or completely removing the Pel enhancer region (Figs. 4 and S8).” ↩ ↩

    4. 7d.

      Section Signatures of selection:  “Patterns of allelic variation showed an excess of derived alleles near the Pel enhancer region of pelvic-reduced populations, as indicated by negative values of Fay & Wu’s H statistic”. ↩

    5. 7e.

      Section Discussion:  “Strikingly, re-introduction of the enhancer and Pitx1 coding region can restore formation of pelvic structures in derived populations that appear to be monomorphic for pelvic reduction.” ↩

    6. 7f.

      Section Discussion:  “Smaller functional lesions might be found in some pelvic-reduced populations, including four populations without obvious deletions. However, three of these populations show unusual morphological features, suggesting that their pelvic loss may have occurred through non-Pitx1-mediated mechanisms (8, 26).” ↩

    7. 7g.

      Section Discussion:  “We hypothesize that sequence features in the Pitx1 locus may predispose the locus to structural changes, possibly explaining the high prevalence of independent deletion mutations fixed in different pelvic-reduced stickleback populations.” ↩

  8. 8.

    Xie, Kathleen T. et al.  “DNA Fragility in the Parallel Evolution of Pelvic Reduction in Stickleback Fish.”  Science 363(6422):81–84.  Jan 4, 2019.  Read from the PMC author manuscript.  Exponents checked in the HTML because text extraction loses superscripts.  Mirror:  PMC

    1. 8a.

      Section Abstract:  “Pitx1 enhancer sequences form alternative DNA structures in vitro, and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by (TG)-dinucleotide repeats.” ↩

    2. 8b.

      Section Main text, paragraph on breakage in yeast artificial chromosomes:  “Constructs without added test regions broke at background rates of 3.37 breaks per 10⁶ divisions (Fig. 2B), consistent with previous reports (14). Chromosomes containing marine Pel broke ~25–50 times more frequently (Fig. 2B), a rate even higher than previously analyzed human fragile sites (14). Pel from freshwater pelvic-reduced populations (but not freshwater pelvic-complete populations (fig. S1)) broke at rates similar to the control (Fig. 2B), suggesting that natural Pel mutations remove breakage-prone regions.” ↩

    3. 8c.

      Section Main text, paragraph on (TG) repeats:  “Three stretches of ~15, ~20, and ~50 (TG)-dinucleotide repeats in marine Pel total ~170 bp”. ↩

    4. 8d.

      Section Main text, paragraph beginning Could elevated mutation rates:  “Population genetic modeling indicates that new mutations occurring at the low rates of typical single nucleotide changes (~10⁻⁹) would rarely arise at a particular locus in postglacial stickleback populations, while mutations occurring at elevated rates (~10⁻⁵ for fragile sites) would arise often.” ↩ ↩ ↩

  9. 9.

    Cleves, Phillip A. et al.  “Evolved Tooth Gain in Sticklebacks Is Associated with a cis-Regulatory Allele of Bmp6.”  PNAS 111(38):13912–13917.  2014.  Read from the PMC copy.  Mirror:  PMC

    1. 9a.

      Section Significance:  “there are few well worked out examples of how morphological features may be gained in natural species (constructive evolution).”  “Genetic mapping and gene expression experiments suggest regulatory changes have occurred in the gene for a bone morphogenetic signaling molecule, leading to increased expression in the freshwater fish that have more teeth.” ↩ ↩

    2. 9b.

      Section Results:  “We sequenced the exons of Bmp6 in marine and benthic fish and found no nonsynonymous coding differences (SI Appendix, Fig. S7).” ↩

  10. 10.

    Ishikawa, Asano et al.  “Copy Number Variation of a Fatty Acid Desaturase Gene Fads2 Associated with Ecological Divergence in Freshwater Stickleback Populations.”  Biology Letters 17(8):20210204.  2021.  Read from the PMC copy.  Used for its summary of Ishikawa et al. 2019, Science 364:886, which was not obtained.  Mirror:  PMC

    1. 10a.

      Section 1. Introduction:  “Marine populations generally possess one Fads2 copy on chromosome 12 and another copy on the X chromosome (chromosome 19), whereas duplication of Fads2 occurred on the X chromosome in several freshwater populations [6]. Accordingly, fish with higher Fads2 copy number better synthesized DHA when fed with DHA-deficient diets, suggesting higher Fads2 copy number is adaptive in DHA-deficient freshwater habitats.” ↩ ↩

  11. 11.

    Colosimo's synonymous difference between the complete and low Eda haplotypes is Ks = 0.01938.5  Two lineages that split t generations ago each collect mutations, so under neutrality Ks ≈ 2μt and t = Ks ÷ 2μ.  With the pedigree rate μ = 5.11 × 10⁻⁹ per base per generation,17 t = 0.01938 ÷ (2 × 5.11 × 10⁻⁹) = 1.90 × 10⁶ generations.  Across Zhang's 95% interval, 4.08 to 6.14 × 10⁻⁹, t runs from 1.58 to 2.38 million generations.  At one generation a year that is 1.9 million years;  at the three years Zhang's team assumed for sea-run fish, 5.7 million.  Caveats:  a genome-wide per-base rate is applied to synonymous sites, which are not all neutral, and the Ks comes from the few genes in the shared region. ↩ ↩

  12. 12.

    For the same Ks to arise in G generations needs μ = Ks ÷ 2G.  With about 4,400 years since the Flood and one generation a year, G = 4,400 and μ = 0.01938 ÷ 8,800 = 2.2 × 10⁻⁶, which is 431 times the measured 5.11 × 10⁻⁹.17  At three years a generation, G ≈ 1,467 and μ = 6.6 × 10⁻⁶, 1,293 times.  Over 6,000 years at one generation a year, μ = 1.6 × 10⁻⁶, 316 times.  New mutations per fish at the measured rate are 2 × 463 × 10⁶ × 5.11 × 10⁻⁹ ≈ 4.7;21 at 431 times, about 2,040.  These figures price a mechanism that no creationist source found proposes.  They are given for scale only. ↩ ↩

  13. 13.

    Barrett, Rowan D. H. et al.  “Natural Selection on a Major Armor Gene in Threespine Stickleback.”  Science 322(5899):255–257.  Oct 10, 2008.  Read from the Schluter lab reprint.  Pages anchored by extracting each PDF page and reading its folio.  The PDF text layer renders the ± sign as "T";  it is printed here as ±.  Mirror:  Zoology.ubc.ca

    1. 13a.

      Page 256:  “We introduced approximately equal numbers of these fish (n = 45 to 46) to each pond in the spring of 2006, initiating replicate freshwater invasions.” ↩

    2. 13b.

      Page 256:  “Genotype frequencies at Eda in the F1 generation were not significantly different from the predicted 1:2:1 ratio (Fig. 2A)”. ↩

    3. 13c.

      Page 256:  “From October 2006 to May 2007, the frequency of the complete allele dropped from 67 to 49%, which reflected the comparatively poor survival of individuals homozygous for the complete allele. We calculated that the selection coefficient (S ) against the complete allele between these dates was 0.52 (± 0.10 SEM) (Fig. 2) (37).” ↩

    4. 13d.

      Page 256:  “In all four ponds, the frequency of the low allele was greater in the first sample of F2 offspring in June 2007 than in all F1 adults sampled in May [June F2: 57.0% (± 4.1% SEM), May F1: 51.6% (± 1.4% SEM)] (Fig. 2A) [one-tailed t test, t(3) = 2.14, P = 0.061].”  The founders were all heterozygotes, so the low allele started at 50%. ↩

    5. 13e.

      Page 257:  “The calculated selection coefficient (S) against the low allele between July and October 2006 was 0.50 (± 0.16 SEM) (Fig. 2C), which offset the gains occurring later in life.” ↩

    6. 13f.

      Page 257:  “This demonstrates not only that countervailing selection pressures diminish the advantage of the low allele over the whole life span but also that the overall fitness effects of Eda do not seem to be determined solely by differences in lateral plate number.” ↩

    7. 13g.

      No printed folios;  page numbers are PDF pages.  Cited here as url https://www.zoology.ubc.ca/~schluter/reprints/barrett%20et%20al%202008%20science%20-%20supporting%20online%20material.pdf.

    8. 13h.

      Page 2:  “We sampled approximately 35,000 fish using unbaited minnow traps. We brought all fish not possessing the full number of lateral plates (n = 354) back to the lab”.  “Of the 354 partial morphs collected from Oyster Lagoon, 182 were heterozygous at Eda. On June 1, 2006, we released the heterozygous fish into 4 artificial ponds”. ↩ ↩

  14. 14.

    Barrett, Rowan D. H. et al.  “Rapid Evolution of Cold Tolerance in Stickleback.”  Proceedings of the Royal Society B 278(1703):233–238.  2011.  Read from the PMC copy.  Published online 4 August 2010.  Mirror:  PMC

    1. 14a.

      Section Abstract:  “Cold tolerance evolved at a rate of 0.63 haldanes to a value 2.5°C lower than that of the ancestral population, matching values found in wild freshwater populations.” ↩

    2. 14b.

      Section 3(c) Evolved populations:  “This rate of phenotypic evolution is among the most rapid to be observed in a natural population [28].” ↩

  15. 15.

    Kitano, Jun et al.  “Reverse Evolution of Armor Plates in the Threespine Stickleback.”  Current Biology 18(10):769–774.  May 20, 2008.  Read from the publisher PDF.  Pages read from each page's folio.

    1. 15a.

      Page 769:  “In 2005, we found that all three lateral-plate morphs were present, with 49% completely plated morphs, 35% partially plated morphs, and 16% low-plated morphs (Figures 1A and 2C). Although a previous study had also shown that all three morphs were present in Lake Washington in 1968–1969, only 6% were classified as completely plated morphs”.  “explains 75.2% of the variance in plate number”. ↩

    2. 15b.

      Page 772:  “By using the empirically estimated values of m, we found that a selection coefficient s (strength of selection for the completely plated morph) of 0.58–0.72 (Table 1) can explain the evolutionary shift from 1969 to 1976 (from 6% completely plated morphs to 40.2% completely plated morphs)”.  A model fit using estimated migration rates, not a direct measurement. ↩

  16. 16.

    Barrett's team trapped about 35,000 fish and found 182 heterozygotes among the 354 with incomplete plates.13  That is at least 182 low alleles among about 70,000 chromosomes, 0.26%.  It is a floor, because fully plated heterozygotes were released untested. ↩

  17. 17.

    Zhang, Chaowei et al.  “Rate of De Novo Mutations in the Three-Spined Stickleback.”  Heredity 134(7):387–395.  2025.  Read from the PMC copy.  CC BY.  Mirror:  PMC

    1. 17a.

      Section Results:  “After applying a series of filters, 115 unique de novo mutations (DNMs) were detected”.  “The overall average per-generation point DNM rate was estimated to be µ = 5.11 × 10^−9/bp/gen (95% confidence interval: 4.08–6.14 × 10^−9/bp/gen, with zero-inflated method being 5.07 × 10^−9/bp/gen)”.  “Assuming a generation time of 3 years for the marine three-spined sticklebacks, the yearly mutation rate would be, on average, 1.70 × 10^−9/bp/year.” ↩ ↩ ↩ ↩ ↩ ↩ ↩

  18. 18.

    Behe, Michael J.  Darwin Devolves:  The New Science About DNA That Challenges Evolution.  HarperOne.  2019.  Read from the EPUB, whose print-page anchors give the page numbers (verified for the sibling article Measuring Biological Information).  The book does not mention sticklebacks.

    1. 18a.

      Page 180:  “A given mutation, then, can either make a new FCT (which I dubbed a gain-of-FCT mutation), destroy an old one (loss-of-FCT mutation), or do something else—either tweak an old FCT in a way that leaves it still working or affect some noncoded feature of a cell (which I called a modification-of-function mutation).” ↩

    2. 18b.

      Page 189:  “In my QRB classification scheme the mutation would be counted as modification-of-function—because no new functional coded element was gained or lost, just copied.”  Said of the citrate duplication in the Lenski experiment.  Applied here to Fads2 by analogy. ↩

  19. 19.

    Catchpoole, David.  “The Stickleback:  Evidence of Evolution?”  Creation Ministries International.  Aug 19, 2009.  Retrieved Sep 27, 2026.  Date from the page metadata.

    1. 19a.

      Section Natural selection at work:  “the so-called “rapid evolution” did not result in the addition of any new features but rather the reduction or even complete loss of existing features—bony plates and pelvic spines.” ↩

    2. 19b.

      Section Natural selection at work:  “For example, when researchers relocated 200 marine sticklebacks to freshwater ponds on the University of British Columbia’s Vancouver campus, they soon observed in the offspring that “natural selection favours reduced armour in freshwater”. From that and other studies, geneticists have now identified that a mutated genetic switch is responsible for the loss of armour. Specifically, it affects expression of the Pitx gene.”  CMI's footnote numbers are left out of the quotation.  Armour plates map to Eda;  Pitx1 controls the pelvis. ↩ ↩

    3. 19c.

      Section Natural selection at work:  “It is certainly no support for fish-to-philosopher evolution, as there is no new genetic information here, rather the mutated (i.e. broken) genetic switch simply prevents expression of a gene.” ↩ ↩

  20. 20.

    Lang, Gregory I., and Amber M. Rice.  “Evolution Unscathed by Behe's Darwin Devolves.”  Evolution 73(4):862–868.  2019.  Read from the first author's early-view PDF, paginated 1–7.  Issue page = early-view page + 861, by arithmetic.

    1. 20a.

      Page 864:  “Counter to Behe’s thesis, none of the observed Eda mutations were predicted to damage the function of the protein (Colosimo et al. 2005). The causative variants are likely cis-regulatory changes that decreased expression of Eda in developing armor, but not in other tissues (O’Brown et al. 2015).”  Early-view page 3. ↩

  21. 21.

    Jones, Felicity C. et al.  “The Genomic Basis of Adaptive Evolution in Threespine Sticklebacks.”  Nature 484:55–61.  Apr 5, 2012.  Read from the publisher PDF.  Pages anchored by extracting each PDF page and reading its running foot.  Mirror:  PMC

    1. 21a.

      Page 55:  “Both coding and regulatory changes occur in the set of loci underlying marine–freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.” ↩

    2. 21b.

      Page 55:  “a total gapped size of 463 Mb, close to previous estimates of 530 Mb (ref. 18).” ↩ ↩ ↩

    3. 21c.

      Page 57:  “SOM/HMM and CSS defined many other loci that also show globally shared marine–freshwater divergence, including 242 regions identified by either method (0.5% of the genome), and 147 regions identified by both (0.2% of the genome).” ↩

    4. 21d.

      Page 59:  “Many of these 64 regions (41%) mapped entirely to non-coding regions of the genome, and presumably contain regulatory changes (Fig. 6a). A smaller fraction contains protein-coding sequences with consistent non-synonymous substitutions between marine and freshwater fish (17%). Finally, a fraction of regions (43%) include both coding and non-coding sequences (including non-coding RNAs), but lack ecotype-specific amino acid substitutions”. ↩

  22. 22.

    Mitchell, Elizabeth.  “Sticklebacks Recycling Ancient Genetic Information.”  Answers in Genesis, News to Know.  Apr 14, 2012.  Retrieved Sep 27, 2026.  On Jones et al. 2012.

    1. 22a.

      “It appears, Kingsley’s team reports in its paper in Nature, that about 80% of the genetic differences arise from alterations in regulatory DNA and about 20% in actual coding DNA.” ↩

    2. 22b.

      “The “ancient genes” the fish seems to switch on and off are the genetic information God provided in order to enable this fish to adapt to various environments, with the ultimate variety in each environment regulated by natural selection and other factors.”  “This variation is not accomplished by acquiring new information but by reshuffling, regulating, and sometimes losing old information.” ↩ ↩

  23. 23.

    Nelson, Thomas C., and William A. Cresko.  “Ancient Genomic Variation Underlies Repeated Ecological Adaptation in Young Stickleback Populations.”  Evolution Letters 2(1):9–21.  2018.  Read from the PMC copy, which has no page numbers, so claims give section headings.  CC BY.  Mirror:  PMC

    1. 23a.

      Section Abstract:  “we find that this habitat‐associated divergent variation averages six million years old, nearly twice the genome‐wide average.” ↩

    2. 23b.

      Section Lineage sorting and time to the most recent common ancestor:  “To convert node ages estimated in BEAST into divergence times, in years, we assumed a 15 million‐year divergence time between threespine and ninespine stickleback at each RAD locus (Aldenhoven et al. 2010).” ↩

    3. 23c.

      Section Results and Discussion, on habitat-associated divergence and species age:  “Nevertheless, in some genomic regions marine and freshwater threespine stickleback are as divergent as threespine and ninespine stickleback, which are classified into separate genera.”  “Genome‐wide variation averaged 4.1 MY old, and TMRCA for the vast majority of RAD loci was under 5 MY old. In contrast, divergence times at habitat‐associated loci averaged 6.4 MYA”. ↩ ↩

    4. 23d.

      Section Results and Discussion, on habitat-associated divergence and species age:  “However, QTL analyses of ecologically relevant armor traits point either to independent genetic architectures in the two species (Shapiro et al. 2009) or to a locus associated with repeated mutational events among threespine stickleback populations (Chan et al. 2010; Shikano et al. 2013).” ↩

    5. 23e.

      Section Adaptive marine‐freshwater sequence divergence involves ancient allelic origins:  “Genome‐wide diversity was similar across populations and habitat types (mean π RS = 0.0032, π BL = 0.0034, π BP = 0.0026, π FW = 0.0038)”. ↩ ↩

  24. 24.

    Nelson and Cresko give 6.4 million years for habitat-linked loci against 4.1 million genome-wide, a ratio of 1.6.23  Roberts Kingman's team gives 1,600 against 700 thousand years, a ratio of 2.3.4  The absolute ages differ about fourfold to sixfold because the papers use different calibrations and different measures:  the age of all alleles at a locus in one, the divergence between two fish in the other.  The ratio does not depend on the calibration. ↩

  25. 25.

    Inputs:  about 3,000 founders in each lake, Cheney in 2009 and Scout in 2011, followed to 2017;4 a genome of 463 Mb;21 μ = 5.11 × 10⁻⁹ per base per generation.17  New mutations per diploid fish = 2 × 463 × 10⁶ × 5.11 × 10⁻⁹ = 4.73.  Fish-generations, if the lakes held only their founding numbers and bred every year:  3,000 × 8 + 3,000 × 6 = 42,000.  Mutations tried:  42,000 × 4.73 ≈ 2.0 × 10⁵.  Gametes:  84,000.  A given base mutates at about μ per gamete, and to one particular other base at no more than about a third of that, so one specific change is expected at most 84,000 × 5.11 × 10⁻⁹ ÷ 3 ≈ 1.4 × 10⁻⁴ times, about 1 in 7,000.  At Xie's fragile-site rate of about 10⁻⁵,8 treated here as a rate per gamete per generation, 84,000 × 10⁻⁵ ≈ 0.84.  Real lake populations grew past their founders, which raises every count;  if the fish bred every second year instead, the generation count halves. ↩

  26. 26.

    The standard neutral formula π = 4Neμ gives Ne = π ÷ 4μ.  With π = 0.0032 in sea-run Rabbit Slough fish23 and μ = 5.11 × 10⁻⁹,17 Ne ≈ 0.0032 ÷ (2.04 × 10⁻⁸) ≈ 157,000.  The formula assumes a population at equilibrium, and the result is a long-term average, not a census. ↩

  27. 27.

    Answers in Genesis.  “Fast Fish Evolutionary Response.”  Answers in Genesis.  Aug 7, 2010.  Retrieved Sep 27, 2026.  Unsigned.  On Barrett et al. 2011, cold tolerance.

    1. 27a.

      “This three-year example of stickleback “evolution” is a soft echo of the adaptation and speciation that would have occurred in the centuries after the Ark landed—and in the millennia since.” ↩

  28. 28.

    Gaskill, Phil B., and Randy J. Guliuzza.  “Engineered Adaptability:  Blockchain-Like Process May Produce Adaptive Traits.”  Institute for Creation Research.  Apr 30, 2019.  Retrieved Sep 27, 2026.  Part of ICR's Continuous Environmental Tracking series.

    1. 28a.

      Section Blockchain Technology and the CET Model:  “When environmental conditions change, sensors trigger an organism’s pre-programmed logic to select pre-existing information designed to respond to the changed condition.” ↩

    2. 28b.

      Section Threespine Sticklebacks:  “This fits expectations if adaptive traits in threespine stickleback are produced “on the fly” by changes to the “instruction code” as hypothesized in the computer program model of DNA.” ↩

  29. 29.

    Answers in Genesis.  “Natural Selection Theory Strengthened By Armored Fish Study.”  Answers in Genesis, News to Know.  Sep 6, 2008.  Retrieved Sep 27, 2026.  Unsigned.  Read from the page bytes, saved in the source library.

    1. 29a.

      “The researchers based their work on the discovery of a mutant gene found in sticklebacks that prohibits the growth of the armor. While the gene is found in less than one percent of marine sticklebacks, it is quite common in the freshwater variety.”  “The team relocated 200 marine sticklebacks that carry the armor-reduction allele to “experimental” freshwater ponds.”  Barrett's supplement gives 182 founders, 45 or 46 per pond.  Colosimo measured 0.2% and 3.8% in two marine samples. ↩

  30. 30.

    Thomas, Brian.  “Amazing Fish Adaptive Design.”  Institute for Creation Research.  May 18, 2012.  Retrieved Sep 27, 2026.  On Jones et al. 2012.

    1. 30a.

      “Both the genes and regulatory sequences were very similar between all the fishes studied, showing that they have been resistant to mutational changes.” ↩

  31. 31.

    Thomas, Brian.  “Tiny Fish Rapidly Adapted to Cold Water.  Was It Evolution?”  Institute for Creation Research.  Aug 17, 2010.  Retrieved Sep 27, 2026.

    1. 31a.

      “By the third generation, descendants remained of those "rare individuals" that already had the capacity to tolerate cold before the experiment began.”  “But in considering this conclusion, it becomes apparent that the reference to "selection" is superfluous.” ↩

  32. 32.

    Zhang, Chaowei et al.  “De Novo Mutation Rates in Sticklebacks.”  Molecular Biology and Evolution 40(9).  2023.  Read from the PMC copy.  CC BY.  Measured in the ninespine stickleback.  Article number msad192.  Mirror:  PMC

    1. 32a.

      Section Abstract:  “Since the estimated µ was 3.1 times smaller than the commonly used substitution rate, recalibration with µ led to substantial increase in estimated divergence times between different stickleback species.” ↩

    2. 32b.

      Section Divergence Time Estimation With DNM Rates:  “the date comparison of scenarios using Guo et al.'s (2013) SSR vs. our DNM rate showed the divergence time of node A shifted from 5.6 Mya (95% highest posterior density [HPD]: 3.3–8.1 Mya) to 17.5 Mya (95% HPD: 10.6–25.3 Mya)”.  Node A is the split between the ninespine and threespine stickleback lineages. ↩

  33. 33.

    Answers in Genesis.  “Reverse Evolution Found In Small Fish Species.”  Answers in Genesis, News to Know.  May 24, 2008.  Retrieved Sep 27, 2026.  Unsigned.  On Kitano et al. 2008.

    1. 33a.

      “Evolution, in its most straightforward biological definition, simply means the change in a biological population over time (ultimately, a change in gene frequencies). Thus, the stickleback population’s loss of armor was evolution, and now the reappearance of the armor is evolution.” ↩