Next Article in Journal
Chromosomally Integrated Human Herpesvirus 6: Models of Viral Genome Release from the Telomere and Impacts on Human Health
Previous Article in Journal
Variability of Emaravirus Species Associated with Sterility Mosaic Disease of Pigeonpea in India Provides Evidence of Segment Reassortment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Natural Resources Institute, University of Greenwich, Central Avenue, Chatham Maritime, Kent ME4 4TB, UK
2
Kebbi State University of Science and Technology Aliero, Sokoto Road, PMB 1144 Birnin Kebbi, Nigeria
3
International Institute of Tropical Agriculture (IITA), Oyo Road, PMB 5320 Ibadan, Nigeria
*
Author to whom correspondence should be addressed.
These authors contributed equally to the work.
Submission received: 28 April 2017 / Revised: 29 June 2017 / Accepted: 4 July 2017 / Published: 11 July 2017
(This article belongs to the Section Viruses of Plants, Fungi and Protozoa)

Abstract

:
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.

Graphical Abstract

1. Introduction

Badnaviruses are plant pararetroviruses (family Caulimoviridae, genus Badnavirus) infecting a broad range of economically important crop plants and have emerged as serious pathogens especially affecting the cultivation of tropical crops, such as banana, black pepper, cacao, citrus, sugarcane, taro and yam [1]. The development of reliable diagnostic tests for badnaviruses is challenging due to high serological and genetic heterogeneity, as experienced in a wide range of crops [2,3,4,5]. The discovery of integrated badnavirus sequences in some host plant genomes of agricultural and horticultural crops complicates the fool-proof diagnosis of episomal infections using nucleic acid-based methods as experienced for banana streak viruses (BSVs) in Musa spp. [6,7,8], and their existence poses further challenges for taxonomy, safe movement of germplasm, and disease management (reviewed by [9]).
Integrated sequences derived from representatives of four out of eight recognized genera of the Caulimoviridae family have been found in 27 species from nine different plant families according to Geering et al. [10]. Termed endogenous pararetroviruses (EPRVs) [11,12,13], integration events have been most studied in banana [14], petunia [15], tobacco [16], rice [17], potato [18] and tomato [19], and the majority of these integrated sequences described so far are fragmented and rearranged [7,20,21,22]. Few of these EPRVs can be activated to give rise to infective episomal forms initiating systemic badnavirus infections de novo [20,22,23]. Four species of the genus Badnavirus described in banana genomes of Musa balbisiana species, namely Banana streak OL virus (BSOLV), Banana streak IM virus (BSIMV), Banana streak MY virus (BSMYV), and Banana streak GF virus (BSGFV) [21,24,25,26], as well as the petuvirus Petunia vein clearing virus (PVCV) in petunia [22], and the solendovirus Tobacco vein-clearing virus (TVCV) in tobacco [16,23] are the only activatable EPRVs discovered to date. Activation is considered in banana to be triggered by tissue culture, hybridization, or temperature differences in newly created banana interspecific hybrids [27,28,29], and a homologous recombination-based model is proposed to explain the release of the BSOLV genome from its endogenous counterpart [30].
Increasing numbers of integrated viral sequences are now identified in next-generation sequencing (NGS) data. This for example includes the recent discovery of the endogenous florendoviruses (family Caulimoviridae, genus Florendovirus) as major components of a large diversity of flowering plant genomes, such as apple, cacao, cassava, citrus, grape, maize, papaya, potato, rice and tomato among others [10]. The analysis of these endogenous sequences is substantially more difficult in plant species still lacking a good quality reference genome as is the case for some staple crops, in particular tropical root and tuber crops. The nature of EPRV sequences can be of concern to breeding programmes and multiplication processes with the goal of distributing large amounts of clean vegetative propagation material. Germplasm-containing activatable EPRVs would need to be removed from such programmes prior to multiplication. In this study, we illustrate the value of PCR-dependent denaturing gradient gel electrophoresis (PCR-DGGE) analysis in unravelling complex mixtures of badnavirus sequences through the example of yam germplasm.
Yam (Dioscorea spp.) is one of the most important staple food crops, playing a major role in food security in West Africa [31,32]. Dioscorea bacilliform viruses (DBVs) (family Caulimoviridae, genus Badnavirus) contain several members that are a concern to the safe movement of the germplasm of this vegetatively-propagated crop [3,33,34]. Sequences representative of the genus Badnavirus have been shown to be prevalent in yam materials tested from all yam growing areas of the world, and are highly heterogeneous [3,34,35,36,37,38]. Integrated badnavirus sequences have been discovered in the genomes of Guinea yam (Dioscorea cayenensis-rotundata complex), termed endogenous DBV (eDBV) sequences [38,39]. The rapid increase in the discovery of new plant pararetroviruses and their integrated sequences over recent years has been the result of improved PCR diagnostics and an increase in research attention [3,21,40,41,42]. Understanding the complexity of badnavirus sequences in plant germplasm is of great importance to virologists and plant breeders, in particular when working with vegetatively-propagated crops with an urgent need for a sustainable supply of virus-free planting material. A pressing need remains in the development of diagnostic tools to differentiate the integrated sequences from those representing episomal particles, such that the potential distribution of virus infections can be assessed and their impact on crop yields determined.
Denaturing gradient gel electrophoresis (DGGE) is a commonly used molecular technique for rapid fingerprint analysis that was first described by Fischer and Lerman [43]. DGGE is capable of separating double-stranded PCR products of similar length but differing sequence composition (reviewed by [44]). The principle of the technique is based on the physicochemical fundamental of DNA base-pairing and the differing mobility of double-stranded and partially denatured DNA when analysed by polyacrylamide gel electrophoresis [44,45]. The addition of a GC clamp to the 5′ end of the primer sequence prevents complete denaturation of DNA molecules during DGGE analysis. The technique detects variation of small DNA fragments (~200–700 bp) that differ by as little as a single base substitution [46,47,48]. DGGE has been used extensively for diversity studies such as microbial biodiversity [45,49,50,51], fungal communities [52], genomes of viral strains [53,54,55], forensic application [56] and plant genome polymorphisms [57,58].
The aim of this study was to determine: (1) whether PCR-DGGE could be used as a rapid technique to screen for differences in complex mixtures of episomal and endogenous badnavirus sequences being present in West African yam landraces and breeding lines; and (2) if this approach would enable the simple identification of integrated badnavirus sequences conserved across germplasm. Historically, this was achieved by cloning and sequencing of PCR products using the generic badnavirus primer pair Badna-forward primer/reverse primer (FP/-RP) [59]. This approach however has the disadvantage of there being a possibility of failing to identify all sequence diversity due to clones selected not being representative of the total diversity present within a tested plant. Restriction fragment length polymorphism (RFLP) analysis has also been used to study yam badnavirus diversity [3,35,60] but this technique shows a lower resolution of diversity than DGGE. However, neither technique (cloning of PCR products and RFLP) allows episomal and endogenous badnavirus sequences to be distinguished.
This study is the first to demonstrate the value of PCR-DGGE (hereafter referred to as DGGE) for unravelling complex sequence mixtures present in badnavirus PCR products amplified from yam breeding lines (Figure 1). This approach resulted in the generation of 114 nucleotide sequences from yam germplasm which have been deposited in the GenBank database under accession numbers KY555456 to KY555569. These DBV sequences fall into six monophyletic species groups and represent several unique DGGE profiles. Interestingly, many bands are conserved across germplasm indicating integrated sequences, and many sequences identified in this study do show high nucleotide identities to eDBV sequences. It is considered that the ability to link DGGE profiles to badnavirus sequences will assist the future identification of badnavirus diversity based on comparing sample profiles to standard profile markers.

2. Materials and Methods

2.1. Plant Material

Yam leaf samples of breeding lines and landraces (47 samples in total) of Dioscorea alata (n = 11), Dioscorea bulbifera (n = 1), Dioscorea cayenensis (n = 1), Dioscorea dumetorum (n = 3), Dioscorea esculenta (n = 2) and Dioscorea rotundata (n = 23) used in this study were collected from screen-houses at the International Institute of Tropical Agriculture (IITA, Ibadan, Nigeria) and placed in polythene bags (14 cm × 14 cm). The breeding lines were a subset of the IITA collection samples (n = 127) consisting of first filial (F1) generations of West African breeding lines of D. rotundata (n = 112) and D. rotundata landraces (n = 15). Details of samples from these collections are given in the Supplementary materials of Seal et al. [38]. Selected yam breeding lines and landraces of D. rotundata (n = 6) were provided by the IITA and tubers were grown in a quarantine aphid-proof glasshouse at the Natural Resources Institute (NRI, Chatham Maritime, UK), as described by Mumford and Seal [61]. Individual leaf samples were collected from each plant in small polythene bags (10 cm × 15 cm) and processed immediately. Dioscorea rotundata accession (TDr) 89/02475A and B, and TDr 1892A and B are clones of the same yam accessions.

2.2. Total Nucleic Acid Extraction from Yam Leaves and PCR Amplification of Badnavirus Sequences

Total nucleic acids were extracted from fresh yam leaf tissue (~100 mg) using a modified cetyltrimethylammonium bromide (CTAB) method as described by Kenyon et al. [3]. Total nucleic acids were screened for the presence of sequences typical of the genus Badnavirus by PCR using the degenerate primer set Badna-forward primer (FP) and Badna-reverse primer (RP) designed by Yang et al. [59]. These Badna-FP/-RP primers amplify a 579-bp region (528 bp excluding primer sequences and representing only complete amino acids) of the reverse transcriptase (RT)-ribonuclease H (RNaseH) gene used for taxonomic assessment of badnaviruses [62]. To confirm the suitability of DNA for PCR amplification all DNA samples were first screened using primers targeting the yam actin gene as described by Silva et al. [63]. PCR amplifications were set up in 25-µL reactions containing 1 µL of template (20 ng), 0.5 µM of each primer, 0.25 mM of each deoxynucleotide triphosphate (dNTP), 1 U DreamTaq DNA polymerase and 1× DreamTaq Green buffer (Thermo Scientific, Loughborough, UK) containing 2 mM MgCl2. The cycle conditions for PCR amplification were 94 °C for 2 min, followed by 40 cycles of 94 °C for 20 s, 55 °C for 30 s, 72 °C for 1 min and a final extension of 72 °C for 10 min. PCR products were analysed by gel electrophoresis through 1.5% (w/v) agarose gels including 1× RedSafe nucleic acid stain (iNtRON Biotechnology, Gyeonggi-do, Korea) in 0.5× Tris-Boric acid-ethylenediaminetetra acetic acid (EDTA) (TBE) buffer. PCR products of the expected size were purified and diluted PCR products (‘nested’) or diluted total nucleic acid extractions (‘direct’) were used as template for the purposes of generating PCR products designated for DGGE analysis. Both Badna-FP and Badna-RP primers were modified by the addition of a GC clamp (5′ CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GAT GCC ITT YGG IIT IAA RAA YGC ICC 3′ and 5′ CGC CCG CCG CGC GCG GCG GGC GGG GCG GGG GCA CGG GGG GCC AYT TRC AIA CIS CIC CCC AIC C 3′, respectively) generating a product of 619 bp. PCR reactions using the Badna-FP/-RP GC-clamp primers used cycle conditions as follows: one cycle at 95 °C for 5 min, then 35 cycles of 94 °C for 1 min, 55 °C for 45 s, 72 °C for 2 min followed by one cycle of extension at 72 °C for 10 min. Prior to DGGE analysis, PCR products were confirmed to be of the correct size by agarose gel electrophoresis. All sequencing in this study was performed by Source BioScience (Nottingham, UK). All primers described were synthesised using Sigma oligo service (Sigma-Aldrich, Irvine, UK) and reSource (Source BioScience, Nottingham, UK), or Sure Clean kit (Bioline, London, UK) purification kits were used to clean PCR products prior to sequencing or cloning.

2.3. Denaturing Gradient Gel Electrophoresis (DGGE)

DGGE was performed using the INGENYphorU-2×2 apparatus (INGENY, Goes, The Netherlands) according to the manufacturers’ instructions and following procedural comments provided in the protocol by Green et al. [44]. Gradient gels containing 6.5% (v/v) polyacrylamide (37.5:1 ratio of acrylamide:bis-acrylamide) (National Diagnostics, Atlanta, USA) were formed using a peristaltic pump (Rietschle Thomas, Schopfheim, Germany) and a gradient maker device (INGENY) with denaturing gradients from 35 to 50% (top to bottom) unless stated otherwise (where 100% is 7 M urea and 40% (v/v) deionized formamide) in 1× Tris-acetate-EDTA (TAE) electrophoresis buffer. Samples (20 µL) were loaded on a stacking gel. Electrophoresis was performed at 80 V at a temperature of 60 °C for 18 h. Gels were stained with 1× SYBR Gold nucleic acid gel stain (Invitrogen, Life Technologies, Paisley, UK) in 1× TAE for 30 min at room temperature and destained in deionized water. Gels were placed on a UV transilluminator (G-box Chemi HR16, Syngene, Cambridge, UK) and visualized. Bands of interest were excised from DGGE gels using a sterile scalpel and DNA eluted by soaking in 100 µL of molecular grade water (Sigma) at 4 °C overnight. Aliquots were diluted 1:10 and re-amplified by PCR using the Badna-FP/-RP primer pair followed by PCR purification. Purified PCR products were cloned prior to sequencing using the pGEM-T Easy vector system (Promega, Southampton, UK) according to the manufacturers’ instructions and standard sequencing primers SP6 and T7. To obtain a consensus sequence and control for cross-contamination, two clones from each excised DGGE band showing a different migration pattern were sequenced.

2.4. Sequence Analysis and Phylogeny

Yam badnavirus partial reverse transcriptase-ribonuclease H (RT-RNaseH) nucleotide sequences generated from plasmid clones were analysed using MEGA version 6.0 [64]. The Badna-FP/-RP and vector sequences were removed and the edited sequences were used for similarity basic local alignment search tool (BLAST) searches in the National Centre for Biotechnology Information (NCBI) GenBank databases [65]. Multiple alignments of the partial RT-RNaseH sequences were performed using the CLUSTALW default settings in MEGA version 6.0, where phylogenetic trees were created using the maximum likelihood method with the Kimura 2-parameter model [66]. The robustness of trees was determined by generating bootstrap consensus trees using 1000 replicates. A nucleotide percent similarity matrix (Table S1) was generated using Multiple Alignment using Fast Fourier Transform [67,68]. Protein sequences were aligned using CLUSTAL OMEGA [69,70] and further processed in BioEdit version 7.2.5 [71]. According to the International Committee on Taxonomy of Viruses (ICTV), sequences of the genus Badnavirus differing in their partial RT-RNaseH coding region by more than 20% meet the species demarcation criteria [62]. Eighty-nine yam badnavirus partial RT-RNaseH sequences (see Figure 4 for accession numbers) and the following virus sequences were obtained from the GenBank and used for comparative analyses: Banana streak OL virus (BSOLV, AJ002234); Cacao swollen shoot virus (CSSV, AJ781003); Commelina yellow mottle virus (ComYMV, NC001343); Rice tungro bacilliform virus (RTBV, X57924); Sugarcane bacilliform MO virus (SCBMOV, M89923); and Taro bacilliform virus (TaBV, AF357836). For consistency, the grouping system reported by Kenyon et al. [3] was adopted in this study. One new group (U12) reported by Umber et al. [39] and three new groups (T13–T15) described by Bömer et al. [33] were also added to the phylogenetic analysis (Figure 4).

3. Results

3.1. DGGE Resolves a Complex Mixture of Badnavirus Sequences Present in Dioscorea Species

To evaluate the potential of DGGE in unravelling the complexity of DBV diversity, yam DNAs (n = 47) were selected based on all those scoring PCR-positive for badnavirus sequences using the generic badnavirus primer pair Badna-FP/-RP. Yang et al. [59] designed this degenerate primer pair based on the consensus sequences of RT and RNaseH coding regions of published badnavirus sequences at the time. These primers are widely used in badnavirus research and also proved to be functional in several DBV diversity studies performed by Bousalem et al. [34], Kenyon et al. [3] and Seal et al. [38] among others. DGGE analysis of Badna-FP/-RP positive PCR products required the addition of a GC clamp to one of the two primers. A GC clamp is usually positioned at the 5′ end of the forward primer [44]. Due to the degeneracy of the generic badnavirus primer set, we decided to test the addition of a GC clamp to both, the Badna-forward as well as the—reverse primer. A subset of eight yam samples grown at the NRI’s quarantine glasshouse was selected and the variability of DBV sequences amplified using the two different Badna GC-clamped primers was compared by DGGE. Distinctive bands that were found to be sharp and intense were excised, cloned and sequenced (Figure 2).
Both Badna GC-clamped primers were able to produce PCR amplifications of the expected 619-bp size with only sample TDr 96/00629 failing in both reactions, the probable cause considered to be PCR inhibitors (Figure 2A,B). ‘Direct’ PCR amplification used the GC-clamped primers with DNA extractions as template, rather than ‘nested’ conditions, where purified Badna-FP/-RP products were used as templates and re-amplified. The ‘direct’ approach was found to improve the resolution of the PCR products analysed by DGGE and resulted in distinct but less intense bands on agarose gels when a lower number of PCR cycles (e.g., 25 cycles) was applied. These relatively sharp and intense DGGE bands could be excised under minimal UV exposure. Conversely, very intense PCR bands as a result of more PCR cycles being applied (35–40 cycles), led to an increase in smearing and less distinct banding patterns during DGGE analysis. Twenty-one bands were excised and successfully processed further (Figure 2).
To obtain phylogenetic information from DGGE analysis, the excised bands were re-amplified, cloned and sequenced. This step lowered the likelihood of obtaining multiple DNA sequences from a single band of interest due to the close migration of different sequences on DGGE gels. With the exception of some samples, two clones per excised DGGE band were sequenced, with both clone sequences being presented unless they were found to be 100% identical to each other. The clone sequence data set (Table 1) was coded as follows: the first two letters stand for the country of origin (NG = Nigeria), ‘b’ represents breeding line samples, ‘l’ represents landrace yam samples, the middle number denotes the position of the excised DGGE band, the next letter denotes the clone (a = clone a and b = clone b) and the last two letters refer to the Dioscorea host species (e.g., Dr = Dioscorea rotundata).
Duplicate clone sequences originating from DGGE band numbers 1–21 resulted in >98% identical sequences, except for band numbers 8, 16 and 18. Only one clone was sequenced successfully for each of the DGGE bands 4 and 5 (Figure 2 and Table 1). Although DGGE bands from different plant samples migrating at the same position in the gel usually resulted in >99% identical sequences (e.g., see bands 10 and 21), sometimes >99% identical sequences migrated at different positions. This was for example the case for bands 12 and 13, that showed 99–100% nucleotide identity to the eDBV 12-clone sequence S1a4Dr (KF829956, [39]) (Table 1). Two of the four clone sequences generated from DGGE bands 12 and 13, NGl12bDr and NGl13aDr are 100% identical across the 528-bp region used for taxonomic assessment of badnaviruses [62]. However, they differ in their primer sequences (determined by examining the cloned excised band sequences), and thus the migration of these sequences at different positions in DGGE is likely a result of the significant degeneracy present in the Badna-FP/-RP primer pair.
The reproducibility of the DGGE technique was tested by running samples from two individual plants for both TDr 89/02475 and TDr 1892 accessions. The DGGE patterns for these biological replicates were found to be identical in both of the primer combinations tested (Figure 2). Nine DGGE bands of PCR amplifications using the Badna-F GC-clamp primer (Figure 2A) included DBV sequences assigned to four different monophyletic groups according to the phylogenetic analysis presented in Figure 4. Sequences clustered into groups K9 (bands 1–3), K8 (all of bands 4–7, one band 8 clone NGb8bDr), U12 (band 9) and T13 (band 8 clone NGb8aDr). Several more bands of low intensity were excised, but could not be processed successfully. The majority of sequences from DGGE bands 10–21 amplified using the Badna-R GC-clamp primer (Figure 1B) clustered into group K8 (Table 1). Band numbers 12 and 13 represent sequences of group U12, showing two distinct positions in the DGGE analysis. The only group K9 sequence was identified for band number 16 (NGb16bDr), however the second clone NGb16aDr represented a K8 sequence. In comparison, these findings suggest that group K9 sequences could be under-represented in PCR amplifications using the Badna-R GC-clamp primer. Hence, we decided to use the Badna-F GC-clamp primer in DGGE analysis hereafter.
Screening of D. rotundata breeding lines (n = 112) and landraces (n = 15) maintained at the IITA generated 100% badnavirus PCR-positive results in a study performed by Seal et al. [38] (details of samples given in Supplementary materials of [38]). In this study, DGGE was used to study the diversity of the badnavirus sequences that exist within a subset of the 100% PCR-positive samples (Table 1). Our DGGE analysis revealed several unique bands but also many bands that are shared in a high proportion of the yam material tested. A total of at least 15 distinct bands were differentiated by DGGE, and following sequencing these could be assigned to different partial RT-RNaseH sequences. The dataset (Table 1) created reflects the complexity of the diversity of badnavirus sequences present in Dioscorea species.
Twenty-one DGGE bands were successfully cloned and sequenced from the DGGE gel presented in Figure 3. Several sequences (DGGE bands) appear to be common to all breeding lines of D. rotundata. This is for example the case for sequences NGl14bDr, NGb30bDr and NGb55aDr, all clustering with species group K8 (Table 1), and being 100% identical to BfA103Dc (AM503393, [34]) and 99% identical to known endogenous sequences of K8, such as S2h9Dr (KF829997, [39]). BLAST similarity matches showing 99–100% identity to BfA103Dc were identified for another 23 DGGE clone sequences (Table 1). Additionally, DGGE patterns can be compared across the samples tested and similarities or differences are easily observed. For example, crossings TDr 89/02475 × TDr 97/00777 and TDr 99/02793 × TDr 1892 showed identical DGGE patterns, suggesting the presence of the same set of DBV sequences.
Eleven cross-breeding lines of D. alata were also analysed in this study (Figure S1). Nine out of these 11 cross-breeding lines contain a double-band made out of sequences NGb60Da and NGb61Da, which were found to be 92% identical to GyJT2Dt (AM503389, [34]) and 93% identical to NG1Da (AM944571, [35]), respectively. Sequence NGb60Da appears to be not only common to the majority of D. alata breeding lines tested in this study, but is also common to most of the D. rotundata material tested (e.g., sequence NGb6aDr corresponding to band 6 in Figure 2). Additionally, eight of the 11 cross-breeding lines analysed contain a second double-band which was not resolved very well, but contained sequence NGb63Da identified as 99% identical to SB42Da (AM072696), which was isolated from a D. alata plant in the Solomon Islands in 2000 [34] and fell into monophyletic group K1 (Table 1). Both double-bands appear to be conserved across most of the D. alata breeding lines analysed in this study and suggest the presence of integrated sequences in D. alata germplasm.

3.2. Phylogenetic Diversity of Dioscorea Badnavirus Sequences

All sequences produced in this study were subjected to similarity BLAST searches in the NCBI GenBank databases and nearest matches as well as percent identities can be found in Table 1. The phylogenetic analysis of the 527–528-bp-long partial RT-RNaseH nucleotide sequences showed that the 114 sequences fall within six monophyletic groups according to the suggested classification of yam badnaviruses [3,33,34,39]. The groups include K1, K5, K8, K9, U12 and T13 (Figure 4 and Table 1). A nucleotide percent similarity matrix was generated using MAFFT and is presented in Table S1. We included all DGGE-derived badnavirus sequences in the phylogenetic analysis, as this approach adds another layer of information by giving an approximate indication about the prevalence of a particular sequence in the material tested.

3.3. Monophyletic Group Assignment of Sequences Identified in This Study

Five DBV sequences, with each of those derived from individual DGGE bands, clustered into monophyletic group K1 described by Kenyon et al. [3]. Two sequences originating from D. dumetorum and D. esculenta samples share 99% nucleotide identity to FJ65bDe (AM072660) and FJ75cDe (AM072663), which were sampled from Fiji in 1999 [34]. Three further sequences were isolated from two D. alata samples and one D. dumetorum sample, showing 99% nucleotide identity to SB42Da (AM072696) isolated from a D. alata plant in the Solomon Islands in 2000 [34]. Six DBV sequences clustered into monophyletic group K5, with five of these originating from D. rotundata samples and sharing 98–99% nucleotide sequence identity to the eDBV5 clone S1g6Dr (KF829974, [39]) and NGl1950Dr (KX008589, [33]). The latter sequence NGl1950Dr was amplified from a D. rotundata plant by rolling circle amplification (RCA) and was considered to therefore most probably be an episomal sequence. Most DGGE-derived sequences (58 out of 114 in total) clustered into monophyletic group K8 and the majority of those originated from D. rotundata samples (Figure 4 and Table 1). Twenty-nine out of 114 DBV sequences in total were assigned to monophyletic group K9. Thirteen DBV sequences originating from D. rotundata samples clustered into monophyletic group U12 described by Umber et al. [39]. Only one sequence identified by DGGE fell into monophyletic group T13 described in [33]. This sequence (NGb8aDr) appears to be comigrating with DGGE sequence NGb8bDr, which is clustering into monophyletic group K8 and common to most of the D. rotundata germplasm tested in this study. A full description of the relationship of sequences is given in Appendix A.

3.4. Conservation of Amino Acid Motifs in Partial RT-RNaseH Badnavirus Sequences

BLAST similarity searches of the partial RT-RNaseH coding region derived from DGGE band sequences showed that most of the sequences had close identity to a number of existing badnavirus sequences (Table 1). The analysis of the deduced amino acid sequences of the badnavirus sequences identified in this study, and published sequences of other badnaviruses have shown some distinctive conserved and semi-conserved regions of the family Caulimoviridae; the regions represent the ‘FIAVYIDDILVFS’ motif [73,74] at position 17–29 of the deduced protein sequence and the ‘LKTTKGLRSWLGILNYAR’ motif [35] at position 95–112 of the deduced protein sequence of the 528-bp-long partial RT-RNaseH. Moreover, the protein alignment presented in Figure S2 allows the simple identification of single amino acid changes or major differences between the phylogenetic groups. Clear patterns specific to every monophyletic species group can be observed, with for example valine at position 98 present in all sequences clustering into group U12. All protein sequences analysed in this study were compared to the reference sequence of Dioscorea bacilliform alata virus (DBALV or DaBVa, X94576-XX94581, [74]).

4. Discussion

4.1. Potential of DGGE in Badnavirus Diversity Studies and Identification of Potential Integrated Sequences

In this study, the variation in complex (e)DBV sequence mixtures across a subset of the 127 IITA D. rotundata samples that scored 100% Badna-positive by PCR in a study by Seal et al. [38] was rapidly and robustly evaluated. The results presented here showed that DGGE is a very useful technique for diversity profiling of the amplified DBV partial RT-RNaseH sequences. This is in agreement with previous reports on the ability of DGGE to differentiate DNA sequences of high similarity in given samples, such as in the diversity studies of viruses, bacteria and phytoplankton [54,75]. A particular strength for the badnavirus sequence diversity is, however, clearly apparent from the similarities in DGGE banding profiles observed between breeding lines of the same Dioscorea species indicating highly conserved sequences that appear to represent eDBVs. For example, crossings TDr 89/02475 × TDr 97/00777 and TDr 99/02793 × TDr 1892 show identical DGGE patterns (Figure 3), suggesting the presence of the same set of DBV sequences, with the majority of those likely to be of endogenous nature and ultimately indicating a very similar genetic background of these samples.
The DGGE banding profiles in this study and the resulting clone sequences thereof, depict the vast complexity of DBV sequences present in yam germplasm. The precise analysis of all DBV sequences present in a given sample using DGGE as a method, however, is not straightforward in the light of this complexity and the difficulty of cloning closely migrating bands. Our efforts focused on cloning clearly defined bands. For example, the majority of breeding lines of the same Dioscorea species share the common DGGE bands 5, 10 and 11 for D. rotundata (Figure 2). Sequences from these bands and others cluster closely together and a total of 22 sequences share 99–100% nucleotide identity to BfA103Dc (AM503393, [34]). It appears probable that these bands may represent integrated sequences as the presence of such high identity episomal viruses in all the different material would not be expected. Episomal badnavirus sequences infecting other crops have been shown to be highly diverse, as illustrated by the genetic diversity of four banana streak virus (BSV) isolates from Australia ranging from 21.8 to 33.6% in a comparison of the amino acid sequences of the ribonuclease H domain in open reading frame 3 (ORF3) [24] and nucleotide diversity of up to 18% of the RT-RNaseH-coding domain of six BSVs from East Africa [41]. Further evidence is that BfA103Dc was shown to be free of episomal viruses by immunosorbent electron microscopy (ISEM) and ELISA [38]. Hybridisation using the common bands as probes in Southern blots would confirm whether DGGE bands 5, 10 and 11 represent common ancient integrated sequence in the host yam genomes rather than an unusually homogeneous episomal viral infection.
The presence of more than one band within every single lane of DGGE analysis reveals that every leaf sample tested contained more than one badnavirus sequence. The presence of sequences representing more than a single badnavirus species in single-leaf samples has been reported before through sequencing of cloned PCR products and generating clones representing several different badnavirus species from a single leaf sample [3,35,60]. This study reveals that in D. rotundata this is in fact the norm, rather than an exception as also indicated in our previous study detecting mixed infections using RCA as a tool to amplify episomal badnavirus sequences [33]. In this context we propose DGGE to be used as complementary method to RCA, enabling the rapid identification of potentially integrated sequences while screening for badnavirus diversity.
Testing of DGGE as a tool in the analysis of complex badnavirus sequence mixtures present in samples of yam breeding lines in this study was based on the assumption that the previous procedure by cloning and sequencing of PCR products is an inefficient method to reveal full sequence diversity. Such techniques require the selection and sequencing of a large number of clones in order to unravel all sequence diversity present in a given sample as demonstrated also for BSV [76]. Equally, RFLP was reported not much a useful tool for yam badnavirus diversity, since sequences with different RFLP patterns were found clustered in the same species group [35]. The above illustrates that the practical advantage of DGGE is that it is a rapid means of detecting sequence diversity in uniformly sized PCR products, eliminating the need for labour-intensive screening of redundant clones [77]. Equally, in a situation when there are a large number of samples to be analysed that contain multiple badnavirus sequences, DGGE is the most appropriate technique in terms of cost and practicality for a diversity study compared to direct sequencing of clones following PCR or RFLP.
Conversely, DGGE suffers from some methodical pitfalls and inherent practical limitations. One concern is the detection limit, with Muyzer et al. [45] suggesting that any target DNA that is less than 1% of the total target pool is unlikely to be detected by DGGE. As such, low badnavirus titers could lead to the unlikely detection of true episomal DBV sequences by DGGE, compared to endogenous sequences potentially existing in higher copy numbers in some yam germplasm. As a related caveat, discrete fingerprint bands may not always be apparent when analysing highly diverse samples, leading to smearing or poorly resolved patterns [44]. Another related limitation is the comigration of DNA molecules of different sequence as reported by Ercolini [49] and also experienced in this study. Although DGGE was able to separate sequences with high similarities, problems included accurate cutting of closely-migrating DGGE bands in particular under suboptimal resolution due to the occurrence of smeared backgrounds as can be seen in several samples presented in Figure 3. Additionally, some technical problems were encountered in carrying out the DGGE technique, namely general handling of fragile DGGE gels and the persistent trouble of background smearing in gels that made band-scoring problematic. Several modifications were performed on the PCR conditions, ranging from altering concentrations of primers, the choice of the Taq polymerase enzyme, the DGGE gradient range and the volume of PCR products loaded per gel lane. In this context, we found the procedural comments on the DGGE technique given by Green et al. [44] useful. Different primers (Badna-FP versus Badna-RP GC-clamped primer) were also tried, but the comparison shown in Figure 2 concluded a limited suitability of the Badna-RP GC-clamp primer in the DBV diversity study because of the potential under-representation of badnavirus sequences clustering in the K9 monophyletic group.
For the future, a worthwhile additional improvement would be the inclusion of a reference DGGE standard in which a profile is created using pooled clones from distinct DGGE bands. Here, clone sequences or PCR re-amplification products would be analysed by DGGE in a lane adjacent to the original sample PCR product from which it was excised. This will assist the typing of unknown populations and eliminate the need for DGGE band cloning and subsequent sequencing. Moreover, the degeneracy of the Badna-FP/-RP primers is expected to contribute to the poor resolution (smearing experienced in this study), potentially creating multiple products per template sequence differing only in their primer binding site. Hence, the design of monophyletic group specific primers should be considered in the future, potentially improving DGGE resolution while reducing the complexity at the same time by focusing the analysis on the most prevalent and interesting DBV species.

4.2. DGGE-Captured Badnavirus Diversity

To investigate the in-depth diversity of badnavirus sequences present in the yam samples analysed in this study, DGGE analysis was applied and representative bands were cloned and sequenced. All sequenced bands were shown to be partial RT-RNaseH sequences of badnaviruses, with 112 DBV sequences analysed clustering into six (K1, K5, K8, K9, U12 and T13) of the 15 putative species groups identified by Bömer et al. [33], Bousalem et al. [34], Kenyon et al. [3] and Umber et al. [39]. A comparison of their amino acid sequences generated from the translated partial RT-RNaseH region sequences together with other members of the family Caulimoviridae revealed conserved and semi-conserved regions similar to the previously published badnavirus sequences [3,33,34,35,38,73,78,79,80]. In particular, a stretch of the conserved ‘FIAVYIDDILVFS’ region of the RT in the C-terminal of the ORF3 polyprotein was observed, as was the semi-conserved ‘LKTTKGLRSWLGILNYAR’ region (Figure S2). This confirms that DNA extracted from Dioscorea spp. samples contained sequences belonging to the genus Badnavirus of the family Caulimoviridae [73,78,79,80,81].
Sequence analyses of the partial RT-RNaseH domain in this work support the classification and diversity of yam badnaviruses proposed by Kenyon et al. [3], and subsequently substantiated upon by Bousalem et al. [34], Umber et al. [39] and Bömer et al. [33]. The maximum variability recorded within the partial RT-RNaseH coding region at the nucleotide level among the 112 yam badnavirus sequences determined in this study was 37%. This degree of variability was within the range of intergroup nucleotide diversity for group K1–K11 (23.1–39.4%) of yam badnaviruses [3]. For badnaviruses present in other plant hosts, high diversity levels have also been reported, with 21.8–33.6% sequence diversity in a comparison of the amino acid sequences of the ribonuclease H domain in ORF3 reported for four BSV isolates from Australia [24], 28% maximum nucleotide diversity for Ugandan BSV isolates [2], 29.4% maximum nucleotide diversity between complete genomes of CSSV isolates [82], and up to 33.5% nucleotide variability within the partial RT-RNaseH sequences for isolates of sugarcane bacilliform viruses (SCBVs) [5].
Two sequences (NGl65De and NGb53Dr) determined in this study shared <60% nucleotide identity with any previously identified yam badnavirus sequence group but were similar to the divergent sequences groups of K12 and K13 reported by Kenyon et al. [3] from South Pacific yams. The sequences clustered closer to RTBV, a more distantly related virus in the family Caulimoviridae, than any member of the genus Badnavirus. It is possible that these sequence groups represent divergent badnaviruses, ancient endogenous pararetrovirus sequences or new genera within the family Caulimoviridae [3]. Further research is needed to characterise the nature of these sequences.

4.3. Endogenous Badnavirus Partial RT-RNaseH Sequences

To date, four of the 15 yam badnavirus species groups identified have been reported to contain eDBV sequences, namely K5, K8, K9 and U12 [38,39]. Phylogenetic analyses of the cloned DGGE bands revealed that all the sequences that originated from DGGE bands common to the majority of D. rotundata (for example NGb30bDr and NGb1aDr) samples clustered into either the K8 or the K9 species groups (Figure 4). These common bands show 99–100% nucleotide identity to several more sequences identified in this study as well as integrated DBV sequences described by Umber et al. [39], and hence, it seems probable that they represent eDBV sequences. A number of DGGE band sequences from D. rotundata were found to both share 99% nucleotide identity to eDBV9 clone S1h2Dr (KF829975, [39]) or eDBV8 clone S2h9Dr (KF829997, [39]) respectively, suggesting that these sequences represent integrated sequences in the samples tested. Several more DGGE band sequences present in D. rotundata were found that clustered with the other two yam badnavirus/eDBV species groups (K5 and U12) (Figure 4). For example, sequences NGl9a/bDr, NGl12a/bDr, and NGl13a/bDr of the U12 species group were from bands that were present in several D. rotundata. They share 99–100% nucleotide identity to eDBV12 clone S1a4Dr (KF829956, [39]), suggesting these sequences could also be integrated sequences in the samples tested. Equally, NGb22a/bDr, NGb23aDr, NGb24aDr, NGb52Dr could also represent integrated sequences as the sequences share 98–99% nucleotide identity to eDBV5 clone S1g6Dr (KF829974, [39]).
In summary, 28 sequences showed eDBV sequences as their nearest match with nucleotide identities above 99%. Twenty-two sequences were also identified that have 99–100% identity to the badnavirus-particle free BfA103Dc (AM503393, [34]) and as a result are likely to represent integrated sequences as described above. Additionally, a further 32 sequences clustering into monophyletic group K8 are highly conserved among themselves (GyJT2DT, AM503389, [34], being the nearest match), suggesting that these sequences are also integrated. These findings together lead to at least 82 of 112 (72%) sequences in this study considered to most probably originate from integrated sequences. It is proposed that DGGE therefore is very effective at revealing eDBV sequences when testing plant germplasm.
In this study, the breeding lines were derived from true seed yam, and this material has been grown over successive cropping seasons at the IITA. The true seed material should have been free of episomal badnavirus particles, as yam badnaviruses have not been reported to be seed-transmissible [62]. However, the origin from true seed yam does not automatically exclude that episomal forms of virus can be present. Vertical transmission of integrated copies via seeds is possible as reported for endogenous PVCV in petunia, where the integrated virus gets activated by direct transcription of integrated PVCV sequences in the form of a tandem array [22]. Vertical transmission of activatable eDBVs cannot be excluded in yams and more research on this potential risk is needed. Nevertheless, the presence of all the common DGGE bands in yam lines analysed in this study suggests that these most likely represent eDBV sequences, whereas additional DGGE bands detected only in some lines potentially represent new infections during the propagation of the material in the field over many cultivation seasons. Further analysis of the DGGE sequences using the Southern blot technique will help to determine if these sequences are diverse eDBVs or represent episomal viruses acquired through screen-house and field propagation.

5. Conclusions

In summary, a workflow combining PCR and DGGE methods for rapid and efficient determination and unravelling of complex mixtures of potentially episomal and endogenous badnavirus sequences has been developed (Figure 1). Here, we used specific PCR primers to amplify a partial sequence of the badnavirus RT-RNaseH coding region to investigate the molecular diversity of (e)DBV sequences in selected West African yam germplasm. PCR products were resolved using DGGE, giving characteristic banding patterns for each yam line examined. This approach is described here for the first time for the assessment of badnavirus sequences present in Dioscorea spp. germplasm and resulted in the identification of complex DGGE profiles representing multiple sequences of DBVs. A total of 112 yam badnavirus sequences were generated (GenBank accession numbers KY555456 to KY555569), falling into six monophyletic species groups. The conserved nature of several DGGE-derived yam badnavirus sequences, as well as high nucleotide identities to eDBV sequences deposited in GenBank suggests that the majority of West African yam germplasm contains a mixture of integrated badnavirus sequences. Our findings highlight DGGE as an extremely useful technique for rapid indication of badnavirus sequence diversity in such samples containing multiple eDBVs, enabling a snapshot of the diversity between genomes to be taken. The approach taken in this study to enable rapid identification of potential candidate integrated badnavirus sequences, indicated by their conserved nature across germplasm, should have wide application for the study of an ever-increasing number of plant species found to contain this previously understudied, but important genus of plant viruses.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4915/9/7/181/s1, Figure S1: Denaturing gradient gel electrophoresis (DGGE) analysis of partial RT-RNaseH badnavirus sequences from eleven samples consisting of cross-breeding lines of D. alata, Figure S2: Protein alignment from deduced amino acid sequences of partial RT-RNaseH nucleotide sequences of 114 yam badnavirus sequences determined in this study together with other members of the family Caulimoviridae, Table S1: Nucleotide percent similarity matrix.

Acknowledgments

The authors gratefully acknowledge the support of this work by the Bill & Melinda Gates Foundation (BMGF) under the “Development of On-farm Robust Diagnostic Toolkits for Yam Viruses” grant to the Natural Resources Institute (NRI) and the “Yam Improvement for Income and Food Security in West Africa (YIIFSWA)” grant to the International Institute of Tropical Agriculture (IITA). Funding to support open access is provided by the Bill & Melinda Gates Foundation.

Author Contributions

Aliyu A. Turaki, Moritz Bömer and Susan E. Seal conceived and designed the experiments; Aliyu A. Turaki and Moritz Bömer performed the experiments; Aliyu A. Turaki, Moritz Bömer and Gonçalo Silva analyzed the data; P. Lava Kumar contributed materials; Aliyu A. Turaki, Moritz Bömer, Gonçalo Silva, P. Lava Kumar and Susan E. Seal wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BFBadna-forward primer
BLASTbasic local alignment search tool
BRBadna-reverse primer
BSOLVBanana streak OL virus
BSVbanana streak virus
ComYMVCommelina yellow mottle virus
CSSVCacao swollen shoot virus
CTABcetyltrimethylammonium bromide
DBALVDioscorea bacilliform alata virus
DBRTVDioscorea bacilliform rotundata (RT) virus
DBVDioscorea bacilliform virus
DGGEdenaturing gradient gel electrophoresis
DOAJDirectory of open access journals
eDBVendogenous Dioscorea bacilliform virus
GCGC clamp
ICTVInternational Committee on Taxonomy of Viruses
IITAInternational Institute of Tropical Agriculture
ISEMimmunosorbent electron microscopy
LDlinear dichroism
MAFFTMultiple Alignment using Fast Fourier Transform
NCBINational Centre for Biotechnology Information
NGNigeria
NGSnext-generation sequencing
NRINatural Resources Institute
o.p.open pollinated
ORFopen reading frame
PCRpolymerase chain reaction
RCArolling circle amplification
RFLPrestriction fragment length polymorphism
RNaseHribonuclease H
RTreverse transcriptase
RTBVRice tungro bacilliform virus
SCBMOVSugarcane bacilliform MO virus
TaBVTaro bacilliform virus
TAETris-acetate-EDTA
TBETris-Boric acid-EDTA
TDaDioscorea alata accession
TDbDioscorea bulbifera accession
TDcDioscorea cayenensis accession
TDeDioscorea esculenta accession
TDrDioscorea rotundata accession

Appendix A

Five DBV sequences, with each of those derived from individual DGGE bands, clustered into monophyletic group K1 described by Kenyon et al. [3]. Sequences NGl71Dd isolated from D. dumetorum and NGl73De from D. esculenta both share 99% nucleotide identity to FJ65bDe (AM072660) and FJ75cDe (AM072663), which were sampled from Fiji in 1999 by Lebas and described by Bousalem et al. [34] as representative sequences of the DeBV-A sub-group A (Dioscorea esculenta bacilliform virus A) monophyletic group. Sequences NGb63Da and NGl72Da isolated from D. alata, as well as NGl70Db from D. dumetorum show 99% nucleotide identity to SB42Da (AM072696) isolated from a D. alata plant in the Solomon Islands in 2000 by Lebas and described by Bousalem et al. [34] as a representative sequence of the DeBV-A sub-group B.
Six DBV sequences clustered into monophyletic group K5, with five of these originating from D. rotundata samples (NGb52Dr, NGb22a/bDr, NGb23aDr and NGb24aDr) and sharing 98–99% nucleotide sequence identity to eDBV5 clone S1g6Dr (KF829974, [39]) and NGl1950Dr (KX008589, [33]). The latter sequence NGl1950Dr was amplified from a D. rotundata plant by rolling circle amplification (RCA) in our previous study [33] and was considered to therefore most probably be an episomal sequence. Another DGGE-derived sequence (NGl66Dd) isolated from D. dumetorum shows 93% nucleotide sequence identity to eDBV5 clone S1un5Dr (KF830000, [39]).
Most DGGE-derived sequences (58 out of 114 in total) clustered into monophyletic group K8 and the majority of those originated from D. rotundata samples (Figure 4 and Table 1). Twenty-three DBV sequences assigned to K8 and isolated from D. rotundata samples cluster very closely together and share 97–100% nucleotide identity to BfA103Dc (AM503393, [34]) and NGb0477Dr (KX008586, [33]). Sequence NGb0477Dr was amplified by RCA in our previous study showing 99% nucleotide identity to BfA103Dc, Gn1633Dr, Gn845Dr and Gn502Dr [33]. Sequence BfA103Dc (AM503393, [34]) was derived from a Pilimpikou yam sample which was reported as viral particle-free, whereas the other three sequences derived from yam samples of Guinea in which badnavirus particles were detected using ISEM in a previous study by Seal et al. [38]. Another K8 sequence (NGl67Dd) originating from D. dumetorum is 98% identical to BN2Da from Benin (AM944584, [35]), while sequence NGl68Dr shares 96% identity with GN2Dr from Ghana (AM944575, [35]), and sequence NGb61Da is 93% identical to NG1Da from Nigeria (AM944571, [35]). Moreover, 32 DGGE-derived sequences share 91–93% nucleotide identity to GyJT2Dt isolated in French Guiana in 1998 (AM503389, [34]).
Twenty-nine out of 114 DBV sequences in total were assigned to monophyletic group K9. Of those, nine DGGE-derived sequences originating from D. rotundata samples share 99% nucleotide sequence identity to BN4Dr (AM944586, [35]) isolated from a D. rotundata plant in Benin. A further four sequences (NGb25aDr, NGb37bDr, NGb54Dr and NGb41bDr) isolated from D. rotundata are 99% identical to eDBV9 clone S1h2Dr (KF829975, [39]), and two more sequences (NGb29bDr and NGb37aDr) from D. rotundata samples are 99% identical to eDBV9 clone S2f8Dr (KF829993, [39]). Moreover, nine DGGE-derived sequences (NGb1bDr, NGb2a/bDr, NGb3a/bDr, NGb35a/bDr, NGb42bDr and NGb16bDr) isolated from D. rotundata samples share 99–100% identity with 432B39Ds (AM503361, [34]) derived from D. sansibarensis isolated in Benin 1999 in which badnavirus particles were detected using ISEM [38] and FJ60aDr (AM072658, [34]) from D. rotundata isolated in Fiji 1999 by Lebas. A further five sequences are 89–100% identical to either eDBV9 clone G1Dr (KF830002, [39]), eDBV9 clone S1h6Dr (KF829977, [39]), eDBV9 clone S1b6Dr (KF829960, [39]), eDBV9 clone S1e3Dr (KF829969, [39]) or eDBV9 clone S1b2Dr (KF829958, [39]), respectively.
Thirteen DBV sequences originating from D. rotundata samples clustered into monophyletic group U12 as described by Umber et al. [39]. Nine of these sequences (NGl9a/bDr, NGl12a/bDr, NGl13a/bDr, NGb41aDr, NGb42aDr and NGb56Dr) are 99–100% identical to eDBV12 clone S1a4Dr (KF829956, [39]) and GHL24Dr isolated in Ghana 1995 (AM072665, [34]). Three further sequences are 99% identical to eDBV12 clone S2a7Dr (KF829978, [39]) and Mt818Dr isolated in Martinique 1998 (AM503378, [34]). The remaining sequence NGb59Dr is 99% identical to eDBV12 clone S2h10Dr (KF829998, [39]).
Only one sequence identified by DGGE fell into the monophyletic group T13 described in our previous study [33]. Sequence NGb8aDr has 99% identity to the episomal Dioscorea bacilliform rotundata (RT) virus 1 (DBRTV1) sequence originating from D. rotundata (KX008596, [33]), which is also available as a full genome sequence (GenBank KX008574). However, the second clone of this DGGE band clustered into group K8, indicating cross-contamination with a closely migrating band.

References

  1. Borah, B.K.; Sharma, S.; Kant, R.; Johnson, A.M.A.; Saigopal, D.V.R.; Dasgupta, I. Bacilliform DNA-containing plant viruses in the tropics: Commonalities within a genetically diverse group. Mol. Plant Pathol. 2013, 14, 759–771. [Google Scholar] [CrossRef] [PubMed]
  2. Harper, G.; Hart, D.; Moult, S.; Hull, R.; Geering, A.; Thomas, J. The diversity of Banana streak virus isolates in Uganda. Arch. Virol. 2005, 150, 2407–2420. [Google Scholar] [CrossRef] [PubMed]
  3. Kenyon, L.; Lebas, B.S.M.; Seal, S.E. Yams (Dioscorea spp.) from the South Pacific Islands contain many novel badnaviruses: Implications for international movement of yam germplasm. Arch. Virol. 2008, 153, 877–889. [Google Scholar] [CrossRef] [PubMed]
  4. Lockhart, B.E.L. Purification and serology of a bacilliform virus associated with Banana streak disease. Phytopathology 1986, 80, 995–999. [Google Scholar] [CrossRef]
  5. Muller, E.; Dupuy, V.; Blondin, L.; Bauffe, F.; Daugrois, J.H.; Laboureau, N.; Iskra-Caruana, M.L. High molecular variability of sugarcane bacilliform viruses in Guadeloupe implying the existence of at least three new species. Virus Res. 2011, 160, 414–419. [Google Scholar] [CrossRef] [PubMed]
  6. Harper, G.; Dahal, G.; Thottappilly, G.; Hull, R. Detection of episomal banana streak badnavirus by IC-PCR. J. Virol. Methods 1999, 79, 1–8. [Google Scholar] [CrossRef]
  7. Ndowora, T.; Dahal, G.; LaFleur, D.; Harper, G.; Hull, R.; Olszewski, N.E.; Lockhart, B. Evidence that badnavirus infection in Musa can originate from integrated pararetroviral sequences. Virology 1999, 255, 214–220. [Google Scholar] [CrossRef] [PubMed]
  8. Le Provost, G.; Iskra-Caruana, M.L.; Acina, I.; Teycheney, P.Y. Improved detection of episomal Banana streak viruses by multiplex immunocapture PCR. J. Virol. Methods 2006, 137, 7–13. [Google Scholar] [CrossRef] [PubMed]
  9. Bhat, A.I.; Hohn, T.; Selvarajan, R. Badnaviruses: The Current Global Scenario. Viruses 2016, 8, 177. [Google Scholar] [CrossRef] [PubMed]
  10. Geering, A.D.W.; Maumus, F.; Copetti, D.; Choisne, N.; Zwickl, D.J.; Zytnicki, M.; McTaggart, A.R.; Scalabrin, S.; Vezzulli, S.; Wing, R.A.; et al. Endogenous florendoviruses are major components of plant genomes and hallmarks of virus evolution. Nat. Commun. 2014, 5, 5269. [Google Scholar] [CrossRef] [PubMed]
  11. Geering, A.D.W.; Scharaschkin, T.; Teycheney, P.Y. The classification and nomenclature of endogenous viruses of the family Caulimoviridae. Arch. Virol. 2010, 155, 123–131. [Google Scholar] [CrossRef] [PubMed]
  12. Mette, M.F.; Kanno, T.; Aufsatz, W.; Jakowitsch, J.; van der Winden, J.; Matzke, M.A.; Matzke, A.J.M. Endogenous viral sequences and their potential contribution to heritable virus resistance in plants. EMBO J. 2002, 21, 461–469. [Google Scholar] [CrossRef] [PubMed]
  13. Staginnus, C.; Iskra-Caruana, M.L.; Lockhart, B.; Hohn, T.; Richert-Pöggeler, K.R. Suggestions for a nomenclature of endogenous pararetroviral sequences in plants. Arch. Virol. 2009, 154, 1189–1193. [Google Scholar] [CrossRef] [PubMed]
  14. LaFleur, D.A.; Lockhart, B.E.L.; Olszewski, N.E. Portions of the banana streak badnavirus genome are integrated in the genome of its host Musa sp. Phytopathology 1996, 86, S100. [Google Scholar]
  15. Richert-Pöggeler, K.R.; Shepherd, R.J. Petunia vein-clearing virus: A plant pararetrovirus with the core sequences for an integrase function. Virology 1997, 236, 137–146. [Google Scholar] [CrossRef] [PubMed]
  16. Jakowitsch, J.; Mette, M.F.; van Der Winden, J.; Matzke, M.A.; Matzke, A.J. Integrated pararetroviral sequences define a unique class of dispersed repetitive DNA in plants. Proc. Natl. Acad. Sci. USA 1999, 96, 13241–13246. [Google Scholar] [CrossRef] [PubMed]
  17. Kunii, M.; Kanda, M.; Nagano, H.; Uyeda, I.; Kishima, Y.; Sano, Y. Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: Implications for integration and evolution. BMC Genom. 2004, 5. [Google Scholar] [CrossRef] [PubMed]
  18. Hansen, C.N.; Harper, G.; Heslop-Harrison, J.S. Characterisation of pararetrovirus-like sequences in the genome of potato (Solanum tuberosum). Cytogenet. Genome Res. 2005, 110, 559–565. [Google Scholar] [CrossRef] [PubMed]
  19. Staginnus, C.; Richert-Pöggeler, K.R. Endogenous pararetroviruses: Two-faced travelers in the plant genome. Trends Plant Sci. 2006, 11, 485–491. [Google Scholar] [CrossRef] [PubMed]
  20. Chabannes, M.; Baurens, F.C.; Duroy, P.O.; Bocs, S.; Vernerey, M.S.; Rodier-Goud, M.; Barbe, V.; Gayral, P.; Iskra-Caruana, M.L. Three infectious viral species lying in wait in the banana genome. J. Virol. 2013, 87, 8624–8637. [Google Scholar] [CrossRef] [PubMed]
  21. Gayral, P.; Noa-Carrazana, J.C.; Lescot, M.; Lheureux, F.; Lockhart, B.E.; Matsumoto, T.; Piffanelli, P.; Iskra-Caruana, M.L. A single Banana streak virus integration event in the banana genome as the origin of infectious endogenous pararetrovirus. J. Virol. 2008, 82, 6697–6710. [Google Scholar] [CrossRef] [PubMed]
  22. Richert-Pöggeler, K.R.; Noreen, F.; Schwarzacher, T.; Harper, G.; Hohn, T. Induction of infectious petunia vein clearing (pararetro) virus from endogenous provirus in petunia. EMBO J. 2003, 22, 4836–4845. [Google Scholar] [CrossRef] [PubMed]
  23. Lockhart, B.E.; Menke, J.; Dahal, G.; Olszewski, N.E. Characterization and genomic analysis of tobacco vein clearing virus, a plant pararetrovirus that is transmitted vertically and related to sequences integrated in the host genome. J. Gen. Virol. 2000, 81, 1579–1585. [Google Scholar] [CrossRef] [PubMed]
  24. Geering, A.D.; McMichael, L.A.; Dietzgen, R.G.; Thomas, J.E. Genetic Diversity Among Banana streak virus Isolates from Australia. Phytopathology 2000, 90, 921–927. [Google Scholar] [CrossRef] [PubMed]
  25. Iskra-Caruana, M.L.; Baurens, F.C.; Gayral, P.; Chabannes, M. A four-partner plant-virus interaction: Enemies can also come from within. Mol. Plant Microbe Interact. 2010, 23, 1394–1402. [Google Scholar] [CrossRef] [PubMed]
  26. Iskra-Caruana, M.L.; Chabannes, M.; Duroy, P.O.; Muller, E. A possible scenario for the evolution of Banana streak virus in banana. Virus Res. 2014, 186, 155–162. [Google Scholar] [CrossRef] [PubMed]
  27. Côte, F.X.; Galzi, S.; Folliot, M.; LamagnÈre, Y.; Teycheney, P.Y.; Iskra-Caruana, M.L. Micropropagation by tissue culture triggers differential expression of infectious endogenous Banana streak virus sequences (eBSV) present in the B genome of natural and synthetic interspecific banana plantains. Mol. Plant Pathol. 2010, 11, 137–144. [Google Scholar] [CrossRef] [PubMed]
  28. Dallot, S.; Acuna, P.; Rivera, C.; Ramirez, P.; Cote, F.; Lockhart, B.E.L.; Caruana, M.L. Evidence that the proliferation stage of micropropagation procedure is determinant in the expression of Banana streak virus integrated into the genome of the FHIA 21 hybrid (Musa AAAB). Arch. Virol. 2001, 146, 2179–2190. [Google Scholar] [CrossRef] [PubMed]
  29. Harper, G.; Hull, R.; Lockhart, B.; Olszewski, N. Viral sequences integrated into plant genomes. Annu. Rev. Phytopathol. 2002, 40, 119–136. [Google Scholar] [CrossRef] [PubMed]
  30. Chabannes, M.; Iskra-Caruana, M.L. Endogenous pararetroviruses—A reservoir of virus infection in plants. Curr. Opin. Virol. 2013, 3, 615–620. [Google Scholar] [CrossRef] [PubMed]
  31. Asiedu, R.; Sartie, A. Crops that feed the World 1. Yams. Food Secur. 2010, 2, 305–315. [Google Scholar] [CrossRef]
  32. Ekanayake, I.J.; Asiedu, R. Problems and Perspectives of Yam-Based Cropping Systems in Africa. J. Crop Prod. 2003, 9, 531–558. [Google Scholar] [CrossRef]
  33. Bömer, M.; Turaki, A.; Silva, G.; Kumar, P.; Seal, S. A Sequence-Independent Strategy for Amplification and Characterisation of Episomal Badnavirus Sequences Reveals Three Previously Uncharacterised Yam Badnaviruses. Viruses 2016, 8, 188. [Google Scholar] [CrossRef] [PubMed]
  34. Bousalem, M.; Durand, O.; Scarcelli, N.; Lebas, B.S.M.; Kenyon, L.; Marchand, J.L.; Lefort, F.; Seal, S.E. Dilemmas caused by endogenous pararetroviruses regarding the taxonomy and diagnosis of yam (Dioscorea spp.) badnaviruses: Analyses to support safe germplasm movement. Arch. Virol. 2009, 154, 297–314. [Google Scholar] [CrossRef] [PubMed]
  35. Eni, A.O.; Hughes, J.D.; Asiedu, R.; Rey, M.E.C. Sequence diversity among badnavirus isolates infecting yam (Dioscorea spp.) in Ghana, Togo, Benin and Nigeria. Arch. Virol. 2008, 153, 2263–2272. [Google Scholar] [CrossRef] [PubMed]
  36. Galzi, S.; Scutt, R.; Prophète, P.; Roumagnac, P.; Filloux, D. Assessment and Characterization of the Genetic Diversity of Viruses Infecting Cultivated Yams (Dioscorea spp.) in Haïti. In Rencontres de Virologie Végétale; Marais Armelle, R.F., Ed.; SFP INRA: Paris, France, 2013; p. 70. [Google Scholar]
  37. Lima, J.S.; Lima, G.S.A.; Micheref, S.J. Variabilidade genética de isolados de badnavírus infectando inhame (Dioscorea spp.) no nordeste do Brasil. Trop. Plant Pathol. 2013, 38, 349–353. [Google Scholar] [CrossRef]
  38. Seal, S.; Turaki, A.; Muller, E.; Kumar, P.L.; Kenyon, L.; Filloux, D.; Galzi, S.; Lopez-Montes, A.; Iskra-Caruana, M.L. The prevalence of badnaviruses in West African yams (Dioscorea cayenensis-rotundata) and evidence of endogenous pararetrovirus sequences in their genomes. Virus Res. 2014, 186, 144–154. [Google Scholar] [CrossRef] [PubMed]
  39. Umber, M.; Filloux, D.; Muller, E.; Laboureau, N.; Galzi, S.; Roumagnac, P.; Iskra-Caruana, M.L.; Pavis, C.; Teycheney, P.Y.; Seal, S.E. The genome of African yam (Dioscorea cayenensis-rotundata complex) hosts endogenous sequences from four distinct badnavirus species. Mol. Plant Pathol. 2014, 15, 790–801. [Google Scholar] [CrossRef] [PubMed]
  40. Gayral, P.; Iskra-Caruana, M.L. Phylogeny of banana streak virus reveals recent and repetitive endogenization in the genome of its banana host (Musa sp.). J. Mol. Evol. 2009, 69, 65–80. [Google Scholar] [CrossRef] [PubMed]
  41. James, A.P.; Geijskes, R.J.; Dale, J.L.; Harding, R.M. Molecular characterisation of six badnavirus species associated with leaf streak disease of banana in East Africa. Ann. Appl. Biol. 2011, 158, 346–353. [Google Scholar] [CrossRef]
  42. Jaufeerally-Fakim, Y.; Khorugdharry, A.; Harper, G. Genetic variants of Banana streak virus in Mauritius. Virus Res. 2006, 115, 91–98. [Google Scholar] [CrossRef] [PubMed]
  43. Fischer, S.G.; Lerman, L.S. Length-independent separation of DNA restriction fragments in two-dimensional gel electrophoresis. Cell 1979, 16, 191–200. [Google Scholar] [CrossRef]
  44. Green, S.J.; Leigh, M.B.; Neufeld, J.D. Denaturing Gradient Gel Electrophoresis (DGGE) for Microbial Community Analysis. In Handbook of Hydrocarbon and Lipid Microbiology; Timmis, K.N., Ed.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 4137–4158. ISBN 978-3-540-77587-4. [Google Scholar]
  45. Muyzer, G.; de Waal, E.C.; Uitterlinden, A.G. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 1993, 59, 695–700. [Google Scholar] [PubMed]
  46. Sheffield, V.C.; Cox, D.R.; Lerman, L.S.; Myers, R.M. Attachment of a 40-base-pair G + C-rich sequence (GC-clamp) to genomic DNA fragments by the polymerase chain reaction results in improved detection of single-base changes. Proc. Natl. Acad. Sci. USA 1989, 86, 232–236. [Google Scholar] [CrossRef] [PubMed]
  47. Short, S.M.; Suttle, C.A. Sequence Analysis of Marine Virus Communities Reveals that Groups of Related Algal Viruses Are Widely Distributed in Nature. Appl. Environ. Microbiol. 2002, 68, 1290–1296. [Google Scholar] [CrossRef] [PubMed]
  48. Top, B. A simple method to attach a universal 50-bp GC-clamp to PCR fragments used for mutation analysis by DGGE. Genome Res. 1992, 2, 83–85. [Google Scholar] [CrossRef]
  49. Ercolini, D. PCR-DGGE fingerprinting: Novel strategies for detection of microbes in food. J. Microbiol. Methods 2004, 56, 297–314. [Google Scholar] [CrossRef] [PubMed]
  50. Hu, X.; Zhang, H.; Lu, H.; Qian, G.; Lv, L.; Zhang, C.; Guo, J.; Jiang, H.; Zheng, B.; Yang, F.; et al. The Effect of Probiotic Treatment on Patients Infected with the H7N9 Influenza Virus. PLoS ONE 2016, 11, e0151976. [Google Scholar] [CrossRef] [PubMed]
  51. Valaskova, V.; Baldrian, P. Denaturing gradient gel electrophoresis as a fingerprinting method for the analysis of soil microbial communities. Plant Soil Environ. 2009, 55, 413–423. [Google Scholar]
  52. Michaelsen, A.; Pinzari, F.; Ripka, K.; Lubitz, W.; Piñar, G. Application of molecular techniques for identification of fungal communities colonising paper material. Int. Biodeterior. Biodegrad. 2006, 58, 133–141. [Google Scholar] [CrossRef]
  53. Carmona, M.; Sepúlveda, D.; Cárdenas, C.; Nilo, L.; Marshall, S.H. Denaturing gradient gel electrophoresis (DGGE) as a powerful novel alternative for differentiation of epizootic ISA virus variants. PLoS ONE 2012, 7, e37353. [Google Scholar] [CrossRef] [PubMed]
  54. Lyttle, D.J.; Orlovich, D.A.; Guy, P.L. Detection and analysis of endogenous badnaviruses in the New Zealand flora. AoB Plants 2011, 11, 1–13. [Google Scholar] [CrossRef] [PubMed]
  55. Short, S.M.; Suttle, C.A. Use of the polymerase chain reaction and denaturing gradient gel electrophoresis to study diversity in natural virus communities. Hydrobiologia 1999, 401, 19–32. [Google Scholar] [CrossRef]
  56. Lerner, A.; Shor, Y.; Vinokurov, A.; Okon, Y.; Jurkevitch, E. Can denaturing gradient gel electrophoresis (DGGE) analysis of amplified 16s rDNA of soil bacterial populations be used in forensic investigations? Soil Biol. Biochem. 2006, 38, 1188–1192. [Google Scholar] [CrossRef]
  57. Liu, B.; Lowes, F. Multiple Divergent ITS1 Copies Were Identified in Single Tomato Genome Using DGGE Analysis. Plant Mol. Biol. Report. 2013, 31, 272–279. [Google Scholar] [CrossRef]
  58. Riedel, G.E.; Swanberg, S.L.; Kuranda, K.D.; Marquette, K.; LaPan, P.; Bledsoe, P.; Kennedy, A.; Lin, B.Y. Denaturing gradient gel electrophoresis identifies genomic DNA polymorphism with high frequency in maize. Theor. Appl. Genet. 1990, 80, 1–10. [Google Scholar] [CrossRef] [PubMed]
  59. Yang, I.C.; Hafner, G.J.; Revill, P.A.; Dale, J.L.; Harding, R.M. Sequence diversity of South Pacific isolates of Taro bacilliform virus and the development of a PCR-based diagnostic test. Arch. Virol. 2003, 148, 1957–1968. [Google Scholar] [CrossRef] [PubMed]
  60. Lebas, B.S.M. Diversity of Viruses Infecting Dioscorea Species in the South Pacific. Ph.D. Thesis, University of Greenwich, Chatham Maritime, UK, 2002. [Google Scholar]
  61. Mumford, R.A.; Seal, S.E. Rapid single-tube immunocapture RT-PCR for the detection of two yam potyviruses. J. Virol. Methods 1997, 69, 73–79. [Google Scholar] [CrossRef]
  62. King, A.M.Q.; Adams, M.J.; Carstens, E.B.; Lefkowitz, E.J. Family—Caulimoviridae. In Virus Taxonomy; Elsevier: San Diego, CA, USA, 2012; pp. 429–443. ISBN 978-0-12-384684-6. [Google Scholar]
  63. Silva, G.; Bömer, M.; Nkere, C.; Lava Kumar, P.; Seal, S.E. Rapid and specific detection of Yam mosaic virus by reverse-transcription recombinase polymerase amplification. J. Virol. Methods 2015, 222, 138–144. [Google Scholar] [CrossRef] [PubMed]
  64. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  65. National Centre for Biotechnology Information (NCBI) GenBank. Available online: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genbank/ (accessed on 10 June 2017).
  66. Kimura, M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 1980, 16, 111–120. [Google Scholar] [CrossRef] [PubMed]
  67. Katoh, K.; Toh, H. Recent developments in the MAFFT multiple sequence alignment program. Brief. Bioinform. 2008, 9, 286–298. [Google Scholar] [CrossRef] [PubMed]
  68. Multiple Alignment Using Fast Fourier Transform (MAFFT). Available online: http://www.ebi.ac.uk/Tools/msa/mafft/ (accessed on 10 June 2017).
  69. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using CLUSTAL OMEGA. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef] [PubMed]
  70. Clustal Omega. Available online: http://www.ebi.ac.uk/Tools/msa/clustalo/ (accessed on 10 June 2017).
  71. Hall, T. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  72. Umber, M.; Gomez, R.-M.; Gélabale, S.; Bonheur, L.; Pavis, C.; Teycheney, P.-Y. 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. Arch. Virol. 2017, 162, 517–521. [Google Scholar] [CrossRef] [PubMed]
  73. Harper, G.; Hull, R. Cloning and sequence analysis of Banana streak virus. Virus Genes 1998, 17, 271–278. [Google Scholar] [CrossRef] [PubMed]
  74. Briddon, R.W.; Phillips, S.; Brunt, A.; Hull, R. Analysis of the sequence of Dioscorea alata bacilliform virus; comparison to other members of the badnavirus group. Virus Genes 1999, 18, 277–283. [Google Scholar] [CrossRef] [PubMed]
  75. Tijdens, M.; Hoogveld, H.L.; Kamst-Van Agterveld, M.P.; Simis, S.G.H.; Baudoux, A.C.; Laanbroek, H.J.; Gons, H.J. Population dynamics and diversity of viruses, bacteria and phytoplankton in a shallow eutrophic lake. Microb. Ecol. 2008, 56, 29–42. [Google Scholar] [CrossRef] [PubMed]
  76. Geering, A.D.W.; Olszewski, N.E.; Harper, G.; Lockhart, B.E.L.; Hull, R.; Thomas, J.E. Banana contains a diverse array of endogenous badnaviruses. J. Gen. Virol. 2005, 86, 511–520. [Google Scholar] [CrossRef] [PubMed]
  77. Ferris, M.J.; Muyzer, G.; Ward, D.M. Denaturing gradient gel electrophoresis profiles of 16S rRNA-defined populations inhabiting a hot spring microbial mat community. Appl. Environ. Microbiol. 1996, 62, 340–346. [Google Scholar] [PubMed]
  78. Bouhida, M.; Lockhartz, B.E.L.; Olszewski, N.E. An analysis of the complete sequence of a sugarcane bacilliform virus genome infectious to banana and rice. J. Gen. Virol. 1993, 74, 15–22. [Google Scholar] [CrossRef] [PubMed]
  79. Medberry, S.L.; Lockhart, B.E.L.; Olszewski, N.L. Properties of Commelina yellow mottle virus’s complete DNA sequence, genomic discontinuities and transcript suggest that it is a pararetrovirus. Nucleic Acids Res. 1990, 18, 5505–5513. [Google Scholar] [CrossRef] [PubMed]
  80. Seal, S.; Muller, E. Molecular analysis of a full-length sequence of a new yam badnavirus from Dioscorea sansibarensis. Arch. Virol. 2007, 152, 819–825. [Google Scholar] [CrossRef] [PubMed]
  81. Geering, A.D.W.; Parry, J.N.; Thomas, J.E. Complete genome sequence of a novel badnavirus, banana streak IM virus. Arch. Virol. 2011, 156, 733–737. [Google Scholar] [CrossRef] [PubMed]
  82. Muller, E.; Sackey, S. Molecular variability analysis of five new complete cacao swollen shoot virus genomic sequences. Arch. Virol. 2005, 150, 53–66. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graphical workflow describing badnavirus diversity studies using denaturing gradient gel electrophoresis (DGGE). PCR: polymerase chain reaction.
Figure 1. Graphical workflow describing badnavirus diversity studies using denaturing gradient gel electrophoresis (DGGE). PCR: polymerase chain reaction.
Viruses 09 00181 g001
Figure 2. Denaturing gradient gel electrophoresis (DGGE) analysis of partial reverse transcriptase-ribonuclease H (RT-RNaseH) badnavirus sequences from seven Dioscorea rotundata breeding lines and one landrace comparing PCR amplifications using the generic badnavirus primer pair Badna-forward primer/reverse primer (FP/-RP) with a GC clamp fused to the forward (A) or reverse (B) primer. PCR conditions were as outlined in the Materials and Methods section, but using undiluted DNA extractions (‘direct’) and only 25 cycles. PCR amplifications (5 µL) were checked for their expected size (619 bp) on 1% w/v agarose in 0.5× Tris-Boric acid-ethylenediaminetetraacetic acid (EDTA) (TBE) gels (black gels on top) before DGGE loading (20 µL). The denaturing gradient was 30–55% and DGGE was performed at 80 V at a temperature of 60 °C for 18 h. Band numbers 1–21 were excised and cloned. The corresponding sequences are presented in Table 1. Dioscorea rotundata accession TDr 89/02475A and B and TDr 1892A and B are clones of the same yam accessions.
Figure 2. Denaturing gradient gel electrophoresis (DGGE) analysis of partial reverse transcriptase-ribonuclease H (RT-RNaseH) badnavirus sequences from seven Dioscorea rotundata breeding lines and one landrace comparing PCR amplifications using the generic badnavirus primer pair Badna-forward primer/reverse primer (FP/-RP) with a GC clamp fused to the forward (A) or reverse (B) primer. PCR conditions were as outlined in the Materials and Methods section, but using undiluted DNA extractions (‘direct’) and only 25 cycles. PCR amplifications (5 µL) were checked for their expected size (619 bp) on 1% w/v agarose in 0.5× Tris-Boric acid-ethylenediaminetetraacetic acid (EDTA) (TBE) gels (black gels on top) before DGGE loading (20 µL). The denaturing gradient was 30–55% and DGGE was performed at 80 V at a temperature of 60 °C for 18 h. Band numbers 1–21 were excised and cloned. The corresponding sequences are presented in Table 1. Dioscorea rotundata accession TDr 89/02475A and B and TDr 1892A and B are clones of the same yam accessions.
Viruses 09 00181 g002
Figure 3. DGGE analysis of partial RT-RNaseH badnavirus sequences from 19 samples consisting of first filial (F1) generations of West African breeding lines (o.p. stands for open pollinated) of D. rotundata comparing patterns of PCR amplifications (20 µL loaded) using the generic badnavirus primer pair Badna-FP/-RP with a GC clamp fused to the forward primer. The denaturing gradient was 35–50% and DGGE was performed at 80 V at a temperature of 60 °C for 18 h. Band numbers 22–42 were excised and cloned. The corresponding sequences are presented in Table 1.
Figure 3. DGGE analysis of partial RT-RNaseH badnavirus sequences from 19 samples consisting of first filial (F1) generations of West African breeding lines (o.p. stands for open pollinated) of D. rotundata comparing patterns of PCR amplifications (20 µL loaded) using the generic badnavirus primer pair Badna-FP/-RP with a GC clamp fused to the forward primer. The denaturing gradient was 35–50% and DGGE was performed at 80 V at a temperature of 60 °C for 18 h. Band numbers 22–42 were excised and cloned. The corresponding sequences are presented in Table 1.
Viruses 09 00181 g003
Figure 4. Bootstrap consensus phylogenetic tree using the maximum likelihood method built from the badnavirus 527–528-bp-long partial RT-RNaseH nucleotide sequences of 114 yam badnavirus sequences determined in this study. Included in the analysis are partial RT-RNaseH sequences with names and accession numbers from GenBank of previously analysed yam samples by Bömer et al. [33], Bousalem et al. [34], Eni et al. [35], Kenyon et al. [3], Seal et al. [38] and Umber et al. [39,72]. Equivalent sequences from Cacao swollen shoot virus (CSSV, AJ781003), Banana streak OL virus (BSOLV, AJ002234), Commelina yellow mottle virus (ComYMV, NC001343), Sugarcane bacilliform MO virus (SCBMOV, M89923), Taro bacilliform virus (TaBV, AF357836) and outgroup Rice tungro bacilliform virus (RTBV, X57924) were added, as well as representative sequences of all monophyletic groups described by Bousalem et al. [34] (where DBV-D: Dioscorea bacilliform virus D), by Umber et al. [39] and by Kenyon et al. [3], denoted by U12 and K1–K11 respectively. Three novel monophyletic groups, T13–15, described by Bömer et al. [33] were also included. Sequences depicted in bold represent partial RT-RNaseH sequences of characterised episomal full-length DBV genomes currently available in GenBank. The phylogenetic tree was divided into sub-groups with groups K8 and U12 presented in (A); K9 and T13 shown in (B) and DBV-D, K1–11, T14, T15 as well as the outgroups shown in (C). The bootstrap analysis of the sequences was 1000 replicates and the cut-off value was 70%. The scale bars show the number of substitutions per base.
Figure 4. Bootstrap consensus phylogenetic tree using the maximum likelihood method built from the badnavirus 527–528-bp-long partial RT-RNaseH nucleotide sequences of 114 yam badnavirus sequences determined in this study. Included in the analysis are partial RT-RNaseH sequences with names and accession numbers from GenBank of previously analysed yam samples by Bömer et al. [33], Bousalem et al. [34], Eni et al. [35], Kenyon et al. [3], Seal et al. [38] and Umber et al. [39,72]. Equivalent sequences from Cacao swollen shoot virus (CSSV, AJ781003), Banana streak OL virus (BSOLV, AJ002234), Commelina yellow mottle virus (ComYMV, NC001343), Sugarcane bacilliform MO virus (SCBMOV, M89923), Taro bacilliform virus (TaBV, AF357836) and outgroup Rice tungro bacilliform virus (RTBV, X57924) were added, as well as representative sequences of all monophyletic groups described by Bousalem et al. [34] (where DBV-D: Dioscorea bacilliform virus D), by Umber et al. [39] and by Kenyon et al. [3], denoted by U12 and K1–K11 respectively. Three novel monophyletic groups, T13–15, described by Bömer et al. [33] were also included. Sequences depicted in bold represent partial RT-RNaseH sequences of characterised episomal full-length DBV genomes currently available in GenBank. The phylogenetic tree was divided into sub-groups with groups K8 and U12 presented in (A); K9 and T13 shown in (B) and DBV-D, K1–11, T14, T15 as well as the outgroups shown in (C). The bootstrap analysis of the sequences was 1000 replicates and the cut-off value was 70%. The scale bars show the number of substitutions per base.
Viruses 09 00181 g004aViruses 09 00181 g004b
Table 1. Basic local alignment search tool (BLAST) analysis of partial RT-RNaseH sequences cloned from DGGE bands.
Table 1. Basic local alignment search tool (BLAST) analysis of partial RT-RNaseH sequences cloned from DGGE bands.
Plant Accession aDGGE Sequence bAccessionPrimers cSize (bp)NCBI Nearest MatchIdentity (%)Species Group d
TDr 89/02475NGb1aDrKY555456BF-GC + BR528eDBV9_S1h6Dr (KF829977)99K9
TDr 89/02475NGb1bDrKY555457BF-GC + BR528432B39Ds (AM503361)99K9
TDr 89/02475NGb2aDrKY555458BF-GC + BR528432B39Ds (AM503361)100K9
TDr 89/02475NGb2bDrKY555459BF-GC + BR528432B39Ds (AM503361)99K9
TDr 89/02475NGb3aDrKY555460BF-GC + BR527432B39Ds (AM503361)99K9
TDr 89/02475NGb3bDrKY555461BF-GC + BR528432B39Ds (AM503361)99K9
TDr 89/02475NGb4DrKY555462BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 89/02475NGb5DrKY555463BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 89/02475NGb6aDrKY555464BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb6bDrKY555465BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 89/02475NGb7aDrKY555466BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb7bDrKY555467BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb8aDrKY555468BF-GC + BR528DBRTV1-[3RT] (KX008598)99T13
TDr 89/02475NGb8bDrKY555469BF-GC + BR528GyJT2Dt (AM503389)92K8
PonaNGl9aDrKY555470BF-GC + BR528eDBV12_S1a4Dr (KF829956)99U12
PonaNGl9bDrKY555471BF-GC + BR528eDBV12_S1a4Dr (KF829956)99U12
PonaNGl10aDrKY555472BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl10bDrKY555473BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl11aDrKY555474BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl11bDrKY555475BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl12aDrKY555476BF + BR-GC528eDBV12_S1a4Dr (KF829956)99U12
PonaNGl12bDrKY555477BF + BR-GC528eDBV12_S1a4Dr (KF829956)100U12
PonaNGl13aDrKY555478BF + BR-GC528eDBV12_S1a4Dr (KF829956)100U12
PonaNGl13bDrKY555479BF + BR-GC528eDBV12_S1a4Dr (KF829956)99U12
PonaNGl14aDrKY555480BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl14bDrKY555481BF + BR-GC528BfA103Dc (AM503393)100K8
PonaNGl15aDrKY555482BF + BR-GC528BfA103Dc (AM503393)99K8
PonaNGl15bDrKY555483BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 89/02475NGb16aDrKY555484BF + BR-GC528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb16bDrKY555485BF + BR-GC527432B39Ds (AM503361)99K9
TDr 89/02475NGb17aDrKY555486BF + BR-GC528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb17bDrKY555487BF + BR-GC528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb18aDrKY555488BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 89/02475NGb18bDrKY555489BF + BR-GC528GyJT2Dt (AM503389)93K8
TDr 89/02475NGb19aDrKY555490BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 89/02475NGb19bDrKY555491BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 89/02475NGb20aDrKY555492BF + BR-GC528GyJT2Dt (AM503389)92K8
TDr 89/02475NGb20bDrKY555493BF + BR-GC528GyJT2Dt (AM503389)93K8
TDr 95/18544NGb21aDrKY555494BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 95/18544NGb21bDrKY555495BF + BR-GC528BfA103Dc (AM503393)99K8
TDr 95/18544 (o.p.)NGb22aDrKY555496BF-GC + BR528eDBV5_S1g6Dr (KF829974)99K5
TDr 95/18544 (o.p.)NGb22bDrKY555497BF-GC + BR528eDBV5_S1g6Dr (KF829974)99K5
TDr 95/18544 (o.p.)NGb23aDrKY555498BF-GC + BR528eDBV5_S1g6Dr (KF829974)99K5
TDr 95/18544 (o.p.)NGb23bDrKY555499BF-GC + BR528BfA103Dc (AM503393)97K8
TDr 95/18544 (o.p.)NGb24aDrKY555500BF-GC + BR528eDBV5_S1g6Dr (KF829974)99K5
TDr 95/18544 (o.p.)NGb24bDrKY555501BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 97/00917 × TDr 99/02607NGb25aDrKY555502BF-GC + BR528eDBV9_S1h2Dr (KF829975)99K9
TDr 97/00917 × TDr 99/02607NGb25bDrKY555503BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb26aDrKY555504BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb26bDrKY555505BF-GC + BR528eDBV9_G1Dr (KF830002)99K9
TDr 97/00917 × TDr 99/02607NGb27aDrKY555506BF-GC + BR528eDBV12_S2a7Dr (KF829978)99U12
TDr 97/00917 × TDr 99/02607NGb27bDrKY555507BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb28aDrKY555508BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb28bDrKY555509BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb29aDrKY555510BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb29bDrKY555511BF-GC + BR528eDBV9_S2f8Dr (KF829993)99K9
TDr 97/00917 × TDr 99/02607NGb30aDrKY555512BF-GC + BR527BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb30bDrKY555513BF-GC + BR528BfA103Dc (AM503393)100K8
TDr 99/02793 × TDr 1892NGb31DrKY555514BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 99/02793 × TDr 1892NGb32DrKY555515BF-GC + BR528eDBV9_S1e3Dr (KF829969)100K9
TDr 04/219 × TDr 98/02677NGb33aDrKY555516BF-GC + BR527GyJT2Dt (AM503389)92K8
TDr 04/219 × TDr 98/02677NGb33bDrKY555517BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 97/00917 (o.p.)NGb34aDrKY555518BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 97/00917 (o.p.)NGb34bDrKY555519BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 96/00629 × TDr 99/02607NGb35aDrKY555520BF-GC + BR528432B39Ds (AM503361)99K9
TDr 96/00629 × TDr 99/02607NGb35bDrKY555521BF-GC + BR528432B39Ds (AM503361)100K9
TDr 97/00917 × TDr 99/02607NGb36aDrKY555522BF-GC + BR528BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb36bDrKY555523BF-GC + BR527BN4Dr (AM944586)99K9
TDr 97/00917 × TDr 99/02607NGb37aDrKY555524BF-GC + BR528eDBV9_S2f8Dr (KF829993)99K9
TDr 97/00917 × TDr 99/02607NGb37bDrKY555525BF-GC + BR528eDBV9_S1h2Dr (KF829975)99K9
TDr 04/219 × TDr 04/219NGb38aDrKY555526BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 04/219 × TDr 04/219NGb38bDrKY555527BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 04/219 × TDr 04/219NGb39aDrKY555528BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 04/219 × TDr 04/219NGb39bDrKY555529BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 04/219 × TDr 04/219NGb40DrKY555530BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 97/00205 ×TDr 1892NGb41aDrKY555531BF-GC + BR528eDBV12_S1a4Dr (KF829956)99U12
TDr 97/00205 ×TDr 1892NGb41bDrKY555532BF-GC + BR528eDBV9_S1h2Dr (KF829975)99K9
TDr 96/00629 × TDr 99/02607NGb42aDrKY555533BF-GC + BR528eDBV12_S1a4Dr (KF829956)99U12
TDr 96/00629 × TDr 99/02607NGb42bDrKY555534BF-GC + BR528432B39Ds (AM503361)99K9
TDr 95/18544NGb43DrKY555535BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 95/18544NGb44DrKY555536BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 95/18544NGb45DrKY555537BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 96/00629NGb46aDrKY555538BF-GC + BR528eDBV12_S2a7Dr (KF829978)99U12
TDr 96/00629NGb46bDrKY555539BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 89/02677NGb47aDrKY555540BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02677NGb47bDrKY555541BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 89/02677NGb48aDrKY555542BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02677NGb48bDrKY555543BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 89/02677NGb49DrKY555544BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 96/00629NGb50DrKY555545BF-GC + BR528eDBV12_S2a7Dr (KF829978)99U12
TDr 89/02677NGb51DrKY555546BF-GC + BR528GyJT2Dt (AM503389)91K8
TDr 97/00917 (o.p.)NGb52DrKY555547BF-GC + BR528eDBV5_S1g6Dr (KF829974)98K5
TDr 97/00917 (o.p.)NGb53DrKY555548BF-GC + BR528NC1 (KJ854414)79K13
TDr 97/00917 × TDr 99/02607NGb54DrKY555549BF-GC + BR528eDBV9_S1h2Dr (KF829975)99K9
TDr 99/02793 (o.p.)NGb55aDrKY555550BF-GC + BR528BfA103Dc (AM503393)100K8
TDr 99/02793 (o.p.)NGb55bDrKY555551BF-GC + BR528BfA103Dc (AM503393)99K8
TDr 96/00621 (o.p.)NGb56DrKY555552BF-GC + BR528eDBV12_S1a4Dr (KF829956)100U12
TDr 97/00917 × TDr 97/00777NGb57DrKY555553BF-GC + BR528GyJT2Dt (AM503389)93K8
TDr 96/00629 × TDr 99/02607NGb58DrKY555554BF-GC + BR528GyJT2Dt (AM503389)92K8
TDr 97/00917 × TDr 99/02607NGb59DrKY555555BF-GC + BR528eDBV12_S2h10Dr (KF829998)99U12
TDa 99/00240 × TDa 95/00310NGb60DaKY555556BF-GC + BR528GyJT2Dt (AM503389)92K8
TDa 99/00240 × TDa 01/00012NGb61DaKY555557BF-GC + BR528NG1Da (AM944571)93K8
TDa 01/00081 × TDa 02/00012NGb62DaKY555558BF-GC + BR528GyJT2Dt (AM503389)93K8
TDa 1/00081×TDa 98/00150NGb63DaKY555559BF-GC + BR528SB42Da (AM072696)99K1
TDa 00/00194 × TDa 98/00150NGb64DaKY555560BF-GC + BR528eDBV9_S1b6Dr (KF829960)89K9
TDe 3049ANGl65DeKY555561BF-GC + BR528WS31aDn (AM421696)73K12
TDd 4118BNGl66DdKY555562BF-GC + BR528eDBV5_S1un5Dr (KF830000)93K5
TDd 4118BNGl67DdKY555563BF-GC + BR528BN2Da (AM944584)98K8
TDr 1950BNGl68DrKY555564BF-GC + BR528GN2Dr (AM944575)96K8
TDc 3808CNGl69DcKY555565BF-GC + BR528eDBV9_S1b2Dr (KF829958)100K9
TDb 3045BNGl70DbKY555566BF-GC + BR528SB42Da (AM072696)99K1
TDd 3778BNGl71DdKY555567BF-GC + BR528FJ65bDe (AM072660)99K1
TDa 1013CNGl72DaKY555568BF-GC + BR528SB42Da (AM072696)99K1
TDe 3028ANGl73DeKY555569BF-GC + BR528FJ65bDe (AM072660)99K1
a The host plants are represented by plant accession. TDa: Dioscorea alata accession; TDb: Dioscorea bulbifera accession; TDc: Dioscorea cayenensis accession; TDd: Dioscorea dumetorum accession; TDe: Dioscorea esculenta accession; TDr: D. rotundata accession; o.p.: open pollinated; b The DGGE clone sequences were coded as follows: the first two letters stand for the country of origin (NG = Nigeria), ‘b’ represents breeding line samples, ‘l’ represents landrace yam samples, the middle number denotes the position of the excised DGGE band, the next letter denotes the clone (a = clone a and b = clone b) and the last two letters refer to the Dioscorea host species (e.g., Dr = Dioscorea rotundata); c BF: Badna FP; BR = Badna-RP; GC: GC clamp (see Methods for details); eDBV: endogenous Dioscorea bacilliform viruses. d According to phylogenetic tree (Figure 4).

Share and Cite

MDPI and ACS Style

Turaki, A.A.; Bömer, M.; Silva, G.; Kumar, P.L.; Seal, S.E. PCR-DGGE Analysis: Unravelling Complex Mixtures of Badnavirus Sequences Present in Yam Germplasm. Viruses 2017, 9, 181. https://0-doi-org.brum.beds.ac.uk/10.3390/v9070181

AMA Style

Turaki AA, Bömer M, Silva G, Kumar PL, Seal SE. PCR-DGGE Analysis: Unravelling Complex Mixtures of Badnavirus Sequences Present in Yam Germplasm. Viruses. 2017; 9(7):181. https://0-doi-org.brum.beds.ac.uk/10.3390/v9070181

Chicago/Turabian Style

Turaki, Aliyu A., Moritz Bömer, Gonçalo Silva, P. Lava Kumar, and Susan E. Seal. 2017. "PCR-DGGE Analysis: Unravelling Complex Mixtures of Badnavirus Sequences Present in Yam Germplasm" Viruses 9, no. 7: 181. https://0-doi-org.brum.beds.ac.uk/10.3390/v9070181

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop