Lesson 14

A

Locked — confirm your name above to begin.

For this whole lesson, we are going to look at some actual empirical data in the same general shape as the simulations you've been running. Our first stop involves the ground finches from Daphne Major, monitored from 1973 to 2012 by the Grants. In the 2000s, genetic sequencing became cheap enough that a lot of the specimens the Grants had been collecting, and especially the new birds they were still monitoring, could be sequenced. The genomes could be analyzed, and the specific alleles being passed on from parent to child could be correlated with the phenotypes that the Grants had been measuring the whole time.

A number of loci were found associated with each of these traits, as you might expect: the characters they were looking at are controlled by a great many loci. However, a few particular loci stuck out. HMGA2 had a strong effect on beak size, depending on which allele you had; ALX1 had a strong effect on beak shape; BMP4 had a strong effect on beak depth; and CALM1 had a strong effect on beak length. Multiple alleles for these and other loci (as, again, each trait was controlled by a large number) were then monitored over time. We're going to track one allele of HMGA2, controlling beak size.

In the simulations below, it's fundamentally the same as what we did before. Rainfall controls how much food is available; how much food is available controls births and deaths. If you don't die, you're a survivor. If you are born, you're a new chick. And the frequency of the big-beak allele is going to be a function of who lives (who survives, who doesn't die), and who is born, and who their parents are.

Where this gets into selection is that beak size is going to have an impact on how rainfall creates food. When the rain falls, seeds are made. All of the birds can eat small seeds, but only certain birds can eat the larger, harder seeds. As a result, how much food is made available to a bird by the falling rain is itself a function of the bird's beak. And so here you're going to run a few simulations to see how the frequency of an allele changes when you can mess with the association of the trait with food acquisition, relative to the amount of births and deaths that are random.

Essentially, in the graph of the model shown below, you're going to be able to influence drift directly, through the role of chance in who lives, who dies and who reproduces; and the role of selection, in terms of how specifically the allele affects food availability as a function of rain, and in what direction. You'll fit this a few times to simulated data, to get a feel for how your manipulations impact it. And then you'll get one chance to try and match the empirical data that was actually seen in the Grants' birds.

  • Simulation based on the medium ground finches from Daphne Major.
  • On the right, a schematic showing loci that vary across several chromosomes in the finches, with particularly strong-effect loci highlighted and labeled. HMGA2, the beak size locus, is the one where we're going to track the two alleles through time.
  • Below, you can control the trait beak size, influenced by HMGA2, and the amount of stochasticity there is in the births and deaths.
  • Try to match the change over time in the frequency of your big-beak allele as a function of rainfall, in terms of both the direction and the certainty.

Predict

Controls

R code


      

To open Stage B

  • Five targets: every one you hit counts. 0 of 5 attempts
  • Try to match actual data. not yet taken
Stage B is open.

B

Solve Stage A to unlock this section.

In Part A, you were predicting the frequency of an allele at a single locus. A diagram was given showing its context in the genome, amongst many other loci, for that chromosome and several others. Meiosis does not pass on individual alleles; rather, in sexually reproducing diploids, you pass on an entire haploid cell. Recombination scrambles some loci, but genomes are still inherited as chunks instead of individual pieces.

Here we're going to run some simulations where we're tracking the frequency of a new allele, but we're going to be concerned not only with its frequency over time, but also its genomic context. The model here will be extremely simple. We are just looking at the advantage of the allele, and the effects of chance (drift, basically), on the rate at which it spreads. The focus here is on how the rate at which it spreads affects the diversity in other places in the genome: places that have no effect whatsoever on the phenotype of interest.

This is based on the FGFR3 receptor, a protein that is related to fibroblast growth factor. If you look at this particular locus in dachshunds versus non-dachshunds, you'll find that dachshunds have almost zero genetic diversity in a region about 20 megabases, or about 20 million base pairs, on either side of FGFR3. All dachshunds are homozygous for all the same alleles. It's a very, very big region of zero diversity compared to normal dogs in the same area, even though leg length has been important to all dog breeds.

The size of the homozygous region is going to be related to how fast the allele went to fixation: how fast the allele drove all other alleles extinct at that locus. The faster this happens, the less time there is for recombination to separate one locus out from its neighbors, producing a situation where a strongly advantageous allele sweeps the area around it free of diversity, just based on the speed of its wiping out variation around it, relative to the speed of recombination.

  • Simulation of one chromosome, 4 million bases long, with 100 marker sites; a new allele arises at a gene in the middle.
  • Top left: 40 chromosomes from the observed data, and 40 from your run. One row is one chromosome; those carrying the new allele are on top; dark is the rarer allele at each marker.
  • Top right: the new allele's frequency, generation by generation. Red: new alleles lost by chance before one spread.
  • Below, you can control the advantage of the new allele, and chance: the thicker the arrow, the fewer individuals breed.
  • Try to match two readings: diversity in the rest of the genome, and the stretch two carrier chromosomes share. Go runs three populations; their average counts.
  • Last, try to match actual data: one shot, no practice. What it is shows after.

Predict

Controls

R code


      

To open Stage C

  • Five targets: every one you hit counts. 0 of 5 attempts
  • Try to match actual data. not yet taken
The selective sweep tells you how strong selection was: a wider area of low variation means selection had to outpace that much recombination. A locus is just a chunk of DNA selection hasn't recombined apart yet. Stage C is open.

C

Solve Stage B to unlock this section.

Do you inherit traits? You don't inherit hair color. You inherit nucleotide sequences that code for proteins that, in the context of the environment and the other proteins you have, produce a trait. Does selection act on the nucleotide sequence, or does it act on the trait? It's the trait. There's a massive disconnect between what gets inherited and what gets selected.

Selection is fast — way faster than mutation, faster than drift in almost every circumstance, sometimes faster than recombination.

  • Stage B's chromosome, now with a trait: the new allele and the environment set it, and survival follows the trait.
  • Each target is observed data based on a real study: the stretch two chromosomes carrying the new allele share, measured when it reached 95%. The selection on the trait is held where the first target puts it.
  • The diagram is the controls. One arrow glows each target; the rest are pale and held. One copy → trait: how much of the allele's effect one copy gives (two copies give all of it). Environment → trait: how much else makes the trait differ. Wet years push the other way. A harmful allele beside it, too close for crossovers to split.
  • Predict: each population is a dot, the generations it took relative to the stretch its carriers share. Go runs six; get their average inside the red band. Dots stay from target to target.
  • Play some runs first; Start the targets when you are ready. Tick practice for a run that does not count.

Predict

Controls

R code


      

To open Stage D

  • Five targets: every one you hit counts. 0 of 5 attempts
No matter how strong selection is on an additive trait, the population only moves by the additive variation in the population each generation. Selection sets direction. Additive variation sets the magnitude. Stage D is open.

D

Solve Stage C to unlock this section.

MHC-I is a billboard. Every cell sticks samples of what it's making on its surface, so killer T cells can check. Normal advertised? Fine. Wrong protein advertised? Killed.

Nothing is ever fixed in time. If scale-eating is rare, do the other fish look out for scale-eaters? No. As scale-eaters become more common, what do the other fish start to do? Watch out for them. Is scale-eating still as good? No. How good scale-eating is depends on how common it is. The fitness function itself changes — partly because environments change, but also because the fitness function changes because of the fitness function.

  • One immune gene with many alleles, like the MHC. It starts with 20 alleles, equally common, and runs 500 generations; one copy in 1,000 turns into a brand-new allele each generation.
  • Pathogens adapt to the host alleles that are common, so a host carrying common alleles is the one they infect. Pathogens → survival: how hard they hit it. Chance: the thicker the arrow, the fewer individuals.
  • A second gene, elsewhere in the same genome, that nothing reads: it starts and mutates the same way.
  • Top left: how common each allele was at each gene, generation by generation. Top right: the alleles a sample of 100 individuals holds at each gene, commonest first, the observed data's and your run's.
  • Predict: alleles at the neutral gene across, alleles at the immune gene up. Go runs three populations; the average of the three must fall inside the red box.
  • Play some runs first; Start the targets when you are ready. Tick practice for a run that does not count.

Predict

Controls

R code


      

To open Stage E

  • Five targets: every one you hit counts. 0 of 5 attempts
The system can sustain a frequency of liars and a frequency of risk-tolerant honest crabs that fluctuates over time. Whichever strategy is rare wins. Stage E is open.

E

Solve Stage D to unlock this section.

If I told you a car was driving at 50 mph and has gone 100 miles, how long ago did it leave? You have a distance and a speed — distance divided by speed gives you time. Number of mutations is your distance. The slope of mutations versus time is your speed. Number of mutations divided by slope gives you time to most recent common ancestor.

  • One gene, 1,200 sites: 300 silent (a change there leaves the protein as it was) and 900 protein-changing. Each site changes in one copy in 100,000 each generation, for 10,000 generations. Go runs ten genes with the same settings.
  • A silent change does nothing. A protein change's effect is drawn from the curve, as in Lesson 12 D: most harm, and the bar at the left is the share that help.
  • Controls: individuals. Drag the curve to move the typical harm (each step to the right is ten times more harmful); drag the bracket for its spread; drag the bar for the share that help. What a target holds is grey.
  • Top left: the first gene's history. Every change that took over, drawn as it spread: silent in grey, protein changes that helped in blue, harmful ones that slipped through in red. Top right: how many of each kind arose, and how many took over.
  • Predict: across, how many silent sites in 1,000 differ between two copies of the gene; up, the protein changes that took over for every silent one, per site. Ten genes make one dot; get it inside the red box. The open diamond is the same ratio among changes still varying.
  • Play some runs first; Start the targets when you are ready. Tick practice for a run that does not count.

Predict

Controls

101
what a new protein change does

R code


      

To finish

  • Five targets: every one you hit counts. 0 of 5 attempts
The more important a gene is, the less variation you'll see at it — not because mutations don't happen there, but because mutations that do happen there never persist.

Good job!

That is all of Lesson 14. Your completion code is below — copy it and hand it in.