Link to this post | posted 05 May, 2017 16:11
welkin	You need to include the amino acid and the anticodon for tRNAs for GenBank ready minimal files. If both programs are undetermined, you can use "tRNA other" and the anticodon. If tRNA scan has made up its mind, and ARagorn hasn't, use the tRNA scan decision.

Link to this post | posted 05 May, 2017 15:33

welkin

Hi Ann, thanks for the data and sorry about the lengthy wait for a reply! for 93: There are too many equally likely hits that do different things— I think the best choice is to pick HTH DNA binding protein. There are antitoxins, transcriptional regulators, and all kinds of these. HHPred is picking up that structural motif, and in this case it could be lots of things. TA systems are really difficult to untangle— the proteins are tiny, so the alignment programs don't work as well. Without bench data, I think we should be cautious. Best, Welkin

Link to this post | posted 05 May, 2017 15:06

welkin

Hi Larry, I moved you to a new thread so it would be easier to find for other people annotating B2 phages. In general, each new function question should be a new thread. For your questions: ORF 2: queuine tRNA-ribosyltransferase, Blastp Phamerator, e=0.0, Glass gp2 yes, go ahead. we should put it on the big list. ORF 6: GTP cyclohydrolase, Blastp Phamerator, e=e-108, LizLemon gp6 yes, go ahead with this one too. ORF 8: peroxide stress protein YaaA, HHPred, prob=99.81 e=2.3e-20, E. coli (HHPred coverage with this match is only 75% of the query sequence.) No. ORF 18: membrane protein, Blastp Phamerator, e=1e-18, Attis gp21 Not certain that this is the same as membrane protein Band-7-like that is in the list for Functional Assignments. The Attis protein was labeled as a membrane protein after it was examined with TMHMM. If you get a similar output, go ahead and call this one a membrane protien too. It is not the same as the Band 7 proteins. ORF 42: tail fiber protein, Blastp NCBI, e=3e-08, Mediterranean phage uvMED (Should this just be called a minor tail protein?) This one is kind of scrawny to be a tail fiber. Let's hold off until we have some bench data. Thanks Larry!

Link to this post | posted 26 Apr, 2017 16:09
welkin	Hi Ann, Thanks for posting the HHPred text output. However i would also need to see how much of your protein aligns to the various matches— probability and e value alone don't get you there, as you can get some pretty high scores for pretty short motif matches. Best, Welkin

Link to this post | posted 24 Apr, 2017 15:34
welkin	5) I like the MRE11 double-strand break endo/exonuclease. We'll add that to the main list. 6). I think we should just stick with HTH DNA binding protein. I know I am wimping out about the membrane part, but i think it is better to stay vague when we don't know what it is doing. Thanks JAckie!

Link to this post | posted 29 Mar, 2017 13:35

welkin

Wow, Larry, you sure don't pick easy ones! I may have to break these up into several posts because it is taking me a while to get through them. First of all, the integrase. This is a tyrosine integrase, rather than a serine (which is the Bxb1 integrase Debbie was referring to). We have several tyrosine integrase starts characterized at the bench, including L5. If I do an ideal alignment using Smith-Waterman at ENI between the two, it looks like the cluster T integrase should be around 1266 in your Nairb. This start lines up better with the Cluster F integrases that have been called on phagesdb and is start "20" on the Starterator report. I'll post others as I get a chance. I am also going to post your .dnam5 file, so others can follow along if they need to. -Welkin 30Kb

Link to this post | posted 20 Mar, 2017 17:54

welkin

Our RNA seq and mass spec data from various clusters supports the guiding principle that the vast majority phage genes and genomes are only transcribed in one direction for any given region. So while I don't have data for this particular instance (and there are always exceptions to every principle), the data that we do have says that you need to choose one or the other when it comes to two genes occupying the same piece of DNA in different frames.

Link to this post | posted 13 Mar, 2017 20:32

welkin

Yes, Greg, from this picture I think you need to delete the tRNA. I can't tell what cluster this is from the picture, but I am guessing there are many other cluster relatives that likely have the same prediction, and the tRNA was deleted. With no bench evidence one way or the other, I am inclined to go for the protein encoding gene that fits into the operon nicely, and to stay consistent with other cluster relatives. Best, Welkin

Link to this post | posted 02 Mar, 2017 15:47

welkin

Hi Jackie, 1) Graham has approved DNA polymerase I, and it is OK to use for those Gordonia phages. 2) I think we want to stay with glucosyl transferase and galactosyl transferase for those two enzymes in the Fs. For more specific I feel like we should demonstrate bench evidence. 3) DNA helicase/methylase is approved as long as you can demonstrate two functional domains. 4) Immunity superinfection protein is too vague based on what we currently know about prophage mediated resistance, and stays off the approved list. I need to go back and figure out how those functions got into those genomes to begin with. 5) Re the MRE11 (what a neat investigation): the HHPred results are pretty confusing, so I can see why this protein has been labeled as both exo- and endo- or nothing at all. What are you proposing? Call all members of this pham "nuclease"? I Only slightly farther down the list are almost equally good alignments to "metallophosphoesterase" and "phage lambda ser/thr kinase". Are you sure that the evidence is good enough to pick one thing out of this list where everything is 98.7 probable or better to a single domain of the protein? I would have to do some serious digging through the literature to figure out which parts of the proteins were aligning to the phage proteins. 6) Excellent use of TMHMM. Which type of DNA binding domain?

Link to this post | posted 01 Mar, 2017 18:56

welkin

Hi Beckie— I love the Gordonia phages! These areas flanking the integrase seem particularly plastic and so it doesn't surprise me that the normal guidelines aren't holding up. Based on your GenemarkS output (which I think is the most compelling evidence you've got here) I want to keep the forwards 43 that is in the middle of that reverse operon. I also think you are right; you make sure that you call the reverse gene to start at 29163 and no longer to accommodate the forwards gene, and delete the other reverse gene that occupies the same space as the forwards gene. I could be persuaded to change my mind though, if you run them all through HHPred and find functions for the reverse gene that you propose to delete. -Welkin