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

Viruses of Yams (Dioscorea spp.): Current Gaps in Knowledge and Future Research Directions to Improve Disease Management

Viruses are a major constraint for yam production worldwide. They hamper the conservation, movement, and exchange of yam germplasm and are a threat to food security in tropical and subtropical areas of Africa and the Pacific where yam is a staple food and a source of income. However, the biology and impact of yam viruses remains largely unknown. This review summarizes current knowledge on yam viruses and emphasizes gaps that exist in the knowledge of the biology of these viruses, their diagnosis, and their impact on production. It provides essential information to inform the implementation of more effective virus control strategies.

Characterization of yam mosaic viruses from Brazil reveals a new phylogenetic group and possible incursion from the African continent

Yam (Dioscorea spp.) is an important crop for smallholder farmers in the Northeast region of Brazil. Wherever yam is grown, diseases caused by yam mosaic virus (YMV) are prevalent. In the present study, the diversity of YMV infecting Dioscorea cayennensis-rotundata was analyzed. In addition, five species of Dioscorea (D. alata, D. altissima, D. bulbifera, D. subhastata, and D. trifida) commonly found in Brazil were analyzed using ELISA and high-throughput sequencing (HTS). YMV was detected only in D. cayennensis-rotundata, of which 66.7% of the samples tested positive in ELISA. Three YMV genome sequences were assembled from HTS and one by Sanger sequencing to group the sequences in a clade phylogenetically distinct from YMV from other origins. Temporal phylogenetic analyses estimated the mean evolutionary rate for the CP gene of YMV as 1.76 × 10^-3 substitutions per site per year, and the time to the most recent common ancestor as 168.68 years (95% Highest Posterior Density, HPD: 48.56-363.28 years), with a most likely geographic origin in the African continent. The data presented in this study contribute to reveal key aspects of the probable epidemiological history of YMV in Brazil.

Homing in on Endogenous Badnaviral Elements: Development of Multiplex PCR-DGGE for Detection and Rapid Identification of Badnavirus Sequences in Yam Germplasm

Viruses of the genus Badnavirus (family Caulimoviridae) are double-stranded DNA-reverse transcribing (dsDNA-RT) plant viruses and have emerged as serious pathogens of tropical and temperate crops globally. Endogenous badnaviral sequences are found integrated in the genomes of several economically important plant species. Infection due to activation of replication-competent integrated copies of the genera Badnavirus, Petuvirus and Cavemovirus has been described. Such endogenous badnaviral elements pose challenges to the development of nucleic acid-based diagnostic methods for episomal virus infections and decisions on health certification for international movement of germplasm and seed. One major food security crop affected is yam (Dioscorea spp.). A diverse range of Dioscorea bacilliform viruses (DBVs), and endogenous DBV (eDBV) sequences have been found to be widespread in yams cultivated in West Africa and other parts of the world. This study outlines the development of multiplex PCR-dependent denaturing gradient gel electrophoresis (PCR-DGGE) to assist in the detection and analysis of eDBVs, through the example of analysing yam germplasm from Nigeria and Ghana. Primers targeting the three most prevalent DBV monophyletic species groups in West Africa were designed to improve DGGE resolution of complex eDBV sequence fingerprints. Multiplex PCR-DGGE with the addition of a tailor-made DGGE sequence marker enables rapid comparison of endogenous badnaviral sequence diversity across germplasm, as illustrated in this study for eDBV diversity in yam.

A Review of Viruses Infecting Yam (Dioscorea spp.)

Yam is an important food staple for millions of people globally, particularly those in the developing countries of West Africa and the Pacific Islands. To sustain the growing population, yam production must be increased amidst the many biotic and abiotic stresses. Plant viruses are among the most detrimental of plant pathogens and have caused great losses of crop yield and quality, including those of yam. Knowledge and understanding of virus biology and ecology are important for the development of diagnostic tools and disease management strategies to combat the spread of yam-infecting viruses. This review aims to highlight current knowledge on key yam-infecting viruses by examining their characteristics, genetic diversity, disease symptoms, diagnostics, and elimination to provide a synopsis for consideration in developing diagnostic strategy and disease management for yam.

Tissue culture and next-generation sequencing: A combined approach for detecting yam (Dioscorea spp.) viruses

In vitro culture offers many advantages for yam germplasm conservation, propagation and international distribution. However, low virus titres in the generated tissues pose a challenge for reliable virus detection, which makes it difficult to ensure that planting material is virus-free. In this study, we evaluated next-generation sequencing (NGS) for virus detection following yam propagation using a robust tissue culture methodology. We detected and assembled the genomes of novel isolates of already characterised viral species of the genera Badnavirus and Potyvirus, confirming the utility of NGS in diagnosing yam viruses and contributing towards the safe distribution of germplasm.

Characterization of a novel member of the family Caulimoviridae infecting Dioscorea nummularia in the Pacific, which may represent a new genus of dsDNA plant viruses

We have characterized the complete genome of a novel circular double-stranded DNA virus, tentatively named Dioscorea nummularia-associated virus (DNUaV), infecting Dioscorea nummularia originating from Samoa. The genome of DNUaV comprised 8139 bp and contained four putative open reading frames (ORFs). ORFs 1 and 2 had no identifiable conserved domains, while ORF 3 had conserved motifs typical of viruses within the family Caulimoviridae including coat protein, movement protein, aspartic protease, reverse transcriptase and ribonuclease H. A transactivator domain, similar to that present in members of several caulimoviridae genera, was also identified in the putative ORF 4. The genome size, organization, and presence of conserved amino acid domains are similar to other viruses in the family Caulimoviridae. However, based on nucleotide sequence similarity and phylogenetic analysis, DNUaV appears to be a distinct novel member of the family and may represent a new genus.

Rolling Circle Amplification to Screen Yam Germplasm for Badnavirus Infections and to Amplify and Characterise Novel Badnavirus Genomes

Since the first discovery of badnaviruses (family Caulimoviridae, genus Badnavirus) in yam (Dioscorea spp.) germplasm in the 1970s (Harrison and Roberts, 1973), several hundred partial badnavirus reverse transcriptase (RT)-ribonuclease H (RNaseH) sequences have been characterised ( Kenyon et al., 2008 ; Bousalem et al., 2009 ), but only a few complete Dioscorea bacilliform virus (DBV) genome sequences have been reported ( Phillips et al., 1999 ; Seal and Muller, 2007; Bömer et al., 2016 and 2017; Sukal et al., 2017 ; Umber et al., 2017 ). We have optimised a workflow involving total nucleic acid extractions and rolling circle amplification (RCA) combined with restriction enzyme analysis for the detection and amplification of DBVs present in yam germplasm. We have employed this approach successfully revealing three novel episomal yam badnaviruses ( Bömer et al., 2016 ). We proposed this to be a complementary method to denaturing gradient gel electrophoresis, which enables a rapid indication of badnavirus diversity as well as the identification of potentially integrated badnavirus sequences in the host genome ( Turaki et al., 2017 ). Here, we describe the step-by-step protocol to screen yam germplasm for badnavirus infections using RCA as an efficient research tool in the amplification and characterization of novel badnavirus genomes.

Complete genome sequence of a new member of the genus Badnavirus, Dioscorea bacilliform RT virus 3, reveals the first evidence of recombination in yam badnaviruses

Yams (Dioscorea spp.) host a diverse range of badnaviruses (genus Badnavirus, family Caulimoviridae). The first complete genome sequence of Dioscorea bacilliform RT virus 3 (DBRTV3), which belongs to the monophyletic species group K5, is described. This virus is most closely related to Dioscorea bacilliform SN virus (DBSNV, group K4) based on a comparison of genome sequences. Recombination analysis identified a unique recombination event in DBRTV3, with DBSNV likely to be the major parent and Dioscorea bacilliform AL virus (DBALV) the minor parent, providing the first evidence for recombination in yam badnaviruses. This has important implications for yam breeding programmes globally.

PCR-DGGE Analysis: Unravelling Complex Mixtures of Badnavirus Sequences Present in Yam Germplasm

Badnaviruses (family Caulimoviridae, genus Badnavirus) have emerged as serious pathogens especially affecting the cultivation of tropical crops. Badnavirus sequences can be integrated in host genomes, complicating the detection of episomal infections and the assessment of viral genetic diversity in samples containing a complex mixture of sequences. Yam (Dioscorea spp.) plants are hosts to a diverse range of badnavirus species, and recent findings have suggested that mixed infections occur frequently in West African yam germplasm. Historically, the determination of the diversity of badnaviruses present in yam breeding lines has been achieved by cloning and sequencing of polymerase chain reaction (PCR) products. In this study, the molecular diversity of partial reverse transcriptase (RT)-ribonuclease H (RNaseH) sequences from yam badnaviruses was analysed using PCR-dependent denaturing gradient gel electrophoresis (PCR-DGGE). This resulted in the identification of complex ‘fingerprints’ composed of multiple sequences of Dioscorea bacilliform viruses (DBVs). Many of these sequences show high nucleotide identities to endogenous DBV (eDBV) sequences deposited in GenBank, and fall into six monophyletic species groups. Our findings highlight PCR-DGGE as a powerful tool in badnavirus diversity studies enabling a rapid indication of sequence diversity as well as potential candidate integrated sequences revealed by their conserved nature across germplasm.

Characterization of badnaviruses infecting Dioscorea spp. in the Pacific reveals two putative novel species and the first report of dioscorea bacilliform RT virus 2

The complete genome sequences of three new badnaviruses associated with yam (Dioscorea spp.) originating from Fiji, Papua New Guinea and Samoa were determined following rolling circle amplification of the virus genomes. The full-length genomes consisted of a single molecule of circular double-stranded DNA of 8106bp for isolate FJ14, 7871bp for isolate PNG10 and 7426bp for isolate SAM01. FJ14 and PNG10 contained three open reading frames while SAM01 had an additional open reading frame which partially overlapped the 3' end of ORF 3. Amino acid sequence analysis of ORF 3 from the three isolates confirmed the presence of conserved motifs typical of other badnaviruses. Phylogenetic analysis revealed the sequences to be closely related to other Dioscorea-infecting badnaviruses. FJ14 and PNG10 appear to be new species, which we have tentatively named dioscorea bacilliform ES virus (DBESV) and dioscorea bacilliform AL virus 2 (DBALV2), respectively, while SAM01 represents a Pacific isolate of the recently published dioscorea bacilliform RT virus 2 and is described as dioscorea bacilliform RT virus 2-[4RT] (DBRTV2-[4RT]).

Characterization of Canna yellow mottle virus in a New Host, Alpinia purpurata, in Hawaii

Canna yellow mottle virus (CaYMV) is an important badnavirus infecting Canna spp. worldwide. This is the first report of CaYMV in flowering ginger (Alpinia purpurata) in Hawaii, where it is associated with yellow mottling and necrosis of leaves, vein streaking, and stunted plants. We have sequenced CaYMV in A. purpurata (CaYMV-Ap) using a combination of next-generation sequencing and traditional Sanger sequencing techniques. The complete genome of CaYMV-Ap was 7,120 bp with an organization typical of other Badnavirus species. Our results indicated that CaYMV-Ap was present in the episomal form in infected flowering ginger. We determined that this virus disease is prevalent in Hawaii and could potentially have significant economic impact on the marketing of A. purpurata as cut flowers. There is a potential concern that the host range of CaYMV-Ap may expand to include other important tropical plants.

The genome sequence of Dioscorea bacilliform TR virus, a member of the genus Badnavirus infecting Dioscorea spp., sheds light on the possible function of endogenous Dioscorea bacilliform viruses

The complete genome sequence of Dioscorea bacilliform TR virus (DBTRV) was determined. The closest relatives of DBTRV are Dioscorea bacilliform AL virus (DBALV) and Dioscorea bacilliform RT virus 1 (DBRTV1). Specific primers were designed and used to determine the prevalence of DBTRV in a yam germplasm collection. It was found that this virus infects Dioscorea alata and D. trifida plants in Guadeloupe and French Guyana. DTRBV was not detected in any of the tested D. cayenensis-rotundata accessions. In silico analysis provided evidence for the presence of DBTRV-like endogenous sequences in the genome of D. cayenensis-rotundata, pointing to a possible role of these sequences in antiviral defense.

A Sequence-Independent Strategy for Amplification and Characterisation of Episomal Badnavirus Sequences Reveals Three Previously Uncharacterised Yam Badnaviruses

Yam (Dioscorea spp.) plants are potentially hosts to a diverse range of badnavirus species (genus Badnavirus, family Caulimoviridae), but their detection is complicated by the existence of integrated badnavirus sequences in some yam genomes. To date, only two badnavirus genomes have been characterised, namely, Dioscorea bacilliform AL virus (DBALV) and Dioscorea bacilliform SN virus (DBSNV). A further 10 tentative species in yam have been described based on their partial reverse transcriptase (RT)-ribonuclease H (RNaseH) sequences, generically referred to here as Dioscorea bacilliform viruses (DBVs). Further characterisation of DBV species is necessary to determine which represent episomal viruses and which are only present as integrated badnavirus sequences in some yam genomes. In this study, a sequence-independent multiply-primed rolling circle amplification (RCA) method was evaluated for selective amplification of episomal DBV genomes. This resulted in the identification and characterisation of nine complete genomic sequences (7.4-7.7 kbp) of existing and previously undescribed DBV phylogenetic groups from Dioscorea alata and Dioscorea rotundata accessions. These new yam badnavirus genomes expand our understanding of the diversity and genomic organisation of DBVs, and assist the development of improved diagnostic tools. Our findings also suggest that mixed badnavirus infections occur relatively often in West African yam germplasm.

Badnaviruses: The Current Global Scenario

Badnaviruses (Family: Caulimoviridae; Genus: Badnavirus) are non-enveloped bacilliform DNA viruses with a monopartite genome containing about 7.2 to 9.2 kb of dsDNA with three to seven open reading frames. They are transmitted by mealybugs and a few species by aphids in a semi-persistent manner. They are one of the most important plant virus groups and have emerged as serious pathogens affecting the cultivation of several horticultural crops in the tropics, especially banana, black pepper, cocoa, citrus, sugarcane, taro, and yam. Some badnaviruses are also known as endogenous viruses integrated into their host genomes and a few such endogenous viruses can be awakened, e.g., through abiotic stress, giving rise to infective episomal forms. The presence of endogenous badnaviruses poses a new challenge for the fool-proof diagnosis, taxonomy, and management of the diseases. The present review aims to highlight emerging disease problems, virus characteristics, transmission, and diagnosis of badnaviruses.

Sequence diversity among badnavirus isolates infecting black pepper and related species in India

The badnavirus, piper yellow mottle virus (PYMoV) is known to infect black pepper (Piper nigrum), betelvine (P. betle) and Indian long pepper (P. longum) in India and other parts of the world. Occurrence of PYMoV or other badnaviruses in other species of Piper and its variability is not reported so far. We have analysed sequence variability in the conserved putative reverse transcriptase (RT)/ribonuclease H (RNase H) coding region of the virus using specific badnavirus primers from 13 virus isolates of black pepper collected from different cultivars and regions and one isolate each from 23 other species of Piper. Of these, four species failed to produce expected amplicon while amplicon from four other species showed more similarities to plant sequences than to badnaviruses. Of the remaining, isolates from black pepper, P. argyrophyllum, P. attenuatum, P. barberi, P. betle, P. colubrinum, P. galeatum, P. longum, P. ornatum, P. sarmentosum and P. trichostachyon showed an identity of >85 % at the nucleotide and >90 % at the amino acid level with PYMoV indicating that they are isolates of PYMoV. On the other hand high sequence variability of 21-43 % at nucleotide and 17-46 % at amino acid level compared to PYMoV was found among isolates infecting P. bababudani, P. chaba, P. peepuloides, P. mullesua and P. thomsonii suggesting the presence of new badnaviruses. Phylogenetic analyses showed close clustering of all PYMoV isolates that were well separated from other known badnaviruses. This is the first report of occurrence of PYMoV in eight Piper spp and likely occurrence of four new species in five Piper spp.

The genome of African yam (Dioscorea cayenensis-rotundata complex) hosts endogenous sequences from four distinct Badnavirus species

Several endogenous viral elements (EVEs) have been identified in plant genomes, including endogenous pararetroviruses (EPRVs). Here, we report the first characterization of EPRV sequences in the genome of African yam of the Dioscorea cayenensis-rotundata complex. We propose that these sequences should be termed 'endogenous Dioscorea bacilliform viruses' (eDBVs). Molecular characterization of eDBVs shows that they constitute sequences originating from various parts of badnavirus genomes, resulting in a mosaic structure that is typical of most EPRVs characterized to date. Using complementary molecular approaches, we show that eDBVs belong to at least four distinct Badnavirus species, indicating multiple, independent, endogenization events. Phylogenetic analyses of eDBVs support and enrich the current taxonomy of yam badnaviruses and lead to the characterization of a new Badnavirus species in yam. The impact of eDBVs on diagnosis, yam germplasm conservation and movement, and breeding is discussed.

The prevalence of badnaviruses in West African yams (Dioscorea cayenensis-rotundata) and evidence of endogenous pararetrovirus sequences in their genomes

Yam (Dioscorea spp.) is an important vegetatively-propagated staple crop in West Africa. Viruses are pervasive in yam worldwide, decreasing growth and yield, as well as hindering the international movement of germplasm. Badnaviruses have been reported to be the most prevalent in yam, and genomes of some other badnaviruses are known to be integrated in their host plant species. However, it was not clear if a similar scenario occurs in Dioscorea yam. This study was conducted to verify the prevalence of badnaviruses, and determine if badnavirus genomes are integrated in the yam genome. Leaf samples (n=58) representing eight species of yam from global yam collections kept at CIRAD, France, and 127 samples of D. rotundata breeding lines (n=112) and landraces (n=15) at IITA, Nigeria, were screened using generic badnavirus PCR primers. Positive amplification of an expected ca. 579bp fragment, corresponding to a partial RT-RNaseH region, was detected in 47 (81%) of 58 samples analysed from CIRAD collections, and 100% of the 127 IITA D. rotundata samples. All the D. cayenensis and D. rotundata samples from the CIRAD and IITA collections tested PCR-positive, and sequencing of a selection of the PCR products confirmed they were typical of the genus Badnavirus. A comparison of serological and nucleic acid techniques was used to investigate whether the PCR-positives were sequences amplified from badnavirus particles or putative endogenous badnavirus sequences in the yam genome. Protein A sandwich-enzyme-linked immunosorbent assay (PAS-ELISA) with badnavirus polyclonal antisera detected cross-reacting viral particles in only 60% (92 of 153) of the CIRAD collection samples analysed, in contrast to the aforementioned 81% by PCR. Immunosorbent electron microscopy (ISEM) of virus preparations of a select set of 16 samples, representing different combinations of positive and negative PCR and PAS-ELISA results, identified bacilliform particles in 11 of these samples. Three PCR-positive yam samples from Burkina Faso (cv. Pilimpikou) were identified in which no viral particles were detected by either PAS-ELISA or ISEM. Southern hybridisation results using a yam badnavirus RT-RNaseH sequence (Gn155Dr) as probe, supported a lack of badnavirus particles in the cv. Pilimpikou and identified their equivalent sequences to be of plant genome origin. Probe Gn155Dr, however, hybridised to viral particles and plant genomic DNA in three D. rotundata samples from Guinea. These results represent the first data demonstrating the presence of integrated sequences of badnaviruses in yam. The implications of this for virus-indexing, breeding and multiplication of seed yams are discussed.

Comparative analysis of endogenous plant pararetroviruses in cultivated and wild Dahlia spp

Two distinct caulimoviruses, Dahlia mosaic virus (DMV) and Dahlia common mosaic virus, and an endogenous plant pararetroviral sequence (DvEPRS) were reported in Dahlia spp. DvEPRS, previously referred to as DMV-D10, was originally identified in the US from the cultivated Dahlia variabilis, and has also been found in New Zealand, Lithuania and Egypt, as well as in wild dahlia species growing in their natural habitats in Mexico. Sequence analysis of three new EPRSs from cultivated dahlias from Lithuania [D10-LT; 7,159 nucleotide level (nt)], New Zealand (D10-NZ, 7,156 nt), and the wild species, Dahlia rupicola, from Mexico (D10-DR, 7,133 nt) is reported in this study. The three EPRSs have the structure and organization typical of a caulimovirus species and showed identities among various open reading frames (ORFs) ranging between 71 and 97 % at the nt when compared to those or the known DvEPRS from the US. Examination of a dataset of seven full-length EPRSs obtained to date from cultivated and wild Dahlia spp. provided clues into genetic diversity of these EPRSs from diverse sources of dahlia. Phylogenetic analyses, mutation frequencies, potential recombination events, selection, and fitness were evaluated as evolutionary evidences for genetic variation. Assessment of all ORFs using phylogenomic and population genetics approaches suggests a wide genetic diversity of EPRSs occurring in dahlias. Phylogenetic analyses show that the EPRSs from various sources form one clade indicating a lack of clustering by geographical origin. Grouping of various EPRSs into two host taxa (cultivated vs. wild) shows little divergence with respect to their origin. Population genetic parameters demonstrate negative selection for all ORFs, with the reverse transcriptase region more variable than other ORFs. Recombination events were found which provide evolutionary evidence for genetic diversity among dahlia-associated EPRSs. This study contributes to an increased understanding of molecular population genetics and evolutionary pathways of these reverse transcribing viral elements.

The common evolutionary history of badnaviruses and banana

Recent plant genome sequencing efforts have revealed myriad viral sequences suggesting a cryptic interaction between both partners. Interestingly, no integration step has ever been reported as an obligatory step in the life cycle of plant viruses. Circular dsDNA viruses belonging to the family Caulimoviridae are the most abundant among integrated plant viral sequences. In this review, we describe how this hitherto hidden interaction could inform the evolutionary history of both partners badnaviruses and banana plants.