Phylogeny inference and data analysis, spring 2011
Lecturer
Sirkka-Liisa Varvio
Scope
4-12 cu.
Time schedule, prerequisites, content
Basic probabilistics and statistical inference are assumed as prerequisites.
The course is tailored for Bioinformatics Master program (MBI) and thus some knowledge on sequence analysis, working with sequence databases, sequence alignment, distance matrix based phylogenies by MEGA-software, are assumed to be known.
The course is, however, open for all students: additional sessions will be arranged for students who are not familiar with these basics (introduced in other MBI-courses).
Part I
III period:
Tuesdays 15-17 and Thursdays 14-16 in B120 and in computer class C128.
Time schedule for extra sessions (for non-MBI -students) in computer class during three first course weeks in negotiable (for example before or after lecture times).
...
Exam (40 points) and homework assignments (20points): 4cr, essay (statistically or biologically focused): 2cr.
50-60 points: grade 5, 30 points: grade 1.
Course is focused on practical working => exam in computer class.
Week 1 program is for non-MBI (Bioinformatics Master program) students. Lectures and practicals in computer class overlap with topics which have in been in Introduction to bioinformatics / Molecules for bioinformatics: Basic phylogeny concepts, basic practical computer working and getting familiar with molecular sequence databases.
18.01., 15-17, computer class C128.
Week 2
Tue 25.01, 15-17, B120, Description of the course, assignments, how to proceed with data collection.
Thu 27.01, 14-16, B120, Alignment with Clustal, editing with Genedoc, MEGA
Week 3
Tue 01.02, 15-16, C128, Practical session with MEGA, for training simple phylo with MEGA.txt, for training simple phylo FASTA-format.txt
Wed 02.02, 12-20, C128, Extra help for those that need, data collections and alignments for datasets 1 and 2.
Thu 03.02, 14-16, D340, Phylogeny books and programs, Distance matrix methods (this is given as a paper copy), Nucleotide substitution modelling.
Week 4
Tue 08.02, 15-16, C128, MrBayes demo, checking that everybody can start with the program. Use this NEXUS-file for training: training.nex.
Wed 09.02, 12-14, C128, Extra help, checking results,dataset 1 should be ready for MrBayes analyses (i.e. you have already done MEGA-analyses) and dataset 2 collected.
Thu 10.02, 14-16, B120, Parsimony, Maximum likelihood.
Week 5
Tue 15.02, 15.00-16, C128, Datasets 3 and 4, Splitstree and Network.
Thu 17.02, 14-16, B120, Lectures on methods continue.
Week 6
Tue 22.02, 15-17, C128, Results from datasets 1 and 2 + lectures on methods (2h session probabaly not enough, we continue until everything is done.)
Thu 24.02, 14-16, B120, Results from dataset 4 + lectures on networks and phylogeny examples (phylodynamics, phylogeography, phylogenomics, phyloprofiling).
Week 7
Thu 03.03, 16-20, C128, EXAM.
Lectures, links, note also the material in "software tools" | Recommended review (and other) papers
|
|---|
Phylogeny books and programs..pdf,
all phylogeny programs from here:
http://evolution.genetics.washington.edu/phylip/software.html
| Before phylogenies, seq collection and alignments:
NCBI BLAST.pdf, Multiple sequence alignment review.pdf,
Collection of tools.pdf, includes links to many aligning,
editing, and phylogeny tree visualization facilities.
|
Distance matrix methods in phylogeny inference
(paper copy in lecture Thu 03.02)
Maximum parsimony phylogeny inference..pdf
Maximum likelihood phylogeny inference..pdf
Bayesian phylogeny inference..pdf
Description of Bayesian inference relevant in MrBayes..pdf
Networks
| Short tutorial article.pdf
.
.
Traditional and Bayesian phylogeny estimation.pdf
Original MrBayes paper.pdf , MCMCMC.pdf
.
.
.
Phylogenetic networks review.pdf , Median-joining networks.pdf
|
Nucleotide substitution modelling..pdf,
Practical tool for model choice:
http://www.hiv.lanl.gov/content/sequence/findmodel/findmodel.html,
(see also Modeltest.pdf)
| Model selection review.pdf
Bootstrap confidence in detail by Efron.pdf
Bootstrap vs posterior prob 1.pdf,Bootstrap vs posterior prob 2.pdf,
Bootstrap vs posterior prob 3.pdf, Posterior prob how meaningful.pdf,
Bayesian inference of character evolution.pdf..
The Bayesian revolution in genetics.pdf |
Examples
|
|---|
Phylogenies from non-sequence material, Worms with legs_Nature240211.pdf, Suppl240211.pdf,Language phylogenies.pdf
http://www.timetree.org/ , http://tolweb.org/tree/, History of molecular clock.pdf , Origins of eukaryotes.pdf, Origins of mitochondria.pdf, The ring of life.pdf
Human origins, A new view of the birth of Homo sapiens.pdf, Finger from Siberia April2010.pdf,Tooth from Siberia December2010.pdf
Phylogenomics.pdf, Pharmacophylogenomics.pdf
Phylodynamics, Phylodynamics review.pdf, Phylodynamics of viruses.pdf, Phylodynamics of influenza viruses.pdf,
Microbiome phylotyping.pdf, http://www.mlst.net/
|
Working in 2-4 student groups. You get 20 points for the total 60 points (40 from exam) by doing these. Note that the course is focused on practicals and you need these skills in exam (reasonable aspects, not sophisticated details).
Dataset 1
Assignment 1 - Instructions for data collection and alignment.pdf , results from neigbor-joining, UPGMA, parsimony phylogenies during weeks (3) and 4 and bayesian results during weeks (4) and 5, whole work ready Tuesday 22.02.
Geneseq not aligned FASTA.txt, Geneseq aligned FASTA.txt, Geneseq in MEGAformat.txt IN CASE YOU WOULD LIKE TO ANALYSE YOUR FLU-VIRUS DATA BY NETWORK-SOFTWARE, SEND YOUR FASTA-FILE TO SIRU TO GET IT IN RDF-FORMAT.
|
Dataset 2
Assignment 2 - Instructions for data collection.pdf, checking data-collection in and preliminary results in weeks (3) and 4. Whole work ready Tuesday 22.02.
Data to be collected from http://www.ncbi.nlm.nih.gov/genomes/FLU/FLU.html
|
Dataset 3
Human gene and its alleles. A change in the plan: The 20 points will be given from assignments 1, 2 and 4. The human examples are illustrated during lecture (you need not do anything in practice).
|
Dataset 4
Data from three bacteria. Select one of them for phylogeny analysis. Compare neighbor-joining (by MEGA) and network phylogenies by Splitstree4 (FASTA-file) and Network4.6 (rdf-file). You will get rdf-files also by email as they are probabaly not ok here.
See your email for instructions. Additional instructions here: Practical instructions for using the Network-software.pdf
Borrelia_ClpAFASTA.fasta, Borrelia_ClpA.rdf
pneumococcusGkiFASTA.txt
StaphylococcusarccFASTA.txt
|
Sequence alignment
The necessary step before phylogeny inference is to align sequence data, inspect the aligned data and correct it. Sequence aligning, like working with sequence databases, are topics in other MBI-courses, and theoretical basics will not be in the program of this phylo-course.
|
Phylogeny software
These program packages have been installed in computer class C128. All these will be used during the course.
Note the links from which you can download the program to your own computer if you prefer working at somewhere else than in C128.
Note also the manuals!
- MEGA4, http://www.megasoftware.net/, Manual: MEGA4 Manual.pdf
- NOTE: MEGA-link now (since 25.Jan 2011) goes to a new version MEGA5 which has been released. As this happened so recently, working in classroom C128 will be with MEGA4, administrator will not be asked to install a new version.
- MrBayes3.1, http://mrbayes.csit.fsu.edu/, Manual: http://mrbayes.csit.fsu.edu/mb3.1_manual.pdf, MrBayes Commands.pdf
- Network4.6, http://www.fluxus-engineering.com/sharenet.htm, Manual: http://www.fluxus-engineering.com/Network4600_user_guide.pdf
- Splitstree4, http://www.splitstree.org/, Splitsree4 Manual.pdf
|
...
Part II
IV period: Practical project work in computer class and at home, time schedule negotiable (4-6 cr), topics to negotiated at the end of part I
The project is based on sequence data collected from databases.
...