<p><strong>Conclusion</strong> The evidence suggests the region you are looking at <strong>Gabi-kat</strong> appears either 27045 downstream or 192250 upstream of syntaxin SNARE (SYP61), but is definitely not within the gene. It is chance homology (e.g. a pseudo-gene) that appears to occur within syntaxin (SYP61). The intergenic region could be much closer to the SYP61 gene if the same reference genome was used for both analyses Genbank LR782542.1 versus Genbank NC_003070.9 (different reference genomes in my analysis). Hence, the University is probably correct.</p>
<p>What I suggest is the genome of <em>A. thaliana</em> is downloaded and Blast analysis is performed locally. This way differences in the reference genome are removed.</p>
<p><em>The analysis is as follows</em>:</p>
<p>The target protein in question is involved in vesicle trafficking (<a href="https://www.ebi.ac.uk/interpro/entry/InterPro/IPR000727/" rel="nofollow noreferrer">here</a>) called SNARE, comprising coiled-coil domains and are highly conserved, for example show 100% conservation between <em>Arabidopsis thaliana</em>, <em>Arabidopsis suecica</em> and   <em>Camelina sativa</em></p>
<blockquote>
<p>MSSAQDPFYIVKEEIQDSIDKLQSTFHKWERISPDMGDQAHVAKELVATCGSIEWQVDELEKAITVAAKDP</p>
</blockquote>
<p>If Gabi-kat was a protein its sequence would be ...</p>
<blockquote>
<p>RISPDMADQAHVAKELVATCGSIEWQVLLFTHGSINPLFIY</p>
</blockquote>
<p>Here's a blast of the whole regions of the Gabi-kat. The region in green is the SNARE proteins ...</p>
<p><a href="https://i.sstatic.net/uBRgA.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/uBRgA.png" alt="enter image description here" /></a></p>
<p>So the issue is what is the terminal 3' protein and Blast that, i.e. the final bit that matches no known plant and the locus you flagged as being anomalous ..</p>
<p><a href="https://i.sstatic.net/6WqUm.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/6WqUm.png" alt="enter image description here" /></a></p>
<p>This is basically zero result via Blast, the protein with low homology are lower metazoans (animals) and a few bacteria ... a long way from plants. In other words this isn't a known protein of natural origin at least on NCBI.</p>
<p>The locus is literally like no protein on earth.</p>
<p>However ... Blasting the nucleotide ...</p>
<ol>
<li><p>Reference <em>Arabidopsis thaliana</em> = syntaxin (SNARE; SYP61)</p>
</li>
<li><p>Gabi-Kat <em>Arabidopsis thaliana</em> top two Blast hits are ...</p>
</li>
</ol>
<ul>
<li>Chromosome 1 at 10039082 - 10039208, or 9819787 to 981991 or</li>
<li>Both value    1e-42</li>
<li>Both dentity is 93%</li>
</ul>
<p>So Gabi-Kat is a real part of the <em>Arabidopsis thaliana</em>, <strong>Gabi-kat is a non-coding region</strong>. <em>Thats the major conclusion</em>.</p>
<p>The question where is the location of SYP61? And that is</p>
<ul>
<li>Chromosome 1 at 10012037 to 10020311</li>
</ul>
<p>Described here from <a href="https://www.ncbi.nlm.nih.gov/gene?Db=gene&amp;Cmd=DetailsSearch&amp;Term=839749" rel="nofollow noreferrer">NCBI's full report</a></p>
<p><a href="https://i.sstatic.net/u46cD.gif" rel="nofollow noreferrer"><img src="https://i.sstatic.net/u46cD.gif" alt="enter image description here" /></a></p>
<p>Thus the mutant is 27045bp downstream or 192250 bp upstream <strong>IF</strong> the same reference genome is used for the blasting (Genome LR782542.1) and the Genbank report NC_003070.9. They could easily be much closer to each other is the same reference genome was used for both sets of analysis.</p>
<hr />
<p>To answer the comments</p>
<ul>
<li>I [the OP] do not follow your [my] blast analysis</li>
<li>My [the OP] blasts gives ... ['hits within the gene']</li>
<li>etc ....</li>
</ul>
<p>It appears to me there is a psuedo-gene phenomenon occurring which is making you think a mutant is located within a gene. However, if you disagree thats perfectly okay and the original version of this answer I had assumed the Gaba-kat team had got this wrong (I always from the basis of supporting the OP).</p>
<p>NCBI was the database repository in the remote analysis, but there's an issue with multiple references ... so to reiterate I suggest you perform this analysis as a local blast to remove multiple reference phenomena.</p>
<p>To conclude I personally don't think there is a basis to believe the University responsible for Gaba-kat has got this wrong <a href="https://www.gabi-kat.de" rel="nofollow noreferrer">https://www.gabi-kat.de</a>. However, if you feel otherwise then I would recommend contacting them with your analysis and any supportive analysis for example from this thread or elsewhere. Please note, I don't have a connection with this university.</p>

CC BY-SA 4.0
2022-12-26T08:21:45.983Z
I do not follow your BLAST of the Gabi-Kat Nucleotide Seq. [My blast](https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Results?tl=5YJb4Lf7RtlEb7qc-21959324-2406285229) gives a 92.9% alignment for the entire stretch (127 bp long) with SYP61 (1:10016601 to 1:10016721). E = 4.6e-45
AvadaMouse
CC BY-SA 4.0
2022-12-26T09:05:03.670Z
Indeed, SYP61 is [located](https://bar.utoronto.ca/thalemine/portal.do?externalids=AT1G28490) from 1:10015856 to 10018047 in NC_003070.9 (GCA_000001735.2; TAIR v. 10.1). Which is inclusive of 1:10016601 to 1:10016721. You have got the specifics wrong.
AvadaMouse
CC BY-SA 4.0
2022-12-26T09:55:49.997Z
Running a nucleotide blast on the underlined part of the NCBI ref-seq (that is the part, which diverged from Gabi-Kat seq; 45 bp) [gives](https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Alignment?tl=d1Rz2jh3iI1XXicd-21960089-2406388817) 100% match to SYP61 (1:10016871 to 1:10016915). E = 3.5e-18.
AvadaMouse
CC BY-SA 4.0
2022-12-26T10:06:34.970Z
To conclude, the 45 bp stretch is **non-contiguous** to the 127 bp stretch in NC_003070.9 (GCA_000001735.2; TAIR v. 10.1).
AvadaMouse
CC BY-SA 4.0
2022-12-25T18:17:44.710Z
2022-12-26T14:05:28.910Z
2022-12-26T14:05:28.910Z
M__
<p>I am tasked with designing primers for a <a href="https://www.gabi-kat.de/db/showseq.php?line=165A12" rel="nofollow noreferrer">particular mutant</a> (Target Gene Locus: At1g28490 i.e. SYP61) of Arabidopsis, obtained from Gabi-Kat:</p>
<ul>
<li>The NCBI reference sequence for the genome <a href="https://www.ncbi.nlm.nih.gov/nuccore/NM_102617.4?report=fasta" rel="nofollow noreferrer">goes</a>:
GGCTTTATAAGATTCTCCAGTAGCCCAAAAATAAACTTTCAGCGAATCTCAAAGCCAACCGCAGGTGATGCGTTTTGCTATGTAAGAGAGACGAATCTCACGCCGGAATCTCTTTTACTTCACCGTCGTTTGTTCTCCGTCGCGTCACCTTTCGATTAATTCCTTCACAGGGTAGCAAAATTTTGTTGATTTCGGG<span class="math-container">$\color{red}{\textrm{[}}$</span>ATGTCTTCAGCTCAAGATCCATTCTACATTGTTAAAGAAGAGATTCAGGATTCTATTGATAAGTTGCAATCTACATTCCACAAATGGG<strong>AGCG</strong>T<strong>ATCTCTCCCGACATGG</strong>G<strong>TGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGT</strong>A͟G͟A͟T͟G͟A͟G͟C͟T͟G͟G͟A͟A͟A͟A͟A͟G͟C͟G͟A͟T͟T͟A͟C͟T͟G͟T͟T͟G͟C͟A͟G͟C͟A͟A͟A͟A͟G͟A͟T͟C͟C͟T͟T͟C͟ATGGTATGGCATTGATGAAGCTGAGCTTGAAAAACGGAGGAGATGGACTAGTAATGCTAGGACACAGGTGCGGAATGTGAAGTCTGGTGTTTTAGCTGGAAAGGTTAGTTCTGGAGCAGGTCATGCAAGTGAAGTGCGTAGAGAACTAATGAGAATGCCAAACTCGGGTGAAGCAAGTAGATATGATCAGTATGGAGGGAGAGACGATGATGGTTTTGTACAGTCAGAATCAGATAGGCAAATGCTGTTGATAAAGCAACAAGACGAAGAATTGGATGAGCTTAGTAAAAGTGTACAGAGAATTGGAGGAGTGGGGCTCACTATACACGACGAACTCGTTGCACAGGAGAGAATAATAGATGAATTGGATACGGAGATGGACAGCACAAAGAACCGTCTTGAATTTGTTCAGAAAAAAGTGGGAATGGTAATGAAGAAAGCTGGAGCAAAAGGTCAAATGATGATGATATGTTTCTTGCTCGTCTTGTTCATCATCCTATTCGTTCTTGTCTTCTTGACCTAA<span class="math-container">$\color{red}{\textrm{]}}$</span>AGAACAGAGCCGCGGGATCTTATAGATGGGAGTTTAATGTCTTACATTATCATGTTTGTGGTGAAACAAATGGAATATTTACATAAGGTTGTCTGTTTTTTCTTCTCCTATCCAATCCCACATATAATTTTGAACTTTGGCTGTGTGATCTCTAGAGAACCTGAAATTAAATACTTTCATAGCTTGATCCCATGGCTAAGTTCCT</li>
<li>Gabi-Kat <a href="https://www.gabi-kat.de/db/showseq.php?line=165A12" rel="nofollow noreferrer">provides</a> the following sequence from the mutant as a confirmatory match with the genome locus:
<strong>AGCG</strong>G<strong>ATCTCTCCCGACATGG</strong>C<strong>TGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGT</strong>C͟C͟T͟C͟C͟T͟A͟T͟T͟C͟A͟C͟C͟C͟A͟C͟G͟G͟T͟T͟C͟T͟A͟T͟T͟A͟A͟T͟C͟C͟C͟C͟T͟C͟T͟T͟T͟A͟T͟A͟T͟A͟T͟A͟A͟</li>
</ul>
<p>I have <strong>bold-faced</strong> the common region to both sequences. What explains the <em>complete mismatch</em> of the trailing part (u͟n͟d͟e͟r͟l͟i͟n͟e͟d͟) of the Gabi-Kat sequence with the ref genome?</p>
<p>It is <strong>not</strong> the case that the mutation happened there; we are in the coding region (demarcated by red square brackets in the NCBI sequence) while Gabi-Kat's mutation is in <a href="https://www.gabi-kat.de/faq.html#faq9" rel="nofollow noreferrer">an untranslated region</a>.</p>
<h2>Update 1:</h2>
<p><a href="https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Results?tl=5YJb4Lf7RtlEb7qc-21959324-2406285229" rel="nofollow noreferrer">Running</a> Ensembl Blast on the GabiKat seq returns a 92.9% identicality with SYP61. Their sequence for the gene (taken from <a href="https://www.ncbi.nlm.nih.gov/nuccore/NC_003070.9?report=fasta" rel="nofollow noreferrer">TAIR10.1</a>) goes:</p>
<p>AGCGGATCTCTCCCGACATGGCTGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGTA͟C͟T͟C͟G͟T͟A͟T͟T͟T͟A͟G͟T͟C͟A͟C͟G͟G͟T͟T͟C͟T͟A͟T͟T͟A͟A͟T͟C͟T͟C͟G͟T͟C͟T͟T͟T͟A͟T͟A͟T͟A͟T͟A͟A͟</p>
<p>Note the difference in the trailing part from NCBI ref; this almost coincides with the Gabi-Kat seq!</p>
<p>So, is the NCBI ref-genome and <a href="https://www.arabidopsis.org/servlets/TairObject?type=sequence&amp;id=2002978274" rel="nofollow noreferrer">the data for the particular TAIR genome</a> inaccurate? Maybe, TAIR updated their whole genome set but forgot to update data for individual genes?</p>
<h2>Update 2:</h2>
<p>My proposed explanation for the anomaly makes no sense. TAIR 7 (2004) has the <a href="https://www.ncbi.nlm.nih.gov/nuccore/NC_003070.5?report=fasta" rel="nofollow noreferrer">same sequence</a>. It looks like I am missing out on something more fundamental.</p>

<gene><genomics><blastn><blastp>
Mismatch between mutant gene and reference gene outside the site of mutation
# Mismatch between mutant gene and reference gene outside the site of mutation

(Asked by: username_0 on 2022-12-25)

I am tasked with designing primers for a [particular mutant](https://www.gabi-kat.de/db/showseq.php?line=165A12) (Target Gene Locus: At1g28490 i.e. SYP61) of Arabidopsis, obtained from Gabi-Kat:

* The NCBI reference sequence for the genome [goes](https://www.ncbi.nlm.nih.gov/nuccore/NM_102617.4?report=fasta):
  GGCTTTATAAGATTCTCCAGTAGCCCAAAAATAAACTTTCAGCGAATCTCAAAGCCAACCGCAGGTGATGCGTTTTGCTATGTAAGAGAGACGAATCTCACGCCGGAATCTCTTTTACTTCACCGTCGTTTGTTCTCCGTCGCGTCACCTTTCGATTAATTCCTTCACAGGGTAGCAAAATTTTGTTGATTTCGGG$\color{red}{\textrm{[}}$ATGTCTTCAGCTCAAGATCCATTCTACATTGTTAAAGAAGAGATTCAGGATTCTATTGATAAGTTGCAATCTACATTCCACAAATGGG**AGCG**T**ATCTCTCCCGACATGG**G**TGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGT**A͟G͟A͟T͟G͟A͟G͟C͟T͟G͟G͟A͟A͟A͟A͟A͟G͟C͟G͟A͟T͟T͟A͟C͟T͟G͟T͟T͟G͟C͟A͟G͟C͟A͟A͟A͟A͟G͟A͟T͟C͟C͟T͟T͟C͟ATGGTATGGCATTGATGAAGCTGAGCTTGAAAAACGGAGGAGATGGACTAGTAATGCTAGGACACAGGTGCGGAATGTGAAGTCTGGTGTTTTAGCTGGAAAGGTTAGTTCTGGAGCAGGTCATGCAAGTGAAGTGCGTAGAGAACTAATGAGAATGCCAAACTCGGGTGAAGCAAGTAGATATGATCAGTATGGAGGGAGAGACGATGATGGTTTTGTACAGTCAGAATCAGATAGGCAAATGCTGTTGATAAAGCAACAAGACGAAGAATTGGATGAGCTTAGTAAAAGTGTACAGAGAATTGGAGGAGTGGGGCTCACTATACACGACGAACTCGTTGCACAGGAGAGAATAATAGATGAATTGGATACGGAGATGGACAGCACAAAGAACCGTCTTGAATTTGTTCAGAAAAAAGTGGGAATGGTAATGAAGAAAGCTGGAGCAAAAGGTCAAATGATGATGATATGTTTCTTGCTCGTCTTGTTCATCATCCTATTCGTTCTTGTCTTCTTGACCTAA$\color{red}{\textrm{]}}$AGAACAGAGCCGCGGGATCTTATAGATGGGAGTTTAATGTCTTACATTATCATGTTTGTGGTGAAACAAATGGAATATTTACATAAGGTTGTCTGTTTTTTCTTCTCCTATCCAATCCCACATATAATTTTGAACTTTGGCTGTGTGATCTCTAGAGAACCTGAAATTAAATACTTTCATAGCTTGATCCCATGGCTAAGTTCCT
* Gabi-Kat [provides](https://www.gabi-kat.de/db/showseq.php?line=165A12) the following sequence from the mutant as a confirmatory match with the genome locus:
  **AGCG**G**ATCTCTCCCGACATGG**C**TGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGT**C͟C͟T͟C͟C͟T͟A͟T͟T͟C͟A͟C͟C͟C͟A͟C͟G͟G͟T͟T͟C͟T͟A͟T͟T͟A͟A͟T͟C͟C͟C͟C͟T͟C͟T͟T͟T͟A͟T͟A͟T͟A͟T͟A͟A͟

I have **bold-faced** the common region to both sequences. What explains the *complete mismatch* of the trailing part (u͟n͟d͟e͟r͟l͟i͟n͟e͟d͟) of the Gabi-Kat sequence with the ref genome?

It is **not** the case that the mutation happened there; we are in the coding region (demarcated by red square brackets in the NCBI sequence) while Gabi-Kat's mutation is in [an untranslated region](https://www.gabi-kat.de/faq.html#faq9).

## Update 1:

[Running](https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Results?tl=5YJb4Lf7RtlEb7qc-21959324-2406285229) Ensembl Blast on the GabiKat seq returns a 92.9% identicality with SYP61. Their sequence for the gene (taken from [TAIR10.1](https://www.ncbi.nlm.nih.gov/nuccore/NC_003070.9?report=fasta)) goes:

AGCGGATCTCTCCCGACATGGCTGATCAAGCACATGTGGCGAAGGAGCTTGTTGCTACTTGTGGAAGCATTGAGTGGCAGGTA͟C͟T͟C͟G͟T͟A͟T͟T͟T͟A͟G͟T͟C͟A͟C͟G͟G͟T͟T͟C͟T͟A͟T͟T͟A͟A͟T͟C͟T͟C͟G͟T͟C͟T͟T͟T͟A͟T͟A͟T͟A͟T͟A͟A͟

Note the difference in the trailing part from NCBI ref; this almost coincides with the Gabi-Kat seq!

So, is the NCBI ref-genome and [the data for the particular TAIR genome](https://www.arabidopsis.org/servlets/TairObject?type=sequence&id=2002978274) inaccurate? Maybe, TAIR updated their whole genome set but forgot to update data for individual genes?

## Update 2:

My proposed explanation for the anomaly makes no sense. TAIR 7 (2004) has the [same sequence](https://www.ncbi.nlm.nih.gov/nuccore/NC_003070.5?report=fasta). It looks like I am missing out on something more fundamental.

**Tags:** gene, genomics, blastn, blastp

---

## 1 answers

### Answer by username_1 on 2022-12-25. Score: 0

**Conclusion** The evidence suggests the region you are looking at **Gabi-kat** appears either 27045 downstream or 192250 upstream of syntaxin SNARE (SYP61), but is definitely not within the gene. It is chance homology (e.g. a pseudo-gene) that appears to occur within syntaxin (SYP61). The intergenic region could be much closer to the SYP61 gene if the same reference genome was used for both analyses Genbank LR782542.1 versus Genbank NC\_003070.9 (different reference genomes in my analysis). Hence, the University is probably correct.

What I suggest is the genome of *A. thaliana* is downloaded and Blast analysis is performed locally. This way differences in the reference genome are removed.

*The analysis is as follows*:

The target protein in question is involved in vesicle trafficking ([here](https://www.ebi.ac.uk/interpro/entry/InterPro/IPR000727/)) called SNARE, comprising coiled-coil domains and are highly conserved, for example show 100% conservation between *Arabidopsis thaliana*, *Arabidopsis suecica* and *Camelina sativa*

> MSSAQDPFYIVKEEIQDSIDKLQSTFHKWERISPDMGDQAHVAKELVATCGSIEWQVDELEKAITVAAKDP

If Gabi-kat was a protein its sequence would be ...

> RISPDMADQAHVAKELVATCGSIEWQVLLFTHGSINPLFIY

Here's a blast of the whole regions of the Gabi-kat. The region in green is the SNARE proteins ...

[![enter image description here](https://i.sstatic.net/uBRgA.png)](https://i.sstatic.net/uBRgA.png)

So the issue is what is the terminal 3' protein and Blast that, i.e. the final bit that matches no known plant and the locus you flagged as being anomalous ..

[![enter image description here](https://i.sstatic.net/6WqUm.png)](https://i.sstatic.net/6WqUm.png)

This is basically zero result via Blast, the protein with low homology are lower metazoans (animals) and a few bacteria ... a long way from plants. In other words this isn't a known protein of natural origin at least on NCBI.

The locus is literally like no protein on earth.

However ... Blasting the nucleotide ...

1. Reference *Arabidopsis thaliana* = syntaxin (SNARE; SYP61)
2. Gabi-Kat *Arabidopsis thaliana* top two Blast hits are ...

* Chromosome 1 at 10039082 - 10039208, or 9819787 to 981991 or
* Both value 1e-42
* Both dentity is 93%

So Gabi-Kat is a real part of the *Arabidopsis thaliana*, **Gabi-kat is a non-coding region**. *Thats the major conclusion*.

The question where is the location of SYP61? And that is

* Chromosome 1 at 10012037 to 10020311

Described here from [NCBI's full report](https://www.ncbi.nlm.nih.gov/gene?Db=gene&Cmd=DetailsSearch&Term=839749)

[![enter image description here](https://i.sstatic.net/u46cD.gif)](https://i.sstatic.net/u46cD.gif)

Thus the mutant is 27045bp downstream or 192250 bp upstream **IF** the same reference genome is used for the blasting (Genome LR782542.1) and the Genbank report NC\_003070.9. They could easily be much closer to each other is the same reference genome was used for both sets of analysis.

---

To answer the comments

* I [the OP] do not follow your [my] blast analysis
* My [the OP] blasts gives ... ['hits within the gene']
* etc ....

It appears to me there is a psuedo-gene phenomenon occurring which is making you think a mutant is located within a gene. However, if you disagree thats perfectly okay and the original version of this answer I had assumed the Gaba-kat team had got this wrong (I always from the basis of supporting the OP).

NCBI was the database repository in the remote analysis, but there's an issue with multiple references ... so to reiterate I suggest you perform this analysis as a local blast to remove multiple reference phenomena.

To conclude I personally don't think there is a basis to believe the University responsible for Gaba-kat has got this wrong <https://www.gabi-kat.de>. However, if you feel otherwise then I would recommend contacting them with your analysis and any supportive analysis for example from this thread or elsewhere. Please note, I don't have a connection with this university.

> **Comments:**
> **username_0**: I do not follow your BLAST of the Gabi-Kat Nucleotide Seq. [My blast](https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Results?tl=5YJb4Lf7RtlEb7qc-21959324-2406285229) gives a 92.9% alignment for the entire stretch (127 bp long) with SYP61 (1:10016601 to 1:10016721). E = 4.6e-45  
> **username_0**: Indeed, SYP61 is [located](https://bar.utoronto.ca/thalemine/portal.do?externalids=AT1G28490) from 1:10015856 to 10018047 in NC_003070.9 (GCA_000001735.2; TAIR v. 10.1). Which is inclusive of 1:10016601 to 1:10016721. You have got the specifics wrong.  
> **username_0**: Running a nucleotide blast on the underlined part of the NCBI ref-seq (that is the part, which diverged from Gabi-Kat seq; 45 bp) [gives](https://plants.ensembl.org/Arabidopsis_thaliana/Tools/Blast/Alignment?tl=d1Rz2jh3iI1XXicd-21960089-2406388817) 100% match to SYP61 (1:10016871 to 1:10016915). E = 3.5e-18.  
> **username_0**: To conclude, the 45 bp stretch is **non-contiguous** to the 127 bp stretch in NC_003070.9 (GCA_000001735.2; TAIR v. 10.1).  

---

