This is a test version of Biostars. For the public version, visit https://www.biostars.org.
100% Identical Transcript Sequences - How Did They Manage To Put Them Into Different Loci?

I stumbled upon two genes:

After more closer look I noticed that they are similar. So similar that I BLAST'ed loci against each other.

enter image description here

Surprise, surprise, they are 100% identical.
Those genes can be grouped as paralogues, but there is still some unclarity for me.

Is such similarity normal? If transcripts are 100% identical, how somebody managed to place them into different loci?

Edit

Three more transcipts from different loci are identical: ENST00000456123.1, ENST00000420149.1, ENST00000436568.1

ensembl
• 4,203 views • link
updated by heboping1215 • heboping1215 • written by PoGibas PoGibas

are their flanking sequences different (utr's?)?

• 0 views • link
written by Pavel Senin Pavel Senin

After extending loci 1000bp to each side I still get 100%: 

blastn \
    -subject <(echo chr15 83087381 83102960 | 
                  awk '{OFS="\t"; print $1,$2-1000,$3+1000}' | 
                  bedtools getfasta -fi hg.19 -bed - -fo -) \
    -query <(echo chr15 82710856 82726435 | 
                awk '{OFS="\t"; print $1,$2-1000,$3+1000}' | 
                bedtools getfasta -fi hg.19 -bed - -fo -) |   
    grep -m1 'Identities'

Identities = 17579/17579 (100%), Gaps = 0/17579 (0%)
• 0 views • link
written by PoGibas PoGibas

interesting indeed. they might have had some long inserts etc., but I agree, that could be a mistake

• 0 views • link
written by Pavel Senin Pavel Senin

Could you tell me how to `BLAST'ed loci against each other`??

• 0 views • link
written by heboping1215 • heboping1215

1 answer

I've asked our friends at Havana to weigh in (either directly or via me) so we'll see what they say.

Follow up: heard back from Havana. Laurens says:

The duplication of this and neighbouring genes is due to a haplotype issue. There is an assembly gap between the two copies and the genomic sequence on the left side is from the RP13 haplotype, and the right side is from RP11 haplotype.

In the next genome assembly (GRCh38 I presume) this region will be replaced with an ungapped single haplotype (RP11) and it will only have one copy of this gene (and of the flanking genes).

• 265 views • link
written by Emily Emily

Emily, thank you for helping!
Can you also look at this question: Gencode V15 Exons < 3Bp (e.g., ENST00000605962 - last exon is 1bp long (same Havana annotation)).

• 0 views • link
written by PoGibas PoGibas

On it (plus some extra characters so BioStar will let me post).

• 0 views • link
written by Emily Emily

I'm replying so that it flags me when there are new answers. I'd like to hear the details.

• 0 views • link
written by Mary Mary

Oh dear. I edited my original post so you might not get a notification. Well this comment will do that.

• 125 views • link
written by Emily Emily

Log in to answer this question.

More posts like this

No similar posts found.