What model works for calculating genetic distances among sequences containing a minisatellite
0
0
Entering edit mode
8.6 years ago
Yongjie Zhang ▴ 110

Hello everyone,

I want to calculate genetic distances between intron allele sequences, but within these sequences is there a minisatellite, where a 19-bp motif may duplicate 2-6 times. I previously used K2P (Kimura-2-parameter) distances for Mantel tests and phylogeny constructions, but a reviewer commented that "since there are many gaps (i.e. indels) in the alignments, K2P can't effectively reflect such variations. In addition, repeat number variation can be faster than SNPs and using K2P or other measures without taking into account that repeat number variation is a single event and not multiple indels could generate very misleading distances and phylogenies. Thus, I believe the authors should modify their evolutionary model of intron sequence variation and reanalyze the data".

Do you know what model I should use to calculate genetic distances? My alignment is attached.

Any comments are welcome.

Yongjie Zhang

>A
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGGGTGGACCCCCGAA--------------------------------------------------------------------------------GCTTGTCCGGGCCACCACTGACGAGACTGGCGCGTTAG
>B
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGGGTGGACCCCCGAA--------------------------------------------------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>C
GTTCGTGCCGTGTGGATCCCCGAAGTGCCGGGTGGACCCCCGAA--------------------------------------------------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>D
GTTCGTGCCGTGTGGATCCCCGAAGTGCCGGGTGGACCCCCGAAGCGCCGTGTGGA-CCCCGAA------------------------------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>E
GTTCGTGCCGGGTGGATCCCCGAAGTGCCGGGTGGATCCCCGAAGCGCCGGCCGGACCCCCGAA------------------------------------------------------------GCTTGCCCGGGCCGCCACTGACGAGACTGGCGTGTTAG
>F
GTTCGTGCCGTGTGGA-CCCCGAAGTTCCGGGTGGACCCCCGAAGTGCCGGGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGAA----------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>G
GTTCGTGCCGTGTGGATCCCCGAAGTGCCGGGTGGACCCCCGAAGCGCCGTGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGAA----------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>H
GTTCGTGCCGTGTGGATCCCCAAAGTGCCGGGTGGACCCCGGAAGTGCCGTGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGAA----------------------------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>I
GTTCGTGCCGTGTGGACCCCCGAAGTGCCGGGTGGATCCCCGAAGTGCCGGGTGGATCCCCGAAGCGCCGGCCGGACCCCCGAA----------------------------------------GCTTGCCCGGGCCGCCACTGACGAGACTGGCGTGTTAG
>J
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGTGTGGA-CCCCGAAGTTCCGGGTGGACCCCCGAAGTGCCGGGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGCA--------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>K
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGTGTGGA-CCCCGAAGTTCCGGGTGGACCCCCGAAGTGCCGGGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGAA--------------------GCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
>L
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGGGTGGACCCCCGAAGTTCCGGGTGGATCCCCAAAGTGCCGGGTGGACCCCCGAAGCGCCGGGTGGACCCCCGAA--------------------GCTTGCCCGGGCCGCCACTGACGAACCTGGCGTGTTAG
>M
GTTCGTGCCGTGTGGA-CCCCGAAGTGCCGTGTGGA-CCCCGAAGTGCCGTGTGGA-CCCCGAAGTTCCGTGTGGACCCCCGAAGTGCCGGGTGGA-CCCCGAAGTGCCGTGTGGACCCCCGAAGCTTGTCCGGGCCGCCACTGACGAGACTGGCGCGTTAG
genetic-distance alignment minisatellite • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 3106 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6