General Genetics/Population Genetics
Population Genetics is the field of genetics which studies allele distributions and genetic variation in populations. Population geneticists study the processes of mutation, migration, natural selection and genetic drift on populations, and in doing so are studying evolution as it occurs.
Foundations of Population Genetics
editTempleton states that the three premises of population genetics are the same premises for population genetics:
- Templates/DNA can replicate
- Templates/DNA can mutate and recombine
- Phenotypes emerge from the interaction of templates/DNA and environment.
Replication of Populations
editThere are three major properties that a population must maintain with replication:
- They are composed of reproducing individuals
- They are distributed over space and time
- They host a population of genes
The first property indicates that individuals of the population must reproduce to keep the population stable. This is necessary because individuals breakdown over time do to the introduction of entropy and inability for an individual to continuously remove the entropy added by the environment. Thus to maintain the population, individuals must pass down their DNA, or organizational encoding, to the next generation. Through continuous reproduction, a population can be maintained over a much longer time than the individuals that comprise it. In addition, the continuous reproduction of over time enables for the population to have properties and components of its own.
The second property is that a population is distributed over a space. Populations can exists as:
- small isolated groups
- a collection of groups with a varied amount of genetic exchange
- a large interbreeding population that exists over a vast space
In general though the population can be divided into a primary group that can be considered as interbreeding and a secondary group that mates occasionally with the primary group. It is this primary group that population geneticists generally study, as it is generally stable. They define the group as a group of interbreeding individuals that share a common system of mating. The secondary group is generally ignored, and treated as noise in the system, unless it is having a major effect on the primary group.
The third property that must be maintained with reproduction is the population's gene pool. The gene pool is the collection of all the genes, organizational templates, in the population that can be used to create new individuals. By studying this gene pool, geneticists can determine the frequency of alleles, and or groups of alleles in the population and how they are changing over time. From the patterns of result that are obtained, geneticists then can start to understand what forces are acting on the population.
Template Mutations
editChange is a requirement of evolution and one method of introducing change is through modification of the templates used sustain the population. In the case of living life, these templates are genes.
Sources of Mutation:
- Insertions
- Deletions
- Single Nucleotide Substitutions (sometimes changing the protein sequences and sometimes not)
- Transpositions
- Duplications
An allele is an alternate form of an template. In the case of biological systems, an allele is a form of a gene. Zooming further out, a version of a region of templates is called a haplotype. Biological systems would call this a sequence of nucleotides, while a in a computer system, this would be a sequence of linked objects.
Modeling Evolution
editInitially we are going to consider populations with genetic architectures of two loci per template.
Hardy-Weinberg Model
editThe genetic architecture of the Hardy Weinberg model is one locus, two allele model (Templeton, p. 35).
Hardy-Weinberg Equilibrium
editIf a population has no forces of evolution acting upon it is in Hardy Weinberg Equilibrium. Quantitatively it says that if the allele proportions of two alleles A and a are denoted p and q then the genotype proportions will be such that the homozygote AA will be of proportion p^{2}, the heterozygote Aa will have proportion 2pq and the homozygote aa will be of proportion q^{2}.
Testing for Hardy-Weinberg
editTo do this test you do what is known as a Chi Square test ( ). Where:
and the degrees of freedom is no. of Genotypes - no. of Alleles if χ ^{2} < table value then we accept that our population is in Hardy-Weinberg.
...note that we don't actually take a square, it is there out of tradition.
Concise Chisquare Table
editDegrees of Freedom | 1 | 2 | 3 |
χ ^{2} Value | 3.84 | 5.99 | 7.83 |
Example
editIf there is a population with genotype proportions AA: 0.1 Aa: 0.4 aa: 0.5
So the allele proportions are
So our expected values are AA: 0.09 Aa: 0.42 aa: 0.49.
Therefore we can accept the hypothesis is in Hardy-Weinberg Equilibrium and that there are no forces of equilibrium on these alleles.
Two Autosomal Loci, Two Allele Model
editExtending the H.W. model to two autosomal locus model with two alleles.
For the purpose of this discussion, the first locus will have alleles A and a and the second locus will have B and b. From this we can get the following gamete types and their frequencies through recombination:
Gamate | Frequency |
---|---|
AB | Freq_{AB} |
Ab | Freq_{Ab} |
aB | Freq_{aB} |
ab | Freq_{ab} |
Sum | 1 |
A population producing the above four gamates can produce the following genotypes:
Gamates | AB | Ab | aB | ab |
---|---|---|---|---|
AB | Freq_{AB}^{ . }Freq_{AB} | Freq_{AB}^{ . }Freq_{Ab} | Freq_{AB}^{ . }Freq_{aB} | Freq_{AB}^{ . }Freq_{ab} |
Ab | Freq_{Ab}^{ . }Freq_{AB} | Freq_{Ab}^{ . }Freq_{Ab} | Freq_{Ab}^{ . }Freq_{aB} | Freq_{Ab}^{ . }Freq_{ab} |
aB | Freq_{aB}^{ . }Freq_{AB} | Freq_{aB}^{ . }Freq_{Ab} | Freq_{aB}^{ . }Freq_{aB} | Freq_{aB}^{ . }Freq_{ab} |
ab | Freq_{ab}^{ . }Freq_{AB} | Freq_{ab}^{ . }Freq_{Ab} | Freq_{ab}^{ . }Freq_{aB} | Freq_{ab}^{ . }Freq_{ab} |
If you analyze the table above, it can be noticed that there are only ten unique combinations. Four combinations correspond to homozygous zygotes and the remaining six are the heterozygous zygotes.
Notice that the sum of Freq_{AB}, Freq_{Ab}, Freq_{aB}, and Freq_{ab} is one. This follows from the earlier model where the sum of p and q equaled one.
Recombination
edit(this section is used to discuss how we get term r)
Homozygous Case
editAB/AB or ab/ab or aB/aB or Ab/Ab
Single Heterozygous Case
editAB/Ab or AB/aB aB/ab or Ab/ab
Double Heterozygous Case
editAB/ab or Ab/aB
Note that recombination is only noticeable in double heterozygotes.
Linkage Disequilibrium
edit(this section is used to discuss how we get term D)
Linkage disequilibrium measure,
editFormally, if we define pairwise LD, we consider indicator variables on alleles at two loci, say . We define the LD parameter as:
Here denote the marginal allele frequencies at the two loci and denotes the haplotype frequency in the joint distribution of both alleles. Various derivatives of this parameter have been developed. In the genetic literature the wording "two alleles are in LD" usually means to imply . Contrariwise, linkage equilibrium, denotes the case .
Pulling Information from the Model
editIs Evolution Occurring
editIf r > 0 and D != 0 then evolution is occurring
If r = 0 or if D equals 0, then no evolution is occurring
Mutation and Disequilibrium - Normalized Linkage Disequilibrium
editMating Systems
edit
Selection
edit
Absolute Fitness
Relative Fitness
Frequency before selection
Frequency after selection
Genetic Drift
editFirst experiments were done in the early 30's. German was a major scientific language before the second world war.
Genetic drift is a function of population size N.
The effects of genetic drift is inversely proportional to population size. This means as the population increases, the deviation from expected allelic frequencies will decrease. (Templeton, p. 84).
Genetic drift is non-directional.
There is no attraction to return ancestral allele frequencies.
Genetic drift is a cumulative function. Changes in allele frequencies from the previous generations are added to the changes that occur in the current generation.
Genetic drift is only occurs when there is variability.
To study genetic drift this section will create a simplified model and expand upon it. The model will just have two states, a) genetic drift is occurring and b) genetic drift is not occurring.
a->b fixation b->a mutations causing differentiable patterns b->a reintroduction of differentiable patterns
Genetic drift can be broken down into two simplified states to start. The f
Example
editIn each population, 2N=32, po=qo=0.5 initially and then did 19 generations of data. (Buri, P. 1956 "Gene frequency in small populations of mutant Drosophila. Evolution 10: 367-402
(See Figure 6.3 from Hedrick, P.W. 2005, Genetics of populations, 3rd edition. Jones and Bartlett, Sudbury, MA)
Note that in this figure, variance is increasing, but mean allele frequency over populations is staying relatively the same.
Simulation of Genetic Drift
editGenetic Drift can be simulated using a Monte Carlo Simulation. (See Figure 6.2 in Hedrick)
The proportion of populations expected to go to fixation for a given allele is equal to the initial frequency of that allele
Only the allele frequencies are changing and the distribution of the allele frequency. The mean allele frequency over multiple replicate populations does not change due to genetic drift.
Heterzygositity or the variance of the allele frequency over the replicate population can be used to understand genetic drift.
We can calculate the number of generations necessary to reduce the reduce the heterozygostity
t=ln(x)*-2*N, where x is how much heterozgostity is left, and N is the population size.
Coalescence Theory
edit
Furthur Reading in Wikipedia
editGlossary
editGene Pool - p. 37 and p. 38 of Templeton's book
Formulas
editHardy Weinberg
editExample 1
editHomozygotes | Heterozygotes | Sum | |
---|---|---|---|
Observed | 85 | 15 | 100 |
Expected | 60 | 40 | 100 |
The probability of observing these numbers is:
Example 2
editGenotype | Count |
---|---|
Homo_{dom} | 20 |
Hetro | 50 |
Homo_{rec} | 30 |
What are the observed frequencies of p and q?
What are the expected number of genotypes?
Do the observed genotypes fit H-W expectation?
If χ^{2} greater than degree of freedom, then the null hypothesis, in this case it fits H.W. can be rejected.
Example 3 - Next Generation
editFor an autosomal locus
Example 4 - Sex Linked (1 locus, 2 allele)
edit
It is possible to solve for the amount of time it will take till an initial frequency is less than or equal to a particular threshold by using the following formula:
Example 5 - Inbreeding - Self Pollination
editGenotype | Count |
---|---|
Homo_{dom} | 20 |
Hetro | 50 |
Homo_{rec} | 30 |
Calculate F Value
Chi-squared then can be calculated by:
Example 6 - Inbreeding (1 Locus / 2 alleles)
editGiven the frequency of the dominant allele and an F value:
Allele | Frequency |
---|---|
p | 0.6 |
q | 0.4 |
What sample size would be necessary to detect an effect of F=0.01 at the S% significant level?
- Look up the X^{2} value for S% significant level (usually S = 5%)
- Use the following formula:
Example 7 - Null Allele
editGenotype | Count |
---|---|
Homo_{dom} | 20 |
Hetro | 50 |
Homo_{rec} | 30 |
Selection Part 1
editExample 1
editGenotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Viability | 0.7 | 1 | 0.6 |
Fertility | 11 | 6 | 3 |
Genotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Viability (v) | 0.7 | 1 | 0.6 |
Fertility (f) | 11 | 6 | 3 |
Absolute Fitness (W) | |||
Relative Fitness (w) |
What are the selection coefficients?
Example 2 - Mean Fitness Calculation
editGenotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Freq | 0.3 | 0.6 | 0.1 |
Rel. Fitness | 1 | 0.6 | 0.1 |
What is the mean fitness of this population?
Example 3 - Natural Selection Analysis
editGenotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Freq | 0.3 | 0.6 | 0.1 |
Rel. Fitness | 1 | 0.6 | 0.1 |
What is the mean fitness of this population?
Selection Part 2
editExample 1
editHomo_{dom} | Heterozygotes | Homo_{rec} | |
---|---|---|---|
Relative Fitness_{females} | 0.6 | 0.8 | 1 |
Relative Fitness_{males} | 1 | 0.8 | 0.6 |
What are the equilibrium frequencies for the pf and qf alleles in females considering that heterozygotes have intermediate fitness?
What are the equilibrium frequencies for the p_{m} and q_{m} alleles in males considering that heterozygotes have intermediate fitness?
What range of s_{f} values will give a stable polymorphism?
What range of s_{m} values will give a stable polymorphism?
Example 2 - Complete Dominance w/ Frequency Dependent Selection
editHomo_{dom} | Heterozygotes | Homo_{rec} | |
---|---|---|---|
Relative Fitness (w) | 1.3 | 1.3 | 1.7 |
What are the selection coefficients s_{homo, dom} and s_{hetero} with q = 0.02?
Calculate the stable equilibrium frequency, q_{e} and p_{e}.
Calculate the mean fitness assuming HW equilibrium and using:
(a) current allele frequencies
(b) stable equilibrium allele frequencies
(c) general formula
Example 3 - Self Incompatibility Locus
editGiven initial frequencies:
(a) What are the frequencies of the S1, S2, and S3 alleles initially, and after one and two generations?
(b) After one generation?
(b) After two generations?
Example 4 - Fitness / Overdominance
editHomo_{dom} | Heterozygotes | Homo_{rec} | |
---|---|---|---|
Relative Fitness | 0.4 | 1.0 | 0.8 |
To find the equilibrium of p:
Example 5 - Selection and Viability
editGenotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Zygotic Frequencies | 0.05 | 0.33 | 0.62 |
Viabilities | 1 | 0.6 | 0.59 |
Genotypes | Homo_{dom} | Heterozygotes | Homo_{rec} |
---|---|---|---|
Zygotic Frequencies (Z) | 0.05 | 0.33 | 0.62 |
Viabilities (V) | 1 | 0.6 | 0.59 |
Adult Freq After Selection |
or
Predicted Genotype Frequencies
Example 6
editHomo_{dom} | Heterozygotes | Homo_{rec} | |
---|---|---|---|
Observed Numbers | 750 | 1000 | 250 |
What are the absolute and relative fitness for each?
Homo_{dom} | Heterozygotes | Homo_{rec} | |
---|---|---|---|
Observed Numbers | 700 | 1000 | 200 |
Predicted Numbers | 500 | 1000 | 500 |