Are all lyssavirus genes equal for phylogenetic analyses?

Individual lyssavirus genes were evaluated for phylogenetic studies from available full genome sequences. The full genome of the ERA rabies virus was sequenced and its accuracy was confirmed through virus recovery by reverse genetics. The full length of the ERA is 11,931 nucleotides (nt), with a leader sequence of 58 nt, the nucleoprotein (N) gene of 1350 nt, phosphoprotein (P) gene of 891 nt, matrix protein (M) gene of 606 nt, glycoprotein (G) gene of 1572 nt, RNA-dependent RNA polymerase (L) gene of 6384 nt, Psi-region (or G–L intergenic region) of 400 nt, and a trailer region of 70 nt. The five mono-cistrons are separated by intergenic regions of 2, 5, 5 and 24 nt, respectively. One obvious difference between the ERA and SAD-B19 rabies virus strains was the putative stop/polyadenylation signal of the G gene, with a poly(A8) tract for ERA, and a poly(A5) for SAD-B19. The TGpoly(A8) sequence tract was identified to be a leaky termination signal in the ERA strain. Through analyses of nt diversity, protein co-variations, structural and functional constraints, and reconstruction of phylogenetic trees from comprehensive datasets, we propose lyssavirus genes probably are of similar value for phylogenetic analyses.