SEA-PHAGES Logo

The official website of the HHMI Science Education Alliance-Phage Hunters Advancing Genomics and Evolutionary Science program.

Welcome to the forums at seaphages.org. Please feel free to ask any questions related to the SEA-PHAGES program. Any logged-in user may post new topics and reply to existing topics. If you'd like to see a new forum created, please contact us using our form or email us at info@seaphages.org.

All posts created by debbie

| posted 21 Feb, 2024 18:20
Hi Beth,
What phage are you annotating? Could you provide the phamerator gene number and its stop codon.
Thanks,
debbie
Posted in: Functional AnnotationMetalloprotease without HEXXH motif?
| posted 17 Feb, 2024 13:54
Fred,
I think you have a good choice!
Best,
debbie
Posted in: Gene or not a GeneIs the current subcluster A4 pham 4236 (as of 2/14/2024) really a gene?
| posted 16 Feb, 2024 17:06
H John,
This is on my plate to fix. They will be in phamerator soon
debbie
Posted in: Cluster FC Annotation TipsAnnotated Cluster FC Phages
| posted 15 Feb, 2024 20:12
I don't think that it matters if they are ATG, GTG, or TTG.
Starts are starts and the amino acids chain start with an fMet.
debbie
Posted in: Choosing Start SitesTwo consecutive ATG start sites
| posted 15 Feb, 2024 15:17
Hi Fred,
The right arm of some of our well studied mycobacteriophages defy common annotation practice.
Here is how I look at it. We have published papers where the right arm has no coding regions because there are RNAs expressed. The paper about mycobacteriophage Giles has the best display of bench data that we have.
I am always cautious to not overcall the right end of a genome. I think these are some of the stats to consider:
This ORF is only found in Cluster A4s.
Phamerator has it as 17/134 genes annotated. Are they all drafts? How many of them have the nucleotide sequence? How conserved is it? Could it be something other than a protein coding sequence?
You mention no coding potential. Nor, is it easy to see how this gene would get easily made.
Check out other cluster A phages like L5 or D29 and note they are not over-calling this region.

Also, what is the risk if you delete it? Do you think a researcher who want to investigate lemur_85 isn't going to find this segment of DNA in Alberto7.

Most important is that you can support whichever call that you make. Looking at all of the data is and not what others have done will get you to a best answer.

Good luck as you make this decision!
debbie
Posted in: Gene or not a GeneIs the current subcluster A4 pham 4236 (as of 2/14/2024) really a gene?
| posted 14 Feb, 2024 20:19
Hi Brandon,
I think that what you are asking me is to decide if the appropriate functional call is Holliday junction resolvase or RuvC-like resolvase, yes?
What do you think it should be? What is your evidence Have you read any papers about them? How are those two functions related? What did your professor think?
debbie
Posted in: Functional AnnotationFunction calling problem
| posted 12 Feb, 2024 19:29
Hi Elizabeth,
The issue is that Cluster FC phages have just recently been found to have direct terminal repeats. Originally they were thought to be circularly permuted, which created a different sequence and end determination. What that means is that all genomes with old (incorrect) sequences at the ends) had to be reposted. That means that all phages look like draft in phamerator.
It also means that the 3 that were final (Atuin, Racecar and Mimi) have to be converted to their new sequence. and revisions sent to GenBank and phamerator. That work is not quite done.

But my challenge to you and your students is that, when the first 3 were done, the comparative data was sparse or non-existent. You have the opportunity to make calls with significant amounts of data and no distractions of what others have interpreted. Currently, with the data for 8 Cluster FC genomes, better information can be gathered from the comparisons - like starterator - to make really good choice about is it a gene, what is its start, and what is its function, without contradicting what others said with less than optimal (or no) comparative data available. Consider embracing this challenge because it is all good!

debbie
Posted in: PhameratorCluster FC in Phamerator
| posted 09 Feb, 2024 21:52
Hi Emily,
Both are good questions.

For the first gene, the coding potential on the GeneMark output is rather robust upstream of 37. How would you evaluate the ability of that gene to start if you are evaluating a start at base 1. What do you want to see upstream of a start to evaluate the start? How would you do that in this case? Why do you say that RBS only favors the second start?

For the second gene, how is that gene going to be made? What favors the best production of that gene?

Keep thinking!
What did you professor have to say?
best,
debbie
Edited 09 Feb, 2024 21:54
Posted in: Singleton Annotation TipsPensacolaC28 Start calling question
| posted 09 Feb, 2024 03:37
Hi Brandon,
More information is needed. What gene (please include coordinates) and genome are your inquiring about?
thanks,
debbie
Posted in: Functional AnnotationFunction calling problem
| posted 09 Feb, 2024 03:33
Hi ACMPhageHunters,
More information is needed to help.
What gene (please provide coordinates) and genome are you inquiring about?
Thanks,
debbie
Posted in: Functional AnnotationBaseplate protein vs. Minor tail protein