Rhabdovirus accessory genes.

The Rhabdoviridae is one of the most ecologically diverse families of RNA viruses with members infecting a wide range of organisms including placental mammals, marsupials, birds, reptiles, fish, insects and plants. The availability of complete nucleotide sequences for an increasing number of rhabdoviruses has revealed that their ecological diversity is reflected in the diversity and complexity of their genomes. The five canonical rhabdovirus structural protein genes (N, P, M, G and L) that are shared by all rhabdoviruses are overprinted, overlapped and interspersed with a multitude of novel and diverse accessory genes. Although not essential for replication in cell culture, several of these genes have been shown to have roles associated with pathogenesis and apoptosis in animals, and cell-to-cell movement in plants. Others appear to be secreted or have the characteristics of membrane-anchored glycoproteins or viroporins. However, most encode proteins of unknown function that are unrelated to any other known proteins. Understanding the roles of these accessory genes and the strategies by which rhabdoviruses use them to engage, divert and re-direct cellular processes will not only present opportunities to develop new anti-viral therapies but may also reveal aspects of cellar function that have broader significance in biology, agriculture and medicine. Crown

Introduction

The Rhabdoviridae is arguably the most diverse family of RNn class="Gene">A viruses. Rhabdoviruses have been isolated from a wide range of organisms including placental mammals, marsupials, birds, reptiles, fish, insects and plants (Kuzmin et al., 2009). Their ecology includes marine, freshwater and terrestrial environments. Transmission can be via insect vectors in which they replicate, by direct horizontal transfer via plant sap (wounds or grafts), aerosol, animal bite or immersion in infected water, or by vertical transovarial passage. The Rhabdoviridae is one of only three virus families that include members that infect plants, vertebrates and invertebrates. Replication can occur in the cytoplasm or the nucleus of infected cells. Infection can be either acute or persistent and rhabdoviruses can cause infections that may be unapparent or result in disease that can vary from being mild to severe or even fatal. Over 200 rhabdoviruses have been identified to date and most of them are poorly described (Dietzgen et al., 2012).

For many years, knowledge of rhabdovirus molecular biology has been based primarily on extensive studies of n class="Species">vesicular stomatitis virus (VSV) and rabies virus (RABV). These rhabdoviruses, share very similar and relatively simple genome organisations with only five structural protein genes (N, P, M, G and L) that are transcribed sequentially from the (–) RNA genome as monocistronic mRNAs. Although some variations are now evident, such as the long non-coding region in the RABV G gene (designated pseudogene ψ) (Tordo et al., 1986) and the small alternative open reading frames (ORFs) in the VSV P gene (C and C′) (Peluso et al., 1996), these viruses have generally been regarded as the rhabdovirus archetypes and as very suitable models for studying virus–host interactions. However, it is becoming evident that the ecological diversity displayed by rhabdoviruses is also apparent in the complexity of their genome organisation and in the increasing array of additional proteins which they encode. Although the functions of many of these proteins are yet unknown, preservation of the five canonical structural protein genes (3′-N-P-M-G-L-5′ in negative polarity) by all known rhabdoviruses suggests that the additional proteins may not be essential for virus replication in culture and are likely to have functions associated with enhancing replication efficiency, blocking host innate immune defences, modulating cellular signalling pathways, re-directing normal cellular functions, and allowing effective cell-to-cell transmission. This expectation has been supported by limited studies to date on some of these proteins. By analogy with other RNA viruses with relatively complex genome organisations (e.g., coronaviruses, lentiviruses, paramyxoviruses), we have therefore adopted the term ‘accessory genes’ to describe the ORFs encoding these novel and potentially informative proteins.

The Rhabdoviridae

Rhabdoviruses are non-segmented (–) an class="Species">ssRNA viruses that, together with paramyxoviruses, filoviruses and bornaviruses, are classified in the order Mononegavirales. The family Rhabdoviridae currently comprises six genera containing 46 species, five species unassigned to a genus and more than 150 viruses that have not been formne">ally classified (Dietzgen et al., 2012). Several likely species and two proposed new genera (an class="Species">Sigmavirus and Tibrovirus) are currently under review for classification, and a number of other viruses have been identified as rhabdoviruses by morphology and/or serology but have not been further characterised.

The genus Lyssavirus includes n class="Species">Rabies virus and the rabies-related viruses Mokola virus, Lagos bat virus, Duvenhage virus, European bat lyssavirus 1, European bat lyssavirus 2 and Australian bat lyssavirus, Irkut virus, Khujand virus, Aravan virus and West Caucasian bat virus as currently designated species. Another unclassified lyssavirus, Shimoni bat virus, has also recently been identified (Kuzmin et al., 2010). Lyssaviruses naturally infect mammals and almost all have been isolated from bats which appear to be the primary natural reservoir (Delmas et al., 2008). Transmission is usually by transfer of saliva in a bite puncture wound. They are neurotropic rhabdoviruses and most have been associated with rabies encephalitis in humans or animals (Rupprecht et al., 2002).

The n class="Disease">genus Vesiculovirus comprises an ecologically diverse but genetically similar group of viruses. Viruses causing vesicular stomatitis in cattle, pigs and horses in the Americas are referred to collectively as serotypes of vesicular stomatitis virus (VSV) and classified as the species Vesicular stomatitis Indiana virus, Vesicular stomatitis New Jersey virus and Vesicular stomatitis Alagoas virus. Other recognised mammalian vesiculovirus species include Cocal virus, Piry virus, Carajas virus and Maraba virus that are endemic in the Americas, Isfahan virus from Asia and Chandipura virus which occurs in both Asia and Africa. Most of these viruses have been isolated from phlebotomine sand flies (Lutzomyia, Phlebotomus and Sergentomyia spp.) or black flies (Simulium spp.) which are considered the primary vectors. However, biting midges (Culicoides spp.), mosquitoes and other insects have also been implicated as potential vectors and transmission can also occur by direct contact between infected vertebrates (Comer and Tesh, 1991, Letchworth et al., 1999, Rodríguez, 2002, Stallknecht et al., 2001). Vesiculoviruses can also infect humans and some have been associated with serious disease (Rao et al., 2004). The genus Vesiculovirus also currently contains the recognised species Spring viraemia of carp virus and several other important pathogens of fish, including pike fry rhabdovirus, trout rhabdovirus, sea trout rhabdovirus, and Siniperca chuatsi rhabdovirus, ulcerative disease rhabdovirus, eel virus American and eel virus European X, also appear to be members of the genus (Dietzgen et al., 2012, Hoffmann et al., 2005). Transmission of fish vesiculoviruses can occur by immersion in infected water and does not appear to involve arthropod vectors (Ahne et al., 2002).

The genus Novirhabdovirus comprises a second group of fish rhabdoviruses that appear to have a very different evolutionary history from the fish vesiculoviruses. Species include n class="Species">Infectious hematopoietic necrosis virus and Viral hemorrhagic septicemia virus which are important pathogens of salmonids, but can also infect a very wide range of other fish species (Walker and Winton, 2010). The species Hirame rhabdovirus and Snakehead rhabdovirus are pathogens of marine fish in temperate and tropical regions of Asia, respectively (Johnson et al., 1999, Kim et al., 2005). As for the fish vesiculoviruses, transmission is directly through mucosal surfaces or the skin by immersion in infected water and does not appear to involve an arthropod vector.

The genus Ephemerovirus comprises arthropod-borne rhabdoviruses that have been isolated from mosquitoes and/or n class="Species">biting midges (Culicoides spp.) and infect primarily ruminants (Walker, 2005). The type species Bovine ephemeral fever virus and the newly recognised species Kotonkan virus each cause a disabling febrile disease in cattle and water buffalo (Kemp et al., 1973, St. George and Standfast, 1988). Other ephemeroviruses, including recognised species Berrimah virus and Adelaide River virus, and likely species Obodhiang virus, Kimberley virus, Malakal virus and Puchong virus, have been isolated only from insects or from healthy cattle and their association with disease is not yet known (Walker et al., in press).

Plant-adapted rhabdoviruses are separated taxonomically into two genera, Cytorhabdovirus and Nucleorhabdovirus, based on their site of virus replication and morphogenesis (Dietzgen et al., 2012). Cytorhabdoviruses replicate in the cytoplasm of n class="Disease">infected cells and virions bud in association with endoplasmic reticulum membranes. The type species of the genus is Lettuce necrotic yellows virus and it includes eight additional recognised species. Nucleorhabdoviruses replicate in the nuclei of infected plant cells. Morphogenesis occurs at the inner nuclear membrane and virions accumulate in perinuclear spaces. The type species is Potato yellow dwarf virus and there are nine additional recognised species in this genus. Viruses in both genera infect monocot and dicot plants and are transmitted by leafhoppers, planthoppers or aphids in which they also replicate (Jackson et al., 2005).

Two other genera are currently under consideration for inclusion in the Rhabdoviridae. The proposed genus n class="Species">Sigmavirus comprises rhabdoviruses infecting only flies (Drosophila spp.) with which they have co-evolved. They are transmitted vertically and confer CO2 sensitivity on the host. There are five proposed species including Drosophila melanogaster sigmavirus (commonly called sigma virus of Drosophila) which is the most extensively studied of the sigmaviruses (Longdon et al., 2010, Longdon et al., 2011). Like the genus Ephemerovirus, the proposed genus Tibrovirus comprises antigenically related rhabdoviruses that primarily infect cattle and water buffalo, and are transmitted by Culicoides species midges. However, tibroviruses and ephemeroviruses have different evolutionary histories as is reflected in their phylogenetic relationships and genome organisations. Tibrogargan virus and Coastal Plains virus, each isolated in Australia, are currently proposed as species in the new genus (Gubala et al., 2011). Bivens Arm virus, isolated from Culicoides insignis in the USA, also infects cattle and water buffalo and may be a third species in the proposed genus (Gibbs et al., 1989).

Complete genome sequences are also available for several other n class="Species">unclassified animne">al rhabdoviruses that are considered in this review. Short sequences in the L and N gene have been used to determine phylogenetic relationships between viruses in the established and proposed genera and mne">any other unclassified animal rhabdoviruses (Bourhy et al., 2005, Kuzmin et al., 2006). One monophyletic cluster identified by this approach includes a geographically and ecologically diverse group of viruses termed the Hart Park group because of the close serological relationship of one of the members, Flanders virus, to Hart Park virus. Both of these viruses have been isolated from mosquitoes and infect birds in North America (Boyd, 1972). The Hart Park phylogenetic group also includes Wongabel virus which was isolated from biting midges (Culicoides spp.) in Australia and appears to infect birds, and Ngaingan virus, also isolated from biting midges in Australia but whose vertebrate hosts include marsupials (Doherty et al., 1973, Gubala et al., 2010). Not included in the Hart Park group but related to each other phylogenetically (but not serologically) are tupaia rhabdovirus isolated from hepatocarcinoma cells of a tree shrew (Tupaia belangeri) imported to the USA from Thailand (Kurz et al., 1986), and Durham virus, which was isolated from the brain of a moribund American coot (Fulica americana) in the USA (Allison et al., 2011). Oak Vale virus, isolated from mosquitoes in Australia and infecting feral pigs, and mandarin fish (Siniperca chuatsi) rhabdovirus (SCRV) are other currently unclassified rhabdoviruses that are considered in this review (Cybinski and Muller, 1990, Zhang et al., 2007).

Rhabdoviruses for which complete genome sequences are available are listed in Table 1 .

Table 1

Rhabdoviruses for which complete genome sequences are currently available.

Genus and species	Abbreviation	GenBank accessiona	Genome length (nt)	ORFs
Vesiculovirus
Vesicular stomatitis Indiana virus	VSINV	NC001560	11,161	6
Vesicular stomatitis Alagoas virus	VSAV	EU373658	11,070	5
Chandipura virus	CHPV	GU212858	11,120	6c
Cocal virus	COCV	EU373657	11,003	6c
Isfahan virus	ISFV	AJ810084	11,088	6c
Maraba virus	MARAV	HQ660076	11,135	6
Spring viraemia of carp virus	SVCV	NC002803	11,019	5
Pike fry rhabdovirus	PFRV	FJ872827	11,097	5
Eel virus European X	EVEX	FN557213	11,806	6
Lyssavirus
Rabies virus	RABV	NC001542	11,932	6c
Mokola virus	MOKV	NC006429	11,940	6c
Duvenhage virus	DUVV	EU293119	11,976	6c
Lagos bat virus	LBV	EU293108	12,016	6c
European bat lyssavirus 1	EBLV1	NC009527	11,966	6c
European bat lyssavirus 2	EBLV2	NC009528	11,930	5
Australian bat lyssavirus	ABLV	NC003243	11,822	6c
Aravan virus	ARAV	EF614259	11,918	5
Irkut virus	IRKV	EF614260	11,980	6c
Khujand virus	KHUV	EF614261	11,903	6c
West Caucasian bat virus	WCBV	EF614258	12,278	6c
Shimoni bat virus	SHIBV	GU170201	12,045	6c
Ephemerovirus
Bovine ephemeral fever virus	BEFV	NC002526	14,900	12
Berrimah virus	BRMV	HM461974	15,024	10
Adelaide River virus	ARV	na	15,119	10
Kotonkon virus	KOTV	HM474855	15,867	11
Obodhiang virus	OBOV	HM856902	14,717	10
Kimberley virus	KIMV	na	15430b	10
Novirhabdovirus
Infectious hematopoietic necrosis virus	IHNV	NC001652	11,131	7c
Viral hemorrhagic septicemia virus	VHSV	NC000855	11,158	7c
Hirame rhabdovirus	HRRV	NC005093	11,034	6
Snakehead rhabdovirus	SHRV	NC000903	11,550	7c
Cytorhabdovirus
Lettuce necrotic yellows virus	LNYV	NC007642	12,807	6
Lettuce yellow mottle virus	LYMoV	NC011532	12,926	6
Northern cereal mosaic virus	NCMV	NC002251	13,222	9
Nucleorhabdovirus
Potato yellow dwarf virus	PYDV	GU734660	12,875	7
Maize fine streak virus	MFSV	NC005974	13,782	7
Maize mosaic virus	MMV	NC005975	12,133	6b
Sonchus yellow net virus	SYNV	NC001615	13,720	6
Rice yellow stunt virus	RYSV	NC003746	14,042	7
Taro vein chlorosis virus	TaVCV	NC006942	12,020	6
Maize Iranian mosaic virus	MIMV	NC011542	12,381	6b
Sigmavirus (proposed genus and species)
Drosophila melanogaster sigma virus	DmSV	GQ375258	12,625	7c
Drosophila obscura sigma virus	DoSV	GQ410979	12,676	7c
Tibrovirus (proposed genus and species)
Tibrogargan virus	TIBV	GQ294472	13,298	8
Coastal Plains virus	CPV	GQ294473	13,203	9
Unassigned species
Flanders virus	FLAV	AH012179	12,965	9
Ngaingan virus	NGAV	NC013955	15,764	15
Wongabel virus	WONV	NC011639	13,196	10
Tupaia virus	TRV	NC007020	11,440	7
Unclassified rhabdoviruses
Durham virus	DURV	FJ952155	11092b	7
Oak Vale virus	OAKV	GQ294474	9414b	8
Siniperca chuatsi rhabdovirus	SCRV	NC008514	11,545	7c
Moussa virus	MOUV	FJ985748	11,526	7c

Where multiple isolates of a virus has been sequenced, the NCBI GenBank reference sequence is presented (denoted by the prefix NC).

Incomplete terminal UTR sequences.

Includes predicted ORF in the P-gene.

Rhabdovirus genome organisation and transcription strategy

The (–) ssRNn class="Gene">A rhabdovirus genome ranges in size from approximne">ately 11 kb to almost 16 kb. The five structural protein genes are arranged in the sne">ame linear order (3′-N-P-M-G-L-5′) and mne">ay be interspersed with one or more additional genes. At each end of the genome, 3′ leader and 5′ trailer sequences of approximately 50–100 nt contain the promoter sequences that initiate replication of the genome and anti-genome, respectively. The 3′ leader and 5′ trailer sequences of plant-adapted rhabdoviruses are considerably longer at 84–206 nt and 145–389 nt, respectively (Jackson et al., 2005). In transcription mode, the RNA-dependent RNA polymerase moves progressively along the genome to transcribe the short leader RNA (l) and then to transcribe, cap and polyadenylate mRNAs from each gene in decreasing molar abundances (Banerjee and Barik, 1992). Critical to this process are transcription initiation (TI) and transcription termination/polyadenylation (TTP) sequences flanking each gene (Barr et al., 2002). These sequences are relatively well conserved in each species and are similar in all species within a genus. The termination of transcription and polyadenylation of mRNAs are critically dependent on a stretch of seven uridine residues that are usually preceded by a guanosine residue. Following polyadenylation, transcription is re-initiated at the next available TI sequence which follows a non-transcribed intergenic region that may be from one to approximately 50 nt in length (Banerjee and Barik, 1992). Corruption of the [U]7 stretch in TPP sequences of animal rhabdoviruses leads to read-through, resulting in polycistronic mRNAs (Barr et al., 1997). In plant-adapted rhabdoviruses, the incorruptible component of the polyadenylation signal is UUMUUUU(U) (Jackson et al., 2005). As discussed below, rhabdovirus genes may contain multiple long open reading frames (ORFs), either overlapping or sequentially in the transcriptional units, and these must be expressed by internal initiation of translation. All evidence suggests that the fundamental processes of genome replication and transcription are common to all rhabdoviruses infecting vertebrates, insects and plants.

Rhabdovirus accessory genes

The five major structural proteins of n class="Species">VSV and RABV are amongst the most intensely studied of all viral proteins. Although there are reported variations in their ancillary functions, the primary functions and structural characteristics of these proteins are conserved across all rhabdoviruses and have been reviewed in detail on many previous occasions (Albertini et al., 2008, Assenberg et al., 2010, Coll, 1995, Jackson et al., 2005, Jayakar et al., 2004, Redinbaugh and Hogenhout, 2005, Roche et al., 2008, Rupprecht et al., 2002). In contrast, the number and locations of accessory genes varies widely both within and between taxa (Fig. 1, Fig. 2 ). In most cases, levels of amino acid sequence identity between cognate accessory proteins are quite low and, with a few exceptions, their specific functions do not appear to be conserved between taxa. As discussed in detail below, some occur as alternative or overlapping ORFs within the major structural protein genes; many occur in the regions between the major structural protein genes as new genes containing single or multiple ORFs. Over 35 rhabdovirus accessory genes have been identified to date and it is likely that many more will emerge as the number of complete genome sequences continues to grow.

Fig. 1

Genome organisation of animal rhabdoviruses (including sigma virus of Drosophila). The positions of ORFs are indicated and those encoding accessory proteins are shaded. Abbreviations of virus names are as listed in Table 1.

Fig. 2

Genome organisation of plant-adapted rhabdoviruses. The positions of ORFs are indicated and those encoding accessory proteins are shaded. Abbreviations of virus names are as listed in Table 1.

Genome organisation of animal rhabdoviruses (including n class="Species">sigma virus of Drosophila). The positions of ORFs are indicated and those encoding accessory proteins are shaded. Abbreviations of virus names are as listed in Table 1.

Genome organisation of plant-adapted rhabdoviruses. The positions of ORFs are indicated and those encoding accessory proteins are shaded. an class="Gene">Abbreviations of virus names are as listed in Table 1.

Vesiculovirus P and M genes

n class="Gene">Although the vesiculovirus genome is the smne">allest and arguably the least complex of all rhabdoviruses, containing only the 5 canonical structural protein genes, additional proteins are expressed from both the P and M genes (Fig. 1). The P genes of most vesiculoviruses encode a smne">all highly basic protein (C) in an alternative ORF. In VSV New Jersey and VSV Indiana, two forms of the C protein (C and C′) are expressed in infected cells by translation initiation at alternative start codons (Peluso et al., 1996, Spiropoulou and Nichol, 1993). There is evidence that the C protein is associated with nucleocapsids and increases the activity and fidelity of viral RNA polymerase activity in vitro (Peluso et al., 1996). However, C and C′ are not essential for replication in cell culture and a VSV infectious clone lacking these proteins has been shown to have identical growth characteristics to wild-type virus (Kretzschmar et al., 1996). The predicted sizes of the putative C proteins in viruses examined to date range from 41 amino acids (Piry virus) to 93 amino acids (Cocal virus). Interestingly, although spring viraemia of carp virus has the potential to encode highly basic proteins of 83 and 57 amino acids in an alternative ORF in the P gene (Hoffmann et al., 2002, Teng et al., 2007), there appears to be no alternative ORFs in the pike fry rhabdovirus P gene (Chen et al., 2009), suggesting its function may not be conserved across all fish vesiculoviruses. However, these data should be viewed cautiously as a single point mutation during adaptation to cell culture would be sufficient and inhibit expression of these accessory proteins by codon modification. It has also been reported that a small (7 kDa) protein is expressed from the VSV P gene from the same ORF as the P protein (Herman, 1986). The presence of alternative ORFs in the P gene of vesiculoviruses mirrors the situation in paramyxoviruses in which the complex expression of the C and V family of proteins in overlapping open reading frame or by RNA editing is very well characterised (Nagai and Kato, 2004). The paramyxovirus C and V proteins have various functions associated with infection and pathogenesis including regulation of viral RNA synthesis, virion assembly and budding, inhibition of apoptosis, and in binding to components of the Jak/STAT signalling pathway to antagonize the cellular interferon response (Fontana et al., 2008, Irie et al., 2010, Kurotani et al., 1998).

The n class="Species">VSV M protein is a multi-functional protein that, as well as being a structural component of virions, is involved during infection in inhibition of viral transcription, assembly and budding, inhibition of host gene expression and perturbation of the cell cytoskeleton (Black and Lyles, 1992, Blondel et al., 1990, Coulon et al., 1990). The M protein also blocks the host anti-viral response and promotes apoptosis both indirectly through inhibition of host gene expression and directly by targeting the extrinsic pathway and promoting expression of pro-apoptotic genes (Kopecky et al., 2001, Pearce and Lyles, 2009). The VSV M gene is also expressed in infected cells from two alternative initiation codons in the same frame, producing the M2 and M3 proteins (Jayakar and Whitt, 2002). Little is known of the function of these proteins but they appear to be involved in cell rounding and cytopathology. Although alternative initiation codons are available in the M proteins of other vesiculoviruses, it is not known if they are utilised.

Lyssavirus ‘pseudogenes’ (ψ) and P genes

The lyssavirus genome also conforms to the basic canonical arrangement of the genes encoding the five major structural proteins but lyssaviruses share the common characteristic of a very long 3′ non-coding region in the G gene (Fig. 1). The region varies in length from approximately 400 to 700 nt in different lyssaviruses and, in one lineage of laboratory-adapted n class="Species">RABV strains, it is truncated by an additional TTP signal 70 nucleotides downstream of the G protein ORF, leading to the suggestion that it mne">ay be the remnant of an ancestral gene or pseudogene (ψ) (Conzelmne">ann et al., 1990, Morimoto et al., 1989, Tordo et al., 1986). However, although the long 3′ non-coding region is preserved, the additional TTP is absent from wild-type rabies viruses isolated directly from infected animals and from other lyssaviruses (Ravkov et al., 1995). In West Caucasian bat virus (WCBV), the region extends 697 nucleotides beyond the G protein ORF and, uniquely, contains an ORF of 180 nucleotides but there is no evidence of expression of the second ORF in WCBV-infected cells (Kuzmin et al., 2008). Universal preservation of the ψ region in wild-type lyssaviruses suggests it has functional significance but a RABV infectious clone lacking the region displayed normal growth characteristics in fibroblasts and neuronal cells in culture and no difference in axonal spread or pathogenicity in mice (Ceccaldi et al., 1998, Schnell et al., 1994). Nevertheless, there is some evidence that, in conjunction with the G protein, the RABV ψ region does contribute to neuroinvasiveness from peripheral sites of infection (Faber et al., 2004). Proposed functions of the region include attenuation of L gene expression, stabilisation of the G mRNA and interaction of the mRNA with host cell proteins (Ravkov et al., 1995).

Like the vesiculoviruses, lyssaviruses also commonly encode small, highly basin class="Gene">c proteins in alternative ORFs in the P gene (Nadin-Dne">avis et al., 2002). The lyssne">avirus ORFs usually occur in the sne">ame region of the P gene as vesiculovirus C ORFs but they vary in size in different isolates and do not appear to be selectively retained by a locally constrained mutation frequency (Nadin-Davis et al., 2002). It is not known if lyssavirus C proteins are expressed and the variability of their detection in different isolates suggests they may have little functional relevance. However, if their role is to facilitate infection in vivo, sequence data on cell culture-adapted isolates should be regarded cautiously as introduced point mutations in the hypervariable P gene may be selected to interrupt functional C ORFs.

In addition to the full-length n class="Gene">P protein, RABV P gene also encodes four smaller products (P2–P5) in the same reading frame as a result of leaky scanning of the ribosome on the P mRNA (Chenik et al., 1995). Due to a nuclear localisation signal (NLS) in the C-terminal domain and a nuclear export signal (NES) in the N-terminal domain, the various forms of the P protein accumulate differentially in the cell such that the P1 (P) and P2 forms are cytoplasmic whereas P3, P4 and P5 are transported to the nucleus (Moseley et al., 2007, Pasdeloup et al., 2005). The P isoforms antagonise cellular interferon responses by several different mechanisms including inhibiting the phosphorylation of cytoplasmic interferon regulatory factor 3 (IRF3), binding to nuclear and cytoplasmic STAT1 and STAT2, and interacting with promyelocytic leukaemia (PML) nuclear bodies (Blondel et al., 2002, Brzozka et al., 2005, Brzozka et al., 2006, Chelbi-Alix et al., 2006). The P3 isoform has also been shown to mediate the sequestration of STAT1 away from the nucleus by a microtubule-associated mechanism, preventing nuclear import and activation of the interferon response (Moseley et al., 2009). The interferon-antagonistic activity of the P protein isoforms are important determinants of the pathogenicity of RABV in mice (Ito et al., 2010, Rieder et al., 2011).

Novirhabdovirus NV genes

In the region of the genome occupied by the lyssavirus ψ non-coding region (i.e., between the G and L genes), the novirhabdoviruses contain an additional gene of approximately 500 nucleotides (Kurath et al., 1985, Kurath et al., 1997, Kurath and Leong, 1985) (Fig. 1). The non-virion (NV) ORF is flanked by TI and TTP sequences, and encodes a neutral to mildly acidin class="Gene">c protein of 110–122 amino acids (13.2–13.6 kDne">a) (Kurath and Leong, 1985, Schutze et al., 1996). The NV proteins of infectious haematopoietic necrosis virus (IHNV) and viral hemorrhagic septicaemia virus (VHSV) are expressed at low levels in infected cells but do not occur in virions (Basurco and Benmansour, 1995, Schutze et al., 1996). The level of amino acid sequence conservation amongst NV proteins is quite low (e.g., VHSV and Hirame rhabdovirus share only 16.5% identity) but the NV gene is present in all four novirhabdovirus species and there is evidence that they are functionally equivalent (Kurath et al., 1997, Thoulouze et al., 2004). The construction of novirhabdovirus infectious clones lacking the NV gene has indicated that it is not essential for virus replication in cell culture but reports vary on the replication efficiency and pathogenicity of the recombinant viruses (Biacchesi, 2011). Deletion of the NV gene from an IHNV infectious clone resulted in reduced replication efficiency in EPC fish epithelial cells and reduced cumulative mortalities following intraperitoneal injection of rainbow trout (Oncorhynchus mykiss) (Biacchesi et al., 2000, Thoulouze et al., 2004) and similar observations have been reported for an NV-knockout infectious clone of VHSV in yellow perch (Perca flavescens) (Ammayappan et al., 2010). However, snakehead rhabdovirus (SHRV) infectious clones either lacking the NV gene or in which the NV ORF was truncated by point mutations were shown to replicate normally in EPC cells, and retained pathogenicity when injected intraperitoneally into zebra fish (Danio rerio) (Alonso et al., 2004, Johnson et al., 2000). It has been suggested that this apparent difference in the biological response to inhibition of NV expression may reflect the natural host range of the viruses which infect cold water fish species (VHSV and IHNV) or warm water fish species (SHRV) and have corresponding temperature ranges for efficient virus replication (Biacchesi, 2011). Recent evidence suggests that the VHSV NV protein has an anti-apoptotic function early in infection in vitro (Ammayappan and Vakharia, 2011). However, as the M protein of IHNV has pro-apoptotic activity (Chiou et al., 2000) and rhabdoviruses generally are known to promote apoptosis, the biological role of the NV protein remains unresolved.

Novirhabdoviruses also encode small basian class="Gene">c proteins in alternative reading frames in the P gene but their functional relevance is unknown.

Sigma virus X gene

n class="Species">Sigmaviruses also contain a single additional gene but it is located between the P and M genes (Teninges and Bras-Herreng, 1987) (Fig. 1). The X gene encodes a polypeptide (PP3) of 298–321 amino acids (Landes-Devauchelle et al., 1995). Although amino acid sequence identity between the PP3 proteins of sigmaviruses infecting Drosophila melanogaster and Drosophila obscura is relatively low (∼11%), sequence similarity is relatively high and each is predicted with high probability to have an N-terminal signal peptide. Each also has three predicted N-glycosylation sites near the N-terminus and so the mature X protein appears to be a secreted glycoprotein. Although original reports suggested that PP3 of Drosophila melanogaster sigmavirus may be distantly related to reverse transcriptase of retrotransposons, the relevant amino acid motifs do not appear to be conserved in the Drosophila obscura sigmavirus PP3 protein (Landes-Devauchelle et al., 1995, Longdon et al., 2010). As sigmaviruses have co-evolved with their hosts and are transmitted only vertically in ecological time frames, further studies of the function of PP3 may reveal interesting aspects of virus–insect interactions and innate anti-viral immunity in insects (Longdon et al., 2011, Tsai et al., 2008).

Accessory genes in plant-adapted rhabdoviruses

The genomes of plant rhabdoviruses are 12–14.5 kb in size and have the same basic organisation as their animal-infecting counterparts, 3′-N-P-M-G-L-5′ (Fig. 2). The unique feature of plant-adapted rhabdoviruses is that they encode a viral movement protein (MP) between P and M genes which facilitates cell-to-cell transport (Jackson et al., 2005). There are no characteristic genus-specific differences in genome size or genome organisation between cytorhabdoviruses and nucleorhabdoviruses. However, nucleorhabdoviruses replicate in the nuclei of n class="Disease">infected plant cells and most viral proteins are directed to the nucleus, while cytorhabdoviruses replicate in the cytoplasm (as all other rhabdoviruses) and viral proteins are mne">ainly localised to the cell periphery and endomembranes (Jackson et al., 2005). Individual members of both genera mne">ay encode up to four additional accessory proteins of unknown function interspersed between N-P, P-M and G-L genes (Fig. 2). The number and location of these additional genes transcends the two genera. Rice yellow stunt virus (RYSV) and rice transitory yellowing virus (RTYV) are strains of the same virus and the nucleotide sequences of their genomes are 98.5% identical (Hiraguri et al., 2010). All other fully sequenced plant rhabdoviruses are only distantly related and have been classified taxonomically as members of distinct virus species.

Genes located between the P and M genes encode 30–35 kDa proteins that appear to be involved in viral intercellular cell-to-cell movement across plasmodesmata, symplastic connections that traverse the plant cell wall of neighbouring cells (Benitez-n class="Gene">Alfonso et al., 2010). ORFs at this position encode the sc4, 4b and Y proteins of sonchus yellow net virus (SYNV), lettuce necrotic yellows virus (LNYV) and potato yellow dwarf virus (PYDV), respectively (Bandyopadhyay et al., 2010, Dietzgen et al., 2006, Scholthof et al., 1994). Other plant-adapted rhabdoviruses also contain accessory genes at this location. Single accessory genes (designated P3) are located at this position in rice yellow stunt virus (RYSV), rice transitory yellowing virus (RTYV), maize mosaic virus (MMV), maize Iranian mosaic virus (MIMV), taro vein chlorosis virus (TaVCV), lettuce yellow mottle virus (LYMoV) and strawberry crinkle virus (SCV) (Heim et al., 2007, Hiraguri et al., 2010, Huang et al., 2005, Massah et al., 2008, Reed et al., 2005, Revill et al., 2005, Schoen et al., 2004). However, maize fine streak virus (MFSV) contains two genes (P3 and P4) and in northern cereal mosaic virus (NCMV) there are four genes (P3, P4, P5 and P6) at this position in the genome (Fig. 2) (Tanno et al., 2000, Tsai et al., 2005).

Secondary structure predictions have shown that the n class="Species">SYNV an class="Gene">sc4, LNYV 4b, LYMoV P3, MMV P3, MFSV P4, RYSV P3 and PYDV Y proteins have some structural similarities in common with those of the ‘30K’ superfamily of plant virus movement proteins (Melcher, 2000). BlastP database searches and detailed secondary structure prediction of the LNYV 4b and LYMoV P3 proteins also revealed close similarities with the movement proteins of viruses in the family Flexiviridae which induce movement-associated tubules in infected cells (Dietzgen et al., 2006, Heim et al., 2007). LYMoV P3 has a 58% amino acid sequence identity with LNYV 4b (Heim et al., 2007). The membrane and cell periphery associations of SYNV sc4 also suggest a role in cell-to-cell movement.

Direct experimental evidence of cell-to-cell movement function exists only for n class="Species">RYSV P3 which is the most characterised protein encoded in the P-M region (Huang et al., 2005). P3 has been shown to trans-complement cell-to-cell movement of a movement-defective heterologous plant virus and binds non-specifically ssRNA in vitro, a common property of MPs of numerous plant viruses (Benitez-Alfonso et al., 2010). Moreover, pull-down assays have shown that P3 interacts with the N protein, providing support for involvement of the nucleocapsid core in cell-to-cell movement (Huang et al., 2005).

When expressed as fusion proteins with green fluorescent protein (GFP), n class="Species">PYDV Y and SYNV sc4 have been shown to localize to the cell periphery (Bandyopadhyay et al., 2010, Goodin et al., 2007) whereas MFSV P4 localized to nuclei (Tsai et al., 2005). The putative MPs of SYNV and PYDV self-interact in bimolecular fluorescence complementation (BiFC) assays (Bandyopadhyay et al., 2010, Min et al., 2010). PYDV Y and M proteins interacted in planta and the complex accumulated inside the nucleus, while SYNV sc4 interacted with a soluble truncated form of G protein and complexes were seen on the nuclear envelope and on punctuate loci along the cell periphery. Sc4 has also been shown to interact with microtubule-associated host proteins and cytoplasm-tethered transcription activators have been implicated in the intercellular movement of SYNV (Min et al., 2010). The available information for these four nucleorhabdoviruses indicates potentially different cell-to-cell movement pathways.

The functions of the other accessory genes, including the four unique n class="Species">NCMV ORFs located between P and M genes, are unknown. Three plant rhabdoviruses, NCMV, SCV and RYSV (and the synonymous RTYV) contain a short ORF between the G and L genes (Fig. 2). The NCMV P9, SCV P6, and RYSV P6 genes are predicted to encode proteins of 52, 68, and 93 amino acids, respectively. They share no identifiable sequence similarities but limited similarity can be identified with other rhabdovirus proteins. The RYSV P6 protein shares 24% sequence identity with the N-terminal 110 amino acids of the SYNV L protein and has limited sequence identity with the NV proteins of novirhabdoviruses (up to 25%) and non-coding regions preceding the L gene in Hendra virus (38.3%) and RABV (36.6%). The SCV P6 and NCMV P9 proteins are moderately basic (pI  ∼ 8.9) but RYSV P6 is highly acidic (pI 3.49). The large negative charge of this protein is reminiscent of the acidic N-terminal domain of the RYSV L protein, although they lack sequence similarity. Of the three small proteins, only RYSV P6 has been characterised experimentally (Huang et al., 2003). P6 is phosphorylated in vitro on both serine and threonine residues. It has also been shown to be associated with virions and may have a structural role but it is present only in low amounts during infection. The function of proteins encoded in this genome position is yet unknown. However, distant similarities of RYSV P6 with the L proteins of RYSV and SYNV may suggest a close evolutionary relationship with L proteins and potentially an ancillary role in rhabdovirus replication.

n class="Species">PYDV is unique amongst the sequenced plant-adapted rhabdoviruses in that it has an ORF X of unknown function located between the N and P genes. The 11 kDne">a X protein has a high content of proline, aspartic and glutamic acids, with a predicted pI of 4.5. These properties are similar to those of MFSV P3 (10 kDa, pI 5.4). However, while MFSV P3-GFP fusions accumulate in punctuate loci in the cytoplasm (Tsai et al., 2005), PYDV X is positive in yeast-based nuclear import assay, indicating its likely nuclear localization (Bandyopadhyay et al., 2010).

Like the animal rhabdoviruses, some cytorhabdoviruses and nucleorhabdoviruses present alternative ORFs in the P gene, but it is unknown if these are expressed in an class="Disease">infected plants or vector insects (Callaghan and Dietzgen, unpublished).

Ephemerovirus accessory genes

Ephemerovirus genomes are amongst the largest and most complex of known rhabdoviruses with multiple accessory genes located in a region that extends up to 4.2 kb between the G and L genes (Fig. 1). Perhaps the most interesting of these is a second glycoprotein (n class="Gene">GNS) gene that immediately follows the G gene (Walker et al., 1992). The GNS gene occurs in all ephemeroviruses and, until recently, has been regarded as a defining characteristic of the genus. The GNS proteins are 60–90 kDa type 1 transmembrane glycoproteins that are highly glycosylated with N-linked glycans (Walker et al., 1992, Wang and Walker, 1993). They share significant amino sequence identity and structural homology with the G proteins of ephemeroviruses and other rhabdoviruses, and appear to have arisen by recombination or gene duplication (see below). The bovine ephemeral fever virus (BEFV) GNS protein has been detected in infected cells but is not present in virions (Walker et al., 1991). Recombinant GNS expressed from a vaccinia virus vector has been detected in association with amorphous structures at the cell surface but not with budding virus particles (Hertig et al., 1996). Despite the structural similarity, the BEFV GNS protein shares no neutralising epitopes with the G protein (Johal et al., 2008) and polyclonal rabbit antibody to the BEFV GNS protein does not neutralise virus produced in either mammalian or insect cell cultures (Hertig et al., 1996). Furthermore, unlike the G protein, BEFV GNS does not induce spontaneous cell fusion at low pH (Johal et al., 2008), suggesting it functions neither in cell attachment nor fusion in vertebrates or invertebrates. In an attenuated vaccine strain of BEFV (Joubert and Walker, unpublished data) and in a highly cell-culture-adapted strain of Adelaide River virus (ARV), corruptions in the TTP sequence at the end of the G gene allow expression of the GNS protein only as a polycistronic mRNA (Wang and Walker, 1993) but this has not been observed in other wild-type or cell culture adapted ephemeroviruses. The function of ephemerovirus GNS proteins is currently unknown.

Immediately downstream of the n class="Gene">GNS gene lies a set of relatively smne">all accessory genes bounded by TI and TTP sequences that vary in number in different ephemeroviruses (McWilliam et al., 1997, Wang et al., 1994). The α gene occurs in all ephemeroviruses and contains two consecutive ORFs (α1 and α2) that are expressed from a single transcript due to the universne">al absence of intervening TTP and TI sequences. The α1 ORF encodes proteins of 10.5–14.5 kDne">a, each with a predicted N-terminal ectodomne">ain containing clusters of aromne">atic residues, a central transmembrane domne">ain and a highly basic C-terminal endodomne">ain (Fig. 3 ). Although sequence identity between the α1 protein is relatively low, the structure is highly conserved and strongly suggests they function as viroporins (Gonzalez and Carrasco, 2003, McWilliam et al., 1997). The α2 ORF encodes proteins that range in size from 10.7 kDa to 14.2 kDa in different ephemeroviruses and have no identifiable motifs that may suggest a particular function. They share identifiable amino acid sequence homology with significantly higher levels of identity (up to 59%) between the more closely related viruses. Conservation of the α2 ORF in all ephemeroviruses and the preservation of sequence identity suggest that it is of functional significance. In most ephemeroviruses the α2 ORF initiation codon is in close proximity to the α1 termination codon suggesting expression may occur by translation re-initiation as has been described for several other RNA viruses (Powell, 2010). In BEFV, ORF α2 also contains a 51 amino-acid ORF (α3) in an alternative reading frame but it does not occur in other ephemeroviruses and its functional significance is unknown (McWilliam et al., 1997).

Fig. 3

An alignment of the sequences of viroporin-like proteins encoded in the genomes of 10 animal rhabdoviruses. The alignment has been structured to emphasise the predicted N-terminal domains, transmembrane domains (shaded), and highly basic C-terminal domains. Aromatic residues in the N-terminal domain basic amino acids in the C-terminal domain are in bold and underlined. Abbreviations of virus names are as listed in Table 1.

n class="Gene">An alignment of the sequences of viroporin-like proteins encoded in the genomes of 10 animne">al rhabdoviruses. The alignment has been structured to emphasise the predicted N-terminal domne">ains, transmembrane domne">ains (shaded), and highly basic C-terminal domne">ains. Aromatic residues in the N-terminal domain basic amino acids in the C-terminal domain are in bold and underlined. Abbreviations of virus names are as listed in Table 1.

n class="Gene">All ephemeroviruses also contain the β gene which is located immediately downstream of the α gene and is expressed as a monocistronic transcript encoding a 16.9–18.4 kDne">a neutral protein (McWilliam et al., 1997, Wang et al., 1994). The ephemerovirus β proteins share significant amino acid sequence identity and have been detected both in infected cells and in virions. However, in highly adapted laboratory strains of BEFV and Kimberley virus (KIMV) the β proteins are severely truncated by point mutations and are not expressed during infection, indicating that they are not essential for growth in vitro and are not essential components of virions (Walker, Joubert and Blasdell, unpublished data). They have no identifiable homology with any other known proteins but, although their function remains unknown, they may well have a key role in pathogenesis and/or blocking the host immune response.

n class="Gene">Although absent from the genomes of ARV and Obodhiang virus (OBOV), BEFV, KIMV, kotonkan virus (KOTV) and Berrimah virus (BRMV) each also contain a γ gene, immediately downstream of the β gene (McWilliam et al., 1997). Proteins encoded in the γ ORF are mildly basic and range in size from 11.8 kDa to 13.6 kDa. They share the highest levels of amino acid sequence identity of any of the ephemerovirus accessory proteins but are not related to any other known proteins and contain no identifiable sequence motifs that suggest a function. The γ proteins have been detected in infected cells and in virions. In wild-type viruses, the γ genes are transcribed as monocistronic mRNAs but in highly cell culture adapted strains, expression of the γ protein is severely attenuated by corruption of the TTP sequence at the end of the β gene, allowing expression only as a bicistronic β–γ transcript. Selective suppression of γ gene expression in culture suggests that its function may be important during infection in vivo.

The 15,870 nt n class="Species">KOTV genome is the largest yet discovered in any rhabdovirus due to an additional accessory gene immediately downstream of the γ gene (Walker and Blasdell, unpublished data). The KOTV δ gene is expressed as a monocistronic mRNA encoding a mildly acidic 12.4 kDa protein. The δ protein has been detected in KOTV-infected cells. It is shares significant amino acid sequence identity with the pleckstrin homology (PH) domain of coactivator-associated arginine methyl transferase (CARM1) which is involved in chromatin re-modelling and transcriptional activation of NF-κB dependent gene expression (Covic et al., 2005, Troffer-Charlier et al., 2007). KOTV δ may mimic this domain, inhibiting activation of the anti-viral innate immune response (Joubert and Walker, unpublished data).

Like most other rhabdoviruses, several ephemeroviruses also encode small highly basian class="Gene">c proteins in alternative ORFs in the P gene but their functional significance is unclear.

The complex genomes of the Hart Park group

Unlike the ephemeroviruses, the complex genomes of the n class="Species">Hart Park group of rhabdoviruses contain multiple accessory genes inserted at various locations (Fig. 1). Wongabel virus (WONV) contains five additional ORFs (U1–U5) interspersed at various locations across the 13,196 nt genome (Gubala et al., 2008). The U1, U2 and U3 ORFs are located between the P and M genes and are each bounded by TI and TTP sequences. ORF U1 and ORF U2 encode acidic proteins of similar size and net charge (21.2 kDa/pI 5.2 and 21.9 kDa/pI 4.5, respectively). They have no remarkable characteristics but do share a significant level of amino acid sequence homology (∼16% identity; 60% similarity) and are more closely related to each other than to any other known protein. The WONV U3 gene follows the U2 gene and also encodes a relatively small acidic protein (16.5 kDa/pI 4.6) that is unrelated to any other known proteins but shares sequence homology with the WONV U1 protein (∼22% identity; 57% similarity), including the common sequence [KSxYDFVWPxxxLxxG] in the central region of each protein (Fig. 4A). Preliminary data suggest that the U3 protein binds to the IAP (inhibition of apoptosis) protein apollon and a component of the chromatin remodelling complex and may be involved in promoting apoptosis (Joubert and Walker, unpublished data). WONV ORF U4 is located within the N gene as a second ORF in the 142 nt region that follows the N protein termination codon. It encodes a putative small acidic protein (5.8 kDa/pI 4.0) of 49 amino acids that could only be translated by internal initiation on a bicistronic (N-U4) mRNA. WONV ORF U5 overlaps the G protein ORF and thus also must be translated by internal initiation on a bicistronic (G-U5) mRNA. Translation from the first initiation codon would result in a 14.9 kDa protein with that is similar in size and structure to the ephemerovirus viroporin-like α1 proteins (Fig. 3). Several low molecular weight proteins (range 20–25 kDa) have been detected in infected cells using WONV hyperimmune antiserum but they are larger than the predicted U1–U5 proteins and may well be degradation products of larger viral proteins (Gubala et al., 2008).

Fig. 4

Clustal X alignments of (A) Wongabel virus (WONV) U1, U2 and U3 proteins and (B) Flanders virus (FLAV) U1 and U2 proteins. A single termination codon interrupting the published sequence of the FLAV U1 protein has been removed to infer the complete amino acid sequence. The FLAV U2 protein is also referred to as the 19K protein. Aligned amino acids are shaded. Fully conserved (*), strongly conserved (:) and weakly conserved (.) amino acids have been assigned by the software.

Clustal X alignments of (n class="Gene">A) Wongabel virus (WONV) U1, U2 and U3 proteins and (B) Flanders virus (FLAV) U1 and U2 proteins. A single termination codon interrupting the published sequence of the FLAV U1 protein has been removed to infer the complete amino acid sequence. The FLAV U2 protein is also referred to as the 19K protein. Aligned amino acids are shaded. Fully conserved (*), strongly conserved (:) and weakly conserved (.) amino acids have been assigned by the software.

The n class="Species">Flanders virus (FLan class="Gene">AV) genome organisation appears to closely resemble that of WONV. The available sequence data (GenBank AH012179) is dated and appears to be corrupted either by sequencing ambiguities or as a result of cell culture adaptation of the virus. However, through alignment with the WONV genome, it is reasonable to infer that FLAV also contains three accessory genes between the P and M genes and a 13.8 kDa viroporin-like protein encoded in a second ORF in the G gene (Fig. 3). The FLAV U1, U2 and U3 ORFs each appear to encode proteins of 18–19 kDa and, as for the corresponding WONV proteins, they share an unusual level of sequence similarity suggesting a common origin and/or functional association. This is most evident between FLAV U1 and U2 (also referred to as the 19K protein) that share ∼20% amino acid sequence identity and ∼58% similarity (Fig. 4B). There is no second consecutive ORF in the FLAV N gene corresponding to WONV U4 but there is a small alternative ORF within the N gene encoding a polypeptide of 66 amino acids (7.9 kDa). In addition to the N, P, M, G and L proteins, three virus-induced proteins have been reported in FLAV infected cells and in virions (Boyd and Whitaker-Dowling, 1988). However, only one of these proteins (19 kDa) corresponds to the predicted size of any FLAV accessory gene product.

The 15,764 nt n class="Species">Ngaingan virus (an class="Species">NGAV) genome is one of the largest and most complex of known rhabdoviruses (Gubala et al., 2010). As in WONV and FLAV, the NGAV genome contains three ORFs between the P and M genes (U1, U2 and U3) encoding neutral to mildly acidic proteins of 15.1–17.1 kDa that share identifiable sequence similarity. However, unlike WONV and FLAV, ORF U1 and ORF U2 are encoded in overlapping reading frames within the same gene and ORF3 is encoded in a separate gene bounded by TI and TTP sequences. Downstream of the U3 gene, the NGAV genome organisation is very complex and varies significantly from WONV and FLAV. Between the M and G genes the NGAV U4 gene encodes a 9.7 kDa protein with a tyrosine-rich C-terminal domain but otherwise has no remarkable features. Downstream of the G gene, however, is a region that resembles the ephemerovirus genomes including a second glycoprotein gene (GNS), followed by three additional ORFs. The GNS gene encodes a type 1 transmembrane glycoprotein that is most closely related to the ephemerovirus GNS proteins and shares sequence homology and structural characteristics with the NGAV G protein and those of all animal rhabdoviruses (Fig. 5 ). The two accessory genes immediately downstream of the NGAV GNS gene also resemble the ephemerovirus α and β genes in format, size and location. The first contains two overlapping ORFs (U5 and U6) that are similar in size to the ephemerovirus α1 and α2 ORFs. However, the 12.6 kDa U5 protein does not have the characteristic structure of the viroporin-like α1 proteins and, unlike the featureless ephemerovirus α2 proteins, the U6 protein has the structure of a small (10–12 kDa) membrane-anchored class 1 glycoprotein with a signal peptide domain, a transmembrane domain and an N-terminal ectodomain containing a single N-glycosylation site. The NGAV U7 gene contains a single ORF encoding an 18.2 kDa protein that is similar in size to ephemerovirus β proteins but shares no identifiable sequence homology and has no remarkable features.

Fig. 5

A minimum evolution phylogenetic tree constructed from a Clustal X multiple sequence alignment of animal rhabdovirus G protein and GNS protein sequences. The tree was constructed using the MEGA program. Viruses assigned to established genera and proposed genera are shaded. Abbreviations of virus names are as listed in Table 1.

n class="Gene">A minimum evolution phylogenetic tree constructed from a Clustal X multiple sequence alignment of animne">al rhabdovirus G protein and GNS protein sequences. The tree was constructed using the MEGA program. Viruses assigned to established genera and proposed genera are shaded. Abbreviations of virus names are as listed in Table 1.

In addition to these independent accessory genes n class="Species">NGAV encodes two smne">all basic proteins (P′1 and P′2) in alternative reading frames within the P gene and a 10.8 kDa protein in an alternative reading frame in the M gene. If all these potential ORFs in translated, NGAV could potentially express a total of 16 functional proteins.

Tibrovirus accessory genes

Whole genome sequences are available for two viruses that have recently been proposed as species in a new genus Tibrovirus. The 13.2–13.3 kb n class="Species">Tibrogargan virus (TIBV) and Coastal Plains virus (CPV) genomes share a similar organisation with 3–4 accessory protein genes (Gubala et al., 2011) (Fig. 1). In each virus, ORFs U1 and U2 are located between the M and G genes. Each is bounded by TI and TTP sequences and encodes a protein of 17–20 kDa. The U1 proteins are acidic and the U2 proteins are basic, and there is extensive amino acid sequence homology between the corresponding proteins. An N-terminal signal peptide is predicted with high probability in the CPV U1 protein, suggesting it may be secreted but it is not predicted in the TIBV U1 protein despite the high level of overall sequence homology and so its significance is uncertain. Features of the U2 proteins are unremarkable. The TIBV and CPV genomes also each contain ORF U3 encoded in an independent gene located between the G and L genes. The U3 proteins have the structural characteristics of viroporin-like proteins that occur in ephemeroviruses as the α1 proteins and in WONV and NGAV as the U5 proteins. In each case, the viroporin ORFs are located between the G and L genes, either directly following the G gene or following the GNS gene in those viruses in which it occurs. In TIBV, U4 is an alternative ORF that overlaps the G protein ORF and could potentially encode a very small basic protein (6.6 kDa/pI 9.6). However, U4 does not occur in CPV and its significance, if any, is unclear. As in many other rhabdoviruses, the CPV P gene also contains an alternative ORF encoding a small basic proteins but none is evident in the TIBV P gene.

The SH genes of other unclassified rhabdoviruses

Several other n class="Species">unclassified rhabdoviruses isolated from insects, mne">ammne">als, birds and fish share a similar genome organisne">ation (Fig. 1). Tupaia rhabdovirus (TUPV), Durham virus (DURV) and Oak Vale virus (OAKV) each encode small hydrophobic (SH) proteins in ORFs flanked by TI and TTP sequences in the region between the M and G genes (Allison et al., 2011, Springfeld et al., 2005). However, although superficially similar, the putative proteins expressed from these genes have different characteristics. The TUPV ORF encodes a small protein with a predicted signal peptide and a central transmembrane domain that would generate an 8.6 kDa non-glycosylated membrane-anchored protein with short ecto- and endo-domains. In DURV, the SH protein is predicted to have similar membrane topology with predicted signal peptidase cleavage generating a 6.9 kDa non-glycosylated membrane-anchored protein. These proteins resemble the small hydrophobic protein encoded in the NGAV U5 gene (see above). The OAKV SH protein is predicted to have an N-terminal signal peptide but no membrane anchor to generate a 4.2 kDa non-glycosylated secreted product (GenBank GQ294474). The mandarin fish (Siniperca chuatsi) rhabdovirus (SCRV) genome does not contain an independent SH gene but does encode a small hydrophobic protein as an alternative reading frame in the upstream M gene (Tao et al., 2008). The predicted structure is similar to that of the OAKV SH protein with a high probability signal peptidase cleavage site generated a secreted product of only 2.9 kDa. Although it has been reported that the putative SCRV SH protein could not be detected in infected fish cells, this may not be unexpected if it is indeed secreted. The X genes of sigmaviruses also encode a protein that is predicted to be secreted (see above) but it is very much larger than these small peptides. Several of these unclassified rhabdoviruses also encode small basic proteins in alternative ORFs in the P gene.

Patterns in the complexity of rhabdovirus genome organisation

n class="Gene">Although characterised by remne">arkable diversity and complexity, a global view of rhabdovirus genomes does reveal some patterns in the nature and location of accessory genes (Fig. 6 ). In various rhabdoviruses, accessory genes occur in each of the junctions between the structural protein genes and as alternative or overlapping ORFs within all structural protein genes except the L gene. Alternative or overlapping ORFs occur least commonly in the N gene (and of course the L gene), possibly reflecting the structural and functional constraints on amino acid sequence variations, and in the G gene that requires sufficient genetic flexibility to evolve under the pressure of neutralising antibody. Alternative ORFs occur most commonly in the P gene in which the preservation of function is far less dependent on conservation of the amino acid sequence, allowing more scope for the evolution and stabilisation of new ORFs in alternative frames (Nadin-Davis et al., 2002, Spiropoulou and Nichol, 1993). The N-P, P-M, M-G, and G-L gene junctions have all been utilised as sites for accessory genes, although in rhabdoviruses examined to date, the P-M and G-L junctions have been utilised more commonly. The primary evolutionary constraint on the occupation of these sites may be the resulting attenuation of transcription in the downstream genes. For example, rare use of the N-P gene junction may be due to the importance of the balance of expression of the N and P proteins in modulating replication and transcription (Banerjee and Barik, 1992).

Fig. 6

Schematic diagram illustrating the use by rhabdoviruses of the five common structural protein genes (N, P, M, G and L) and the intergenic regions as sites for insertion of ORFs encoding putative accessory proteins genes. The genes are defined as transcriptional units bounded by TI and TTP sequences.

n class="Gene">Alternative ORFs in the P gene occur very commonly and usually encode smne">all basic proteins. Although their function has not been examined in any detail, by analogy with paramyxoviruses, they may have multiple roles in replication, virion assembly and in modulating host responses to infection (Fontana et al., 2008, Irie et al., 2010, Kurotani et al., 1998). The P-M gene junction has been utilised as a site of accessory genes in many rhabdoviruses including genes encoding movement proteins in cytorhabdoviruses and nucleorhabdoviruses and the apparently secreted and glycosylated PP3 protein of sigmaviruses. The P-M junction is also the location of the set of three related 15–22 kDa proteins in the Hart Park group viruses that may be involved in the promotion of apoptosis. The M-G gene junction is occupied less commonly but genes encoding small proteins that are either secreted or membrane-anchored do occur here, including the TIBV U1 protein, and the small hydrophobic (SH) proteins of the unclassified viruses TUPV, DURV and OAKV (Allison et al., 2011, Springfeld et al., 2005). Interestingly, in the mandarin fish rhabdovirus (SCRV), and SH protein is encoded just upstream as an alternative ORF in the M gene (Tao et al., 2008).

The G-L gene junction is the most commonly occupied site of accessory genes in animal rhabdoviruses. Genes encoding the non-structural glycoprotein (n class="Gene">GNS) invariably occur immediately downstream of the G gene and share a common genetic lineage (Gubala et al., 2010, Walker et al., 1992, Wang and Walker, 1993). Genes encoding viroporin-like proteins also usually occur in this region, immediately following GNS gene in ephemeroviruses (McWilliam et al., 1997, Wang et al., 1994) and following the G gene in FLAV and tibroviruses which lack GNS (Gubala et al., 2011). Interestingly, the WONV viroporin-like protein is encoded immediately upstream of the G-L junction in an ORF that overlaps the end of the G gene. The viroporin-like proteins all have a similar structure with a central hydrophobic domain and a highly basic C-terminal tail, and so may have evolved from a common ancestor (Fig. 3). Many other accessory genes are located in the G-L intergenic region including the novirhabdovirus NV gene, and genes coding proteins of unknown function in ephemeroviruses, NGAV and some plant-adapted rhabdoviruses. The G-L intergenic region is also the site of the lyssavirus ψ region (Tordo et al., 1986). Sequence insertion at this site may be most favoured as it would primarily affect expression of the L protein which is required only in catalytic amounts.

The evolution of rhabdovirus genome complexity: gene duplication or recombination?

The ecological diversity of rhabdoviruses suggests they have an ancient evolutionary history, perhaps originating as insect viruses that have adapted to replication in plants or animals (Hogenhout et al., 2003). There is recent evidence that endogenous viral elements (EVEs) derived from rhabdoviruses occur commonly in insect genomes with which they have co-evolved over many millions of years (Katzourakis and Gifford, 2010). Rhabdovirus genomic diversity will also therefore have evolved over geological time frames, possibly involving rare events that may appear highly improbable in a contemporary ecological context. We have argued previously that the evolution of rhabdovirus accessory genes may have occurred by a copy-choice mechanism involving polymerase-jumping (Wang and Walker, 1993). n class="Gene">According to this model, upstream relocation of the polymerase during replication would generate a sequence duplication in a process mirroring the downstream relocation that has been proposed for the generation of defective-interfering (DI) genomes (Lazzarini et al., 1981, Perrault, 1981). There seems no primne">a face reason to suggest that upstream relocation and downstream relocation should not occur with similar frequency but purifying selection would enrich for rapidly replicating DI genomes and quickly remove larger genomes unless the duplication event has provided a specific selective advantage. Although survival of such larger genomes may be exceedingly rare, preservation of the ψ region in lyssaviruses does indicate that the selective advantage of maintaining long stretches of apparently non-functional sequence may sometimes be cryptic (Kuzmin et al., 2008, Ravkov et al., 1995) and sequence similarities in adjacent genes of some rhabdoviruses do suggest that duplication may have contributed to the evolution of accessory genes. Such sequence similarities are evident in the GNS genes in ephemeroviruses and NGAV (Gubala et al., 2010, Walker et al., 1992, Wang and Walker, 1993), the related U1, U2 and U3 genes in the P-M region of Hart Park group viruses (see above), and the P6 gene in rice yellow stunt virus (Huang et al., 2003).

Other possible mechanisms for the acquisition of accessory genes include homologous genetic recombination and lateral gene transfer (LGT). Life-long persistent n class="Disease">infections occur commonly in insects and mixed infections with different viruses have been reported (Chen et al., 2004, Mourya et al., 2001, Thavara et al., 2006). Although genetic recombination is believed to occur rarely in (–) ssRNA viruses (Holmes, 2009), possible recombination events have been detected and may have contributed to their long term evolution (Chare et al., 2003, Han et al., 2008, Lukashev, 2005, Qin et al., 2008, Sibold et al., 1999). Unlike gene duplication, in which the new gene must be preserved while evolving a new function, recombination could offer immediate selective advantage in larger more complex genomes by increasing tissue tropism, host range or replication efficiency or more effective blocking of host defensive mechanisms. Homologous genetic recombination must therefore be considered an equally plausible explanation for the generation of genomic structural diversity in rhabdoviruses. LGT also offers a mechanism by which accessory genes could be obtained from the host or an unrelated virus to provide immediate selective advantage (Holmes, 2009). There is evidence that the coronavirus haemagglutinin-esterase glycoprotein was acquired from influenza C virus (Luytjes et al., 1988) and that the envelope glycoprotein of the tick-borne orthomyxovirus Thogoto may have been derived from a baculovirus (Morse et al., 1992). Several other examples of possible LGT involving transfer of host genes to viruses have been reported but these events appear to be very rare and none has yet been reported to our knowledge in rhabdoviruses or other unsegmented (–) ssRNA viruses (Holmes, 2009).

It can also be argued that purifying selection should drive evolution towards smaller and perhaps more efficient genomes rather than large complex genomes. Indeed, the generation of DI particles, which occurs with high frequency during n class="Disease">infection in vitro, would appear to provide the mechanism for rapid elimination of redundant sequences (Lazzarini et al., 1981, Perrault, 1981) and phylogenetic analyses do suggest that accessory genes mne">ay have been lost during the evolution of some lineages. As shown in Fig. 5, phylogenetic analysis of available rhabdovirus G protein sequences places NGAV with WONV and FLAV to form the Hart Park group and analyses of the L and N protein sequences supports these relationships (Bourhy et al., 2005, Gubala et al., 2010). However, the NGAV GNS protein clusters with the ephemerovirus GNS proteins (Fig. 5), suggesting that it was obtained by a common ancestral virus and subsequently lost from WONV and FLAV. Alternatively, NGAV may have acquired GNS independently in a later event involving recombination with a virus in the ephemerovirus lineage.

an class="Gene">All arguments considered, perhne">aps the most congruent hypothesis is thne">at the evolution of rhne">abdovirus genomic diversity is a dynamic process in which accessory genes hne">ave been gne">ained and lost during the course of continual adaptation and purifying selection.

Future research

Despite the importance of rhabdoviruses as pathogens of n class="Species">humans, livestock, fish and plants, it is evident from this review that very little is known of the functions of the proteins encoded in their mne">any accessory genes. In some respects, this is very surprising as five mne">ajor rhabdovirus structural proteins have been the subject of intensive investigation and have provided some of the most important models for understanding fundamental processes during viral infection. On the basis of available evidence, and by analogy with other viruses that are rich in accessory functions (e.g., HIV, SARS coronavirus), we can anticipate that many of the rhabdovirus accessory proteins may be important determinants of virulence and pathogenesis and function to increase the efficiency of viral replication, block host innate and/or adaptive immune responses, regulate apoptosis or induce cytopathology. Understanding the strategies employed by rhabdoviruses to engage, divert and re-direct cellular processes will not only present opportunities to develop new anti-viral therapies but may also reveal aspects that have broader significance in biology, agriculture and medicine.

One of the most attractive targets for future study is the role of rhabdovirus accessory genes during n class="Disease">infection in insects and in transmission to mne">ammne">als or plants. Despite their importance as vectors of disease, and over 100 years of immunology, the anti-viral defences of insects remne">ain poorly understood (Huszar and Imler, 2008). RNAi has been identified as an important anti-viral mechanism and there is evidence of involvement of the Toll and Imd pathways that function primarily in innate immune responses to bacteria and parasites (Kemp and Imler, 2009). In Drosophila melanogaster, the innate immune response to sigmavirus infection has been shown to involve upregulated peptidoglycan reporter and antimicrobial peptide gene expression (Tsai et al., 2008). However, innate anti-viral sensor and effector mechanisms in insects are largely unexplored and may be critical determinants of their susceptibility to infection and ability to act as efficient vectors. As most rhabdoviruses naturally infect insects, their accessory genes may well be involved in engaging defensive responses and so provide useful tools for exploring invertebrate anti-viral pathways.

There has been a significant increase in the number of fully sequenced rhabdovirus genomes in recent years and the availability of next-generation sequencing platforms, rapid sample preparation methods and new multiplexing technologies will ensure that GenBank depositions continue to grow rapidly. Indeed, it can be anticipated that most or all known rhabdoviruses will be sequenced in the foreseeable future. This will reveal much about the phylogenetic relationships and evolutionary history of rhabdoviruses. It will also present the opportunity to characterise and exploit an increasing array of accessory genes to understand the fundamental processes of transmission, an class="Disease">infection and immunity.