GS 561: Molecular Population Genetics and Evolution. T R 9:00- 10:20. Foege S110.

Instructor: James H. Thomas (Jim), Office: Foege S340B, Email: jht@uw.edu

http://faculty.washington.edu/jht/GS561_2017/

The aim of this course is to cover topics in molecular evolution and population genetics. The format will be a discussion of a general topic focused around a specific paper on Tuesdays, followed on Thursday by a computer lab on the same topic.

For the first week, we will have a discussion led by the instructor of Maximum-likelihood and Bayesian tree methods.

In subsequent weeks, on Tuesday everyone will turn in a one-paragraph summary of the primary assigned paper for the week. Class will be a discussion of the paper.

On Thursday we will use software related to each topic to analyze sample data.

I will not hold formal office hours, but you are welcome to drop by my office any time (I am usually there). Foege S340B.




Week 1    Course summary, Maximum-likelihood, MCMC and Bayesian phylogeny (Thomas).

        Papers:

                Lewis (2001). Phylogenetic systematics turns over a new leaf. Trends in Ecol. and Evol. 16: 30-37.
                Huelsenbeck et al. (2001). Bayesian inference of phylogeny and its impact on evolutionary biology. Science 294: 2310-2314.

        Maximum-likelihood tree slides

         Bayesian tree slides



Week 2    Genome-scale evolutionary conservation

          Paper guide: Let's focus this week on the biological results a little more than the method. For the methodology, understand Figure 1 (in detail, what all the symbols mean etc., you may have to do some background reading on Hidden Markov Models). For the biology, focus (as the paper does) on the relative amounts of conserved DNA in different eukaryotic groups, and also think about how this probably extends to bacteria, archaea, and viruses (not covered in the paper). A second focus should be the fraction of conserved non-coding DNA in different groups and what it means. A hard question to consider: how should we define and measure organismal complexity? This is a term used often in biology, rarely with real meaning (often it means roughly "the closer it is to human the more complex it is"). Things NOT to focus on: 1) the data for nematodes (as the authors several times point out, the two taxa available at that time are too divergent to give reasonable results; unfortunately they include the results anyway), 2) the section on secondary structure in non-coding HCEs.

        Primary Paper:

                Siepel et al. (2007). Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Research 15: 1034-1050.

        Additional Papers:

                Bejerano et al. (2004). Ultraconserved elements elements in the human genome. Science 304: 1321-1325.

                   Lowe et al. (2011). Three periods of regulatory innovation during vertebrate evolution. Science 333: 1019-1024.



Week 2 Lab - data sets and scripts can be gotten from /net/thomas/vol1/home/jht/GS561 (this is not publicly available - you will have to copy the files from a GS cluster login)

You won't be able to run phastCons in situ because you will need write permission (unless you edit the scripts to make output to your directories, which you can't do directly because the scripts are in my directory - Catch-22).

The simplest approach is just to copy everything (nothing is too huge). From your own directory (or somewhere you have write permission):
cp /net/thomas/vol1/home/jht/GS561/* .


The file NOTES.txt has a description of what the various files are and what you can do with them. You will have to run phastCons on the GS cluster, but the post-processing Python programs should run any where Python 2.x is installed.

The basic lab will be to run phastCons using different rho, mu, and nu parameters and analyze the output for its general properties.
rho is the multiplier that phastCons uses to get the C (conserved) tree starting from the U (unconserved) tree (which I have provided, made from 4D sites).
mu is the transition probability for C to U
nu is the transition probably for U to C


                   


Week 3    Positive Selection

         Paper guide: The primary paper uses the methods for dN/dS analysis developed in a series of papers by Yang and Nielsen (and often Goldman and a few other authors). In the Swanson paper, read the statistical analyses Methods section carefully and relate it to the Yang et al. paper (see below). In the main text and Table 1 you can discount the M3 vs. M0 tests (these are not very realistic models) and concentrate on the M7 vs. M8 tests, which are the most commonly used models. It is rather complicated, but try to understand all the numbers in Table 1, especially the Likelihood Ratio Tests (2 deltaL) and Parameter estimates. As for the biology, pay close attention to the Introduction and to the interpretations given in Discussion (you may have to read other papers or sources to understand fully why these processes might give rise to positive selection). Also think about other processes that might drive positive selection, and strengths and weakness of the dN/dS approach to detecting them.

         The Yang et al. paper is long and very complex. You can ignore all of the tests on real data (the Swanson paper is more interesting) and most of the various model combinations the authors discuss. Concentrate on the Theory section (up to Statistical distributions, which you can skim), especially the M7 and M8 models and the Discussion. We won't cover how specific sites are identified (it is complex and the essence is simple - it looks for sites where the maximum-likelihood dN/dS value is high).

        Primary Paper:

               Swanson et al. (2001). Positive Darwinian selection drives the evolution of several female reproductive proteins in mammals. PNAS 98: 2509-2514.

        Additional Paper:

               Yang et al. (2000). Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155: 431-449.


Week 3 Lab - All of this week's lab can be done on your laptop or on the GS compute cluster (but you will have to work out the details on your own). BEFORE CLASS, to be sure everything is working, download and unpack PAML software (Win paml4.9d.tgz or Mac/Linux paml4.8a.macosx.tgz). Do NOT use PAML-X. I don't think the PAML version will matter but let me know if you have problems. Download input files: codeml_test.ctl, codons_test.paln, and tree_test.dnd. Put all three input files and the PAML program codeml (or codeml.exe on Win, found in the bin directory after PAML is unpacked) together in a new directory (or add the codeml directory to your path). From a command prompt, navigate to your directory and type "codeml codeml_test.ctl" to make sure you don't get any errors (it should take a few minutes to run and will write output to codeml_test_out.txt). If you do have trouble, come see me. Don't worry about what all the stuff that writes to the screen - just make sure codeml runs without errors. BTW, you can get rid of the screen progress monitor but adding a line "noisy = 0" to the .ctl file.

Week 3 Lab




Week 4  Deleterious and neutral mutations

          Paper guide: There is great interest in predicting which alleles segregating in human populations are deleterious (i.e. cause "disease", or more broadly cause phenotypes). Most such alleles are thought to be neutral or nearly neutral (have little or no phenotypic consequence). The Kircher paper takes a new and different approach to this issue by using fixed nucleotide changes on the human lineage as a sample of nearly neutral changes and comparing them to all possible changes (simulated), which will be a mixture of neutral and deleterious alleles. Focus your thinking on the general issues - how this approach can inform analysis of any genome population - and don't focus on human disease per se. The paper spends most of its time on additional efforts to combine approaches and datasets and using machine-learning to evaluate the combined data to predict human disease alleles. Try to understand this part of the paper at a high level, but don't focus too much time there. Class discussion will be focused more on general population genetic issues - where mutations come from, how they become fixed in populations, etc.

          Unless you have a strong background in population genetics and molecular evolution, I strongly recommend reading the Duret overview of some of the key issues.

        Primary Paper:

                Kircher et al. (2014). A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46: 310-315.

                Duret L. 2008, a short overview of neutral theory.


Week 4 No Lab



Week 5  Population Structure

        Paper guide: Because of a long history of migration of small founder populations, followed by relative genetic isolation of geographically restricted groups, there are substantial genetic differences among various human populations (though the differences are restricted to a small subsest of all variation). Concentrate on the Rosenberg paper, but the Pritchard paper is a natural companion because the Structure software is published there. Don't get too caught up trying to figure out all the names of different indigenous groups, it is sufficient to think of each as a distinct genetic population. Pay close attention to the two beautiful figures, including all the details of mixed inference and what the different K values mean. (You will need to read part of the Pritchard paper to fully understand the K values.) 

       Here are a few questions to think about: how can you explain the preponderance of diversity shared by all groups (remember there were sometimes multiple rounds of small founder groups)? Similarly where do you think the population-private alleles come from? How can you explain the groups that prominently have two colors in Figure 1 (e.g. Mozabite and Hazara)?

        Do take a pass at the Pritchard paper, in particular to understand the K values and the methods used to cluster populations, but the method is complex so don't get too bogged down in the details. [Ignore the font changes here - my HTML editor has a bug that is preventing me from changing it even when I edit the HTML source code.]


      Primary Papers:

                Rosenberg et al. (2002). Genetic structure of human populations. Science 298: 2381-2385.

                Pritchard et al. (2000). Inference of population structure using multilocus genotype data. Genetics 155: 945-959.


        Additional Papers:

                Rasmussen et al. (2011). An aboriginal Australian genome reveals separate human dispersals into Asia. Science 334: 94-98.


Week 5 Lab - software download and test BEFORE class please.

Week 5 Lab - Structure