maclandrol / profileNJ

Gene tree correction using species tree and NJ

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

AttributeError: 'NoneType' object has no attribute 'get_child_at'

lychen83 opened this issue · comments

Hello Emmanuel,

I am using profileNJ to fix gene tree.
My gene tree is as following:
(red_Nyctaginaceae_Anulocaulis_annulatus@27853:1.00000050003e-06,((red_Nyctaginaceae_Boerhavia_ciliata@13258:0.0750229228896,red_Nyctaginaceae_Boerhavia_ciliata@47639:0.393322292936)81:0.0503980792703,((red_Nyctaginaceae_Commicarpus_scandens@37240:0.0853400496541,(red_Nyctaginaceae_Commicarpus_scandens@37241:0.0693061575783,red_Nyctaginaceae_Commicarpus_scandens@37242:1.00000050003e-06)52:0.00188295368716)99:0.118423458122,(red_Nyctaginaceae_Anulocaulis_annulatus@27850:1.00000050003e-06,((((((red_Sarcobataceae_Sarcobatus_vermiculatus@10866:0.122532346073,red_Phytolaccaceae_Rivina_humilis@6133:0.0771373097868)77:0.0457652835781,red_Phytolaccaceae_Phytolacca_americana@19681:0.0957523478508)24:0.0232077280165,((((red_Gisekiaceae_Gisekia_pharnaceoides@65620:0.0162456785951,(red_Gisekiaceae_Gisekia_pharnaceoides@65621:0.00681970420165,red_Gisekiaceae_Gisekia_pharnaceoides@65622:1.00000050003e-06)93:0.0137781609124)95:0.132874653756,(red_Aizoaceae_Mesembryanthemum_crystallinum@8529:0.223075138192,(blue_Kewaceae_Kewa_bowkeriana@31410:0.228851579754,red_Aizoaceae_Delosperma_echinatum@7254:0.141714511709)29:0.0396085895611)17:0.0154021255408)20:0.0367093001471,((((red_Anacampserotaceae_Anacampseros_filamentosa@14002:0.0826118742079,red_Anacampserotaceae_Anacampseros_filamentosa@14001:0.003197010308)100:0.109493084342,(red_Talinaceae_Talinum_fruticosum@66847:0.0676274464062,red_Talinaceae_Talinum_fruticosum@66849:0.0139847411967)74:0.0511618653015)70:0.091045419994,((red_Didieraceae_Decarya_madagascariensis@23199:0.0852317909327,((red_Amaranthaceae_Alternanthera_philoxeroides@17136:0.454697936118,red_Basellaceae_Basella_alba@23542:0.0199662285352)23:0.046662612322,red_Anacampserotaceae_Anacampseros_filamentosa@49652:0.595301597061)4:0.0337060396382)4:0.0429593180478,((red_Anacampserotaceae_Anacampseros_filamentosa@14000:0.0841031764793,blue_Portulacaceae_Portulaca_oleracea_SRA@46368:0.109269449235)26:0.0268324229235,(red_Cactaceae_Opuntia_arenaria@39809:0.0315454355484,red_Cactaceae_Opuntia_arenaria@39810:0.00343597093487)100:0.09104635865)39:0.0515253411516)2:0.037503228393)5:0.0635440257922,(unknown_Microteaceae_Microtea_debilis@17540:0.145860437869,red_Montiaceae_Claytonia_nevadensis@56211:0.187699045965)20:0.0705402302503)4:0.0742239111562)0:0.0184869798936,(((blue_Molluginaceae_Glinus_lotoides_var_virens@36790:0.146012741869,(((blue_Molluginaceae_Glinus_lotoides_var_virens@14738:0.00983571271816,blue_Molluginaceae_Glinus_lotoides_var_virens@14739:0.0269370027329)99:0.195303617112,blue_Molluginaceae_Mollugo_verticillata@40080:0.150193781851)26:0.0320365128929,blue_Molluginaceae_Mollugo_pentaphylla@44374:0.218825470856)24:0.0361546632693)19:0.022154212428,(((blue_Macarthuriaceae_Macarthuria_australis@11352:0.00512290112574,blue_Macarthuriaceae_Macarthuria_australis@11351:0.0125256548542)96:0.142732271888,(((unknow_Achatocarpaceae_Achatocarpus_gracilis@77644:0.0105323245142,unknow_Achatocarpaceae_Phaulothamnus_spinescens@17223:0.0137385511835)99:0.174326161,((blue_Limeaceae_Limeum_aethiopicum@23901:0.00347291975739,(blue_Limeaceae_Limeum_aethiopicum@23900:0.020879247556,blue_Limeaceae_Limeum_aethiopicum@23899:0.0204887166203)69:0.00440202539874)100:0.257186692901,(blue_Molluginaceae_Mollugo_nudicaulis@169477:0.0411641859481,(blue_Molluginaceae_Pharnaceum_exiguum@2865:0.111007925512,(blue_Molluginaceae_Suessenguthiella_caespitosa@36653:0.0710636432298,blue_Molluginaceae_Mollugo_cerviana@40369:0.0791069642423)44:0.0211564098103)86:0.0655998388691)51:0.0241342277802)33:0.018360023148)12:0.0257142317318,((((blue_Caryophyllaceae_Paronychia_jamesii@203596:1.00000050003e-06,blue_Caryophyllaceae_Paronychia_jamesii@59891:1.00000050003e-06)100:0.116801488685,blue_Caryophyllaceae_Paronychia_jamesii@87181:0.0367037061832)100:0.154892753682,((blue_Caryophyllaceae_Polycarpon_tetraphyllum@79315:0.00943712161928,blue_Caryophyllaceae_Polycarpon_tetraphyllum@79314:1.00000050003e-06)97:0.135688581605,(((blue_Caryophyllaceae_Arenaria_serpyllifolia@39076:0.140821562494,blue_Caryophyllaceae_Dianthus_caryophyllus@29848_1:0.235153719659)28:0.0207193309636,(blue_Caryophyllaceae_Spergularia_marina@25124:0.00857472979863,(blue_Caryophyllaceae_Spergularia_marina@25123:1.00000050003e-06,blue_Caryophyllaceae_Spergularia_marina@25125:0.032413342875)96:0.0360590845676)100:0.171391427605)16:0.00742687466695,(blue_Caryophyllaceae_Silene_latifoliasra@61371:0.608803217036,((blue_Caryophyllaceae_Cerastium_arvense@29092:0.176471004778,(blue_Caryophyllaceae_Eremogone_hookeri_ssp_desertorum@24385:1.00000050003e-06,blue_Caryophyllaceae_Eremogone_hookeri_ssp_desertorum@24384:0.028997709389)100:0.102532078937)63:0.0615824913541,(blue_Caryophyllaceae_Silene_latifoliasra@33251:0.0860763478519,blue_Caryophyllaceae_Silene_latifoliasra@33249:0.0473470645264)100:0.142510307827)28:0.0086055003488)23:0.0298212089246)57:0.0785513976465)36:0.0397157106442)46:0.0858246940584,(((((unknow_Achatocarpaceae_Achatocarpus_gracilis@15738:0.0751052162584,(unknow_Achatocarpaceae_Achatocarpus_gracilis@10733:1.00000050003e-06,unknow_Achatocarpaceae_Achatocarpus_gracilis@10732:1.00000050003e-06)97:0.137939787409)100:2.52918046416,(blue_Caryophyllaceae_Paronychia_jamesii@124808:1.28939472946,unknow_Achatocarpaceae_Achatocarpus_gracilis@54893:1.13110559451)48:0.263587566679)98:0.656578144061,blue_Caryophyllaceae_Paronychia_jamesii@104139:0.0146226875442)13:0.126328798414,((blue_Caryophyllaceae_Paronychia_jamesii@104138:1.00000050003e-06,blue_Caryophyllaceae_Paronychia_jamesii@38862:0.184517130543)79:0.0632214730562,blue_Caryophyllaceae_Paronychia_jamesii@45552:0.192720512027)37:0.0562810045063)15:0.16001611266,((red_Chenopodiaceae_Halogeton_glomeratus@13825:0.119182644484,red_Chenopodiaceae_Haloxylon_ammodendron@77005:0.0505640087097)95:0.0679183037396,((red_Amaranthaceae_Alternanthera_philoxeroides@28652:0.22807223126,red_Amaranthaceae_Beta_vulgaris@Bv4_077820_hmhd_t1:0.100321540569)60:0.0668036619595,red_Chenopodiaceae_Suaeda_glauca@37135:0.213103441475)19:0.0136174880121)66:0.0782098971661)5:0.025204635845)0:0.0372880400021)2:0.0296039839173)1:0.0616609262282,(((((blue_Frankeniaceae_Frankenia_laevis@4934:1.00000050003e-06,blue_Frankeniaceae_Frankenia_laevis@4935:0.0108193420069)100:0.113372259493,blue_Tamaricaceae_Tamarix_hispida@8044:0.181599665814)53:0.0660229906552,(((blue_Ancistrocladaceae_Ancistrocladus_robertsoniorum@15866:0.137899700795,(blue_Nepenthaceae_Nepenthes_alata@66492:0.219953678896,(blue_Drosophyllaceae_Drosophyllum_lusitanicum@6817:0.0323063765371,blue_Drosophyllaceae_Drosophyllum_lusitanicum@6812:1.00000050003e-06)100:0.322394276114)42:0.0277593744386)42:0.0423170412953,blue_Nepenthaceae_Nepenthes_alata@85433:0.216188745092)13:0.0161457887812,(blue_Droseraceae_Dionaea_muscipula@48108:0.258854402013,(((outgroup_Amborella_trichopoda@evm_27_model_AmTr_v1_0_scaffold00071_156:0.426480249177,((outgroup_Aquilegia_coerulea@Aqcoe6G010400_1:0.140132688544,outgroup_Aquilegia_coerulea@Aqcoe6G010200_1:0.249208415282)61:0.0724924230133,((outgroup_Zea_mays@Zm00008a013145_P01:0.0522620018993,outgroup_Zea_mays@Zm00008a032782_P01:0.101009498335)70:0.083846286344,(outgroup_Oryza_sativa@Os01t0878900_01:0.139702605473,(outgroup_Oryza_sativa@Os01t0878800_02:1.00000050003e-06,outgroup_Oryza_sativa@Os01t0878800_01:1.00000050003e-06)100:0.117502646818)100:0.223505781583)98:0.368763377311)50:0.0422207746664)59:0.0690541524215,(((outgroup_Populus_trichocarpa@Potri_004G135300_1:0.208653785333,(outgroup_Populus_trichocarpa@Potri_004G135200_4:0.0764610072188,(outgroup_Populus_trichocarpa@Potri_018G136800_1:0.00936388850823,outgroup_Populus_trichocarpa@Potri_T045400_1:1.00000050003e-06)99:0.0742228558494)99:0.153808836863)95:0.136055379112,((outgroup_Glycine_max@Glyma_20G214400_1:0.256716064194,(outgroup_Glycine_max@Glyma_09G156900_1:0.0621670151985,outgroup_Glycine_max@Glyma_16G207700_1:0.043820860631)94:0.16033721009)94:0.171766116703,(outgroup_Glycine_max@Glyma_09G157100_1:0.46175941102,(outgroup_Glycine_max@Glyma_16G207800_1:0.028439010553,outgroup_Glycine_max@Glyma_09G157000_1:0.183539879894)79:0.0463315056194)79:0.10699724804)71:0.0747566305492)55:0.0760353820259,((((outgroup_Vitis_vinifera@GSVIVT01031190001:0.0256653007821,outgroup_Vitis_vinifera@GSVIVT01006088001:1.00000050003e-06)94:1.653372061,outgroup_Vitis_vinifera@GSVIVT01018544001:0.213136249472)63:0.0554542474739,outgroup_Vitis_vinifera@GSVIVT01033293001:0.0200492076484)62:0.0348137824811,outgroup_Vitis_vinifera@GSVIVT01000251001:0.0838410031877)74:0.0737423107332)25:0.0276842030798)29:0.044197619284,((outgroup_Solanum_lycopersicum@Solyc01g090180_3_1:0.277606741204,outgroup_Coffea_arabica@Cc07_g20640:0.343008889097)42:0.0841144654861,((outgroup_Solanum_lycopersicum@Solyc11g072550_2_1:0.302338404267,outgroup_Arabidopsis_thalinana@AT4G15093_1:0.436816722819)24:0.0188066298664,outgroup_Coffea_arabica@Cc11_g03920:0.479777401031)20:0.0337535007977)20:0.0561779037806)26:0.07451027086)13:0.0475201197687)6:0.0237981169035)6:0.0238729787111,((unknow_Simmondsiaceae_Simmondsia_chinensis@100924:0.189600963278,(blue_Plumbagianceae_Limonium_californicum@45428:0.232131193752,(blue_Polygonaceae_Rumex_hastatulus@16758:0.344598423995,(((blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0071717_1_g000001_aua_1:1.00000050003e-06,(blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0004016_1_g000002_aua_1:1.00000050003e-06,(blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0004016_1_g000012_aua_1:0.117004146291,blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0112001_1_g000001_aua_1:0.0190768407961)99:0.121485347836)60:0.0253442949593)99:0.106331605743,((blue_Polygonaceae_Persicaria_virginiana@54779:0.0857392960927,blue_Polygonaceae_Persicaria_virginiana@54783:1.00000050003e-06)91:0.0513962488511,blue_Polygonaceae_Persicaria_virginiana@54780:1.00000050003e-06)97:0.0985315130729)71:0.0516459100847,(blue_Polygonaceae_Polygonum_aviculare@57297:1.00000050003e-06,(blue_Polygonaceae_Polygonum_aviculare@57296:0.0152913818203,blue_Polygonaceae_Polygonum_aviculare@57295:0.0157085187162)85:0.019380761212)98:0.15237605624)56:0.0209571040719)94:0.118636776846)85:0.128781974648)32:0.0478493621649,unknown_Physenaceae_Physena_madagascariensis@83959:0.180516013582)23:0.0562197202921)6:0.0416687112748,(red_Stegnospermataceae_Stegnosperma_halimifolium@50374:0.648279810586,(((blue_Macarthuriaceae_Macarthuria_australis@44248:0.0124010328353,blue_Macarthuriaceae_Macarthuria_australis@44247:0.0111550655233)98:0.168312638375,((unknow_Achatocarpaceae_Achatocarpus_gracilis@8287:1.00000050003e-06,(unknow_Achatocarpaceae_Phaulothamnus_spinescens@22347:1.00000050003e-06,unknow_Achatocarpaceae_Phaulothamnus_spinescens@30700:1.00000050003e-06)60:0.00582980293387)99:0.122423853759,(((((red_Amaranthaceae_Alternanthera_philoxeroides@86019:0.0650918391032,(red_Amaranthaceae_Beta_vulgaris@Bv4_077830_epmy_t1:0.0195146423759,(red_Amaranthaceae_Beta_vulgaris@Bv4_077740_snqf_t1:0.0274875683039,(red_Amaranthaceae_Beta_vulgaris@Bv4_077190_ccnz_t1:0.702250620046,(red_Amaranthaceae_Beta_vulgaris@Bv4_077770_spoh_t1:1.00000050003e-06,red_Amaranthaceae_Beta_vulgaris@Bv4_077220_kucj_t1:0.0139877857738)73:0.0237639603752)59:1.00000050003e-06)94:0.0632171196311)96:0.0803466894744)86:0.0609606246414,(red_Amaranthaceae_Alternanthera_philoxeroides@86014:0.0232819599774,red_Amaranthaceae_Alternanthera_philoxeroides@86017:0.0149069679789)94:0.0274363110871)68:0.0476858150693,(red_Chenopodiaceae_Spinacia_oleracea@Sp_189620_irkq_t1:0.14034902607,((red_Chenopodiaceae_Haloxylon_ammodendron@38175:0.0959680984296,red_Chenopodiaceae_Suaeda_glauca@52504:0.1181330888)31:0.0145171355125,((red_Chenopodiaceae_Chenopodium_quinoa@AUR62000602_RA:0.0222443983653,red_Chenopodiaceae_Chenopodium_quinoa@AUR62006953_RA:0.0238207330162)86:0.0425558886521,(red_Amaranthaceae_Amaranthus_hypochondriacus@AH010349_RA:0.167903535242,red_Chenopodiaceae_Halogeton_glomeratus@8530:0.324573525072)21:0.0337687859831)6:0.0159677437698)11:0.011131701743)2:0.0076572672178)33:0.0416230550603,(((red_Chenopodiaceae_Spinacia_oleracea@Sp_021190_czop_t1:0.142668019173,((red_Chenopodiaceae_Chenopodium_quinoa@AUR62012347_RA:0.0440833377229,red_Chenopodiaceae_Chenopodium_quinoa@AUR62022994_RA:0.0362098257302)97:0.0622828015524,(red_Amaranthaceae_Amaranthus_hypochondriacus@AH023177_RA:0.226611785068,(red_Chenopodiaceae_Chenopodium_quinoa@AUR62022991_RA:0.0372520148039,red_Chenopodiaceae_Chenopodium_quinoa@AUR62012345_RA:0.0296406692536)98:0.0404952312027)33:0.0360037621528)30:0.042973830363)30:0.0331461789366,((red_Amaranthaceae_Alternanthera_philoxeroides@59746:0.109727385325,((red_Chenopodiaceae_Suaeda_glauca@59843:0.158546287598,(red_Chenopodiaceae_Suaeda_glauca@59851:0.0827794173069,red_Chenopodiaceae_Suaeda_glauca@59838:1.00000050003e-06)100:0.104450982226)100:0.288492016851,(red_Chenopodiaceae_Haloxylon_ammodendron@27698:1.00000050003e-06,red_Chenopodiaceae_Haloxylon_ammodendron@84557:0.0401013036465)94:0.152118633137)77:0.0525355882634)35:0.0457146142091,red_Amaranthaceae_Beta_vulgaris@Bv2_029920_ptjc_t1:0.140531018989)13:0.0202192576867)99:0.222562245513,((red_Amaranthaceae_Amaranthus_hypochondriacus@AH010348_RA:0.418013579455,red_Amaranthaceae_Beta_vulgaris@Bv4_077860_jqdf_t1:0.140404750424)69:0.0550583265208,(red_Chenopodiaceae_Spinacia_oleracea@Sp_187560_ynow_t1:0.169284483664,(red_Chenopodiaceae_Chenopodium_quinoa@AUR62000600_RA:0.127035308782,red_Chenopodiaceae_Chenopodium_quinoa@AUR62006948_RA:0.208273277319)48:0.0607782665212)61:0.0436433396201)98:0.13433121035)97:0.117263252707)45:0.028761209398,(((blue_Caryophyllaceae_Spergularia_marina@12586:1.00000050003e-06,blue_Caryophyllaceae_Spergularia_marina@12590:0.0708983599054)85:0.0141765525137,blue_Caryophyllaceae_Spergularia_marina@12589:0.0423836284827)100:0.140513459984,(blue_Caryophyllaceae_Polycarpon_tetraphyllum@72832:0.0208665162873,(blue_Caryophyllaceae_Polycarpon_tetraphyllum@72830:1.00000050003e-06,((blue_Caryophyllaceae_Polycarpon_tetraphyllum@72834:0.0620786731196,blue_Caryophyllaceae_Polycarpon_tetraphyllum@72835:1.00000050003e-06)33:0.013312815008,(blue_Caryophyllaceae_Polycarpon_tetraphyllum@72838:0.0630690141182,blue_Caryophyllaceae_Polycarpon_tetraphyllum@72831:1.00000050003e-06)46:0.0497470084722)8:0.0102520339444)24:0.0361707344392)100:0.323204470543)97:0.108163013927)71:0.0551810315231)80:0.0460231906404)53:0.0327926923682,((((blue_Limeaceae_Limeum_aethiopicum@68164:0.703913656537,((((red_Didieraceae_Decarya_madagascariensis@29408:0.00504786656397,red_Didieraceae_Decarya_madagascariensis@29404:0.00453455164026)98:0.0147650367185,(((red_Anacampserotaceae_Anacampseros_filamentosa@32996:1.00000050003e-06,red_Anacampserotaceae_Anacampseros_filamentosa@32995:0.144633259033)100:0.160393197818,((red_Talinaceae_Talinum_fruticosum@30305:0.0513541873124,(red_Talinaceae_Talinum_fruticosum@30304:0.0300932696031,red_Talinaceae_Talinum_fruticosum@30303:1.00000050003e-06)27:0.0057272875112)97:0.0486290976427,(blue_Portulacaceae_Portulaca_oleracea_SRA@9226:0.142598063477,(red_Cactaceae_Opuntia_arenaria@90271:0.00286346607762,red_Cactaceae_Opuntia_arenaria@90268:0.0572947065624)99:0.0406036303321)64:0.0147009068631)40:0.0237005664637)80:0.0387982693978,red_Didieraceae_Decarya_madagascariensis@29406:0.0195913782293)47:0.00520175879545)81:0.050562959147,(red_Basellaceae_Basella_alba@24974:0.108787644441,(red_Montiaceae_Claytonia_nevadensis@33794:1.00000050003e-06,red_Montiaceae_Claytonia_nevadensis@33795:0.0489141182157)100:0.156769343229)59:0.037280707043)100:0.331012798304,(red_Aizoaceae_Mesembryanthemum_crystallinum@15579:0.314878530955,((red_Phytolaccaceae_Phytolacca_americana@24205:0.192587068548,((red_Nyctaginaceae_Boerhavia_ciliata@78055:0.0712257077249,(red_Nyctaginaceae_Commicarpus_scandens@9477:1.00000050003e-06,red_Nyctaginaceae_Commicarpus_scandens@9476:0.206796475121)100:0.0589447652263)100:0.258388890044,red_Phytolaccaceae_Rivina_humilis@15447:0.233259111355)34:0.0590215451687)47:0.0468064498957,red_Gisekiaceae_Gisekia_pharnaceoides@439:0.314077141929)72:0.0501799496883)99:0.184799085537)36:0.0310736343068)6:0.0579244740088,((blue_Kewaceae_Kewa_bowkeriana@59573:0.171226204832,((red_Phytolaccaceae_Rivina_humilis@4997:0.278071557072,(((red_Nyctaginaceae_Boerhavia_ciliata@53349:0.0640395197814,(red_Nyctaginaceae_Anulocaulis_annulatus@4970:1.00000050003e-06,red_Nyctaginaceae_Anulocaulis_annulatus@4971:0.0177672426299)90:0.0365636252802)57:0.01003317354,red_Nyctaginaceae_Boerhavia_ciliata@50164:0.150016967513)99:0.152452547943,red_Phytolaccaceae_Phytolacca_americana@40:0.127683522769)84:0.0687452805015)84:0.0995787224804,(red_Phytolaccaceae_Phytolacca_americana@14690:0.0988624870575,((red_Phytolaccaceae_Rivina_humilis@421:0.0292166250364,red_Phytolaccaceae_Rivina_humilis@420:1.00000050003e-06)99:0.0877269308588,(red_Sarcobataceae_Sarcobatus_vermiculatus@70860:0.114869901882,(red_Nyctaginaceae_Boerhavia_ciliata@16248:0.0439063917821,(red_Nyctaginaceae_Anulocaulis_annulatus@1049:0.00895819714886,red_Nyctaginaceae_Commicarpus_scandens@67003:0.113415314455)58:0.0243939569324)96:0.144960080333)35:0.0365676636161)34:0.0359884729771)48:0.0228386693822)38:0.0606384720903)33:0.0413098376907,(((((red_Aizoaceae_Delosperma_echinatum@65143:0.113343002194,red_Aizoaceae_Mesembryanthemum_crystallinum@37003:0.0800643730405)65:0.0447476440776,blue_Kewaceae_Kewa_bowkeriana@59574:1.27221620025)36:0.0189329291169,red_Gisekiaceae_Gisekia_pharnaceoides@69646:0.163433973089)46:0.0381454867639,(red_Phytolaccaceae_Rivina_humilis@4998:0.602461193837,red_Gisekiaceae_Gisekia_pharnaceoides@20109:0.0610639931338)57:0.23752221865)6:0.0238600429415,(((blue_Portulacaceae_Portulaca_oleracea_SRA@31090:0.23951117607,((red_Anacampserotaceae_Anacampseros_filamentosa@49054:0.0563600999188,red_Anacampserotaceae_Anacampseros_filamentosa@49055:1.00000050003e-06)100:0.617623275802,red_Didieraceae_Decarya_madagascariensis@7459:0.129782108119)64:0.0507264239777)40:0.0982751355127,((((red_Didieraceae_Decarya_madagascariensis@820:1.00000050003e-06,red_Didieraceae_Decarya_madagascariensis@819:1.00000050003e-06)100:0.119619573996,(red_Didieraceae_Decarya_madagascariensis@37748:0.223839920574,red_Basellaceae_Basella_alba@10747:0.0710524210091)41:0.0345237329096)32:0.0341565452711,((red_Cactaceae_Opuntia_arenaria@90869:0.54711050801,red_Cactaceae_Opuntia_arenaria@90870:1.00000050003e-06)88:0.0548605611695,((blue_Portulacaceae_Portulaca_oleracea_SRA@58254:0.0322522260078,blue_Portulacaceae_Portulaca_oleracea_SRA@58252:0.186624845484)93:0.113927472006,(blue_Portulacaceae_Portulaca_oleracea_SRA@58255:1.00000050003e-06,blue_Portulacaceae_Portulaca_oleracea_SRA@58253:0.0559260279786)69:0.0413434156758)65:0.0607221321249)21:0.0089626285605)10:0.0176964143689,red_Montiaceae_Claytonia_nevadensis@33142:0.148116943125)5:0.0318106957008)3:0.0097697518217,blue_Portulacaceae_Portulaca_oleracea_SRA@58257:0.149733368299)8:0.00433431258927)2:0.0307428363496)19:0.0396773744486)15:0.0457702284728,(unknown_Microteaceae_Microtea_debilis@174616:0.291433884905,blue_Limeaceae_Limeum_aethiopicum@68165:0.270121843777)13:0.0195031292037)8:0.0180102469356,red_Stegnospermataceae_Stegnosperma_halimifolium@1183:0.195477808334)6:0.0061178864135)51:0.155768673644)46:0.158957479493)2:0.0472965681448)0:0.00987529304223)2:0.0348702553549,(red_Stegnospermataceae_Stegnosperma_halimifolium@56270:1.00000050003e-06,red_Stegnospermataceae_Stegnosperma_halimifolium@56269:0.00565862058507)100:0.209823739846)8:0.0414095351517)3:0.0328220368533)3:0.06645549128,(blue_Polygonaceae_Bistorta_bistortoides@72793:0.204146176305,(blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0000031_1_g000051_aua_1:0.0823949985446,blue_Polygonaceae_Fagopyrum_esculentum@Fes_sc0000031_1_g000052_aua_1:0.346863879859)99:0.213262225963)100:0.647882120908)96:0.230009731892,red_Nyctaginaceae_Commicarpus_scandens@40276:0.129934563369)49:0.0572269717227,(red_Nyctaginaceae_Boerhavia_ciliata@69613:0.159106231322,red_Nyctaginaceae_Anulocaulis_annulatus@27851:0.0151942477554)47:0.0459910420282)96:0.136964073569)100:0.47270153196)98:0.0817617987936)100:0.020620013969,red_Nyctaginaceae_Anulocaulis_annulatus@27852:1.00000050003e-06):0;

The species tree is as :
(red_Aizoaceae_Mesembryanthemum_crystallinum:0.0450604,red_Aizoaceae_Delosperma_echinatum:0.0661721,((((((((red_Talinaceae_Talinum_fruticosum:0.04799,((blue_Portulacaceae_Portulaca_oleracea_SRA:0.086149,red_Anacampserotaceae_Anacampseros_filamentosa:0.0710509)0.227561:0.0243853,red_Cactaceae_Opuntia_arenaria:0.0592674)0.204741:0.0204919)0.192849:0.0215703,red_Didieraceae_Decarya_madagascariensis:0.0756626)0.158779:0.0255401,red_Basellaceae_Basella_alba:0.0977187)0.172439:0.0252723,red_Montiaceae_Claytonia_nevadensis:0.142093)0.489594:0.0515856,(((blue_Molluginaceae_Mollugo_verticillata:0.0762864,blue_Molluginaceae_Glinus_lotoides_var_virens:0.0952005)0.517477:0.0502438,(((blue_Molluginaceae_Suessenguthiella_caespitosa:0.0280913,blue_Molluginaceae_Pharnaceum_exiguum:0.0355129)0.899156:0.0675433,blue_Molluginaceae_Mollugo_cerviana:0.0546466)0.763359:0.0706458,blue_Molluginaceae_Mollugo_nudicaulis:0.0725936)0.241623:0.0323627)0.139654:0.0280663,blue_Molluginaceae_Mollugo_pentaphylla:0.0869379)0.349136:0.0419757)0.0882282:0.0184875,((unknown_Microteaceae_Microtea_debilis:0.115314,((((((((outgroup_Zea_mays:0.123303,outgroup_Oryza_sativa:0.103723)0.98409:0.389963,outgroup_Amborella_trichopoda:0.333062)0.54552:0.092093,outgroup_Aquilegia_coerulea:0.235061)0.46131:0.0676768,(outgroup_Vitis_vinifera:0.143998,((outgroup_Populus_trichocarpa:0.172474,outgroup_Glycine_max:0.253492)0.263158:0.0539216,outgroup_Arabidopsis_thalinana:0.328932)0.292166:0.0451055)0.114825:0.0309429)0.126476:0.0274854,(outgroup_Solanum_lycopersicum:0.205785,outgroup_Coffea_arabica:0.176291)0.696264:0.0778749)0.646766:0.0867224,(((blue_Tamaricaceae_Tamarix_hispida:0.123722,blue_Frankeniaceae_Frankenia_laevis:0.130271)0.792527:0.108798,((((blue_Polygonaceae_Polygonum_aviculare:0.0996086,blue_Polygonaceae_Rumex_hastatulus:0.114415)0.490317:0.043164,blue_Polygonaceae_Fagopyrum_esculentum:0.115399)0.246203:0.0315209,(blue_Polygonaceae_Persicaria_virginiana:0.065746,blue_Polygonaceae_Bistorta_bistortoides:0.0667565)0.569225:0.0429856)0.821776:0.120827,blue_Plumbagianceae_Limonium_californicum:0.230353)0.629088:0.0671343)0.211169:0.0295152,((blue_Nepenthaceae_Nepenthes_alata:0.143197,blue_Droseraceae_Dionaea_muscipula:0.210585)0.289514:0.0445673,(blue_Drosophyllaceae_Drosophyllum_lusitanicum:0.12284,blue_Ancistrocladaceae_Ancistrocladus_robertsoniorum:0.092231)0.577822:0.0480156)0.467979:0.0461271)0.166734:0.0249106)0.367376:0.0394082,unknow_Simmondsiaceae_Simmondsia_chinensis:0.146691)0.307513:0.0399512,unknown_Physenaceae_Physena_madagascariensis:0.130179)0.44781:0.0583514)0.0903978:0.0304107,(((red_Stegnospermataceae_Stegnosperma_halimifolium:0.128971,blue_Macarthuriaceae_Macarthuria_australis:0.0993202)0.128646:0.026511,blue_Limeaceae_Limeum_aethiopicum:0.107646)0.0405785:0.0221891,(((unknow_Achatocarpaceae_Phaulothamnus_spinescens:0.0109613,unknow_Achatocarpaceae_Achatocarpus_gracilis:0.0125008)0.906147:0.0845274,((((red_Chenopodiaceae_Chenopodium_quinoa:0.084165,red_Chenopodiaceae_Spinacia_oleracea:0.0576629)0.696344:0.0359273,(red_Chenopodiaceae_Suaeda_glauca:0.0880383,(red_Chenopodiaceae_Haloxylon_ammodendron:0.0228602,red_Chenopodiaceae_Halogeton_glomeratus:0.0253254)0.904701:0.0704752)0.623865:0.0342216)0.188027:0.0168836,red_Amaranthaceae_Beta_vulgaris:0.0647222)0.212696:0.0173415,(red_Amaranthaceae_Amaranthus_hypochondriacus:0.116166,red_Amaranthaceae_Alternanthera_philoxeroides:0.139112)0.446766:0.0346947)0.627079:0.0499217)0.311209:0.0232434,(((((blue_Caryophyllaceae_Silene_latifoliasra:0.076559,blue_Caryophyllaceae_Eremogone_hookeri_ssp_desertorum:0.0918624)0.247248:0.0267116,((blue_Caryophyllaceae_Arenaria_serpyllifolia:0.0866925,blue_Caryophyllaceae_Cerastium_arvense:0.103289)0.306308:0.0289678,blue_Caryophyllaceae_Dianthus_caryophyllus:0.116041)0.1227:0.0224426)0.434873:0.0324654,blue_Caryophyllaceae_Spergularia_marina:0.107852)0.548895:0.0430384,blue_Caryophyllaceae_Polycarpon_tetraphyllum:0.15344)0.41366:0.0363426,blue_Caryophyllaceae_Paronychia_jamesii:0.158148)0.607955:0.0537153)0.289675:0.0271112)0.0171957:0.0133588)0.05681:0.0157166)0.119807:0.0196167,((red_Phytolaccaceae_Phytolacca_americana:0.055412,(((red_Nyctaginaceae_Commicarpus_scandens:0.0309045,(red_Nyctaginaceae_Boerhavia_ciliata:0.0291848,red_Nyctaginaceae_Anulocaulis_annulatus:0.0191104)0.597589:0.0282621)0.923423:0.131052,red_Phytolaccaceae_Rivina_humilis:0.123754)0.213821:0.0288509,red_Sarcobataceae_Sarcobatus_vermiculatus:0.159495)0.1272:0.0240484)0.146726:0.0212431,red_Gisekiaceae_Gisekia_pharnaceoides:0.0978699)0.168582:0.0211642)0.100763:0.0189328,blue_Kewaceae_Kewa_bowkeriana:0.0993581)0.751145:0.0619381);

I have confirmed that the species names in gene tree and species tree are completely identifcal.
However, when I run the command:

profileNJ -s species_tree.tre -g gene_tree.tre --sep '@' --spos 'prefix' -r best

I got error:
Traceback (most recent call last):
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/bin/profileNJ", line 4, in
import('pkg_resources').run_script('ProfileNJ==20161018', 'profileNJ')
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/pkg_resources/init.py", line 765, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/pkg_resources/init.py", line 1545, in run_script
exec(code, namespace, namespace)
File "/panfs/roc/groups/0/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/ProfileNJ-20161018-py2.7.egg/EGG-INFO/scripts/profileNJ", line 285, in
verbose=args.verbose, sol_limit=args.sol_limit, method=args.cluster, path_limit=args.path_limit)
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/ProfileNJ-20161018-py2.7.egg/profileNJ/PolytomySolver/Multipolysolver.py", line 660, in solvePolytomy
path_limit) + method), ptree, sptree, matrice, order, path_limit, cluster_method=method, verbose=verbose)
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/ProfileNJ-20161018-py2.7.egg/profileNJ/TreeLib/memorize.py", line 26, in call
output = self.function(*args, **kwargs)
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/ProfileNJ-20161018-py2.7.egg/profileNJ/PolytomySolver/Multipolysolver.py", line 140, in polySolver
:], verbose=verbose, method=cluster_method, cost=cost_table[xsize - 1, 0]))
File "/home/yangya/clingyun/biosofts/python2.7.6_installed/lib/python2.7/site-packages/ProfileNJ-20161018-py2.7.egg/profileNJ/PolytomySolver/Multipolysolver.py", line 454, in constructFromPath
[dup_tree.get_child_at(0).name, dup_tree.get_child_at(1).name])
AttributeError: 'NoneType' object has no attribute 'get_child_at'

The example data in your software package can run. Could you kindly help me resolve this prolem?

Thank you so much

Best

Lingyun

Hi Lingyun,

ProfileNJ only accepts rooted species tree as input, since its a requirement of the gene tree polytomy resolution algorithm.

But yes, the error is not informative at all, and I will fix it as soon as I have some time.

If you root your species tree, it should work. I have noticed, there are outgroups in the species tree. Try rooting the tree with the corresponding clade (or maybe just midpoint rooting).
Also please note that, as correction relies on the species tree, different rooting will produce different results. Although the method is somewhat robust to that, you should still keep that in mind.

Best

Emmanuel

Hi Emmanuel,
Thank you for help. After rooting the species tree, the analysis working.
I still have one issue. I did not see any change on the topology of my gene tree, although there is obviously some strange topology (bootstrap value is lower than 30 for them) comparing to the species tree. Is it possible to set the profileNJ more aggressive? Can I use the bootstrap rather than distance matrix?
Thank you so much.
Best regards,
Lingyun

Hello Lingyun,

You just need to define a bootstrap threshold to contract edges with low support.

(Since the threshold is not defined and you input genetree is not rooted profileNJ only test the best rooting with the lowest reconciliation score, therefore, RF = 0 compared to your original tree).

If you set the threshold to 90 for example (any subtree with parental branch below that threshold will be modified), the correct command line is:

profileNJ -s species_tree.tre -g gene_tree.tre --sep '@' --spos 'prefix' -r best --seuil 90

For more options see: profileNJ -h

As a side note, you could test several threshold (and even different dup/loss cost), then use consel to evaluate the statistical significance of corrected trees against your input, that I assume is the, maximum likelihood tree. This will allow you to somewhat explore the distribution of trees that are equivalent to the best ML tree, but are in agreement with the species phylogeny.

Best,

Emmanuel

Thanks Emmanuel,
After used the command '--seuil 90', it working!
Best
Lingyun