is there a list of profiles or pairwise distances in profiles ?
> This combination leaves no doubt that the 3 BP deletion was acquired via recombination.
I think most people have lots of doubt here.
all H5N1 in Egypt is Qinghai, how should Shanxi-polymorphisms come to Egypt
I'm currently not so happy with the grouping algorithm
here is another grouping, which seems to be better for small files:
Code:
00000000000000000000000000000000000000000000000000000011111111111111111111111111111111
00000112222223344444555555555666666666777777888888999900000011222233333444444445555666
23499371224462905669000145678001566789113567246678338901115639355848899003456792679015
77503856250911979892018292838128136679162184292795362964680577323716958781054741359548
-codon-position------------------------------ 1 1 2 1 2 12 1 12 1 2 1 1 1 1 21 12 1 1 2
---Index-------------------------------------ACAGGAGTTACTCCCGCATATCGACCGGTGCCGACGATGGCTCGGTGCGACTGCAAGGATGCTAAGCGCCCACGTGCGAGAAGTAT
1 >A/Ck/Egypt/06553-NLQP,2006/03/21 ...................................................................................... 0 0 , 0 0 1:>A/Ck/Egypt/06553-NLQP,2006/03/21
2 >A/Ck/Egypt/092-NLQP,2009/01 ......A.C..C.T.A....C......A...T..............A..G...............................G.... 11 11 , 11 11 2:>A/Ck/Egypt/092-NLQP,2009/01
3 >A/Gs/Egypt/0910sm-NLQP,2009/01 ......A.C..C.T.A....C......A...T..............A..G...............................G.... 12 11 , 12 11 3:>A/Gs/Egypt/0910sm-NLQP,2009/01
4 >A/Dk/Egypt/0821S-NLQP,2008/08/25 ......A.C..C.T.A....C......A...T..............A..G...............................G.... 17 11 , 12 11 4:>A/Dk/Egypt/0821S-NLQP,2008/08/25
5 >A/Dk/Egypt/09332sm-NLQP,2009/05 ....A.A.CC.C.T.A....C......A...T.................G.........................A.....G...- 23 14 , 15 13 5:>A/Dk/Egypt/09332sm-NLQP,2009/05
6 >A/Ck/Egypt/0945-NLQP,2009/03 ......A.C..C.T.A....C......A...T.................G...........T......T.........G..G.... 13 13 , 13 13 6:>A/Ck/Egypt/0945-NLQP,2009/03
7 >A/Egypt/N03228,2009/03/30 ..G...A.C..C.T.A.G..C...A..A...T......AA.......T.G.C................T..........T.G.... 19 19 , 19 19 7:>A/Egypt/N03228,2009/03/30
8 >A/Egypt/N03272,2009/04/01 ..G...A.C..C.T.A.G..C...A..A...T......A........T.G.C................T..........T.G.... 19 18 , 19 18 8:>A/Egypt/N03272,2009/04/01
9 >A/Egypt/N02407,2009/03/09 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 16 13 , 16 13 9:>A/Egypt/N02407,2009/03/09
10 >A/Egypt/N04395,2009/05/13 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 13 13 , 13 13 10:>A/Egypt/N04395,2009/05/13
11 >A/Egypt/N04396,2009/05/17 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 16 13 , 16 13 11:>A/Egypt/N04396,2009/05/17
12 >A/Egypt/N03439,2009/04/18 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 13 13 , 13 13 12:>A/Egypt/N03439,2009/04/18
13 >A/Gs/Egypt/09324sm-NLQP,2009/04 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 14 13 , 14 13 13:>A/Gs/Egypt/09324sm-NLQP,2009/04
14 >A/Egypt/N02039,2009/03/03 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 14 13 , 14 13 14:>A/Egypt/N02039,2009/03/03
15 >A/Dk/Egypt/0918-NLQP,2009/02 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 13 13 , 13 13 15:>A/Dk/Egypt/0918-NLQP,2009/02
16 >A/Egypt/N01310,2009/02/28 ..G...A.C..C.T.A....C......A...T......A..........G..................T............G.... 14 13 , 14 13 16:>A/Egypt/N01310,2009/02/28
17 >A/Ck/Egypt/0987alb-NLQP,2009/05/12 G.G...A.C..C.T.A....C......A...T......A..........G....G.............T............G.... 18 15 , 15 15 17:>A/Ck/Egypt/0987alb-NLQP,2009/05/12
18 >A/Ck/Egypt/0987-NLQP,2009/05/12 ..G...A.C..C.T.A....C......A...T......A..........G....G.............T............G.... 16 14 , 16 14 18:>A/Ck/Egypt/0987-NLQP,2009/05/12
19 >A/Dk/Egypt/09274sm-NLQP,2009/04 ..G...A.C..C.T.A....C......A...T......A..........G..................T....A.......G.... 23 14 , 16 14 19:>A/Dk/Egypt/09274sm-NLQP,2009/04
20 >A/Ck/Egypt/09485s-NLQP,2009/05 ..G...A.C..C.T.A....C.A....A...T......A..........G..................T............G...- 22 15 , 14 14 20:>A/Ck/Egypt/09485s-NLQP,2009/05
21 >A/Egypt/N00605,2009/02/03 ......A.C..C.T.A....C.A....A...T.............C...G..................T............G.... 14 13 , 14 13 21:>A/Egypt/N00605,2009/02/03
22 >A/Egypt/N03438,2009/04/16 ......A.C..C.T.A....C......A...T.................G..................T............G.... 15 11 , 15 11 22:>A/Egypt/N03438,2009/04/16
23 >A/Ck/Egypt/0894-NLQP,2008/11/27 ......A.C..C.T.A....C......A...T.................G..................T............G.... 11 11 , 11 11 23:>A/Ck/Egypt/0894-NLQP,2008/11/27
24 >A/Dk/Egypt/0972-NLQP,2009/04 ......A.C..C.T.A....C......A...T.................G..................T............G.... 13 11 , 13 11 24:>A/Dk/Egypt/0972-NLQP,2009/04
25 >A/Dk/Egypt/093-NLQP,2009/01 ......A.C..CTT.A..C.C......A...T.................G...............................G...- 20 13 , 12 12 25:>A/Dk/Egypt/093-NLQP,2009/01
26 >A/Egypt/N00001,2009/01/11 ......A.C..CTT.A..C.C......A...T.................G...............................G.... 14 12 , 14 12 26:>A/Egypt/N00001,2009/01/11
27 >A/Gs/Egypt/0929-NLQP,2009/02 ...A..A.C..C.T......C......A...T.................G...............................G.... 11 10 , 11 10 27:>A/Gs/Egypt/0929-NLQP,2009/02
28 >A/Dk/Egypt/0967smf-NLQP,2009/02 ......A.C..C.T.A....C......A...T.................G...............................G.... 16 10 , 13 10 28:>A/Dk/Egypt/0967smf-NLQP,2009/02
29 >A/Gs/Egypt/09102smf-NLQP,2009/03 -.....A.C..C.T.A....C......A...T................AG...............................G.... 21 12 , 16 11 29:>A/Gs/Egypt/09102smf-NLQP,2009/03
30 >A/Ck/Egypt/0896-NLQP,2008/12/24 ......ACC..C.T.A....C.....AA...T................AG........G......................G.... 15 14 , 15 14 30:>A/Ck/Egypt/0896-NLQP,2008/12/24
31 >A/Dk/Egypt/09315s-NLQP,2009/04 --....ACC..C.T.AT...C.....AA...T.................G...T...............T...........G.... 25 17 , 19 15 31:>A/Dk/Egypt/09315s-NLQP,2009/04
32 >A/Dk/Egypt/0926-NLQP,2009/02 ......A....C.T.A....C......A...T...........A.....G...............................G...C 14 11 , 14 11 32:>A/Dk/Egypt/0926-NLQP,2009/02
33 >A/Egypt/2629-NAMRU3,2007/ ---...A......T.A....C......A...T.................G...............................G.... 17 11 , 8 8 33:>A/Egypt/2629-NAMRU3,2007/
34 >A/Egypt/2621-NAMRU3,2007/ ......A......T.A....C......A...T.................G...............................G.... 8 8 , 8 8 34:>A/Egypt/2621-NAMRU3,2007/
35 >A/Egypt/0636-NAMRU3,2007/01 ......A......T.A....C......A.....................G......................T........G.... 10 8 , 10 8 35:>A/Egypt/0636-NAMRU3,2007/01
36 >A/Ck/Egypt/9400NAMRU3-CLEVB211,2007/ ....A.A......T.A....C.....................T......G...............................G.... 15 8 , 15 8 36:>A/Ck/Egypt/9400NAMRU3-CLEVB211,2007/
37 >A/Egypt/1394-NAMRU3,2007/ ...A..A......T.A....C..................R.........G...............................G...C 12 9 , 12 9 37:>A/Egypt/1394-NAMRU3,2007/
38 >A/Egypt/2289-NAMRU3,2008/03/02 ---...A...T..T.A....C..................T.........G.....G...........A........TA...G.... 22 15 , 13 12 38:>A/Egypt/2289-NAMRU3,2008/03/02
39 >A/Egypt/2546-NAMRU3,2008/03/08 ---...A...T..T.A....C..................A.........G.....G...........A........T...GG.... 24 15 , 15 12 39:>A/Egypt/2546-NAMRU3,2008/03/08
40 >A/Ck/Egypt/0847-NLQP,2008/01/18 ......A......T.A....C............................G............................G..G.... 15 7 , 15 7 40:>A/Ck/Egypt/0847-NLQP,2008/01/18
41 >A/Egypt/2631-NAMRU3,2007/ ---...A.....TT.A....C............................G....................T..........G.... 18 11 , 9 8 41:>A/Egypt/2631-NAMRU3,2007/
42 >A/Ck/Egypt/R2,2007/01/01 ......A......T.A....C....T..................A....G....................T..........G.... 18 9 , 18 9 42:>A/Ck/Egypt/R2,2007/01/01
43 >A/Dk/Egypt/R5,2007/02/20 ......A......T.A....C............G...............G....................T..........G.... 11 8 , 11 8 43:>A/Dk/Egypt/R5,2007/02/20
44 >A/Gs/Egypt/R4,2007/02/01 ......A......T.A...GC...................T........GT......A....A.......T.........GG.... 15 13 , 15 13 44:>A/Gs/Egypt/R4,2007/02/01
45 >A/Ck/Egypt/R3,2007/02/01 ......A......T.A...GC............G......T........GT...........A..................G.... 12 11 , 12 11 45:>A/Ck/Egypt/R3,2007/02/01
46 >A/Ck/Egypt/R6,2007/02/26 .G....A......T.A....C..........T.GT.........A....GT......A............................ 24 12 , 24 12 46:>A/Ck/Egypt/R6,2007/02/26
47 >A/Dk/Egypt/09349S-NLQP,2009/04 -----.A....C.T.A....C......A...TA...T....C.......G..A............A...........A...G.... 39 20 , 18 15 47:>A/Dk/Egypt/09349S-NLQP,2009/04
48 >A/Egypt/N04394,2009/05/11 ......A....C.T.A....C......A...TA...T....C.A.....G..A............A...............G.... 19 15 , 19 15 48:>A/Egypt/N04394,2009/05/11
49 >A/Egypt/N02752,2009/03/24 ......A....C.T.A....C......A...T...A.....CT......G...A.....CT..........GT........G.... 19 17 , 19 17 49:>A/Egypt/N02752,2009/03/24
50 >A/Egypt/N00606,2009/02/07 ......A......T.A.G..C......A...T...A.............G...A.....CT..........GT........GT... 18 16 , 18 16 50:>A/Egypt/N00606,2009/02/07
51 >A/Dk/Egypt/093L-NLQP,2009/02 ......A......T.AT...C....T.A...T...A.............G...A.....CT..........GTA.......GT.G. 20 19 , 20 19 51:>A/Dk/Egypt/093L-NLQP,2009/02
52 >A/Ck/Egypt/092L-NLQP,2009/02 ......A..G...T.A....C....T.A...T...A.............G...A.....CT..........GTA.......GT.G. 19 19 , 19 19 52:>A/Ck/Egypt/092L-NLQP,2009/02
53 >A/Dk/Egypt/0912-NLQP,2009/02 ......A..G...T.A....C....T.A...T...A.............G...A.....CT..........GTA.......GT.G. 19 19 , 19 19 53:>A/Dk/Egypt/0912-NLQP,2009/02
54 >A/Dk/Egypt/0930smL-NLQP,2009/02 ......A......T.A....C......A...T...A.............G...A.....CT........T.GT........G.... 20 15 , 20 15 54:>A/Dk/Egypt/0930smL-NLQP,2009/02
55 >A/Egypt/N05056,2009/06/06 .....GA......T.A....C......A...T...A.C...........G...A.....CT..........GT........G.... 22 16 , 22 16 55:>A/Egypt/N05056,2009/06/06
56 >A/Dk/Egypt/0970-NLQP,2009/04 .....GA......T.A....C.A....A...T...A.C.......C...G...A.....CT..........GT........G.... 20 18 , 20 18 56:>A/Dk/Egypt/0970-NLQP,2009/04
57 >A/Egypt/N03450,2009/04/19 ......A......T.A....CAA....A.A.T...A.............G.C.A..A..CT..........GT........G.... 21 19 , 21 19 57:>A/Egypt/N03450,2009/04/19
58 >A/Egypt/N04527,2009/05/18 ......A......T.A....CA.....A.A.T...A.............G.C.A..A..CT..........GT........G.... 19 18 , 19 18 58:>A/Egypt/N04527,2009/05/18
59 >A/Dk/Egypt/0969smL-NLQP,2009/05 ......A......T.A....CA.....ACA.T...A.............G...A..A..CTT.........GT........G.... 19 19 , 19 19 59:>A/Dk/Egypt/0969smL-NLQP,2009/05
60 >A/equine/Egypt/av1,2009/03/24 ......A.C..C.T.A....C......ACA.T.GT..............G................T..............G.... 21 15 , 21 15 60:>A/equine/Egypt/av1,2009/03/24
61 >A/Dk/Egypt/0945smf-NLQP,2009/02 -.....A.C..C.T.A....C......A.A.T.................G................T..............G...- 25 14 , 13 12 61:>A/Dk/Egypt/0945smf-NLQP,2009/02
62 >A/Ck/Egypt/0910-NLQP,2009/01 G.....A....C.T.A....C......A...T.................GT..............................G.C.. 14 12 , 14 12 62:>A/Ck/Egypt/0910-NLQP,2009/01
63 >A/Ck/Egypt/0915-NLQP,2009/02 G.....A....C.T.A....C......A...T.G...............GT.......G................A.....G.C.. 17 15 , 17 15 63:>A/Ck/Egypt/0915-NLQP,2009/02
64 >A/Gs/Egypt/0912smg-NLQP,2009/02 G.....A....C.T.A....C......A...T.................GT........................A.....G.C.. 17 13 , 17 13 64:>A/Gs/Egypt/0912smg-NLQP,2009/02
65 >A/Ck/Egypt/093smg-NLQP,2009/01 G.....A....C.T.A....C......A...T.................GT.................T......A.....G.C.. 17 14 , 17 14 65:>A/Ck/Egypt/093smg-NLQP,2009/01
66 >A/Egypt/N00585,2009/01/23 .A....A....C.TTA....C.A....A...T.................G...............................G.... 14 12 , 14 12 66:>A/Egypt/N00585,2009/01/23
67 >A/Ck/Egypt/0920-NLQP,2009/02 .A....A....C.TTA....C.A....A...T.................G...............................G.... 15 12 , 15 12 67:>A/Ck/Egypt/0920-NLQP,2009/02
68 >A/Dk/Egypt/09339S-NLQP,2009/04 --....A.C..C.T.A....C.A....A..AT.................G.C...........GC.........C......G.... 23 18 , 17 16 68:>A/Dk/Egypt/09339S-NLQP,2009/04
69 >A/Egypt/N04316,2009/05/09 ......A.C..C.T.A....C.A....A..AT.................G.............GC.........C......G.... 16 15 , 16 15 69:>A/Egypt/N04316,2009/05/09
70 >A/Egypt/N02563,2009/03/14 ......A.C..C.T.A....C.A....A..AT.................G...........T.GC.........C......G.... 17 16 , 17 16 70:>A/Egypt/N02563,2009/03/14
71 >A/Egypt/N04823,2009/05/25 ......A.C..C.T.A....C.AG...A...T..............A..G...............................G.... 15 13 , 15 13 71:>A/Egypt/N04823,2009/05/25
72 >A/Dk/Egypt/0990SM-NLQP,2009/02 ......A.C..C.T.A....C.AG...A...T..............A..G...............................G.... 14 13 , 14 13 72:>A/Dk/Egypt/0990SM-NLQP,2009/02
73 >A/Dk/Egypt/09238sm-NLQP,2009/03 ......A.C..C.T.A....C....T.A...T..............A..G...............................G.... 15 12 , 15 12 73:>A/Dk/Egypt/09238sm-NLQP,2009/03
---Index-------------------------------------ACAGGAGTTACTCCCGCATATCGACCGGTGCCGACGATGGCTCGGTGCGACTGCAAGGATGCTAAGCGCCCACGTGCGAGAAGTAT
-codon-position------------------------------ 1 1 2 1 2 12 1 12 1 2 1 1 1 1 21 12 1 1 2
00000000000000000000000000000000000000000000000000000011111111111111111111111111111111
00000112222223344444555555555666666666777777888888999900000011222233333444444445555666
23499371224462905669000145678001566789113567246678338901115639355848899003456792679015
77503856250911979892018292838128136679162184292795362964680577323716958781054741359548
Since this thread is on the H5N1 donkey sequence, it is worth going through that sequence to demonstrate additional examples of recombination.
The above rearranged table now has the donkey sequence on row 60. As noted yesterday, the doney sequence has a “triple” that is also found in the duck sequence which is now on line 59. The “triple” is really a “quadruple”, because 4 polymorphisms in a row match in rows 59 and 60, even though the two sequences are quite distinct. The quadruple is G573A, T588C, G601A, and C618T. The outside markers, G573A and C618T are fairly widespread, but T588C and G601A are quite rare.
T588C is only found in the two sequences on lines 59 and 60. Not only isn’t it in other H5N1 sequences with the 3BP deletion, but it is not in any other H5N1 in Egypt, or even any clade 2.2 sequence at Genbank. It is in only one other sequence at Genbank, a 2004 H5N1 isolate from Guangdong, China. Because it is on two different backgrounds in Egypt, it arose by independent events. It either jumped from one sequence to the other via recombination, or is a coincidental copy error. The fact that it is so rare, yet is found in the same 3 BP deletion sub-clade in isolates collected on March 24 and May in 2009 at relatively close locations (Beni Suef and Sohag), makes the coincidence argument EXTREMELY unlikely.
However, when combined with the adjacent change , G601A, the coincidental random error explanation moves into the CDC/WHO/GSGS "hopes and dreams" category because G601A is also exceedingly rare and is only found in the 5 isolates on lines 57-62, but clearly falls into two distinct groups. One group includes the duck sequence on line 59 as well as the two human sequences on lines 57 and 58. The two human sequences where both from young children (4M) who lived in Sohag. The collections were on Apr 19, 2009 and May 18, 2009. The duck isolate was also from Sohag and also collected in May, 2009.
The other group includes the donkey sequence from Beni Suef, as well as a duck sequence collected in Feb in Fayoum, which is adjacent to Beni Suef,
Thus, once again a VERY rare polymorphism is found in 5 sequences which represent two distinct backgrounds collected over the same time frame (Feb-May 2009) in the same general area (Fayoum, Beni Suef, Sohag), which are do to two independent events, recombination between co-circulating H5N1’s which include the two sequences positive for T588C, which are also positive for G610A, or random error coincidences.
The matching of these two adjacent polymorphisms in the same two sub-clades is well into the “too many coincidences” category to support the “random mutation” fantasy.
Moreover, a bit further down is a “triple” (C618T, C651G, A663T) which is in the donkey sequence on line 60, as well as the chicken sequence on line 46, which is another background, but also from Beni Suef, providing more support for recombination and less support for more random coincidences.