CSB352编程代写、代写Python，Java程序

日期：2022-02-17 11:15

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