CSB352: Lab Report 1
This assignment asks you to investigate a gene of unknown function from Arabidopsis
thaliana, by identifying orthologous sequences from other organisms and performing a
multiple sequence alignment of your discovered orthologs. In a nutshell, your lab report
will introduce your gene (see arabidopsis.org to as described in the Gene Assignment
section of the Lab Reports part of Portal; details on the gene/protein are also available
from this database), give explicit details of the analysis you performed (BLAST and
MSA), describe the particulars of what you found and comment on what you think it
means in a biological context – what could the function of this gene be? You will also
perform a phylogenetic analysis on a subset of BLAST, generating the best-resolved (i.e.
greatest number of nodes that are well supported based on their bootstrap scores)
phylogenetic tree you can. Is there evidence of a gene duplication event in a recent
ancestor of Arabidopsis thaliana and the other species whose sequences are homologous?
The report should be no longer than 6 pages, double-spaced (excluding figures and
references)! You should cite all primary sources plus all tools/websites/data sources used,
using the journal Nature’s reference format.
Things you need to keep in mind
about your query:
-
are you going to BLAST a nucleotide or protein sequence (or both)? Genomic
DNA or mRNA? Is your sequence long enough?about your search database:
-
do you want an all-encompassing database or a well-curated one?
about your BLAST parameters:
- which BLAST search algorithm is best for finding orthologs in other species?
- are you going to use blastn, blastp, blastx? Why?
-
which scoring matrix should you use?
about your alignment parameters:
- which alignment program are you going to use?
- how rigorous are your gap penalty parameters?
-
does your alignment“make sense”?
about your phylogenetic analysis:
- which method(s) are you going to use?
- include the closest A. thaliana paralog in the analysis (if there is one)
- note that duplicated genes are typically not identical and that the“nonredundant”
database does actually have redundancies in it.
CSB352H Bioinformatic Methods Winter 2022 – Lab Report 1
Your report should consist of the following sections:
Introduction: (4 marks)
A short, two-paragraph description of your gene and the analyses performed.
Be sure to include the following information:
- organism
- sequence source (i.e GI or AGI ID number, sequence type)
- physical location (if available)
- functional product (if known) – pI, MW, and other parameters as appropriate
Summarize the genes and taxa included in your analysis. You can make note of any
pertinent alignment patterns you detected in your MSA (ie: conserved regions, range of
sequence identity, etc.), and state what you are trying to ascertain with your phylogenetic
analyses (do not copy the above objective!).
Methods: (8 marks)
A concise explanation of the parameters, data sources, web sites used and their citations,
and settings for each in silico experiment performed in the final analysis you’re
presenting in this report. Note: The marker must be able to replicate your results solely
from the information in this section!
Be sure to detail the BLAST search(es) you performed along with any parameter
adjustments you made and justify them. Include the Blast version and databases version.
Describe which alignment program you chose as well as all parameters used in the
analysis. Explain the design of your phylogenetic analysis in prose form, including the
software used, the alignment and phylogenetic tree-building algorithm(s)
used (justify your choices), the parameters you chose for each algorithm. DO NOT
INCLUDE RESULTS IN THIS SECTION (see below.)
Tip: The literature accepted E-value cutoff for accepting a sequence as orthologous to
the query sequence is 1e-20 as per Lecture 2. If your assigned gene returns few or no
orthologs in the search, then contact abraham.yang@utoronto.ca to request a new gene.
Results: (12 marks)
Your results will include a short discussion on what was found and anything you did that
was“out of the ordinary”and why. Also it will include 3 figures, complete with proper
captions. Each figure should include a caption that describes the figure. The caption should
be detailed enough that the readers do not need to read the text in order to understand the
figure:
1) a table of your orthologous sequences (and the most similar A. thaliana paralog if
CSB352H Bioinformatic Methods Winter 2022 – Lab Report 1
there is one) and their associated BLAST statistical values. Be sure to justify your
choice of sequences you selected for alignment (are all sequences truly
“nonredundant”?), and include only the sequence retrieval (do not include pairwise
alignments, please).
2) your best multiple sequence alignment of the top 15 orthologs from your BLAST
search – include only one other A. thaliana sequence that is not the same as your
assigned gene, that is, include the closest Arabidopsis thaliana paralog, if one exists.
Make sure the MSA is formatted in such a way that conserved regions, residues, or
nucleotides are highlighted. Trim it accordingly (i.e.: don’t include uneven tail-ends,
as they are uninformative.)
3) A write-up in which the results of your phylogenetic analysis are displayed and
described (but not analyzed). Include within this page your preferred phylogenetic
tree with branch lengths and bootstrap values. In the text, explain any pertinent
trends, divergences, or similarities seen in your tree.
If your alignment or tree viewer program doesn’t generate image files to import into your
report document, use screen shots with the PrintScrn button then paste followed by
cropping.
Ensure that the gene name and organism name (abbreviated if need be) is readable for
each sequence in the alignment. Your query sequence should be at the top, followed the
closest paralog and orthologs sorted from most to least similar.
Discussion: (6 marks)
Interpret your results in a biological context. Can you suggest a possible function of your
gene? Do its relationships with other taxa make sense to you? Is this gene widespread in
the tree of life (or at least in plants)? Can you thus infer whether it is relatively old or
new? Is it highly conserved across taxa? Are some regions highly conserved and others
not? What might this tell you about its function? Do any species have duplicated versions
of this gene? When did the duplication occur?
DUE ON FEBRUARY 18, 2022
SUBMIT ONLINE