Answers.Body: <p>If you <a href="http://www.ensembl.org/Homo_sapiens/Transcript/Summary?g=ENSG00000012048;r=17:43093586-43125364;t=ENST00000477152" rel="nofollow noreferrer">look at the original data at Ensembl</a> you'll notice that most of these are labeled "CDS 3' incomplete" <s>and have a TSL (transcript support level) of 1, which is as low as it goes</s>. It seems likely that that is simply an incomplete annotation. I'm not surprised that there are a bunch of these for BRCA1, since there was a long time when there were a LOT of legal issues surrounding working on it.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-15T08:57:14.750Z
Answers.Text: @DevonRyan TSL1 is the most reliable, TSL5 is the least. Also, OP, there are some times of day when you might get an answer instantly, but mostly I would assume a full working day is a reasonable amount of time to wait before impatiently crossposting everywhere. You don't know where in the world the one or two experts who can answer your question are, and what time they'll get onto the web.
Answers.UserDisplayName: Emily_Ensembl
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-15T09:10:34.290Z
Answers.Text: @EmilyPerry Gah! Thanks for letting me know about my TSL mistake. How the heck did an annotated "incomplete" transcript get a TSL of 1?
Answers.UserDisplayName: Devon Ryan
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-15T10:58:57.443Z
Answers.Text: Still won't let me comment. Devon, TSL and incomplete are independent. It's TSL1 because there is at least one non-suspect piece of evidence that covers the whole thing. If all we've got is evidence of a fragment, then that's all we've got, but the evidence of that fragment is strong.
Answers.UserDisplayName: Emily_Ensembl
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-15T11:07:15.280Z
Answers.Text: @EmilyPerry Sorry, that's an annoying stackexchange thing, you need 50 rep to comment. Anyway, thanks for clarifying that TSL and "3' incomplete" are independent annotations!
Answers.UserDisplayName: Devon Ryan
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-16T00:33:45.990Z
Answers.Text: @EmilyPerry I am not a bioinformatician or a biologist, so I don't really understand what a "non-suspect piece of evidence that covers the whole thing" means. Does it mean that this sequence with no asterisk in the end was found as a protein synthesized by the gene? Probably not since there is no stop codon, right? Or does it mean that an mRNA that matches this peptide sequence was measured? The mRNA without a stop codon would be degraded, though.. Sorry if this is a dumb question, but I'd really appreciate some clarification on that.
Answers.UserDisplayName: user5054
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-15T07:25:59.873Z
Answers.LastActivityDate: 2017-08-15T09:11:38.693Z
Answers.LastEditDate: 2017-08-15T09:11:38.693Z
Answers.UserDisplayName: Devon Ryan
Answers.Body: <p>If it's 3' incomplete that means the evidence used to create it was a fragment. <a href="http://www.ensembl.org/Homo_sapiens/Transcript/SupportingEvidence?db=core;g=ENSG00000012048;r=17:43093586-43125364;t=ENST00000477152" rel="nofollow noreferrer">Here</a>'s the evidence used to construct BRCA1-214 ENST00000477152.5, a 3' incomplete.</p>
  
  <p>You can see that <a href="http://www.ebi.ac.uk/ena/data/view/AK307553.1" rel="nofollow noreferrer">there</a>'s a full length cDNA from EMBL, AK307553.1, which was used to create this model.</p>
  
  <p>The sequence was mapped against the genome sequence to create the transcript model. When searching for the translation, it was found that the open reading frame was started, but there was no stop codon. This suggests that the cDNA AK307553.1 is actually a fragment, that the mRNA was broken or cleaved before it was reverse-transcribed sequenced. We display it in Ensembl in the hope that this will lead someone to identify the full length transcript that it represents.</p>
  
  <p>Many people choose not to work with these incomplete transcripts. You can do this by filtering by Gencode Basic in biomaRt.</p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-16T08:49:48.813Z
Answers.LastActivityDate: 2018-01-10T07:56:14.690Z
Answers.LastEditDate: 2018-01-10T07:56:14.690Z
Answers.UserDisplayName: Emily_Ensembl
Body: <p><em>Note: this question has also been asked on <a href="https://www.biostars.org/p/267545/" rel="nofollow noreferrer">BioStars</a></em></p>
  
  <p>I am getting peptides using biomaRt library in R for the gene 'BRCA1'. 27 different aminoacid sequences are returned, and 12 of those have an asterisk (*) in the end, while 15 do not. What does that mean? I know that asterisk is used to indicate stop codon. So, does that mean the peptides without an asterisk do not have a stop codon? Is that possible at all? Are all of those valid peptides or not?</p>
  
  <p>Below is the R code I used: </p>
  
  <pre><code>library(biomaRt);
  mart = useMart('ensembl', dataset='hsapiens_gene_ensembl');
  seq = getSequence(id='BRCA1', type='hgnc_symbol', seqType='peptide', mart = mart);
  </code></pre>
  
  <p><code>MDLSALRVEEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKVLLCCPSWSTVVRS*</code> is an example sequence with an asterisk in the end. <code>MDLSALRVEEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKSLQESTRFSQLVEEL‌​LKIICAFQLDTGLEYANSYN‌​FAKKENNSPEHLKDEVSI</code> is an example without an asterisk.</p>
  
Tags: <gene><proteins><biomart>
Title: Why are some of the gene peptides returned by biomaRt missing an asterisk in the end?
ThreadText: # Why are some of the gene peptides returned by biomaRt missing an asterisk in the end?
  
  (Asked by: username_0 on 2017-08-15)
  
  *Note: this question has also been asked on [BioStars](https://www.biostars.org/p/267545/)*
  
  I am getting peptides using biomaRt library in R for the gene 'BRCA1'. 27 different aminoacid sequences are returned, and 12 of those have an asterisk (\*) in the end, while 15 do not. What does that mean? I know that asterisk is used to indicate stop codon. So, does that mean the peptides without an asterisk do not have a stop codon? Is that possible at all? Are all of those valid peptides or not?
  
  Below is the R code I used:
  
  ```
  library(biomaRt);
  mart = useMart('ensembl', dataset='hsapiens_gene_ensembl');
  seq = getSequence(id='BRCA1', type='hgnc_symbol', seqType='peptide', mart = mart);
  
  ```
  
  `MDLSALRVEEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKVLLCCPSWSTVVRS*` is an example sequence with an asterisk in the end. `MDLSALRVEEVQNVINAMQKILECPICLELIKEPVSTKCDHIFCKSLQESTRFSQLVEEL‌​LKIICAFQLDTGLEYANSYN‌​FAKKENNSPEHLKDEVSI` is an example without an asterisk.
  
  **Tags:** gene, proteins, biomart
  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_2 on 2017-08-16. Score: 8 | Accepted answer
  
  If it's 3' incomplete that means the evidence used to create it was a fragment. [Here](http://www.ensembl.org/Homo_sapiens/Transcript/SupportingEvidence?db=core;g=ENSG00000012048;r=17:43093586-43125364;t=ENST00000477152)'s the evidence used to construct BRCA1-214 ENST00000477152.5, a 3' incomplete.
  
  You can see that [there](http://www.ebi.ac.uk/ena/data/view/AK307553.1)'s a full length cDNA from EMBL, AK307553.1, which was used to create this model.
  
  The sequence was mapped against the genome sequence to create the transcript model. When searching for the translation, it was found that the open reading frame was started, but there was no stop codon. This suggests that the cDNA AK307553.1 is actually a fragment, that the mRNA was broken or cleaved before it was reverse-transcribed sequenced. We display it in Ensembl in the hope that this will lead someone to identify the full length transcript that it represents.
  
  Many people choose not to work with these incomplete transcripts. You can do this by filtering by Gencode Basic in biomaRt.
  
  
  ---
  
  ### Answer by username_1 on 2017-08-15. Score: 4
  
  If you [look at the original data at Ensembl](http://www.ensembl.org/Homo_sapiens/Transcript/Summary?g=ENSG00000012048;r=17:43093586-43125364;t=ENST00000477152) you'll notice that most of these are labeled "CDS 3' incomplete" ~~and have a TSL (transcript support level) of 1, which is as low as it goes~~. It seems likely that that is simply an incomplete annotation. I'm not surprised that there are a bunch of these for BRCA1, since there was a long time when there were a LOT of legal issues surrounding working on it.
  
  > **Comments:**
  > **username_2**: @username_1 TSL1 is the most reliable, TSL5 is the least. Also, OP, there are some times of day when you might get an answer instantly, but mostly I would assume a full working day is a reasonable amount of time to wait before impatiently crossposting everywhere. You don't know where in the world the one or two experts who can answer your question are, and what time they'll get onto the web.  
  > **username_1**: @EmilyPerry Gah! Thanks for letting me know about my TSL mistake. How the heck did an annotated "incomplete" transcript get a TSL of 1?  
  > **username_2**: Still won't let me comment. username_1, TSL and incomplete are independent. It's TSL1 because there is at least one non-suspect piece of evidence that covers the whole thing. If all we've got is evidence of a fragment, then that's all we've got, but the evidence of that fragment is strong.  
  > **username_1**: @EmilyPerry Sorry, that's an annoying stackexchange thing, you need 50 rep to comment. Anyway, thanks for clarifying that TSL and "3' incomplete" are independent annotations!  
  > **username_0**: @EmilyPerry I am not a bioinformatician or a biologist, so I don't really understand what a "non-suspect piece of evidence that covers the whole thing" means. Does it mean that this sequence with no asterisk in the end was found as a protein synthesized by the gene? Probably not since there is no stop codon, right? Or does it mean that an mRNA that matches this peptide sequence was measured? The mRNA without a stop codon would be degraded, though.. Sorry if this is a dumb question, but I'd really appreciate some clarification on that.  
  
  ---
  
  