Skip to main content

The proposed new species, cacao red vein virus, and three previously recognized badnavirus species are associated with cacao swollen shoot disease

Abstract

Background

Cacao swollen shoot virus (CSSV), Cacao swollen shoot CD virus (CSSCDV), and Cacao swollen shoot Togo A virus (CSSTAV) cause cacao swollen shoot disease (CSSD) in West Africa. During 2000–2003, leaf and shoot-swelling symptoms and rapid tree death were observed in cacao in Cote d’Ivoire and Ghana. Molecular tests showed positive infection in only ~50–60% of symptomatic trees, suggesting the possible emergence of an unknown badnavirus.

Methods

The DNA virome was determined from symptomatic cacao samples using Illumina-Hi Seq, and sequence accuracy was verified by Sanger sequencing. The resultant 14, and seven previously known, full-length badnaviral genomic and RT-RNase H sequences were analyzed by pairwise distance analysis to resolve species relationships, and by Maximum likelihood (ML) to reconstruct phylogenetic relationships. The viral coding and non-coding sequences, genome organization, and predicted conserved protein domains (CPDs) were identified and characterized at the species level.

Results

The 21 CSSD-badnaviral genomes and RT-RNase H sequences shared 70–100% and 72–100% identity, respectively. The RT-RNase H analysis predicted four species, based on an ≥80% species cutoff. The ML genome sequence tree resolved three well-supported clades, with ≥70% bootstrap, whereas, the RT-RNase H phylogeny was poorly resolved, however, both trees grouped CSSD isolates within one large clade, including the newly discovered Cacao red vein virus (CRVV) proposed species. The genome arrangement of the four species consists of four, five, or six predicted open reading frames (ORFs), and the CPDs have similar architectures. By comparison, two New World cacao-infecting badnaviruses encode four ORFs, and harbor CPDs like the West African species.

Conclusions

Three previously recognized West African cacao-infecting badnaviral species were identified, and a fourth, previously unidentified species, CRVV, is described for the first time. The CRVV is a suspect causal agent of the rapid decline phenotype, however Koch’s Postulates have not been proven. To reconcile viral evolutionary with epidemiology considerations, more detailed information about CSSD-genomic variability is essential. Also, the functional basis for the multiple genome arrangements and subtly distinct CPD architectures among cacao-infecting badnaviruses is poorly understood. New knowledge about functional relationships may help explain the diverse symptomatologies observed in affected cacao trees.

Background

The Theobroma cacao (L.) tree, or cacao, is the source of the chocolate bean. The species is endemic to the Amazon Basin in South America and was introduced into West Africa during the 1880’s [1], where 70% of the world’s bulk cocoa supply is now produced. Soon after commercial cacao plantations were established, virus-like symptoms were reported in cacao trees consisting of foliar vein-banding, swellings on vegetative shoots (chupons), and reduced pod size, pod number, and quality of beans [2, 3]. Three to five years post-symptom development, infected trees decline and eventually die [4, 5]. Graft transmission was reported by Steven in 1936 [3], leading to the hypothesis that the swollen shoot disease was caused by a plant virus, endemic to West Africa, referred to as Cacao swollen shoot virus (CSSV) [4, 6]. The virus is transmitted in a semi-persistent manner by 14 mealybug species [7], and is not known to be seed-transmitted in cacao [8]. Efforts to develop virus-resistant cacao genotypes have been largely unsuccessful, and so over time, crop losses have been substantial in Cote d’Ivoire [9, 10], Ghana [11,12,13], Nigeria [14, 15], Sierra Leone [16], and Togo [17, 18]. Until now, CSSV and two other cacao-infecting virus species, Cacao swollen shoot CD virus (CSSCDV), and Cacao swollen shoot Togo A virus (CSSTAV), have been associated with cacao swollen shoot disease (CSSD) [19], and are members of the genus Badnavirus (family, Caulimoviridae) [20,21,22,23], also referred to as pararetroviruses.

Badnaviruses have a circular, double-stranded DNA genome 7.0–9.2 kilobase pairs (kbp) in size, encapsidated in a non-enveloped bacilliform particle [22, 23]. Replication occurs by a viral-encoded reverse transcriptase (RT), and proceeds through an RNA intermediate [24]. The three CSSD-associated genomes encode four to six open reading frames (ORFs), referred to as ORFs 1–4, X, and Y [10, 25]. The precise function of the badnaviral 16 kDa ORF1 protein is not known. The ORF2 encodes a 15 kDa protein with DNA and RNA binding activity [26], and ORF3 encodes a ~212 kDa polyprotein containing several domains that are cleaved to release a functional movement protein (MP), coat protein (CP), aspartic protease (AP), and the RT and ribonuclease H (RNase H) proteins. The ORFs 4, X, and Y, which overlap ORF3, encode 95, 13, and 14 kDa proteins, respectively, all of unknown function [10, 27].

During the early 2000’s to the present, characteristic leaf and shoot swelling CSSD symptoms, accompanied by rapid tree decline and death about one year after the occurrence of initial foliar symptoms, have been observed in a portion of the trees growing in commercial plantations located in western Ghana and eastern Cote d’Ivoire [25; authors’ personal observations]. Concurrently, previously reliable serological and several polymerase chain reaction (PCR) amplifications failed to detect virus in 40–60% of samples, despite presence of CSSD-like symptoms [2, 18, 28, 29]. Also, primers designed to amplify a hypervariable region of the viral MP (located on ORF3), referred to as the ‘ORF3A’ primers, reported to distinguish nine CSSD-badnavirus variants, collectively failed to detect virus in 25–50% of samples from Cote d’Ivoire and Ghana [10, 18, 28, 30]. Even though this MP locus is not phylogenetically informative at the species level [20], the region is informative of extensive intraspecific variability. Similarly, in a recent study of 124 field isolates from Cote d’Ivoire, only half of the samples were positive for CSSD badnavirus infection by PCR amplification [31]. The primers were designed based on the seven CSSD-associated genome sequences from Cote d’Ivoire, Ghana and Togo available in GenBank, to direct amplification of CP, MP, RNase H, RT, and non-coding region fragments. While one or more primer pairs amplified badnavirus amplicon(s), confirmed as badnavirus by DNA sequencing, overall ~50% of symptomatic leaf and shoot samples yielded no PCR product, indicating greater than expected genomic variability among CSSD badnaviruses [28]. Until now, there has been no systematic study of the extent of genomic variability of CSSD badnaviruses in West Africa.

The lack of sufficient representative full-length CSSD-associated genome sequences has precluded a reconciliation of evolutionary and epidemiological information required to inform disease management practices. The objective of this study was to better understand the extent of genomic variability of CSSD-badnaviruses in Cote d’Ivoire and Ghana associated with swollen shoot symptoms, including the recently observed ‘rapid decline’ phenotype. Here, a genomic pathology approach was taken utilizing Illumina Hi-Seq DNA sequencing for ‘virome discovery’, and validation by PCR amplification with virus-specific primers and Sanger DNA sequencing. The 14 apparently full-length genomes were characterized with respect to pairwise nucleotide (nt) identity, phylogenetic relationships, genome organization, and conserved protein domain architecture, in relation to the seven previously reported CSSD-badnaviral genome sequences from Ghana, Cote d’Ivoire and Togo.

Methods

Plant samples, DNA isolation, and rolling circle amplification

Sixty-five leaf or shoot samples were collected from symptomatic cacao trees in commercial plantations, or at experimental plots located at the Centre National de Recherche Agronomique in Cote d’Ivoire during 2012. An additional thirty-five samples were collected in 2015 from commercial plantations or the historical collection maintained at the Cocoa Research Institute of Ghana. Each sample was placed in a 50-ml screwcap tube containing 100% glycerol, transported to the University of Arizona, and held at 4 °C. Total DNA was isolated from 100 mg of plant tissue from all samples using the cetyltrimethyl ammonium bromide (CTAB) DNA extraction method of Doyle and Doyle [32]. Purified DNA (2 μl) was enriched for circular DNA using the Templiphi rolling circle amplification (RCA) kit (GE Healthcare Bio-Sciences), according to manufacturer’s instructions, with modifications [33, 34].

Illumina sequencing and assembly of paired-end reads

Illumina paired-end libraries were constructed from purified DNA or from RCA products using the TruSeq PE cluster kit, to obtain a 350 bp mean insert size. The Illumina HiSeq 2500 sequencing run was carried out at the University of Arizona Genetics Core facility (Tucson, AZ). The reads were de-multiplexed, and quality was assessed using FASTQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). The adapters were removed, and reads were trimmed using TRIMMOMATIC v0.32 [35] and assembled de novo using SeqMan NGen software v12 (DNASTAR), with and without plant host sequence filtering. For filtering, the T. cacao nuclear genome (host plant) contigs were downloaded from GenBank Accessions CM001879.1, CM001888.1, FR7222157.1, and KE132922.1 [36, 37]. Also used for filtering were the cacao chloroplast genome sequence, NC014676.2 and cotton Gossypium hirsutum (L.) (cacao relative) mitochondria genome, JX065074.1.

Annotation of badnavirus-like genomes

The assembled apparently full-length genome sequences were first annotated using BLAST2GO [38] to confirm their identity as badna-like viruses, and then subjected to BLASTn [39] analysis to determine the top hits (similarity scores) based on all sequences available in the Genbank database. Based on established convention for the genus Badnavirus, the first nt of the predicted tRNAmet priming site was assigned as the first nt coordinate [40, 41] for each full-length CSSD-badnaviral genome sequence. The predicted ORFs were identified for CSSD genome sequences using the NCBI ORF Finder tool, and the conserved protein domains (CPD) were predicted using the NCBI Conserved Domain Database (CDD) [42]. The predicted CPDs were aligned for comparisons using MUSCLE [43], implemented in CLC Sequence Viewer 7.5.

Pairwise and phylogenetic analyses of coding regions and complete genome sequences

Genome sequences for 122 badnavirus isolates representing 37 species were downloaded from GenBank and the genome was considered intact or based on the RT-RNase H region used for taxonomic assignment to species. A haplotype search implemented in FaBox v1.41 was used to identify and remove sequences sharing 100% nt identity [44], leaving 87 full-length genome and RT-RNase H region sequences for the analyses, respectively. Pairwise distances were determined using Standard Demarcation Tool (SDTv1.2) software [43], for 87 full-length genomes or RT-RNase H regions (1230 bp). The sequence alignments were carried out in MUSCLE [43], implemented in CLC Sequence Viewer 7.5.

Phylogenetic analysis of genomic and RT-RNase H sequences was carried out with the Maximum Likelihood (ML) algorithm available in MEGA6 [45]. Tree reconstructions used the General Time Reversible substitution model, determined by the highest Bayesian Information Criterion score obtained, with gamma distribution for invariable sites, and 1000 bootstrap iterations. Confidence values were placed at major nodes having a ≥ 70% bootstrap.

Validation of CSSD-associated badnaviral genome sequences by sanger sequencing

To validate the de novo-assembled Illumina full-length genome sequences, six representative isolates, three each from Cote d’Ivoire and Ghana, representing Cacao vein virus (CRVV) and the established CSSV species group (four and two sequences, respectively) (Table 1), were selected for PCR amplification, cloning, and bi-directional capillary (Sanger) DNA sequencing. Purified total DNA was enriched for circular DNA by RCA. The partially enriched badnaviral DNA was used as template for PCR amplification of full-length viral genomes with sequence-specific abutting primers (Additional file 1: Table S1). Amplicons were cloned using the CloneAmp™ HiFi PCR Premix (Clonetech), according to the manufacturer’s instructions. The reaction contained 1X CloneAmp™ HiFi PCR Premix, 0.2 μM each of reverse and forward primers, 2 μL of RCA product, and nuclease-free water to a final volume of 50 μL. The PCR conditions were: denaturation at 98 °C for 2 min, followed by 40 cycles of denaturation at 98 °C for 20 s, annealing at 55 °C for 15 s, extension at 72 °C for 8 min, and final extension at 72 °C for 10 min. Amplicons were separated by agarose gel (0.8%) electrophoresis in TAE buffer, pH 8.0. Bands ~7 kbp in size were excised and gel-purified using the illustra GFX PCR and DNA Gel Band Purification kit (GE Healthcare Bio-Sciences), according to manufacturer’s instructions. The concentration of DNA was determined using the Nanodrop 2000 UV–vis spectrophotometer (Thermo Scientific). Amplicons were cloned into the pGEM5 plasmid vector (Promega), previously linearized by digestion with Not I (New England Biolabs). Transformation was carried out using the In-Fusion HD Cloning kit (Clontech), according to the manufacturer’s instructions. Insert size was confirmed following purification of the plasmid vector with the GeneJET Plasmid Miniprep Kit (ThermoFisher Scientific), and by restriction with Not I. The DNA sequence was determined for plasmids bearing a cloned insert of ≥7 kbp in size by bi-directional, capillary Sanger DNA sequencing, and primer walking (Eton Bioscience, San Diego, CA). Sequences were assembled using SeqMan Pro v.12 (DNASTAR, Madison, WI). Viral ORFs, non-coding regions, and predicted, CPDs were identified for comparison to the assembled Illumina sequences using ORF Finder, CDD, and pairwise distance analysis (SDT), as described above.

Table 1 Predicted conserved protein domains identified in cacao-infecting badnavirus genome sequences

Results

Illumina and sanger sequences

The Illumina reads were assembled into 14 full-length badnaviral genomes ranging from 6920 to 7172 bp in size. Also, BLASTn analysis against the GenBank database identified a large number of partial contigs with significant matches to CSSD-badnaviral sequences. However, in this report, only complete badnaviral genome sequences were considered, and they are derived from 410,662–31,417,994 reads with a depth of coverage of 36–6600.

PCR amplification was carried out using sequence-specific primers (Additional file 1: Table S1) and Sanger sequencing of six full-length genome sequences, representing two CSSV and three CRVV full-length genomes, and one partial genome assembled from CRVV isolate CI286. Illumina and Sanger sequence comparisons indicated that the cloned viral amplicons were nearly identical in length and sequence, at 99.0–99.9% nt identity, between the respective sequences. BLASTn analysis for the 14 apparently full-length CSSD-badnaviral genome sequences indicated the top hits in the GenBank database were consistent to previously reported CSSD-badnaviral sequences based on similarity score, sequence coverage, and e-value of zero. A comparison of the resultant Next-gen Illumina and Sanger DNA sequences using the above criteria, as well as: location of the first nt coordinate, number and location of predicted ORFs ≥10 kDa, and presence of the hallmark-badnaviral non-coding sequence region, indicated that both approaches yielded the same genome sequence. The Illumina sequences have been deposited in GenBank as the Accession numbers, KX592571-KX592584.

Genome characterization

Identification of predicted coding regions (>50 aa) on the CSSD-badnaviral sense strand for the 14 CSSD-genomes revealed three types of genome arrangements, arbitrarily herein referred to as I, II, and III (Fig. 1). Samples CI275, CI286, GH64, and GH67, representing isolates from Cote d’Ivoire and from Ghana, had a type I arrangement, which consisted of four predicted ORFs, 1–3 and Y (Fig. 1a). Seven isolates, CI44, CI134, CI135, CI215, CIS2, CIS3, and CIT5, all from Cote d’Ivoire, had a type II genome arrangement, comprising five predicted ORFs, 1–3, X, and Y (Fig. 1b). And, CI301, CI311, and GH75, which represent isolates from Cote d’Ivoire and Ghana, had a type III arrangement, or six ORFs 1–4, X, and Y (Fig. 1c). The CSSD-associated genome sequences each contained a predicted tRNAmet binding site, located at nt coordinates 1–18, at which reverse transcription is initiated [24, 46].

Fig. 1
figure 1

The three types of genome arrangements found among the cacao swollen shoot- badnaviral genome sequences. The predicted coding regions are indicated by filled arrows as: blue = ORF1, orange = ORF2, green = ORF3, pink = ORF4, purple = ORFX, and black = ORFY. a The four ORF-genome arrangement (type I), shown for isolate GH67, was also observed for isolates CI275, GH64 and CI286. b The five ORF-genome arrangement (type II), shown for isolate CIS2, was present among isolates CIS3, CI134, CI135, CIT5, CI44 and CI215. c The six ORF genome arrangement (type III), represented by isolate GH75, was present in isolates CI311 and CI301. The small black arrowhead at the top of each genome map indicates the first nucleotide coordinate, which also is the first 5′-nt of the tRNAmet binding site. Country name and abbreviation: CI (Cote d’Ivoire), GH (Ghana), and TG (Togo) indicating country of sample collection

Conserved protein domains

Predicted CPD searches for the 21 genomes indicated a number of similarities, but also discernable differences, among CSSD-species and well-studied non-cacao-associated badnaviruses [10, 25, 47]. With respect to all of the CSSD-badnavirus species, CPDs were identified on ORFs 1 and 3. The ORF1 consistently harbored a domain of unknown function (DUF) 1319 [42], classified as a member of a badnavirus-specific, uncharacterized protein superfamily (Fig. 2a, Table 1). The ORF3 CPDs were identified as domains of zinc knuckle finger (Zn), pepsin-like aspartate protease (Pepsin), RT, and RNase H proteins (Table 1) present in all genomes except for CSSCDV [28], which lacked the Zn domain. Domains that were shared among all cacao-infecting badnaviruses were DUF1319 in ORF1, and Zn, Pepsin, RT and RNase H domains on ORF3. No CPDs were identified for CSSD-badnaviral species in the ORFs 2, 4, X, or Y, except for CSSCDV ORFY, which harbored a DUF3187 domain.

Fig. 2
figure 2

Predicted, conserved protein domains (CPDs) and corresponding amino acid sequence alignments. a Linear maps (not to scale) show the three known genome arrangements for the West African cacao-infecting badnaviruses, with either four (blue), five (black), or six (green) open reading frames (ORFs), compared with cacao-infecting badnaviruses from Trinidad, Cacao mild mosaic virus (CaMMV) (gray box), and Cacao yellow-vein banding virus (CYVBV) (red box), each with four ORFs. Abbreviation with virus name: CSSV = Cacao swollen shoot virus; CSSCDV = Cacao swollen shoot CD virus; CSSTAV = Cacao swollen shoot Togo A virus; CRVV = Cacao red vein virus; DUF = domain of unknown function; Zn = zinc knuckle finger; Pepsin = pepsin-like aspartate protease; RT = reverse transcriptase; RNase H = ribonuclease H; PHD = plant homeodomain finger; and Trim = trimeric-dUTPase. b Amino acid alignment of the CPDs for four representative CSSD-badnaviral genomes from West Africa, and CaMMV and CYVBV from Trinidad. Amino acid residues that differed among most or all of the cacao-infecting badnaviral isolates are shaded. The ‘consensus’ indicates the amino acids (aa) that are 100% conserved in CPDs, across the six badnaviral species, and excludes non-conserved residues, indicated by the “–“. The black vertical bars indicate the level of aa residue conservation, ranging from 0 to 100%, for the six representative genome arrangements with CPD architecture

A survey of CPDs for the two New World cacao-infecting species, Cacao mild mosaic virus (CaMMV) and Cacao yellow vein-banding virus (CYVBV) [48], indicated that certain CPDs are conserved in ORFs 1 and 3 by all six species, while several domains were found only on ORF3 of CaMMV or CYVBV, and absent in CSSD-genomes (Table 1, Fig. 2a). In particular, DUF4200, which is unique to CaMMV ORF3, was located between the Zn and Pepsin domains. In CYVBV ORF3, the plant homeodomain finger (PHD) was identified adjacent to the Zn domain, and the trimeric-dUTPase (Trim) domain was located between the zinc finger and the aspartate protease domains (Table 1, Fig. 2a).

Pairwise nucleotide identity for species demarcation

The SDT analysis of the RT-RNase H region identified four species based on ≥80% badnavirus species cutoff established by the International Committee on the Taxonomy of Viruses (ICTV): (i) CRVV, (ii) CSSV, (iii) Cacao swollen shoot CD (Cote d’Ivoire) virus (CSSCDV), and (iv) Cacao swollen shoot Togo A virus (CSSTAV) (Table 2). The best-represented species is CSSV, at 10 isolates, including those reported here and GenBank, from Ghana (AJ608931, AJ609019, and AJ609020) and Togo (AJ534983 and L14546), making it the most widely distributed species thus far. CSSCDV and CSSTAV are reported so far from Cote d’Ivoire (JN606110) and Togo (AJ781003), respectively. And, the herein proposed species, CRVV, represented by isolates from Cote d’Ivoire (CI275, CI286) and Ghana (GH64 and GH67), had the greatest extent of intraspecific divergence among the known CSSD-badnaviruses, at 28–30% (70–72% identity) (Table 3).

Table 2 Pairwise nucleotide identity of the RT-RNase H sequence for cacao-infecting badnavirus isolates
Table 3 Pairwise nucleotide identity for the complete genome sequence of cacao-infecting badnavirus isolates determined in this study

By pairwise distance analysis of the RT-RNase H region, the 14 newly determined and seven GenBank CSSD-associated genome sequences shared 72–100% nt identity (Table 2). In contrast, the corresponding RT-RNase H region for the 21 isolates and 80 other badnaviral RT-RNase H sequences (GenBank) was 58–70% identical. Notably, CaMMV and CYVBV are as divergent from each other (40%), at 60% shared nt identity, as they are from West African CSSD-badnaviruses. Although SDT distinguished four species (≥80%) irrespective of RT-RNase H or complete viral genome sequences, species group composition was inconsistent, in particular, for the CSSD isolates. For example, SDT analysis for complete CRVV genomes (≥80% threshold) delimited CI286, and the group containing CI275, GH64, and GH67, as two groups. In contrast, the RT-RNase H analysis identified the same four isolates as one species (Tables 2 and 3). Also, the RT-RNase H analysis placed TG_AJ781003 with CSSTAV species, while by complete genome analysis, the isolate grouped with CSSV species (Tables 2 and 3). Further, the CSSTAV complete genome shared 80% nt identity with Togo isolates L14546 and AJ534983, compared to 75% for the corresponding RT-RNase H region, resulting in a conflicted species status. Finally, only CSSCDV (CI_JN606110) was delimited as the same species by both the RT-RNase H and the full-length genome sequence analysis.

Phylogenetic analyses

The phylogenetic relationships were determined, based on the RT-RNase H region of 87 GenBank accessions and the 14 CSSD-badnaviral sequences. Despite the accepted use of this region for badnavirus species demarcation, the ML tree was unresolved, e.g. as a polytomy, at ≥70% bootstrap support. Within the badnavirus polytomy, the West African CSSD isolates grouped separately from the other badnaviruses, with 99% bootstrap support, indicating a close phylogenetic relationship among cacao-infecting viruses in West Africa (Fig. 3a). In contrast, ML analysis of the complete genome sequences for the same isolates (CSSD and 87 badnaviral genome sequences from GenBank), resolved three clades, each well supported by ≥70% bootstrap values. For discussion purposes here, the clades have been arbitrarily labeled as I, II and III (Fig. 3b). Clade I (70% bootstrap) contained all known (21) CSSD-associated badnaviral genomes (Fig. 3b). Clade I also contained badnaviruses associated with diverse host species, including one recently reported cacao-infecting virus from Trinidad, CaMMV [48]. A second cacao-infecting badnavirus from Trinidad, CYVBV, grouped with seven non-cacao-infecting badnaviruses in Clade II (Fig. 3b). Clade III contained 18 badnaviruses from non-cacao hosts, including the banana streak virus species group (Fig. 3b).

Fig. 3
figure 3

Phylogenetic analysis using Maximum likelihood (ML), implemented in MEGA 6 [45], for 87 badnaviral genomes (a) RT-RNase H region (1230 bp; ~coordinates 5196–6425) and (b) complete genome sequences. Horizontal branch lengths are proportional to genetic distance, and the values placed at major nodes indicate ≥70% bootstrap support (1000 iterations). Each genome sequence is indicated by the code for the country from which the sample was collected, and sample number. CSSD-badnavirus reference sequences are indicated by Genbank Accession number, and country code. The blue, black and green letters for the CSSD subclade isolates indicate number of coding regions per genome as: four, five and six open reading frames, respectively. Country codes and names are: CI = Cote d’Ivoire, GH = Ghana, and TG = Togo. The names of viruses are abbreviated as: PYMaV-Pagoda yellow mosaic associated virus; TBV-Taro bacilliform virus; SPBVA-Sweetpotato badnavirus A; SPBVB-Sweetpotato badnavirus B; MBV1-Mulberry badnavirus 1; DMV-Dracaena mottle virus; CoYMV-Commelina yellow mottle virus; BsCVBV-Bougainvillea spectabilis chlorotic vein-banding virus; SCBIMV-Sugarcane bacilliform IM virus; SCBMorV-Sugarcane bacilliform Mor virus; BSUMV-Banana streak UM virus; BSUIV-Banana streak UI virus; BSULV-Banana streak UL virus; CSSV-Cacao swollen shoot virus; CSSTAV-Cacao swollen shoot Togo A virus; CSSCDV-Cacao swollen shoot CD virus; CRVV-Cacao red vein virus; PYMoV-Piper yellow mottle virus; GVBV-Gooseberry vein banding virus; RYNV-Rubus yellow net virus; HBV-Hibiscus bacilliform virus; FBV1-Fig badnavirus 1; CiYMV-Citrus yellow mosaic virus; PVBV-Pelargonium vein banding virus; BVF-Blackberry virus F; PBCoV-Pineapple bacilliform comosus virus; GVCV-Grapevine vein clearing virus; BSGFV-Banana streak GF virus; BSMYV-Banana streak MY virus; DBV-Dioscorea bacilliform virus; SCBV-Sugarcane bacilliform virus; KTSV-Kalanchoe top spotting virus; BSUAV-Banana streak UA virus; BSOLV-Banana streak OL virus; BSCAV-Banana streak CA virus; BSIMV-Banana streak IM virus; BSVAcV-Banana streak Vietnam Acuminata virus, CaMMV-Cacao mild mosaic virus; and CYVBV-Cacao yellow vein-banding virus

Discussion

In this study, 14 full-length badnavirus genome sequences were determined by Illumina ‘virome discovery’, and confirmed by Sanger DNA sequencing of cloned, PCR amplicons. In all instances, the badnaviral genome sequences were obtained from leaf samples collected from cacao trees exhibiting foliar discoloration and swollen shoot symptoms, whereas, several of the trees also showed atypical symptoms consisting of accelerated tree decline, and rapid death. Although speculative, the latter symptoms are reminiscent of rapid necrosis and death, often associated with hypersensitive-like responses to pathogen infection.

Based on genome size and arrangement, and nt and aa sequence comparisons the genomic sequences determined here were most closely related to other previously reported badnaviruses, with the closest relatives being the three known cacao-infecting badnaviruses from West Africa: CSSV, CSSCDV, and CSSTAV [10, 25] (Additional file 2: Table S2). In addition, all of the CSSD-genomes contained the badnavirus hallmark tRNAmet binding site, which is reminiscent of tRNA sequences present in plant genomes [49], leading to speculation that badnaviral-associated tRNA sequences are host-derived.

The CSSD-genomes were variable with respect to predicted coding regions and genome arrangement, but were consistent among other badnaviruses, which are known to have three, four, or five predicted ORFs, variously arranged (Fig. 1). This observation has been reported previously for other CSSD isolates, which have been recently recognized as three distinct species [10]. The proposed new badnavirus species of cacao, CRVV, identified for the first time in this study, is highly divergent from CSSV, CSSDV, and CSSTSV, with which it shares only 70–72% nt identity. And, the previously known species encode five or six ORFs, compared to four ORFs predicted for the CRVV genome, indicating that CRVV is unique among the West African CSSD-badnaviruses identified so far (Additional file 2: Table S2). Compared to the four CSSD species, the New World cacao-infecting badnaviruses, CaMMV and CYVBV, described thus far encode four ORFs [48], whereas, the genomes of all known non-cacao infecting badnaviruses (genus-wide) have from three to seven ORFs.

The predicted conserved domains among the 14 genomes showed similarities to those reported in the genus Badnavirus, including the DUF (ORF1), and Zn, Pep, RT, and RNase H, all found in ORF3 (Table 1, Fig. 3). The only exception was the CSSCDV genome [28], which lacked a detectable Zn domain, and had an additional, unique DUF3187 domain in ORFY that is annotated as an ‘outer membrane hypothetical protein’ in certain Proteobacteria [42]. The absence of the Zn domain is difficult to reconcile because it has been reported to be an essential coat protein motif [25]. Interestingly, no domains have been detected in all of the CSSD-associated genomes in ORFs 4, Y, and X, including the 14 reported here (Table 1, Fig. 2). The inability to detect predicted conserved domain(s) may suggest that ‘badnavirus-novel’ domains are not necessarily absent, but that they may not be discoverable using available CDD tools [42].

A comparison between CSSD-associated genomes and the only two other cacao-infecting badnaviruses, CaMMV and CYVBV, showed interesting differences and some similarities in the CPDs. The DUF4200, which is unique to CaMMV ORF3, is annotated as a coiled-coil domain of unknown function (eukaryotic) [42]. In pararetroviruses, coiled-coil domains have been identified in virion-associated proteins, implicated in aphid transmission and for Cauliflower mosaic virus cell-to-cell movement [50, 51]. The PHD, unique to CYVBV ORF3, has been associated with transcriptional regulation and chromatin-associated functions [42], and the Trim also found in CYVBV ORF3, has predicted dUTPase activity e.g. catalyzes hydrolysis of dUTP-Mg complexes to dUMP and pyrophosphate [42, 52]. Although the Trim domain is unstudied for CYVBV, its presence in CYVBV and in two divergent badnaviruses, Piper yellow mottle virus [53] and Dioscorea bacilliform virus [54], suggestive of a possibly conserved, genus-wide function. The presence of the three unique domains in the two Trinidad viral genomes, and their absence in the four CSSD-associated species, make them important targets for future study.

High genomic variability was discovered among the 21 genomes from cacao, and four distinct badnavirus species were identified based on the ICTV-established ≥80% nt identity threshold on the RT-RNase H region (Table 2). Here, the previously unreported CRVV is proposed to constitute a new badnavirus species, and it is partially characterized for the first time. Until recently, all CSSD-badnaviruses were referred to as CSSV. In 2015, CSSV became recognized by the ICTV as the type species, and CSSV, CSSCDV, and CSSTAV were formally designated as distinct badnaviral species [19]. Based on these results, CRVV is the fourth cacao-infecting badnaviral species known to be endemic to West Africa. Previously, a single causal badnavirus was associated with cacao plants exhibiting different swollen shoot disease symptoms. However, the new evidence presented here shows that numerous badnaviral species and strains, evident from extensive within and between clade genomic variability, are associated with disease symptoms.

The concept of species group has not been frequently applied to plant viruses. However, a closely related ‘group of species’, banana streak virus complex, is recognized within the badnavirus genus that is reminiscent of CSSD-badnaviruses. Both comprise a group of closely related badnaviral species that are more closely related to one another than to other known badnavirus, with apparently restricted host-associations [55,56,57]. Similarly, the Sugarcane mosaic virus (Potyviridae) group contains closely-related, species and strains that are divergent from other known potyviruses [58]. And, within the genus, Begomovirus (Geminiviridae), five or more closely related species and multiple strains, identified thus far, incite leaf curl disease of malvaceous hosts, including the cotton crop, widespread on the Indian Subcontinent [59, 60].

Phylogenetic comparisons of all available badnaviral genome sequences, including the genome sequences determined in this study, indicated that CSSD-badnaviral genomes are more closely related to each other than to other known badnaviruses, based on the tight affiliation as a ‘sister clade’, in relation to other badnaviruses. This suggests that all CSSD-associated badnavirus species known thus far to cause swollen shoot disease share a common ancestor(s), providing robust support for West African endemism, and for the region as the center of CSSD- badnaviral diversification.

The first outbreak of CSSD was reported in Ghana [3], followed by Cote d’Ivoire ten years later. Among the available CSSV genome sequences, eleven originate from Cote d’Ivoire and three are from Ghana (Fig. 3b). The presence of CSSV in these neighboring countries is possibly suggestive of phylogeographical distribution e.g. regional CSSV endemism, however, additional sequences are required to clarify the centers of endemism of CSSV as well as the other CSSD species. Irrespective of RT-RNase H region or complete genome comparisons, CSSTAV and CSSCDV, are highly divergent from each other by 24%, and from CSSV and CRVV, by 20–25% and 28–30%, respectively (Table 3). This degree of genomic variability is suggestive of long-standing separation, and reminiscent of phylogeographical and/or host associated co-evolution. The CSSD-badnaviruses are considered endemic viruses to West Africa, and have been documented to infect many endemic tropical tree and shrub genera, including T. cacao relatives at the family-level, making it likely that additional undiscovered badnaviruses infect wild hosts in West Africa, with potential to undergo host shifting under opportune conditions.

Although genomic sequences are available for a relatively small number of isolates overall, the hypothesis that CRVV may have emerged in cacao as the result of a host jump from its wild reservoir(s), may be supportable based on its low shared nt identity, at ~70–72%, with other CSSD-badnavirus species extant in cacao, and could also explain why it is most closely phylogenetically-related to them. Evidence for CSSD-badnavirus infection of endemic tropical tree species is based nearly entirely on biological studies, e.g. grafting and/or mealybug transmission tests, whose results have not been verified using molecular methods. Based on these studies, a large number of CSSD hosts have been reported, including Adansonia digitata L., Ceiba pentandra L., Cola chlamydantha K. Schum., Cola gigantean A. Chev., and Sterculia tragacantha Lindl [61,62,63], and are among ~90 species in 30 plant families used as shade for cacao and other crops [64]. Molecular confirmation of suspect CSSD host-infection and accurate badnavirus identification when found, are important first steps to enabling knowledge of CSSD-badnavirus evolution and origin to be reconciled with specific epidemiological factors leading to outbreaks, to inform short-term management approaches and CSSD breeding strategies to enable sustainable production of the crop in the long term.

Historically, CSSD isolates have been characterized as ‘mild’ or ‘severe’, with mild isolates causing mild foliar symptoms that persisted for only several days, such as CSSV-N1A (AJ609020), and the characteristically prevalent, severe isolates, such as CSSV-New Juaben (AJ608931), that cause persistent foliar and shoot symptoms, decline, and tree death tree in 3–5 years [2, 65]. Although the N1A and New Juaben ‘strains’ are members of the CSSV species (Table 2), and same subclade (Fig. 3a, b), reported differences in pathogenicity are confounding. There is no definitive link between the rapid decline phenotype strains of recent, and strains previously recognized, as ‘severe’. The genetic basis for differences in pathogenicity e.g. virulence, among the CSSD-badnaviruses has not been investigated, however important clues may reside in the diverse genome arrangements and CPD architectures (Fig. 2b), particularly those associated with predicted functions in pathogenicity.

Despite nearly 100 years of CSSD research, there is no definitive understanding of the connections between origins or pathways of disease spread, or of the latter relationships with genomic diversification of cacao-infecting badnaviruses. The ‘unique to West Africa-genome type’ of the proposed CRVV species (Table 2; Fig. 3b), and its only recent association with the atypical, rapid decline and death phenotype observed in cacao trees in western Ghana during 2000, and next in eastern and then western Cote d’Ivoire by 2003, points to a possible CRVV-origin near the border of the two countries, with subsequent spread into plantations recently established in Western Cote d’Ivoire. Because CRVV-like MP sequences were obtained from symptomatic samples collected during surveys carried out in Cote d’Ivoire and Ghana in 2000–2003, [28] albeit, unknowingly, until this report, CRVV and perhaps other unknown variants may possibly have already spread to other cacao-producing countries, including Nigeria and Togo. These observations underscore the urgent need to identify the causal agent of the ‘rapid decline’ phenotype, and circumvent its further spread in cacao. This goal can only be achieved through the coordinated development and use of reliable molecular diagnostic tools, and well-supported surveillance efforts to track this dynamic badnavirus complex, and carry out epidemiological studies on a regional level in cacao plantations and in nearby suspect, endemic plant hosts of the CSSD-complex.

Conclusions

This is the first discovery of CRVV, a previously undescribed badnavirus species infecting cacao trees in West Africa, represented thus far, by two isolates each from Ghana and Cote d’Ivoire. The CRVV-isolates were collected from locations where the first swollen shoot ‘rapid decline and death’ outbreaks were observed during 2000–2003. This is suggestive of a potential link between CRVV distribution, and the locations where the onset of rapid tree decline and death were observed within one year after foliar and swollen shoot symptom development. Although CRVV is the suspect causal agent of the rapid decline and death syndrome, Koch’s Postulates have not been proven, thus, urgent studies are now required to test this hypothesis. In addition, the phylogenetic analysis of all available cacao-associated genome sequences indicated that certain CSSV isolates in Ghana, which have been historically referred to as ‘severe strains’ of CSSV because of the debilitating symptoms they cause, are distinct from CRVV, albeit, it is important to note that CRVV and CSSV occupy an overlapping geographical range. This knowledge, together with additional, intensive CSSD surveillance using Next Generation DNA sequencing to facilitate genomic pathology studies of West African cacao plantations and suspect endemic host reservoirs, is expected to provide a data-rich resource with far-reaching applications, including the development and use of field and laboratory diagnostics to gain new insights into epidemiological and viral-host evolutionary dynamics, and provide vital assistance to breeding programs developing durable badnavirus-resistant cacao for long-term management.

Abbreviations

CaMMV:

Cacao mild mosaic virus

CDD:

Conserved domain database

CP:

Coat protein

CPD:

Conserved protein domain

CRVV:

Cacao red vein virus

CSSCDV:

Cacao swollen shoot CD (Cote d’Ivoire) virus

CSSD:

Cacao swollen shoot disease

CSSTAV:

Cacao swollen shoot Togo A virus

CSSV:

Cacao swollen shoot virus

CYVBV:

Cacao yellow vein-banding virus

DUF:

Domain of unknown function

ICTV:

International committee on taxonomy of viruses

ML:

Maximum Likelihood

MP:

Movement protein

ORF:

Open reading frame

PHD:

Plant homeodomain

RCA:

Rolling circle amplification

RNase H:

Ribonuclease H

RT:

Reverse transcriptase

SDT:

Sequence demarcation tool

References

  1. Bartley BGD. The genetic diversity of cacao and its utilization. Wallingford, Oxfordshire, UK: CABI Publishing; 2005.

    Book  Google Scholar 

  2. Ollennu L, Owusu GK. Field evaluation of the protective capacity of CSSV mild strain N1 against severe strain new Juaben (1A) isolate. Ghana J Agric Sci. 2003;36:3–12.

    Google Scholar 

  3. Steven WF. Swollen shoot and die-back—a new disease of cocoa. Gold Coast Farmer. 1936;V:144.

  4. Posnette AF. Virus diseases of cacao in West Africa. V. Alternative host plants. Ann. Appl. Biol. 1950;34:388–402.

    Google Scholar 

  5. Adegbola MOK. The reaction of four Nigerian isolates of CSSV on the germination and growth of cacao seedlings under different environmental conditions. Proc. 5th Int. Cocoa Res. Conf. Ibadan, Nigeria; 1975. p. 338–343.

  6. Thresh JM. The spread of virus disease in cacao. West African Cocoa Res Inst Tech Bull. 1958;5:36.

    Google Scholar 

  7. Dzahini-Obiatey H, Fox RTV. Early signs of infection in cacao swollen shoot virus (CSSV) inoculated cocoa seeds and the discovery of the cotyledons of the resultant plants as rich sources of CSSV. African J Biotechnol. 2010;9:593–603.

    Article  CAS  Google Scholar 

  8. Ameyaw GA, Wetten A, Dzahini-Obiatey H, Domfeh O, Allainguillaume J. Investigation on Cacao swollen shoot virus (CSSV) pollen transmission through cross-pollination. Plant Pathol. [Internet]. 2013 [cited 2014 Sep 6];62:421–7. Available from: http://0-doi-wiley-com.brum.beds.ac.uk/10.1111/j.1365-3059.2012.02640.x

  9. Mangenot G, Alibert H, Basset A. Sur les caracteres du swollen shoot en Cote d’Ivoire. Rev Int Bot Appl Paris. 1946;283:13.

    Google Scholar 

  10. Muller E, Sackey S. Molecular variability analysis of five new complete cacao swollen shoot virus genomic sequences. Arch Virol. 2005;150:53–66.

    Article  CAS  PubMed  Google Scholar 

  11. Domfeh O, Dzahini-Obiatey H, Ameyaw GA, Abaka-Ewusie K, Opoku G. Cocoa swollen shoot virus disease situation in Ghana : a review of current trends. African J Agric Res. 2011;6:5033–9.

    Google Scholar 

  12. Dzahini-Obiatey H, Ameyaw GA, Ollennu LA. Control of cocoa swollen shoot disease by eradicating infected trees in Ghana: A survey of treated and replanted areas. Crop Prot. [Internet]. 2006 [cited 2014 Sep 12];25:647–52. Available from: http://0-linkinghub-elsevier-com.brum.beds.ac.uk/retrieve/pii/S0261219405002504

  13. Dzahini-Obiatey H, Domfeh O, Amoah FM. Over seventy years of a viral diseaseof cocoa in Ghana: from researcher’s perspective. African J. Agric. Res. 2010;5:476–85.

    Google Scholar 

  14. Thresh JM, Tinsley TW. The viruses of cocoa. West African cocoa res. Inst. Tech. Bull. 1959;7:1–32.

    Google Scholar 

  15. Dongo LN, Orisajo SB. Status of cocoa swollen shoot virus disease in Nigeria. African J. Biotechnol. 2007;6(17):2054–61.

    Article  Google Scholar 

  16. Attafuah A, Blencowe JW, Brunt AA. Swollen shoot disease of cocoa in the Sierra Leone. Trop Agric Trinidad. 1963;40:229–32.

    Google Scholar 

  17. Partiot M, Amefia YK, Djiekpor EK, Bakar KA. Le swollen shoot du cacaoyer au Togo. Investaire preliminaire et premiere estimation des pertes causees par la maladie. Cafe-cacao Thea. 1978;22:217–28.

    Google Scholar 

  18. Oro F, Mississo E, Okassa M, Guilhaumon C, Fenouillet C, Cilas C, et al. Geographical differentiation of the molecular diversity of cacao swollen shoot virus in Togo. Arch. Virol. [Internet]. 2012 [cited 2014 Sep 6];157:509–14. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/22160623.

  19. Muller E. Creation of four new species in the genus Badnavirus [internet]. 2015. Available from: https://talk.ictvonline.org/ICTV/proposals/2015.026aP.A.v2.Badnavirus_4sp.pdf​.

  20. King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. Badnavirus in: virus taxonomy. Ninth report of the international committee on taxonomy of viruses. London: Elsevier Academic Press; 2012.

    Google Scholar 

  21. Lot H, Djiekpor E, Jacquemond M. Characterization of the genome of cacao swollen shoot virus. J Gen Virol. 1991;72:1735–9.

    Article  CAS  PubMed  Google Scholar 

  22. Bhat AI, Hohn T, Selvarajan R. Badnaviruses : the current global scenario. Viruses. 2016;8:e177.

    Article  PubMed  Google Scholar 

  23. Borah BK, Sharma S, Kant R, Johnson AMA, Saigopal DVR, Dasgupta I. Bacilliform DNA-containing plant viruses in the tropics: commonalities within a genetically diverse group. Mol Plant Pathol. 2013;14:759–71.

    Article  CAS  PubMed  Google Scholar 

  24. Geering ADW. Caulimoviridae (Plant pararetroviruses). eLS [Internet]. 2014; Available from: http://0-doi-wiley-com.brum.beds.ac.uk/10.1002/9780470015902.a0000746.pub3

  25. Hagen LS, Jacquemond M, Lepingle A, Lot H, Tepfer M. Nucleotide sequence and genomic organization of cacao swollen shoot virus. Virology. 1993;196:619–28.

    Article  CAS  PubMed  Google Scholar 

  26. Jacquot E, Hagen LS, Jacquemond M, Yot P. The open reading frame 2 product of cacao swollen shoot badnavirus is a nucleic acid-binding protein. Virology. 1996;225:191–5.

    Article  CAS  PubMed  Google Scholar 

  27. Jacquot E, Hagen LS, Michler P, Rohfritsch O, Stussi-Garaud C, Keller M, et al. In situ localization of cacao swollen shoot virus in agroinfected Theobroma Cacao. Arch. Virol. [internet]. 1999;144:259–71. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/10470252.

  28. Kouakou K, Kébé BI, Kouassi N, Aké S, Cilas C, Muller E. Geographical distribution of cacao swollen shoot virus molecular variability in Côte d ’ Ivoire. Plant Dis. 2012;96:1445–50.

    Article  Google Scholar 

  29. Sagemann W, Lesemann D-E, Paul HL, Adomako D, Owusu GK. Detection and some comparison of Ghanaian isolates of cacao swollen shoot virus (CSSV) by enzyme-linked immunosorbent assay (ELISA) and immuno-electron microscopy (IEM) using an antiserum to CSSV strain 1A. J Phytopathol. 1985;114:79–89.

    Article  Google Scholar 

  30. Abrokwah F, Dzahini-Obiatey H, Galyuon I, Osae-Awuku F, Muller E. Geographical distribution of cacao swollen shoot virus molecular variability in Ghana. Plant Dis. 2016;100:2011–7.

    Article  Google Scholar 

  31. Chingandu N, Kouakou K, Aka R, Gutierrez OA, Brown JK. Unexpected genome variability at multiple loci suggests Cacao swollen shoot virus comprises multiple, divergent molecular variants. J. Emerg. Dis. Virol. [Internet]. 2017;3:1–10. Available from: https://www.sciforschenonline.org/journals/virology/article-data/JEDV-3-128/JEDV-3-128.pdf

  32. Doyle JJ, Doyle JL. Isolation of plant DNA from fresh tissue. Focus (Madison). 1990;12:13–5.

    Google Scholar 

  33. Stevens H, Rector A, Van Ranst M. Multiply primed rolling-circle amplification method for the amplification of circular DNA viruses. Cold spring Harb. Protoc. [internet]. 2010 [cited 2014 Sep 19];2010:pdb.prot5415. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/20360369.

  34. Rector A, Tachezy R, Van Ranst M. A sequence-independent strategy for detection and cloning of circular DNA virus genomes by using multiply primed rolling-circle amplification. J. Virol. [Internet]. 2004;78:4993–8. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=400362&tool=pmcentrez&rendertype=abstract

  35. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Argout X, Salse J, Aury J, Guiltinan MJ, Droc G, et al. The genome of Theobroma Cacao. Nat Genet. 2011;43:101–8.

    Article  CAS  PubMed  Google Scholar 

  37. Motamayor JC, Mockaitis K, Schmutz J, Haiminen N, Iii DL, Cornejo O, et al. The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color. Genome Biol. [Internet]. BioMed Central Ltd; 2013;14:r53. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=4053823&tool=pmcentrez&rendertype=abstract

  38. Conesa A, Gotz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:1–12.

    Article  Google Scholar 

  39. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  CAS  PubMed  Google Scholar 

  40. Hull R, Covey SN. Does cauliflower mosaic virus replicate by reverse transcription? Trends Biochem Sci. 1983;8:119–21.

    Article  CAS  Google Scholar 

  41. Dixon LK, Hohn T. Initiation of translation of the cauliflower mosaic virus genome from a polycistronic mRNA: evidence from deletion mutagenesis. EMBO J. 1984;3:2731–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, et al. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2015;43:D222–6.

    Article  CAS  PubMed  Google Scholar 

  43. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Villesen P. FaBox: an online toolbox for FASTA sequences. Mol Ecol Notes. 2007;7:965–8.

    Article  CAS  Google Scholar 

  45. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evolution. 2013;30:2725–9.

    CAS  Google Scholar 

  46. Medberry SL, Lockhart BEL, Olszewski NL. 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–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Hohn T, Fütterer J, Hull R. The proteins and functions of plant pararetroviruses : Knowns and unknowns. CRC Crit Rev Plant Sci. 1997;16:133–61.

    Article  CAS  Google Scholar 

  48. Chingandu N, Zia-ur-rehman M, Sreenivasan TN, Surujdeo-Maharaj S, Umaharan P, Gutierrez OA, et al. Molecular characterization of previously elusive badnaviruses associated with symptomatic cacao in the New World. Arch. Virol. [Internet]. 2017;1–9. Available from: https://0-doi-org.brum.beds.ac.uk/10.1007/s00705-017-3235-2

  49. Ghosh P, Ghosh K, Simsekt M, Rajbhandaryt UL. Nucleotide sequence of wheat germ cytoplasmic initiator methionine transfer ribonucleic acid. Nucleic Acids Res. 1982;10:3241–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Stavolone L, Herzog E, Leclerc D, Hohn T. Tetramerization is a conserved feature of the virion- associated protein in plant pararetroviruses. J Virol. 2001;75:7739–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Stavolone L, Villani ME, Leclerc D, Hohn T. A coiled-coil interaction mediates cauliflower mosaic virus cell-to-cell movement. Proc. Natl. Acad. Sci. U. S. A. [Internet]. 2005;102:6219–24. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=15837934

  52. Tormo-Mas MA, Donderis J, Garcia-Caballer M, Alt A, Mir-Sanchis I, Marina A, et al. Phage dUTPases control transfer of virulence genes by a proto-oncogenic G protein-like mechanism. Mol Cell. 2013;49:947–58.

    Article  CAS  PubMed  Google Scholar 

  53. Hany U, Adams IP, Glover R, Bhat AI, Boonham N. The complete genome sequence of piper yellow mottle virus (PYMoV). Arch Virol. 2014;159:385–8.

    Article  CAS  PubMed  Google Scholar 

  54. Briddon RW, 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–83.

    Article  CAS  PubMed  Google Scholar 

  55. 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–20.

    Article  CAS  PubMed  Google Scholar 

  56. Geering AD, McMichael LA, Dietzgen RG, Thomas JE. Genetic diversity among banana streak virus isolates from Australia. Phytopathology. 2000;90:921–7.

    Article  CAS  PubMed  Google Scholar 

  57. Lockhart BEL, Olszewski NE. Serological and heterogeneity of Banana streak badnavirus: Implications for virus detection in Musa germplasm. In: Ganry J, editor. Breed. Banan. plantain Resist. to Dis. pests. Montepellier, France: CIRAD, INIBAP; 1993. p. 105–13.

  58. Yang ZN, Mirkov TE. Sequence and relationships of sugarcane mosaic and sorghum mosaic virus strains and development of RT-PCR-based RFLPs for strain discrimination. Phytopathology. 1997;87:932–9.

    Article  CAS  PubMed  Google Scholar 

  59. Brown JK, Zerbini FM, Navas-Castillo J, Moriones E, Ramos-Sobrinho R, Silva JCF, et al. Revision of Begomovirus taxonomy based on pairwise sequence comparisons. Arch Virol. 2015;160:1593–619.

    Article  CAS  PubMed  Google Scholar 

  60. Briddon RW. Cotton leaf curl disease, a multicomponent begomovirus complex. Mol Plant Pathol. 2003;4:427–34.

    Article  CAS  PubMed  Google Scholar 

  61. Posnette AF, Robertson NF, Todd JM. Virus diseases of cacao in West Africa. V. Alternative host plants. Ann. Appl. Biol. 1950;37:229–40.

    Google Scholar 

  62. Tinsley TW. The ecology of cacao viruses. I. The role of wild hosts in the incidence of swollen shoot virus in West Africa. J Appl Ecol. 1971;8:491–5.

    Article  Google Scholar 

  63. Todd JM. An indigenous source of swollen shoot disease of cacao. Nature. 1951;167:952–3.

    Article  CAS  PubMed  Google Scholar 

  64. Richard A, Ræbild A. Tree diversity and canopy cover in cocoa systems in Ghana. New For Springer Netherlands. 2016;47:287–302.

    Google Scholar 

  65. Ameyaw GA, Domfeh O, Dzahini-Obiatey HK, Ollennu LAA, Owusu GK. Appraisal of cocoa swollen shoot virus (CSSV) mild isolates for cross protection of cocoa against severe strains in Ghana. Plant Dis. 2016;100:810–5.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We gratefully acknowledge funding for this project from the USDA-ARS through the Specific Cooperative Agreement #6038-21000-023-07 titled: Development and Optimization of Molecular Diagnostics Method for Qualitative and Quantitative Detection of Cacao Swollen Shoot Virus with MARS Inc. through Trust Agreement #58-6631-6-123, titled: Genetic Improvement of Theobroma cacao, and from the Borlaug Fellowship Program, funded by the USDA-Foreign Agricultural Services with support from the World Cocoa Foundation.

Funding

Funding was provided by USDA-ARS through a Specific Cooperative Agreement, Project # 6038–21,000–023-07, supported in part by MARS, Inc. Trust Agreement #58–6631–6-123, and by the World Cocoa Foundation Foundation Borlaug Fellowship Program for KK and RA.

Availability of data and materials

The complete genome sequences generated during the current study are available in the NCBI GenBank repository, and they were assigned the Accession numbers KX592571- KX592584.

Author information

Authors and Affiliations

Authors

Contributions

NC and JKB designed the experiments. NC performed the experiments, including primer design, PCR and cloning, sequence assemblies, and data analysis. NC and JKB interpreted the results. NC wrote the first draft, and NC and JKB edited the manuscript. KK, RA, and GA provided plant samples, and provided constructive comments and edits for the manuscript. Plant sample collections were made under authorization of the latter (KK, RA, and GA) institutions in Cote d’Ivoire (CNRA) and Ghana (CRIG), which conduct research and extension in cacao and so have permission to visit farms and public lands to collect the cacao samples used for virus analysis. The plant samples were received under USDA-APHIS permit number P526P-15-02089 (exp. 04.21.2018) issued to JKB for diagnostic purposes. HWH cooperated in the initial assembly and editing of Illumina sequence data, and provided constructive comments on the submitted manuscript. OAG contributed to the design of the study and manuscript editing. JKB developed the research plan, provided oversight and coordination for the project, and was the recipient of the research grant. OAG collaborated on the work plan, and provided guidance in cacao genetics, as well as financial support through his program at USDA-ARS. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Nomatter Chingandu or Judith K. Brown.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

The USDA-ARS review has been completed. No other consents are required.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1: Table S1.

Primers for polymerase chain reaction amplification of representative full-length badnaviral genomes associated with cacao swollen shoot disease in West Africa. The viral-specific sequences for each primer are underlined, and the first 8 or 9 bases indicate the plasmid vector sequence to which a Not I restriction site (indicated in italics) was included to facilitate cloning. The abbreviations indicating the country of sample collection, as GH = Ghana and CI = Cote d’Ivoire. (DOCX 13 kb)

Additional file 2: Table S2.

The open reading frames (ORFs) predicted on the plus-strand of cacao-associated badnaviral genomes determined here, using the BLAST and ORF Finder tool algorithms (NCBI-GenBank database). BLASTn and BLASTp scores were considered for the top twenty Genbank hits. The abbreviations indicate the country of sample collection, as CI = Cote d’Ivoire, GH = Ghana, and TG = Togo. The notation of ‘ND’ denotes where no ORF was predicted. (DOCX 17 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chingandu, N., Kouakou, K., Aka, R. et al. The proposed new species, cacao red vein virus, and three previously recognized badnavirus species are associated with cacao swollen shoot disease. Virol J 14, 199 (2017). https://0-doi-org.brum.beds.ac.uk/10.1186/s12985-017-0866-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12985-017-0866-6

Keywords