Variation in the genome of rice tungro bacilliform virus: molecular characterization of six isolates

Analysis of codon usage patterns and influencing factors in rice tungro bacilliform virus

Rice tungro bacilliform virus (RTBV) belongs to genus Tungrovirus within the family Caulimoviridae harbors circular double-stranded DNA (dsDNA). Rice tungro disease (RTD) caused by RTBV, responsible for severe rice yield losses in South and Southeast Asia. Here, we performed a systematic evolutionary and codon usage bias (CUB) analysis of RTBV genome sequences. We analysed different bioinformatics techniques to calculate the nucleotide compositions, the relative synonymous codon usage (RSCU), and other indices. The results indicated slightly or low codon usage bias in RTBV isolates. Mutation and natural selection pressures have equally contributed to this low codon usage bias. Additionally, multiple factors such as host, geographical distribution also affect codon usage patterns in RTBV genomes. RSCU analysis revealed that RTBV shows mutation bias and prefers A and U ended codons to code amino acids. Codon usage patterns of RTBV were also found to be influenced by its host. This indicates that RTBV have evolved codon usage patterns that are specific to its host. The findings from this study are expected to increase our understanding of factors leading to viral evolution and fitness with respect to hosts and the environment.

Unusual genomic features of a badnavirus infecting mulberry

Mulberry badnavirus 1 (MBV1) has been characterized as the aetiological agent of a disease observed on a mulberry tree in Lebanon (accession L34). A small RNA next-generation sequencing library was prepared and analysed from L34 extract, and these data together with genome walking experiments have been used to obtain the full-length virus sequence. Uniquely among badnaviruses, the MBV1 sequence encodes a single ORF containing all the conserved pararetrovirus motifs. Two genome sizes (6 kb and 7 kb) were found to be encapsidated in infected plants, the shortest of which shares 98.95 % sequence identity with the full L34 genome. In the less-than-full-length deleted genome, the translational frame for the replication domains was conserved, but the particle morphology, observed under electron microscopy, was somehow altered. Southern blot hybridization confirmed the coexistence of the two genomic forms in the original L34 accession, as well as the absence of cointegration in the plant genome. Both long and deleted genomes were cloned and proved to be infectious in mulberry. Differently from other similar nuclear-replicating viruses or viroids, the characterization of the MBV1-derived small RNAs showed a reduced amount of the 24-mer class size.

Replication of Plant Viruses

Viruses replicate using both their own genetic information and host cell components and machinery. The different genome types have different replication pathways which contain controls on linking the process with translation and movement around the cell as well as not compromising the infected cell. This chapter discusses the replication mechanisms, faults in replication and replication of viruses co-infecting cells.

Viruses replicate using both their own genetic information and host cell components and machinery. The different genome types have different replication pathways which contain controls on linking the process with translation and movement around the cell as well as not compromising the infected cell. This chapter discusses the replication mechanisms, faults in replication and replication of viruses coinfecting cells.

Bacilliform DNA-containing plant viruses in the tropics: commonalities within a genetically diverse group

Plant viruses, possessing a bacilliform shape and containing double-stranded DNA, are emerging as important pathogens in a number of agricultural and horticultural crops in the tropics. They have been reported from a large number of countries in African and Asian continents, as well as from islands from the Pacific region. The viruses, belonging to two genera, Badnavirus and Tungrovirus, within the family Caulimoviridae, have genomes displaying a common plan, yet are highly variable, sometimes even between isolates of the same virus. In this article, we summarize the current knowledge with a view to revealing the common features embedded within the genetic diversity of this group of viruses.

TAXONOMY

Virus; order Unassigned; family Caulimoviridae; genera Badnavirus and Tungrovirus; species Banana streak viruses, Bougainvillea spectabilis chlorotic vein banding virus, Cacao swollen shoot virus, Citrus yellow mosaic badnavirus, Dioscorea bacilliform viruses, Rice tungro bacilliform virus, Sugarcane bacilliform viruses and Taro bacilliform virus.

MICROBIOLOGICAL PROPERTIES

Bacilliform in shape; length, 60-900 nm; width, 35-50 nm; circular double-stranded DNA of approximately 7.5 kbp with one or more single-stranded discontinuities.

HOST RANGE

Each virus generally limited to its own host, including banana, bougainvillea, black pepper, cacao, citrus species, Dioscorea alata, rice, sugarcane and taro.

DISEASE SYMPTOMS

Foliar streaking in banana and sugarcane, swelling of shoots in cacao, yellow mosaic in leaves and stems in citrus, brown spot in the tubers in yam and yellow-orange discoloration and stunting in rice.

USEFUL WEBSITES

http://www.dpvweb.net.

The complete genome sequence of a south Indian isolate of Rice tungro spherical virus reveals evidence of genetic recombination between distinct isolates

In this study, complete genome of a south Indian isolate of Rice tungro spherical virus (RTSV) from Andhra Pradesh (AP) was sequenced, and the predicted amino acid sequence was analysed. The RTSV RNA genome consists of 12,171 nt without the poly(A) tail, encoding a putative typical polyprotein of 3,470 amino acids. Furthermore, cleavage sites and sequence motifs of the polyprotein were predicted. Multiple alignment with other RTSV isolates showed a nucleotide sequence identity of 95% to east Indian isolates and 90% to Philippines isolates. A phylogenetic tree based on complete genome sequence showed that Indian isolates clustered together, while Vt6 and PhilA isolates of Philippines formed two separate clusters. Twelve recombination events were detected in RNA genome of RTSV using the Recombination Detection Program version 3. Recombination analysis suggested significant role of 5' end and central region of genome in virus evolution. Further, AP and Odisha isolates appeared as important RTSV isolates involved in diversification of this virus in India through recombination phenomenon. The new addition of complete genome of first south Indian isolate provided an opportunity to establish the molecular evolution of RTSV through recombination analysis and phylogenetic relationship.

Further support of genetic conservation in Indian isolates of Rice tungro bacilliform virus by sequence analysis of an isolate from North-Western India

The genomic sequence of an isolate of Rice tungro bacilliform virus (RTBV), collected from the state of Punjab (Pb), a non-endemic tungro region from North-Western India was determined. In silico comparison of the 7931-bp sequence with isolates from Southeast Asia and the three previously characterized Indian isolates, revealed not only similar genome size to other Indian isolates but also high degree of homology both at nucleotide (>93 %) and amino acid (>96 %) levels among them. On the other hand, like the other Indian isolates, RTBV-Pb showed much lower nucleotide (<87 %) and amino acid (<90 % in most of the open reading frames) identities with the Southeast Asian isolates owing to several nucleotide substitutions and indels. In-depth annotation comparisons reinforce the hypothesis that Indian isolates of RTBV have diverged sufficiently from the Southeast Asian ones to form a separate group.

The molecular diversity and evolution of Rice tungro bacilliform virus from Indian perspective

Rice tungro disease is caused by a combination of two viruses: Rice tungro spherical virus and Rice tungro bacilliform virus (RTBV). This study was performed with the objective to decipher the molecular variability and evolution of RTBV isolates present in the tungro-affected states of Indian subcontinent. Phylogenetic analysis based on ORF-I, ORF-II, and ORF-IV sequences showed distinct divergence of Indian RTBV isolates into two groups; one consisted isolates from Hyderabad (Andhra Pradesh), Cuttack (Orissa), and Puducherry and another from West Bengal, Chinsura West Bengal, and Kanyakumari (Tamil Nadu). The results obtained from phylogenetic analysis were further supported with the single nucleotide polymorphisms (SNPs), insertion and deletion (INDELs) and evolutionary distance analysis. In addition, sequence difference count matrix revealed a maximum of 56 (ORF-I), 13 (ORF-II) and 73 (ORF-IV) nucleotides differences among all the Indian RTBV isolates taken in this study. However, at the protein level these differences were not significant as revealed by K (a)/K (s) ratio calculation. Sequence identity at nucleotide and amino acid level was 92-100 % (ORF-I), 96-100 % (ORF-II), 94-100 % (ORF-IV) and 86-100 % (ORF-I), 98-100 % (ORF-II) and 95-100 % (ORF-IV), respectively, among Indian isolates of RTBV. The divergence of RTBV isolates into two independent clusters of Indian and non-Indian was shown with the help of the data obtained from phylogeny, SNPs, and INDELs, evolutionary distance analysis, and conserved motifs analysis. The important role of ORF-I and ORF-IV in RTBV diversification and adaptation to different rice growing regions is also discussed.

Genetic variation of coat protein gene among the isolates of Rice tungro spherical virus from tungro-endemic states of the India

Rice tungro disease, one of the major constraints to rice production in South and Southeast Asia, is caused by a combination of two viruses: Rice tungro spherical virus (RTSV) and Rice tungro bacilliform virus (RTBV). The present study was undertaken to determine the genetic variation of RTSV population present in tungro endemic states of Indian subcontinent. Phylogenetic analysis based on coat protein sequences showed distinct divergence of Indian RTSV isolates into two groups; one consisted isolates from Hyderabad (Andhra Pradesh), Cuttack (Orissa), and Puducherry and another from West Bengal, Coimbatore (Tamil Nadu), and Kanyakumari (Tamil Nadu). The results obtained from phylogenetic study were further supported with the SNPs (single nucleotide polymorphism), INDELs (insertion and deletion) and evolutionary distance analysis. In addition, sequence difference count matrix revealed 2-68 nucleotides differences among all the Indian RTSV isolates taken in this study. However, at the protein level these differences were not significant as revealed by Ka/Ks ratio calculation. Sequence identity at nucleotide and amino acid level was 92-100% and 97-100%, respectively, among Indian isolates of RTSV. Understanding of the population structure of RTSV from tungro endemic regions of India would potentially provide insights into the molecular diversification of this virus.

The large intergenic region of Rice tungro bacilliform virus evolved differentially among geographically distinguished isolates

Rice tungro bacilliform virus (RTBV) is a plant pararetrovirus. The large intergenic region (LIGR) of RTBV having a single transcriptional promoter produces more than genome length pregenomic RNA (pgRNA) which directs synthesis of circular double-stranded viral DNA and serves as a polycistronic mRNA. By computer-aided analysis of LIGR, the 11 RTBV isolates sequenced so far were compared with respect to structural organization of promoter and pgRNA 5'-leader. The results revealed only 74.90% identity at LIGR between 'Southeast Asian' (SEA) and 'South Asian' (SA) isolates of RTBV indicating considerable variation between two groups which was also reflected during analysis of promoter and leader sequence. The predicted promoter region of SA isolates exhibited major variations in terms of transcription start site and consensus sequences of cis motifs expecting further exploitation of promoter region of SA isolates. The reduced length of leader sequence along with less numbers and different arrangements of small open reading frames (sORFs) in case of SA isolates might have some alterations in the control of expression of ORF II and III between the two groups. In spite of these variations, the leader sequence of both SEA and SA type isolates showed formation of stable secondary or stem-loop structure having identical features for efficient translation. The conservation of sORF1 at seven nucleotides upstream of stable stem-loop, CU-rich sequence following the sORF1 stop codon and AU-rich shunt landing sequence immediately downstream of the secondary structure suggested conservation of ribosomal shunt mechanism in all RTBV isolates irrespective of their geographical distribution.

Analysis of the complete DNA sequence of rice tungro bacilliform virus from southern India indicates it to be a product of recombination

The complete nucleotide sequence of an isolate of rice tungro bacilliform virus (RTBV), collected from Kanyakumari, India, where RTBV was reported recently for the first time, has been analyzed. Sequence comparison revealed that the RTBV isolate from Kanyakumari (RTBV-KK) has a high degree of identity to the two previously reported RTBV sequences from India, RTBV-AP and RTBV-WB, which had been collected from field locations about 10 years ago and 1000-2000 km away from the collection site of RTBV-KK. Most of the sequence domains reported previously in other RTBV isolates were found to be conserved in RTBV-KK. Closer inspection revealed RTBV-KK to be a possible recombinant between RTBV-AP and RTBV-WB in the genomic region encompassing the coat protein gene.

Phylogenetic analysis of Rice tungro bacilliform virus ORFs revealed strong correlation between evolution and geographical distribution

A new isolate of Rice tungro bacilliform virus (RTBV) was collected from Chinsura, West Bengal, India. The full genome was sequenced and deposited to GenBank designating the new one as Chinsura isolate. The four open reading frames (ORFs) of the new isolate were compared with those of previously reported 'South-east Asian' (SEA) and 'South Asian' (SA) isolates emphasizing the ORF3, which is the largest and functionally most important gene of RTBV. In the ORFs, Chinsura isolate shared 90.0-100.0% identity at amino acid level with SA isolates, but only 58.76-88.63% identity with SEA isolates for the same. Similarly, the amino acid identity of ORFs between SEA and SA isolates ranged from 58.77 to 88.64, whereas within each group the corresponding value was >96.0%. The phylogenetic analysis based on nucleotide and amino acid sequences of each ORF made two broad clusters of SEA- and SA-types including Chinsura isolate within SA cluster. Moreover, the relative positions and length of functional domains corresponding to movement protein (MP), coat protein (CP), aspartate protease (PR) and reverse transcriptase/ribonuclease H (RT/RNase H) of ORF3 of Chinsura isolate were completely identical with SA isolates. The clustering pattern indicated strong influence of geographical habitat on genomic evolution. Comparison of ORF3 among all the isolates revealed major variations at non-functional regions in between the functional domains and at the hypervariable 3'-terminal end of ORF3, while PR appeared to have evolved differentially in SA isolates expecting further characterization.

Transgenic expression of coat protein gene of Rice tungro bacilliform virus in rice reduces the accumulation of viral DNA in inoculated plants

Rice tungro, a devastating disease of rice in south and southeast Asia, is caused by the joint infection of Rice tungro bacilliform virus (RTBV) and Rice tungro spherical virus (RTSV). In order to obtain transgenic resistance against RTBV, indica rice cultivar Pusa Basmati-1 was transformed to express the coat protein (CP) gene of an Indian isolate of RTBV. Rice plants containing the transgene integrated in low copy numbers were obtained, in which the CP was shown to accumulate in the leaf tissue. The progenies representing three independent transformation events were challenged with Indian isolates of RTBV using viruliferous Green leafhoppers, and the viral titers in the inoculated plants were monitored using DNA dot-blot hybridization. As compared to non-transgenic controls, two independent transgenic lines showed significantly low levels of RTBV DNA, especially towards later stages of infection and a concomitant reduction of tungro symptoms.

Dilemmas caused by endogenous pararetroviruses regarding the taxonomy and diagnosis of yam (Dioscorea spp.) badnaviruses: analyses to support safe germplasm movement

The discovery of endogenous pararetroviral sequences (EPRVs) has had a deep impact on the approaches needed for diagnosis, taxonomy, safe movement of germplasm and management of diseases caused by pararetroviruses. In this article, we illustrate this through the example of yam (Dioscorea spp.) badnaviruses. To enable progress, it is first necessary to clarify the taxonomical status of yam badnavirus sequences. Phylogeny and pairwise sequence comparison of 121 yam partial reverse transcriptase sequences provided strong support for the identification of 12 yam badnavirus species, of which ten have not been previously named. Virus prevalence data were obtained, and they support the presence of EPRVs in D. rotundata, but not in D. praehensilis, D. abyssinica, D. alata or D. trifida. Five yam badnavirus species characterised by a wide host range seem to be of African origin. Seven other yam badnavirus species with a limited host range are probably of Asian-Pacific origin. Recombination under natural circumstances appears to be rare. Average values of nucleotide intra-species genetic distances are comparable to data obtained for other RNA and DNA virus families. The dispersion scenarios proposed here, combined with the fact that host-switching events appear common for some yam badnaviruses, suggest that the risks linked to introduction via international plant material exchanges are high.

A variable region of the sugarcane bacilliform virus (SCBV) genome can be used to generate promoters for transgene expression in sugarcane

Four promoters derived from sugarcane bacilliform virus (SCBV) were compared and characterised. Three were obtained by PCR amplification of purified virion DNA extracted from three sugarcane cultivars. The fourth promoter was obtained by subcloning from an almost genome-length clone of SCBV. All promoters were able to drive stable expression of beta-glucuronidase in sugarcane. The PCR-derived promoter sequences shared more DNA homology with banana streak virus than to the subcloned SCBV. The subcloned promoter was the strongest expressing and was able to drive reporter gene expression in vitro and in the leaves, meristems and roots of glasshouse-grown sugarcane. Expression levels were at least equal to or higher than those measured for the maize polyubiquitin promoter.

Detection of molecular variability in rice tungro bacilliform viruses from India using polymerase chain reaction-restriction fragment length polymorphism

Rice tungro bacilliform virus (RTBV) with rice tungro spherical virus (RTSV) causes the destructive tungro disease of rice. In order to ascertain the molecular variability of RTBV in India, primers were designed to amplify a polymorphic DNA fragment of the virus. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis on a number of field isolates indicated mixed infections and molecular heterogeneity in the viral genome.

Early detection of cacao swollen shoot virus using the polymerase chain reaction

A polymerase chain reaction assay was developed which allows early detection of cacao swollen shoot virus (CSSV) in DNA extracts from cacao plantlets agroinoculated with the Togolese isolate Agou 1. The primers used were derived from badnavirus conserved sequences and nucleic acid was extracted with the Plant DNeasy extraction kit (Qiagen). CSSV genome was detectable between 6 and 20 days after inoculation. The first leaf symptoms appeared after 4 weeks and the first shoot swelling symptoms after 8 weeks.

Sequence changes in six variants of rice tungro bacilliform virus and their phylogenetic relationships

The DNA of three biological variants, G1, Ic and G2, which originated from the same greenhouse isolate of rice tungro bacilliform virus (RTBV) at the International Rice Research Institute (IRRI), was cloned and sequenced. Comparison of the sequences revealed small differences in genome sizes. The variants were between 95 and 99% identical at the nucleotide and amino acid levels. Alignment of the three genome sequences with those of three published RTBV sequences (Phi-1, Phi-2 and Phi-3) revealed numerous nucleotide substitutions and some insertions and deletions. The published RTBV sequences originated from the same greenhouse isolate at IRRI 20, 11 and 9 years ago. All open reading frames (ORFs) and known functional domains were conserved across the six variants. The cysteine-rich region of ORF3 showed the greatest variation. When the six DNA sequences from IRRI were compared with that of an isolate from Malaysia (Serdang), similar changes were observed in the cysteine-rich region in addition to other nucleotide substitutions and deletions across the genome. The aligned nucleotide sequences of the IRRI variants and Serdang were used to analyse phylogenetic relationships by the bootstrapped parsimony, distance and maximum-likelihood methods. The isolates clustered in three groups: Serdang alone; Ic and G1; and Phi-1, Phi-2, Phi-3 and G2. The distribution of phylogenetically informative residues in the IRRI sequences shared with the Serdang sequence and the differing tree topologies for segments of the genome suggested that recombination, as well as substitutions and insertions or deletions, has played a role in the evolution of RTBV variants. The significance and implications of these evolutionary forces are discussed in comparison with badnaviruses and caulimoviruses.

Differentiation of rice tungro bacilliform virus strains by restriction analysis and DNA hybridization

Rice tungro bacilliform virus (RTBV) is one of the two viruses that cause tungro disease. Four RTBV strains maintained in the greenhouse for 4 years, G1, G2, Ic, and L, were differentiated by restriction fragment length polymorphism (RFLP) analysis of the native viral DNA. Although strains G1 and Ic had identical restriction patterns when cleaved with Pst1, BamHI, EcoRI, and EcoRV, they can be differentiated from strains G2 and L by EcoRI and EcoRV digestion. These same endonucleases also differentiate strain G2 from strain L. When total DNA extracts from infected plants were used instead of viral DNA, and digested with EcoRV, identical restriction patterns for each strain (G2 and L) were obtained from roots, leaves, and leaf sheaths of infected plants. The restriction patterns were consistent from plant to plant, in different varieties, and at different times after inoculation. This technique can be used to differentiate RTBV strains and determine the variability of a large number of field samples.

Genetic variation of rice tungro bacilliform virus in the Philippines

Restriction fragment length polymorphisms (RFLPs) were found in 27 full genome length clones from a glasshouse isolate of rice tungro bacilliform pararetrovirus (RTBV) from the International Rice Research Institute (IRRI), the Philippines and from clones from 5 field isolates from different parts of the Philippines. There was much less variation between the IRRI clones than between the field isolate clones. The RFLPs were due to single base changes and represented about 10% of the potential sites. Sequencing across the region between nt 7772 and 7989 confirmed that the field isolates differed from the published sequence more than did the IRRI clones. The most common substitutions were G > A, A > G and T > C. Sequence heterogeneity was also noted in PCR products from RTBV DNA from the isolates. These observations are discussed in relation to the quasispecies population concept of viruses.

Sequences of the three coat protein genes of a Malaysian isolate of rice tungro spherical virus reveal a close relationship to the Philippine isolate

Coat protein genes CP1, CP2 and CP3 of an isolate (MaP1) of rice tungro spherical virus (RTSV) from Malaysia were isolated, cloned and sequenced. Comparative analysis indicated that MaP1 isolate is closely related to the Philippine isolate.

Rice tungro spherical virus: nucleotide sequence of the 3' genomic half and studies on the two small 3' open reading frames

Rice tungro spherical virus (RTSV) consists of a single-stranded RNA genome of about 12 kilobases that contains one large open reading frame, ORF 1 and two small ORFs 2 and 3 at its 3' end (Shen et al., 1993, Virology 193:621-630); it was suggested that ORF 2 was expressed via a frameshift. To study the genomic information of RTSV and the variation between different RTSV isolates, the 3' half of a Philippine isolate and parts of a Thai and an Indian isolate were cloned and sequenced. Significant sequence differences were found in ORF 2 and in the 3' non-translated region. Additional stop codons have been revealed in the previously described ORF 2 in several independent clones from the three different virus isolates, the most conserved stop codon in the middle of ORF 2 being confirmed by direct RNA sequencing. These results suggest that ORF 2 could only express a peptide of about 5 kDa instead of 12 kDa as proposed earlier. Polyclonal antisera were raised against ORF 2 and 3 proteins as fusions with glutathione-S-transferase. Using these antisera we failed to detect any virus-specific peptides in extracts from infected rice plants and in virus preparations. The nucleotide sequence of the 3' end of our RTSV isolates contains several small ORFs and does not contain a repeat of 256 nucleotides found in the published sequence. These results indicate that RTSV could contain an unusually long 3' non-coding region of 1240 nucleotides in length.

Molecular biology of rice tungro viruses

Rice tungro, the most important virus disease of rice in South and Southeast Asia, is caused by a complex of two viruses, rice tungro bacilliform virus (RTBV) and rice tungro spherical virus (RTSV). RTBV is a plant pararetrovirus with bacilliform particles, the structure of which is based on T = 3 icosahedral symmetry cut across the threefold axis.The particles encapsidate a circular double-stranded DNA of 8 kbp that encodes four proteins. The current information on the properties, functions, and expression of these proteins is discussed, as is the evidence for replication by reverse transcription. Two major strains of RTBV have been recognized, one from the Indian subcontinent and the other from Southeast Asia. RTSV particles contain a single-stranded RNA genome of 12 kb that encodes a large polyprotein and possibly one or two smaller proteins. The properties and processing of the polyprotein are described and the resemblance to picornaviruses noted.