Link to this post | posted 05 Aug, 2022 14:56
kcornely	Karen, this is very helpful, thank you. An example of a RecB-like protein would be very helpful. I understand your thought process and your logic as identifying this protein as a Cas4 endonuclease. But why do so many other P1 cluster phages assign the function as RecB-like exonuclease/helicase? Shouldn't proteins in the same pham all have the same function? Thanks again for your reply, Kathleen

Link to this post | posted 03 Aug, 2022 18:48

kcornely

Greetings, all! I'd like to come back to my previous post regarding Starterator. We are encountering quite a few genes in Frankenweenie in which the suggested Starterator start is not the same as the most annotated start. I give an example above; here's another example: For gene 30, the original Glimmer call @bp 15143 has strength 3.13 ** not called by GeneMark; ST: start @15143 is the suggested start, but start @15146 has five most annotated. If we choose the start at 15146, which we believe is correct, we obtain q1:s1 data for all five BM cluster phages, with an alignment of 98.8%. In addition, the 15146 start produces a -1 bp overlap and expresses a gene length of 249 bp. All of the other five BM cluster phages that express the gene of this pham (7493) are also 249 bp. The Glimmer/Starterator start at 15143 provides a -4 bp overlap. So this gene has both -4 and -1 bp overlap options and the advice we got from the forum is to go with the -1 bp overlap. So we believe that the most annotated start, 15146, is correct. But we don't understand why Starterator didn't call the most annotated start. As always, thanks for helping me to understand Starterator better! Best wishes, Kathleen

Link to this post | posted 03 Aug, 2022 18:32

kcornely

Greetings all! We generated our R2I for Langerak and were pleased to see that we had assigned all functional assignments correctly except for one: Gene 46 (31633-32562). This gene belongs to pham 37230. On Phamerator, all of the other most closely related P1 cluster phages (Donovan, FirstPlacePfu, HUHilltop, Jebeks and Techage) all have this same pham, and the assigned function is RecB-like exonuclease/helicase. On the SEA PHAGES function list, we are told the following: RecB-like exonuclease/helicase: If both a helicase and nuclease domain are present, the RecB label should be used. Cas4-family exonuclease: This family of exonucleases is similar to the exonuclease domain of RecB. The Cas4 label should be used if the gene includes only the exonuclease region. IF the gene also includes a helicase domain, the RecB label should be used. Cas4 family nucleases tend to have alignments to the crystal structure 4R5Q_A, 41C1_A and to the PD-(D/E)XK nuclease superfamily (PF12705.7, among others) On HHpred for Langerak, we did get hits to the Cas4 family nucleases mentioned in the above paragraph. But we also got hits to helicases: 3U4Q_B, 4CEI_A, 6PPU_B, 1W36_E, and others. We decided to assign the function as the RecB-like exonuclease/helicase because (1) we thought we had fulfilled the requirement stated on the SEA PHAGES official function list that we should have hits to both helicases and exonucleases and (2) All of the other P1 cluster phages with the same pham used this functional assignment on Phamerator. But we were marked wrong. I'd like to get this right going forward, so could someone please help me out? As always, thanks for helping me "review to improve"! Best wishes, Kathleen

Link to this post | posted 11 Jul, 2022 20:36

kcornely

Greetings, all! We are annotating Frankenweenie, a BM cluster phage and had some questions regarding the interpretation of Starterator data. Here's an example: The start of gene 144 is called by Glimmer and GeneMark at 89678. This is also the suggested start on Starterator. But the most annotated start on Starterator is 89666 for all of the other five BM cluster phages. It also results in q1:s1 BLAST data for all of the other five BM cluster phages and as an added bonus, produces a -4 bp overlap and a gene length more in line with the other five BM cluster phages. So our question is: why are we seeing so many "suggested starts" on Phamerator for this phage that are not the same as the "most annotated" starts? And why, when we do choose the most annotated start (rather than the SS), it seems the better fit? Thanks in advance for helping me to understand Starterator better!

Link to this post | posted 11 Jul, 2022 18:38

kcornely

Greetings, all! I have another comment about the 4 bp overlap. We are annotating Frankenweenie, a BM cluster phage that has only six members total. We are occasionally encountering situations where it doesn't make sense to call the start at the position that provides a 4 bp overlap. For example, for gene 95, Glimmer and GeneMark call the start at 65479, which provides a 4 bp overlap. But to us, it makes sense to call the start at 65470. The 65470 start (1) is not the SS, but is MA on Starterator, (2) it provides us with q1:s1 BLAST his to the other five BM cluster phages, and (3) the product is similar in size to the other five BM cluster phages. I'd love to hear your insights on this observation, and thanks in advance!

Link to this post | posted 17 Jun, 2022 02:00
kcornely	Thanks Chris and Debbie! This is incredibly helpful.

Link to this post | posted 16 Jun, 2022 15:46
kcornely	Thank you Debbie! This is very helpful.

Link to this post | posted 16 Jun, 2022 15:45

kcornely

Greetings, all! We are annotating Frankenweenie, a Streptomyces BM cluster phage. We haven't annotated a BM cluster phage before. We were wondering if gp65 and gp66 were tail assembly chaperones and if gp66 was translated as a frameshift. There are six BM cluster phages, and the only phage in which the frameshift is called is Satis. It is not called for the other BM cluster phages. A frameshift seems likely for Frankenweenie because there is considerable overlap between gp65 and gp66. Thanks in advance for your help! Kathleen

Link to this post | posted 15 Jun, 2022 21:20

kcornely

Greetings, all! We are annotating a BM cluster phage, Frankenweenie. I have a team of four students working on this annotation project this summer. This is the first time that we've annotated a phage other than a mycobacteriophage, so we will be posting a lot of questions to the forum this summer! Our first question: Do Streptomyces phages have the 4 base pair overlap that is common in mycobacteriophages? I'm looking at gp30 now, and the BLAST data and the Starterator report fit better with a -1 bp overlap start. Should we go for the -4 bp start anyway? Thanks in advance for your help! Kathleen

Link to this post | posted 14 May, 2022 17:35
kcornely	Thanks for the speedy reply, Mitch!