Molecular characterisation of six badnavirus species associated with leaf streak disease of banana in East Africa

Banana leaf streak disease, caused by several species of Banana streak virus (BSV), is widespread in East Africa. We surveyed for this disease in Uganda and Kenya, and used rolling-circle amplification (RCA) to detect the presence of BSV in banana. Six distinct badnavirus sequences, three from Uganda and three from Kenya, were amplified for which only partial sequences were previously available. The complete genomes were sequenced and characterised. The size and organisation of all six sequences was characteristic of other badnaviruses, including conserved functional domains present in the putative polyprotein encoded by open reading frame (ORF) 3. Based on nucleotide sequence analysis within the reverse transcriptase/ribonuclease H-coding region of open reading frame 3, we propose that these sequences be recognised as six new species and be designated as Banana streak UA virus, Banana streak UI virus, Banana streak UL virus, Banana streak UM virus, Banana streak CA virus and Banana streak IM virus. Using PCR and species-specific primers to test for the presence of integrated sequences, we demonstrated that sequences with high similarity to BSIMV only were present in several banana cultivars which had tested negative for episomal BSV sequences.

Viruses of banana in East Africa

Bananas are one of the world's most important food crops, providing sustenance and income for millions of people in developing countries and supporting large export industries. Viruses are considered major constraints to banana production, germplasm multiplication and exchange, and to genetic improvement of banana through traditional breeding. In Africa, the two most important virus diseases are bunchy top, caused by Banana bunchy top virus (BBTV), and banana streak disease, caused by Banana streak virus (BSV). BBTV is a serious production constraint in a number of countries within/bordering East Africa, such as Burundi, Democratic Republic of Congo, Malawi, Mozambique, Rwanda and Zambia, but is not present in Kenya, Tanzania and Uganda. Additionally, epidemics of banana streak disease are occurring in Kenya and Uganda. The rapidly growing tissue culture (TC) industry within East Africa, aiming to provide planting material to banana farmers, has stimulated discussion about the need for virus indexing to certify planting material as virus-free. Diagnostic methods for BBTV and BSV have been reported and, for BBTV, PCR-based assays are reliable and relatively straightforward. However for BSV, high levels of serological and genetic variability and the presence of endogenous virus sequences within the banana genome complicate diagnosis. Uganda has been shown to contain the greatest diversity in BSV isolates found anywhere in the world. A broad-spectrum diagnostic test for BSV detection, which can discriminate between endogenous and episomal BSV sequences, is a priority. This PhD project aimed to establish diagnostic methods for banana viruses, with a particular focus on the development of novel methods for BSV detection, and to use these diagnostic methods for the detection and characterisation of banana viruses in East Africa. A novel rolling-circle amplification (RCA) method was developed for the detection of BSV. Using samples of Banana streak MY virus (BSMYV) and Banana streak OL virus (BSOLV) from Australia, this method was shown to distinguish between endogenous and episomal BSV sequences in banana plants. The RCA assay was used to screen a collection of 56 banana samples from south-west Uganda for BSV. RCA detected at least five distinct BSV isolates in these samples, including BSOLV and Banana streak GF virus (BSGFV) as well as three BSV isolates (Banana streak Uganda-I, -L and -M virus) for which only partial sequences had been previously reported. These latter three BSV had only been detected using immuno-capture (IC)-PCR and thus were possible endogenous sequences. In addition to its ability to detect BSV, the RCA protocol was also demonstrated to detect other viruses within the family Caulimoviridae, including Sugar cane bacilliform virus, and Cauliflower mosaic virus. Using the novel RCA method, three distinct BSV isolates from both Kenya and Uganda were identified and characterised. The complete genome of these isolates was sequenced and annotated. All six isolates were shown to have a characteristic badnavirus genome organisation with three open reading frames (ORFs) and the large polyprotein encoded by ORF 3 was shown to contain conserved amino acid motifs for movement, aspartic protease, reverse transcriptase and ribonuclease H activities. As well, several sequences important for expression and replication of the virus genome were identified including the conserved tRNAmet primer binding site present in the intergenic region of all badnaviruses. Based on the International Committee on Taxonomy of Viruses (ICTV) guidelines for species demarcation in the genus Badnavirus, these six isolates were proposed as distinct species, and named Banana streak UA virus (BSUAV), Banana streak UI virus (BSUIV), Banana streak UL virus (BSULV), Banana streak UM virus (BSUMV), Banana streak CA virus (BSCAV) and Banana streak IM virus (BSIMV). Using PCR with species-specific primers designed to each isolate, a genotypically diverse collection of 12 virus-free banana cultivars were tested for the presence of endogenous sequences. For five of the BSV no amplification was observed in any cultivar tested, while for BSIMV, four positive samples were identified in cultivars with a B-genome component. During field visits to Kenya, Tanzania and Uganda, 143 samples were collected and assayed for BSV. PCR using nine sets of species-specific primers, and RCA, were compared for BSV detection. For five BSV species with no known endogenous counterpart (namely BSCAV, BSUAV, BSUIV, BSULV and BSUMV), PCR was used to detect 30 infections from the 143 samples. Using RCA, 96.4% of these samples were considered positive, with one additional sample detected using RCA which was not positive using PCR. For these five BSV, PCR and RCA were both useful for identifying infected samples, irrespective of the host cultivar genotype (Musa A- or B-genome components). For four additional BSV with known endogenous counterparts in the M. balbisiana genome (BSOLV, BSGFV, BSMYV and BSIMV), PCR was shown to detect 75 infections from the 143 samples. In 30 samples from cultivars with an A-only genome component there was 96.3% agreement between PCR positive samples and detection using RCA, again demonstrating either PCR or RCA are suitable methods for detection. However, in 45 samples from cultivars with some B-genome component, the level of agreement between PCR positive samples and RCA positive samples was 70.5%. This suggests that, in cultivars with some B-genome component, many infections were detected using PCR which were the result of amplification of endogenous sequences. In these latter cases, RCA or another method which discriminates between endogenous and episomal sequences, such as immuno-capture PCR, is needed to diagnose episomal BSV infection. Field visits were made to Malawi and Rwanda to collect local isolates of BBTV for validation of a PCR-based diagnostic assay. The presence of BBTV in samples of bananas with bunchy top disease was confirmed in 28 out of 39 samples from Malawi and all nine samples collected in Rwanda, using PCR and RCA. For three isolates, one from Malawi and two from Rwanda, the complete nucleotide sequences were determined and shown to have a similar genome organisation to previously published BBTV isolates. The two isolates from Rwanda had at least 98.1% nucleotide sequence identity between each of the six DNA components, while the similarity between isolates from Rwanda and Malawi was between 96.2% and 99.4% depending on the DNA component. At the amino acid level, similarities in the putative proteins encoded by DNA-R, -S, -M, - C and -N were found to range between 98.8% to 100%. In a phylogenetic analysis, the three East African isolates clustered together within the South Pacific subgroup of BBTV isolates. Nucleotide sequence comparison to isolates of BBTV from outside Africa identified India as the possible origin of East African isolates of BBTV.

Global genetic diversity and geographical and host-species distribution of beak and feather disease virus isolates

Psittacine beak and feather disease (PBFD) has a broad host range and is widespread in wild and captive psittacine populations in Asia, Africa, the Americas, Europe and Australasia. Beak and feather disease circovirus (BFDV) is the causative agent. BFDV has an ~2 kb single stranded circular DNA genome encoding just two proteins (Rep and CP). In this study we provide support for demarcation of BFDV strains by phylogenetic analysis of 65 complete genomes from databases and 22 new BFDV sequences isolated from infected psittacines in South Africa. We propose 94% genome-wide sequence identity as a strain demarcation threshold, with isolates sharing > 94% identity belonging to the same strain, and strain subtypes sharing> 98% identity. Currently, BFDV diversity falls within 14 strains, with five highly divergent isolates from budgerigars probably representing a new species of circovirus with three strains (budgerigar circovirus; BCV-A, -B and -C). The geographical distribution of BFDV and BCV strains is strongly linked to the international trade in exotic birds; strains with more than one host are generally located in the same geographical area. Lastly, we examined BFDV and BCV sequences for evidence of recombination, and determined that recombination had occurred in most BFDV and BCV strains. We established that there were two globally significant recombination hotspots in the viral genome: the first is along the entire intergenic region and the second is in the C-terminal portion of the CP ORF. The implications of our results for the taxonomy and classification of circoviruses are discussed. © 2011 SGM.

Sequence and organization of barley yellow dwarf virus genomic RNA

The nucleotide sequence of the genomic RNA of barley yellow dwarf virus, PAV serotype was determined except for the 5′-terminal base, and its genome organization deduced. The 5,677 nucleotide genome contains five large open reading frames (ORFs). The genes for the coat protein (1) and the putative viral RNA-dependent RNA polymerase were identified. The latter shows a striking degree of similarity to that of carnation mottle virus (CarMV). By comparison with corona- and retrovirus RNAs, it is proposed that a translational frameshift is involved in expression of the polymerase. An ORF encoding an Mr 49,797 protein (50K ORF) may be translated by in-frame readthrough of the coat protein stop codon. The coat protein, an overlapping 17K ORF, and a 3′ 6.7K ORF are likely to be expressed via subgenomic mRNAs. © 1988 IRL Press Limited.

Rice ragged stunt oryzavirus genome segments S7 and S10 encode non-structural proteins of M(r) 68 025 (Pns7) and M(r) 32364 (Pns10) : brief report

The nucleotide sequences of genome segments S7 and S10 of a Thai-isolate of rice ragged stunt virus (RRSV) were determined. The 1938 bp S7 sequence contains a single large open reading frame (ORF) spanning nucleotides 20 to 1 843 that is predicted to encode a protein of M(r) 68 025. The 1 162 bp S10 sequence has a major ORF spanning nucleotides 142 to 1 032 that is predicted to encode a protein of M(r) 32364. This S10 ORF is preceded by a small ORF (nt 20-55) which is probably a minicistron. Coupled in vitro transcription-translation from the two major ORFs gave protein products of the expected sizes. However, no protein was visualised from S10 when the small ORF sequence was included. Proteins were expressed in Escherichia coli from the full length ORF of S7 (P7) and from a segment of the S10 ORF (P10) fused to the ORF of glutathione S-transferase (GST). Neither fusion protein was recognised by polyclonal antibodies raised against RRSV particles. Furthermore, polyclonal antibodies raised against GST-P7 fusion protein did not recognise any virion structural polypeptides. These data strongly suggest that the proteins P7 and P10 do not form part of RRSV particle. This is further supported by observed sequence homology (though very weak) of predicted.

The M(r) 43K major capsid protein of rice ragged stunt oryzavirus is a post-translationally processed product of a M(r) 67 348 polypeptide encoded by genome segment 8

The nucleotide sequence of DNA complementary to rice ragged stunt oryzavirus (RRSV) genome segment 8 (S8) of an isolate from Thailand was determined. RRSV S8 is 1 914 bp in size and contains a single large open reading frame (ORF) spanning nucleotides 23 to 1 810 which is capable of encoding a protein of M(r) 67 348. The N-terminal amino acid sequence of a ~43K virion polypeptide matched to that inferred for an internal region of the S8 coding sequence. These data suggest that the 43K protein is encoded by S8 and is derived by a proteolytic cleavage. Predicted polypeptide sizes from this possible cleavage of S8 protein are 26K and 42K. Polyclonal antibodies raised against a maltose binding protein (MBP)-S8 fusion polypeptide (expressed in Escherichia coli) recognised four RRSV particle associated polypeptides of M(r) 67K, 46K, 43K and 26K and all except the 26K polypeptide were also highly immunoreactive to polyclonal antibodies raised against purified RRSV particles. Cleavage of the MBP-S8 fusion polypeptide with protease Factor X produced the expected 40K MBP and two polypeptides of apparent M(r) 46K and 26K. Antibodies to purified RRSV particles reacted strongly with the intact fusion protein and the 46K cleavage product but weakly to the 26K product. Furthermore, in vitro transcription and translation of the S8 coding region revealed a post-translational self cleavage of the 67K polypeptide to 46K and 26K products. These data indicate that S8 encodes a structural polypeptide, the majority of which is auto- catalytically cleaved to 26K and 46K proteins. The data also suggest that the 26K protein is the self cleaving protease and that the 46K product is further processed or undergoes stable conformational changes to a ~43K major capsid protein.

Nucleotide sequence of coat protein gene for GPV isolate of barley yellow dwarf virus and construction of expression plasmid for plant

GPV is a Chinese serotype isolate of barley yellow dwarf virus (BYDV) that has no reaction with antiserum of MAV, PAV, SGV, RPV and RMV The sequence of the coat protein (CP) of GPV isolate of BYDV was identified and its amino acid sequence was deduced. The coding region for the putative GPV CP is 603 bases nucleotides and encodes a Mr 22 218 (22 ku) protein. The same as MAV, PAV and RPV, GPV contained a second ORF within the coat protein coding region. This protein of 17 024 Mr (17 ku) is thought to correspond to the Virion protein genome linked (Vpg). Sequence comparisons of the CP coding region between the GPV isolate of BYDV and other isolates of BYDV have been done. The nucleotide and amino acid sequence homology of GPV has a greater identity to the sequence of RPV than those of PAV and MAV. The GPV CP sequence stored 83.7% of nucleotide similarity and 77.5% of deduced amino acid similarity, whereas that of the PAV and MAV shared 56.9%, 53.2% and 44.1%, 43.8% respectively. According to BYDV-GPV CP sequence, two primers were designed. The cDNA of CP was produced by RT-PCR. Full-length cDNA of CP was inserted into plasmid to construct expression plasmids named pPPI1, pPPI2 and pPPI5 based on different promoters. The recombinant plasmids were identified by using α-32P-dATP labelled CP probe, α-32P-ATP labelled GPV RNA probe and sequencing to confirm real GPV CP gene cDNA in plasmids.

Characterisation of the subgenomic RNAs of an Australian isolate of barley yellow dwarf luteovirus

We have characterised the subgenomic RNAs of an Australian isolate of BYDV-PAV. Northern blot analyses of infected plants and protoplasts have shown that this isolate synthesises three subgenomic RNAs. Precise mapping of the transcription start sites of all three subgenomic RNAs and translational analyses of subgenomic RNA 2 and 3 have revealed a number of features. First, the transcription start site of subgenomic RNA 1 in this isolate differs markedly from the start site determined for an Illinois isolate of BYDV-PAV. Second, the start sites of subgenomic RNA 1 and 2 occur at a sequence that closely resembles the 5' end sequence of the genomic RNA (5'AGUGAAGA). Third, subgenomic RNA 2 appears to express ORF 6 of BYDV-PAV but the gene product is truncated due to the appearance of a new stop codon in the sequence. Last, subgenomic RNA 3, which is abundantly transcribed and encapsidated by the virus particle, appears to have no coding ability. We postulate that this novel subgenomic RNA has a regulatory function.

Infectious in vitro transcripts from a cloned cDNA of barley yellow dwarf virus

A full-length cDNA clone of barley yellow dwarf virus (BYDV-PAV serotype) has been constructed and fused to the bacteriophage T7 RNA polymerase promoter. RNA transcripts produced in vitro, either capped or uncapped, were infectious in Triticum monococcum protoplasts. Protoplasts inoculated with in vitro-transcribed BYDV RNA accumulated coat protein, synthesized new viral RNAs, and produced virus particles. Aphid feeding on extracts from protoplasts inoculated with in vitro RNA transcripts can be used to transfer the virus progeny to whole plants. Introduction of mutations which interrupt specific BYDV-PAV open reading frames (ORFs) V and VI eliminated infectivity while an ORF I mutant remained infectious. Infectious RNA transcripts derived from BYDV cDNA clones will facilitate analysis of the molecular aspects of BYDV infection and further enhance our understanding of this economically important virus.

Cytochrome P450 aromatase in grey mullet: cDNA and promoter isolation; brain, pituitary and ovarian expression during puberty.

In a study towards elucidating the role of aromatases during puberty in female grey mullet, the cDNAs of the brain (muCyp19b) and ovarian (muCyp19a) aromatase were isolated by RT-PCR and their relative expression levels were determined by quantitative real-time RT-PCR. The muCyp19a ORF of 1515 bp encoded 505 predicted amino acid residues, while that of muCyp19b was 1485 bp and encoded 495 predicted amino acid residues. The expression level of muCyp19b significantly increased in the brain as puberty advanced; however, its expression level in the pituitary increased only slightly with pubertal development. In the ovary, the muCyp19a expression level markedly increased as puberty progressed. The promoter regions of the two genes were also isolated and their functionality evaluated in vitro using luciferase as the reporter gene. The muCyp19a promoter sequence (650 bp) contained a consensus TATA box and putative transcription factor binding sites, including two half EREs, an SF-1, an AhR/Arnt, a PR and two GATA-3s. The muCyp19b promoter sequence (2500 bp) showed consensus TATA and CCAAT boxes and putative transcription binding sites, namely: a PR, an ERE, a half ERE, a SP-1, two GATA-binding factor, one half GATA-1, two C/EBPs, a GRE, a NFkappaB, three STATs, a PPAR/RXR, an Ahr/Arnt and a CRE. Basal activity of serially deleted promoter constructs transiently transfected into COS-7, [alpha]T3 and TE671 cells demonstrated the enhancing and silencing roles of the putative transcription factor binding sites. Quinpirole, a dopamine agonist, significantly reduced the promoter activity of muCyp19b in TE671. The results suggest tissue-specific regulation of the muCyp19 genes and a putative alternative promoter for muCyp19b.

Abaca bunchy top virus, a new member of the genus Babuvirus (family Nanoviridae)

Two isolates of a novel babuvirus causing "bunchy top" symptoms were characterised, one from abaca (Musa textilis) from the Philippines and one from banana (Musa sp.) from Sarawak (Malaysia). The name abacá bunchy top virus (ABTV) is proposed. Both isolates have a genome of six circular DNA components, each ca. 1.0-1.1 kb, analogous to those of isolates of Banana bunchy top virus (BBTV). However, unlike BBTV, both ABTV isolates lack an internal ORF in DNA-R, and the ORF in DNA-U3 found in some BBTV isolates is also absent. In all phylogenetic analyses of nanovirid isolates, ABTV and BBTV fall in the same clade, but on separate branches. However, ABTV and BBTV isolates shared only 79-81% amino acid sequence identity for the putative coat protein and 54-76% overall nucleotide sequence identity across all components. Stem-loop and major common regions were present in ABTV, but there was less than 60% identity with the major common region of BBTV. ABTV and BBTV were also shown to be serologically distinct, with only two out of ten BBTV-specific monoclonal antibodies reacting with ABTV. The two ABTV isolates may represent distinct strains of the species as they are less closely related to each other than are isolates of the two geographic subgroups (Asian and South Pacific) of BBTV.

p53 and little brother p53/47: linking IRES activities with protein functions

The tumor suppressor p53 represents a paradigm for gene regulation. Its rapid induction in response to DNA damage conditions has been attributed to both increased half-life of p53 protein and also increased translation of p53 mRNA. Recent advances in our understanding of the post-transcriptional regulation of p53 include the discovery of internal ribosome entry sites (IRESs) within the p53 mRNA. These IRES elements regulate the translation of the full length as well as the N-terminally truncated isoform, p53/47. The p53/47 isoform is generated by alternative initiation at an internal AUG codon present within the p53 ORF. The aim of this review is to summarize the role of translational control mechanisms in regulating p53 functions. We discuss here in detail how diverse cellular stress pathways trigger alterations in the cap-dependent and cap-independent translation of p53 mRNA and how changes in the relative expression levels of p53 isoforms result in more differentiated p53 activity.

Lost in Translation : Translation Mechanisms in Production of Cocksfoot Mottle Virus Proteins

The present study focuses on the translational strategies of Cocksfoot mottle virus (CfMV, genus Sobemovirus), which infects monocotyledonous plants. CfMV RNA lacks the 5'cap and the 3'poly(A) tail that ensure efficient translation of cellular messenger RNAs (mRNAs). Instead, CfMV RNA is covalently linked to a viral protein VPg (viral protein, genome-linked). This indicates that the viral untranslated regions (UTRs) must functionally compensate for the lack of the cap and poly(A) tail. We examined the efficacy of translation initiation in CfMV by comparing it to well-studied viral translational enhancers. Although insertion of the CfMV 5'UTR (CfMVe) into plant expression vectors improved gene expression in barley more than the other translational enhancers examined, studies at the RNA level showed that CfMVe alone or in combination with the CfMV 3'UTR did not provide the RNAs translational advantage. Mutation analysis revealed that translation initiation from CfMVe involved scanning. Interestingly, CfMVe also promoted translation initiation from an intercistronic position of dicistronic mRNAs in vitro. Furthermore, internal initiation occurred with similar efficacy in translation lysates that had reduced concentrations of eukaryotic initiation factor (eIF) 4E, suggesting that initiation was independent of the eIF4E. In contrast, reduced translation in the eIF4G-depleted lysates indicated that translation from internally positioned CfMVe was eIF4G-dependent. After successful translation initiation, leaky scanning brings the ribosomes to the second open reading frame (ORF). The CfMV polyprotein is produced from this and the following overlapping ORF via programmed -1 ribosomal frameshift (-1 PRF). Two signals in the mRNA at the beginning of the overlap program approximately every fifth ribosome to slip one nucleotide backwards and continue translation in the new -1 frame. This leads to the production of C-terminally extended polyprotein, which encodes the viral RNA-dependent RNA polymerase (RdRp). The -1 PRF event in CfMV was very efficient, even though it was programmed by a simple stem-loop structure instead of a pseudoknot, which is usually required for high -1 PRF frequencies. Interestingly, regions surrounding the -1 PRF signals improved the -1 PRF frequencies. Viral protein P27 inhibited the -1 PRF event in vivo, putatively by binding to the -1 PRF site. This suggested that P27 could regulate the occurrence of -1 PRF. Initiation of viral replication requires that viral proteins are released from the polyprotein. This is catalyzed by viral serine protease, which is also encoded from the polyprotein. N-terminal amino acid sequencing of CfMV VPg revealed that the junction of the protease and VPg was cleaved between glutamate (E) and asparagine (N) residues. This suggested that the processing sites used in CfMV differ from the glutamate and serine (S) or threonine (T) sites utilized in other sobemoviruses. However, further analysis revealed that the E/S and E/T sites may be used to cleave out some of the CfMV proteins.

First complete genome sequence of a capsicum chlorosis tospovirus isolate from Australia with an unusually large S RNA intergenic region

The first complete genome sequence of capsicum chlorosis virus (CaCV) from Australia was determined using a combination of Illumina HiSeq RNA and Sanger sequencing technologies. Australian CaCV had a tripartite genome structure like other CaCV isolates. The large (L) RNA was 8913 nucleotides (nt) in length and contained a single open reading frame (ORF) of 8634 nt encoding a predicted RNA-dependent RNA polymerase (RdRp) in the viral-complementary (vc) sense. The medium (M) and small (S) RNA segments were 4846 and 3944 nt in length, respectively, each containing two non-overlapping ORFs in ambisense orientation, separated by intergenic regions (IGR). The M segment contained ORFs encoding the predicted non-structural movement protein (NSm; 927 nt) and precursor of glycoproteins (GP; 3366 nt) in the viral sense (v) and vc strand, respectively, separated by a 449-nt IGR. The S segment coded for the predicted nucleocapsid (N) protein (828 nt) and non-structural suppressor of silencing protein (NSs; 1320 nt) in the vc and v strand, respectively. The S RNA contained an IGR of 1663 nt, being the largest IGR of all CaCV isolates sequenced so far. Comparison of the Australian CaCV genome with complete CaCV genome sequences from other geographic regions showed highest sequence identity with a Taiwanese isolate. Genome sequence comparisons and phylogeny of all available CaCV isolates provided evidence for at least two highly diverged groups of CaCV isolates that may warrant re-classification of AIT-Thailand and CP-China isolates as unique tospoviruses, separate from CaCV.

Nucleotide sequence of the 3' terminal region of belladonna mottle virus-Iowa (renamed Physalis mottle virus) RNA and an analysis of the relationships of tymoviral coat proteins

The 3prime terminal 1255nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3prime terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.