TIPS IN SEED
Defining gene names
- Role1 @ Role2
A monofunctional enzyme that can catalyze two different reactions, both role 1 or role 2
- Role1 / Role2
A bifunctional enzyme that does two unrelated roles, role 1 and role 2
- Role1 ; Role2
Does one of either role 1 or role 2 but we don't know which
There are several occasions where you are sure what something is not but you don't know what it is. For example, there is a protein in your spreadsheet that is a duplicate, and you are sure that it is the wrong one because it is not a cluster. To remove it, you can add "But Not" to the end of the annotation. This is very useful as you can see what things have been checked and shown not to be something!
- Suppose that you wish to rewrite all assignments of the form homoserine kinase with Homoserine kinase (EC 184.108.40.206).
To do this, you first go to some fields in the initial page that have a button just beneath them called "Generate Assignments via Translation". You fill out three fields:
- About 1.5 inches up there is a field called "Save as user". You put your user name (without master: in front of it) in that box
- The From box gets "homoserine kinase"
- The To box get "Homoserine kinase (EC 220.127.116.11)"
You click on the "Generate Assignments via Translation". This builds something called an "assignments set" under your user name. Note that no assignments have yet been made.
- Then, you need to process the assignments set. You do this by
going to the search form about 1.5 inches below where you were to were
it says "Process Saved Assignment Sets". You fill in your user name (WITH the master:) and click on "Process Assignment Sets".
- You will see the assignment sets that have been saved, you select one, and you "process it". Processing an assignment set amounts to
- checking it to make sure that you want to make those
- deleting the ones you don't want to make
- making the assignments
- deleting the assignment set
A very useful tool for managing the variants of your subsystems.
Close your subsystem before using this tool to change the variants codes.
Here is how to make a tree and use it to edit annotations, as described by Ross.
- Pick a column in one of your subsystems that contains, say, 20-50 PEGs.
For myself, I am choosing TcuB: works with TcuA to oxidize tricarballylate to cis-aconitate in Tricarballylate Utilization.
- Go down to For sequences in a column (i.e., role): and pick the appropriate column. Then click on Align Sequences in Column
This may run for a while, since the program must try to construct one or mor multiple sequence alignments.
- When it comes back, go to the very end of the subsystem stuff, and you should see one or more trees shown. They all have numbers. Mine is "2.1" where the "2" means "role 2" and the ".1" means the first alignment (of possibly several) constructed from the sequences in the column.
- Now go up to Realign subgroup within a column (adding homologs):
- Type in your tree number (in my case, 2.1)
- In the box labeled Include homologs that pass the following threshhold: and type in "1.0e-20" (which means look for sequences that are similar to those in the column at a cutoff of 1.0e-20)
- Now click on Realign Sequences in Column
- When the page comes back, go to the bottom, select the tree you made (in my case "2.1"), and click on use_tree.
- You are now in an environment where you can look at the phylogenetic tree, you can reroot it, or you can make assignments using it.
Here are some things you can do with the tree
- If you go and change some annotations (by clicking on the leaves), you may wish to just redraw the tree with the annotations updated. That is what Redraw tree does.
- To reroot the tree, check a single node, and click on Reroot tree. This will reroot the tree between the checked node and its parent.
This is very useful as you try to make the tree "coherent"
- You can delete subtrees by checking their root nodes and clicking on Delete subtrees.
- You can assign to all nodes in a subtree by checking the node, pasting in the desired assignment (at the bottom) and clicking on assign.
- Alternatively, you can check a node (not a leaf) and a single leaf. Then, clicking assign will assign the function from the leaf to every leaf in the subtree below the checked nonleaf.
- You can ask that homogeneous subtrees be collapsed (this is useful as you progress in a big tree to making things consistent)