# Novel Ampeloviruses Infecting Cassava in Central Africa and the South-West Indian Ocean Islands

Yves Kwibuka, Espoir Bisimwa, Arnaud G Blouin, Claude Bragard, Thierry T. Candresse, Chantal Faure, Denis Filloux, Jean-Michel Lett, François Maclot, Armelle Marais, et al.

► To cite this version:

Yves Kwibuka, Espoir Bisimwa, Arnaud G Blouin, Claude Bragard, Thierry T. Candresse, et al. Novel Ampeloviruses Infecting Cassava in Central Africa and the South-West Indian Ocean Islands. Viruses, 2021, 13 (6), pp.1030. 10.3390/v13061030. hal-03244473

HAL Id: hal-03244473

https://hal.inrae.fr/hal-03244473

Submitted on 1 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Novel Ampeloviruses Infecting Cassava in Central Africa and the South-West Indian Ocean Islands

Yves Kwibuka<fnref n="1" />,<fnref n="2" />,* ID, Espoir Bisimwa<fnref n="2" />, Arnaud G. Blouin<fnref n="1" />, Claude Bragard<fnref n="3" /> ID, Thierry Candresse <fnref n="4" /> ID,
Chantal Faure<fnref n="4" />, Denis Filloux<fnref n="5" />,<fnref n="6" />, Jean-Michel Lett <fnref n="7" /> ID, François Maclot<fnref n="1" />, Armelle Marais<fnref n="4" />,
Santatra Ravelomanantsoa<fnref n="8" /> ID, Sara Shakir <fnref n="9" /> ID, Hervé Vanderschuren <fnref n="9" />,<fnref n="10" /> and Sébastien Massart <fnref n="1" />,*

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Kwibuka, Y.; Bisimwa, E.; Blouin, A.G.; Bragard, C.; Candresse, T.; Faure, C.; Filloux, D.; Lett, J.-M.; Maclot, F.; Marais, A.; et al. Novel Ampeloviruses Infecting Cassava in Central Africa and the South-West Indian Ocean Islands. *Viruses* 2021, 13, 1030. https://doi.org/10.3390/v13061030

Academic Editors: Jan Frederik Kreuze and Adrian Fox

Received: 28 April 2021

Accepted: 25 May 2021

Published: 29 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Cassava is one of the most important staple crops in Africa and its production is seriously damaged by viral diseases. In this study, we identify for the first time and characterize the genome organization of novel ampeloviruses infecting cassava plants in diverse geographical locations using three high-throughput sequencing protocols [Virion-Associated Nucleotide Acid (VANA), dsRNA and total RNA], and we provide a first analysis of the diversity of these agents and of the evolutionary forces acting on them. Thirteen new *Closteroviridae* isolates were characterized in field-grown cassava plants from the Democratic Republic of Congo (DR Congo), Madagascar, Mayotte, and Reunion islands. The analysis of the sequences of the corresponding contigs (ranging between 10,417 and 13,752 nucleotides in length) revealed seven open reading frames. The replication-associated proteins have three expected functional domains: methyltransferase, helicase, and RNA-dependent RNA polymerase (*RdRp*). Additional open reading frames code for a small transmembrane protein, a heat-shock protein 70 homolog (*HSP70h*), a heat shock protein 90 homolog (*HSP90h*), and a major and a minor coat protein (*CP* and *CPd* respectively). Defective genomic variants were also identified in some cassava accessions originating from Madagascar and Reunion. The isolates were found to belong to two species tentatively named Manihot esculenta-associated virus 1 and 2 (MEaV-1 and MEaV-2). Phylogenetic analyses showed that MEaV-1 and MEaV-2 belong to the genus *Ampelovirus*, in particular to its subgroup II. MEaV-1 was found in all of the countries of study, while MEaV-2 was only detected in Madagascar and Mayotte. Recombination analysis provided evidence of intraspecies recombination occurring between the isolates from Madagascar and Mayotte. No clear association with visual symptoms in the cassava host could be identified.

**Keywords:** Manihot esculenta; Ampelovirus; high-throughput sequencing; Central Africa; Indian Ocean islands; *Closteroviridae*

<footnote n="1">¹ Plant Pathology Laboratory, TERRA-Gembloux Agro-Bio Tech, University of Liège, Passage des Déportés, 2, 5030 Gembloux, Belgium; arnaud.blouin@uliege.be (A.G.B.); Francois.maclot@chuliege.be (F.M.)</footnote>

<footnote n="2">² Faculté des Sciences Agronomiques, Université Catholique de Bukavu, BP 285 Bukavu, Democratic Republic of the Congo; ebisimwa@yahoo.com</footnote>

<footnote n="3">³ Earth and Life Institute, Applied Microbiology-Phytopathology, UCLouvain, 1348 Louvain-la-Neuve, Belgium; claude.bragard@uclouvain.be</footnote>

<footnote n="4">⁴ Université Bordeaux, INRAE, UMR BFP, CS20032, CEDEX, 33882 Villenave d'Ornon, France; thierry.candresse@inrae.fr (T.C.); chantal.faure@inrae.fr (C.F.); armelle.marais-colombel@inrae.fr (A.M.)</footnote>

<footnote n="5">⁵ CIRAD, UMR PHIM, 34090 Montpellier, France; denis.filloux@cirad.fr</footnote>

<footnote n="6">⁶ PHIM Plant Health Institute, Université Montpellier, CIRAD, INRAE, Institut Agro, IRD, 34000 Montpellier, France</footnote>

<footnote n="7">⁷ CIRAD, UMR PVBMT, Pôle de Protection des Plantes, Saint-Pierre, F-97410 Ile de la Reunion, France; jean-michel.lett@cirad.fr</footnote>

<footnote n="8">⁸ FOFIFA-CENRADERU, Laboratoire de Pathologie Végétale, BP 1444 Ambatobe, Madagascar; srakotoarisoaravel@gmail.com</footnote>

<footnote n="9">⁹ Plant Genetics Laboratory, TERRA-Gembloux Agro-Bio Tech, University of Liège, Passage des Déportés, 2, 5030 Gembloux, Belgium; shakir.sara@yahoo.com (S.S.); herve.vanderschuren@uliege.be (H.V.)</footnote>

<footnote n="10">¹⁰ Laboratory of Tropical Crop Improvement, Division of Crop Biotechnics, Biosystems Department, KU Leuven, 3000 Leuven, Belgium</footnote>

<footnote>* Correspondence: bisimwayves@gmail.com (Y.K.); sebastien.massart@uliege.be (S.M.)</footnote>
# 1. Introduction

Approximately 40% of the world's area dedicated to the cultivation of roots and tubers is planted with cassava. In Africa, 56% of the total production of roots and tubers comes from cassava, 51% of which are from Nigeria and the Democratic Republic of Congo (DR Congo) [1]. Cassava is one of the most important staple crops, as it ranks fourth as a source of calories for human consumption [1]. Many of the poorest farmers and most undernourished households in Africa depend on cassava as a principal nutrition source. It also constitutes an important income source in rural and marginal areas and has multiple uses, most notably as a food security and regular food crop [2].

Due to its remarkable adaptability to a wide range of soil and environmental conditions, cassava is bestowed with resilience to global warming and climate change, and with a potential for better return under adverse soil and weather conditions [3]. The exceptional soil carbon sequestration properties of cassava makes it a potential crop for improving the green revolution fatigue [4].

Despite the large contribution of Africa to global cassava production (60%), its performance in terms of yield is the lowest (9 tons/ha on average) [1]. In most African cassava-producing areas, the yield is far below the potential [5]. In DR Congo, yields of approximately 8 tons/ha can be obtained under farmers' fields conditions, but yields are highly variable depending on the local pedo-climatic conditions.

Cassava suffers from many pests and diseases, including viruses, which can seriously affect the quality and quantity of the harvest as well as the quality of the planting materials. Two viral diseases are of major economic importance in sub-Saharan Africa, namely, cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) [6]. CBSD is associated with two *Ipomovirus* species, collectively named cassava brown streak viruses (CBSVs), and CMD is associated with nine *Begomovirus* species, collectively named cassava mosaic geminiviruses (CMGs) [6]. Yield reduction due to CMD may be severe and losses up to 82% have been reported, especially in cassava plants dually infected with African cassava mosaic virus (ACMV) and the Ugandan strain of East African cassava mosaic virus (EACMV-UG) [7].

The family *Closteroviridae* is a large and diverse group of filamentous plant viruses (particles of 650–2200 nm in length) with single-stranded RNA genomes. Some species can be transmitted semi-persistently by aphids, whiteflies, mealybugs or soft scales insects [8,9]. The members of this family are known to affect several crops of major economic importance, such as sugar beet, citrus, tomato, lettuce, potato, sweet potato, grapevine, pineapple, cherry, and some ornamentals [9–13].

To date, 56 virus species have been classified as definitive or tentative members of the family *Closteroviridae* [14]. They are grouped into the following four genera: *Ampelovirus* (monopartite genome, mealybugs and soft scale insect vectors), *Closterovirus* (monopartite genome, aphid vectors), *Velarivirus* (monopartite genome, no known vectors) and *Crinivirus* (bipartite genome, whitefly vectors) [14,15].

The accepted demarcation criteria between the species in the *Closteroviridae* family [9,16] include the following aspects: (i) particle size, (ii) size of CP, as determined by the deduced amino acid sequence data, (iii) genome structure and organization (number and relative location of the open reading frames (ORFs)), (iv) amino acid sequence of the relevant gene products (CP, RdRp, HSP70h), differing by more than 25%, (v) vector species and specificity, (vi) magnitude and specificity of the natural and experimental host range, and (vii) cytopathological features (aspect of inclusion bodies and origin of cytoplasmic vesicles) [14].

Specifically, the genus *Ampelovirus* comprises species with linear particles of 1400–2000 nm long. Their genome is a monopartite positive sense, single-stranded RNA of 13.7–18.5 kb with a number of ORFs varying between 7 and 12 [9]. The members of this genus are divided into two subgroups accommodating, respectively, seven species with large (15,000 to over 18,000 nt) and complex (9 to 12 ORFs) genomes (subgroup I), and five species with smaller (13,000–14,000 nt) and simpler (7 ORFs) genomes (subgroup II) [9,16].
The majority of the *Ampelovirus* species have been identified from woody hosts (grapevine, *Prunus* sp., figs) and pineapple. Although their pathogenicity sometimes remains unclear, several members are reported to induce a diverse range of symptoms, while others are reported to have no association with symptoms. Natural vectors are *Pseudococcidae* mealybugs and soft scale insects, which transmit viruses in a semipersistent manner [9,12]. None of the known ampeloviruses are transmitted through seeds or mechanically. All ampeloviruses persist in plant parts used for vegetative propagation and are disseminated with them over long distances. The geographical distribution is therefore usually wide [9].

In this study, we identify for the first time and characterize the genome organization of novel ampeloviruses infecting cassava plants originating from diverse geographical locations in central Africa and the southwestern Indian Ocean islands. We also provide a first analysis of the diversity of these agents and of the evolutionary forces acting on them.

## 2. Materials and Methods

### 2.1. *Origin of the Analyzed Cassava Samples and High-Throughput Sequencing*

Field surveys were conducted from March to May 2016 in the South Kivu province (DR Congo). Cuttings were collected on cassava landraces plants showing foliar CBSD-like symptoms. These cuttings did not express the foliar symptoms anymore after growing in the greenhouse for six months. In addition, leaf samples from plants of different landraces were collected from the germplasm collections of the Conseil Départemental (DARTM in Mayotte), of CIRAD (Reunion) and of FOFIFA (Madagascar) between 2015 and 2017 (Table 1).

<table><caption>**Table 1.** Identity and origin of cassava landraces used and indication on the high-throughput sequencing (HTS) protocol performed on these samples.</caption><thead><tr><th>Cassava Landraces</th><th>Country of Origin</th><th>Collection Date DD/MM/YY</th><th>Symptoms</th><th>Sequencing Strategy and Target</th><th>Isolates Detected</th><th>Other Virus Detected</th></tr></thead><tbody><tr><td>Menatana</td><td>Madagascar</td><td>18/10/2017</td><td>Asymptomatic in greenhouse</td><td>VANA</td><td>MG-Men-241, MG-Men-9</td><td>-</td></tr><tr><td>Miandrazaka</td><td>Madagascar</td><td>18/10/2017</td><td>Asymptomatic in greenhouse</td><td>VANA</td><td>MG-Mia-10, MG-Mia-362, MG-Mia-403, MG-Mia-2-2</td><td>-</td></tr><tr><td>Long Java</td><td>Reunion</td><td>25/02/2015</td><td>Asymptomatic in field</td><td>dsRNAs</td><td>RE-Ljv</td><td>-</td></tr><tr><td>Reunion</td><td>Mayotte</td><td>26/03/2015</td><td>Symptoms of CMD and CBSD</td><td>dsRNAs</td><td>MY-Ren-5, MY-Ren-8</td><td>EACMV+UCBSV</td></tr><tr><td>6 Mois Blanc</td><td>Mayotte</td><td>26/03/2015</td><td>Symptoms of CMD</td><td>dsRNAs</td><td>MY-6mb-4, MY-6mb-6</td><td>EACMV</td></tr><tr><td>Nambiyombiyo</td><td>DR Congo</td><td>13/01/2016</td><td>Asymptomatic in greenhouse</td><td>Total RNA</td><td>CG-Nmb</td><td>CBSV</td></tr><tr><td>Kahunde</td><td>DR Congo</td><td>13/01/2016</td><td>Asymptomatic in greenhouse</td><td>Total RNA</td><td>CG-Kah</td><td>-</td></tr></tbody></table>

VANA: virion-associated nucleic acid; dsRNAs: double-stranded RNAs; EACMV: East African cassava mosaic virus; UCBSV: Ugandan cassava brown streak virus; CBSV: cassava brown streak virus; CMD: cassava mosaic disease; CBSD: cassava brown streak disease.

For the two samples from DR Congo, total RNA was extracted from cassava leaves using the RNeasy Plus Plant Mini Kit (QIAGEN®, Hilden, Germany) according to the manufacturer's instruction and quality tested using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). DNase treatment was applied using Amplification grade DNASE I (Life Technologies, California, USA) according to the manufacturer's instructions. Ribosomal RNAs were depleted using a RibZero plant leaf kit for RNA-seq (Life Technologies Limited, Paisley, UK), and libraries were prepared following the manufacturer's instructions using a TrueSeq stranded total RNA kit (Illumina,
New York, NY, USA). The RNA libraries were sequenced on a Nextseq 500 sequencing machine at the University of Liege (Liege, Belgium), with a read length of 2 × 75 nt.

For the three samples from Mayotte and Reunion, double-stranded RNAs (dsRNAs) were extracted from leaf samples using the procedure of Marais et al. [17]. Purified dsRNAs were converted to cDNA and amplified using a random whole-genome amplification procedure [17] and finally sequenced using Illumina Myseq technology (2 × 250 nt paired reads) on the Genotoul INRAE platform (Toulouse, France).

For the two samples from Madagascar, HTS was performed using a virion-associated nucleic acid (VANA)-based metagenomics approach as described by Palanga et al. [18].

## 2.2. *Bioinformatic Analyses*

Samples from Madagascar, Mayotte and Reunion: Following demultiplexing and quality trimming, reads were assembled into contig using CLC Genomic Workbench 8.5.1 and following versions. Contigs were annotated by BLASTx analysis against the GenBank protein database. Contigs were then manually assembled into scaffolds, which were consolidated and extended by several rounds of mapping of reads to yield the finalized scaffolds sequences. For two of the isolates (RE-Ljv and MY-Ren), reads corresponding to deletion events generating defective molecules were identified during this assembly process and scaffolds corresponding to these DI RNAs were therefore also reconstructed.

Samples from DR Congo: Sequence analyses were done using the Geneious 11.1.3 environment (www.geneious.com and embedded plugins). Data were submitted to a pre-processing step consisting of setting paired reads, and trimming adaptor sequences and low-quality reads (by BBDUK). Unique sequences were generated as clean reads by merging paired reads (by BBMerge) and removing duplicate reads (Dedupe). SPADES was used for de novo assembly. Reconstructed contigs were screened against the Refseq viral database retrieved from NCBI (12 October 2018) using BLASTn and BLASTx searches and annotated through tBLASTx. Viral contigs were further analyzed directly on the NCBI site using a BLASTn and BLASTx search with standard parameters.

Identification of functional gene domains for all samples was done by submitting predicted open reading frames (ORFs) directly on the NCBI's conserved domain search tool (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) (accessed from October to December 2018).

Multiple sequence alignments were built with ClustalW embedded in Mega X. Maximum likelihood phylogenetic trees were reconstructed in Mega X using the GTR+GI model for the alignment of nucleotide sequences and a Poisson model with uniform distribution for amino acid sequence alignments. Bootstrap analysis (100 replicates) was performed in order to evaluate the stability and significance of branches. Phylogenetic trees of amino acid sequences were reconstructed using only isolates for which the amino acid sequence of the analyzed protein was complete or nearly complete.

Putative recombination events were detected and evaluated using the RDP4 program (version 4.99) [19] and a ClustalW-built multiple alignment of the complete genome sequences.

## 2.3. *Confirmatory RT-PCR and Sequencing*

The presence of the detected viral contigs of MEaV-1 (isolates CG-Nmb and CG-Kah) in the corresponding plants was confirmed by RT-PCR using specific primer pairs (Table 2) designed according to the consensus sequence of the reconstructed viral contigs. This also allowed to confirm the sequence of the contigs of isolates mentioned.
<table><caption>**Table 2.** RT-PCR primer pairs used for the confirmation of the viral sequences in the samples.</caption><thead><tr><td>PRIMER</td><td>TARGET ORF</td><td>SEQUENCES (5'==> 3')</td><td>Tm</td></tr></thead><tbody><tr><td>12.495 F cp</td><td>CP ORF5</td><td>AATTTGGGAGGAGTGCGACC</td><td>60.3</td></tr><tr><td>12.715 R cp</td><td>CP ORF5</td><td>AGACCGACTTGTGCTACTCTTG</td><td>60.0</td></tr><tr><td>11.735 F hsp90</td><td>HSP90 ORF4</td><td>GTCCCGGCCTACATGCAATT</td><td>60.8</td></tr><tr><td>11.959 R hsp90</td><td>HSP90 ORF4</td><td>ACCGCCAGTCAACTCTCGTA</td><td>60.9</td></tr><tr><td>9.021 F hsp70</td><td>HSP70 ORF4</td><td>TATGGTTGTTGCTCGCGACT</td><td>60.0</td></tr><tr><td>9.292 R hsp70</td><td>HSP70 ORF4</td><td>CTGACAAACCAGCAGCAGTTG</td><td>60.3</td></tr><tr><td>7.244 F rdrp</td><td>RdRp ORF1B</td><td>GGACAACCTCCGAAACCGTAT</td><td>60.1</td></tr><tr><td>7.544 R rdrp</td><td>RdRp ORF1B</td><td>CTTTCGCTGCCATTGGTGTC</td><td>60.1</td></tr></tbody></table>

## 2.4. Characterization of Defective RNA (D-RNA Molecules)

### 2.4.1. RNA Isolation, cDNA Synthesis and Amplification of D-RNAs

Total RNA extracted from cassava accessions in which defective RNA molecules have been identified (RE-Ljv and MG-Mena-242) was purified and used for cDNA synthesis using the ProtoScript II First Strand cDNA Synthesis Kit (NEB) following the manufacturer's instructions. Briefly, 0.5 µg of cDNA was mixed with 2 µL of random hexamers primers (60 µM), heated at 65 °C for 5 min, and immediately cooled on ice for 5 min. The cooled solution was mixed with 10 µL of 2x ProtoScript II reaction mix and 2 µL of 10x ProtoScript II enzyme mix, mixed well and centrifuged for 5 s. This reaction mixture was incubated at 25 °C for 5 min and 42 °C for 1 h. The reaction was later terminated by incubation at 80 °C for 5 min.

Specific primer pairs designed on genome regions flanking the deletion in each of the four defective RNA molecules (Table 3) were used to amplify a region spanning the deletion of the D-RNAs. PCR amplification was performed in 20 µL total volume with 10 µL of Q5 High-Fidelity Master Mix (NEB), 2 µL of cDNA template, 0.5 µL of 10 µM of each primer and nuclease-free water for the remaining volume. The amplification program was set using 98 °C for 3 min of initial denaturation, followed by 35 cycles of 98 °C for 10 s, 57 °C for 10 s, 72 °C for 30 s and a final extension at 72 °C for 5 min. The amplified PCR products were analyzed on 1% agarose gel to validate the presence of amplicons of the expected size. The amplicons were gel-purified using the Monarch DNA Gel Extraction Kit (NEB) and cloned into pJET1.2 (Thermo Fisher Scientific).

<table><caption>**Table 3.** RT-PCR primer pairs used for the characterization of the defective RNA molecules.</caption><thead><tr><td>Primers Name</td><td>Sequences (5'==>3')</td><td>Expected Amplicon Size</td></tr></thead><tbody><tr><td>MEaV-1-RE-Ljv-D-RNA1-F</td><td>ACATCTAAATGCTAACGAACGAAGAG</td><td rowspan="2">380 bp</td></tr><tr><td>MEaV-1-RE-Ljv-D-RNA1-R</td><td>CAACGCCAGAATCTTCGTACA</td></tr><tr><td>MEaV-1-RE-Ljv-D-RNA2-F</td><td>AGGCTTTCGACAGTGAAGAAGTG</td><td rowspan="2">330 bp</td></tr><tr><td>MEaV-1-RE-Ljv-D-RNA2-R</td><td>CGTAGCCATACTGAAGGATAGCA</td></tr><tr><td>MEaV-2-RE-Ljv-D-RNA3-F</td><td>TGGAAGCCGCTGGTAAACTACA</td><td rowspan="2">400 bp</td></tr><tr><td>MEaV-2-RE-Ljv-D-RNA3-R</td><td>CAAGCACGTTCAATATTAGGAATAGTAC</td></tr><tr><td>MEaV-2-MG-Mena-D-RNA4-F</td><td>AGACATATGAAAGAGTTGCATTGGTG</td><td rowspan="2">500 bp</td></tr><tr><td>MEaV-2-MG-Mena-D-RNA4-R</td><td>ACCTACAAATAATTTCGCTCGTCTG</td></tr></tbody></table>

### 2.4.2. Sequence Analysis and Alignments

The amplicons clones were Sanger-sequenced using pJET1.2 forward and reverse primers. Two clones were sequenced for each defective RNA molecule and no sequence variation was observed among them. The obtained sequences were aligned to reference
parental full-length segments (RE-Ljv, MG-Mena-9) using the Clustal Omega multiple sequence alignment tool (online version, October 2020), allowing to confirm the existence of molecules containing the expected deletions.

## 3. Results

### 3.1. Identification of Closteroviridae Members in Cassava

Thirteen contigs or scaffolds of more than 10 kilobases of *Closteroviridae* isolates could be reconstructed using the sequencing reads from the seven analyzed landraces (Table 4). The lengths of the contigs, which represent large parts of the corresponding genomes, ranged from 10,417 to 13,770 nucleotides (nt). Some of the analyzed plants contained multiple contigs corresponding to different isolates (Table 4). Thorough screening of these reconstructed contigs against viral reference databases revealed nucleotide and amino acid sequence homologies with members of the family *Closteroviridae* and, in particular, with the genus *Ampelovirus*. In addition to these previously unknown viruses, the results further indicated the presence of CBSVs and/or CMGs in three cassava samples (Table 1). Cassava brown streak virus (CBSV, isolate KOR6, GU563327) was detected in one of the plants from the DR Congo, while Uganda cassava brown streak virus (UCBSV, isolate KM) [20] and/or EACMV were detected in samples from Mayotte.

<table><caption>**Table 4.** Lengths, number of reads integrated, percent of total reads and average coverage for *Closteroviridae* contigs from the various analyzed cassava samples.</caption><thead><tr><th>Cassava Landraces</th><th>Isolate Name <sup>1</sup></th><th>Accession Number (Genbank)</th><th>Virus Name</th><th>Contig Length [Nt]</th><th>Number of Reads Mapped</th><th>Percent Of Total Reads</th><th>Average Genome Coverage</th><th>Internal Un-determined Nucleotides</th><th>Internal Gaps</th></tr></thead><tbody><tr><td rowspan="2">Menatana</td><td>MG-Men-241</td><td>MT773584</td><td>MEaV-1</td><td>13,528</td><td>38,475</td><td>12.1</td><td>506.8X</td><td>0</td><td>0</td></tr><tr><td>MG-Men-9</td><td>MT773591</td><td>MEaV-2</td><td>13,752</td><td>10,945</td><td>4.10</td><td>153.6X</td><td>0</td><td>0</td></tr><tr><td rowspan="3">Miandrazaka</td><td>MG-Mia-10</td><td>MT773592</td><td>MEaV-2</td><td>13,770</td><td>26,987</td><td>4.8</td><td>368X</td><td>377</td><td>1</td></tr><tr><td>MG-Mia-362</td><td>MT773594</td><td>MEaV-2</td><td>12,007</td><td>19,421</td><td>8.24</td><td>322X</td><td>0</td><td>0</td></tr><tr><td>MG-Mia-403</td><td>MT773590</td><td>MEaV-1</td><td>10,518</td><td>26,917</td><td>11.4</td><td>480.9X</td><td>0</td><td>0</td></tr><tr><td rowspan="2">Long java</td><td>MG-Mia-2-2</td><td>MT773596</td><td>MEaV-2</td><td>10,417</td><td>7676</td><td>1.4</td><td>136.4X</td><td>843</td><td>1</td></tr><tr><td>RE-Ljv</td><td>MT773586</td><td>MEaV-1</td><td>13,640</td><td>21,294</td><td>18.5</td><td>312.5X</td><td>0</td><td>0</td></tr><tr><td rowspan="2">Reunion</td><td>MY-Ren-5</td><td>MT773589</td><td>MEaV-1</td><td>13,172</td><td>1976</td><td>3.4</td><td>29.6X</td><td>72</td><td>1</td></tr><tr><td>MY-Ren-8</td><td>MT773595</td><td>MEaV-2</td><td>13,762</td><td>377</td><td>0.7</td><td>5.4X</td><td>2889</td><td>14</td></tr><tr><td rowspan="2">6 Mois blanc</td><td>MY-6mb-4</td><td>MT773585</td><td>MEaV-1</td><td>13,615</td><td>4992</td><td>12.7</td><td>73.3X</td><td>288</td><td>4</td></tr><tr><td>MY-6mb-6</td><td>MT773593</td><td>MEaV-2</td><td>13,764</td><td>2631</td><td>6.7</td><td>39.2X</td><td>990</td><td>9</td></tr><tr><td>Nambiyombiyo</td><td>CG-Nmb</td><td>MT773587</td><td>MEaV-1</td><td>13,612</td><td>7502</td><td>0.12</td><td>42X</td><td>0</td><td>0</td></tr><tr><td>Kahunde</td><td>CG-Kah</td><td>MT773588</td><td>MEaV-1</td><td>13,616</td><td>15,351</td><td>0.23</td><td>75X</td><td>11</td><td>1</td></tr></tbody></table>

The isolates were named considering two letters of the code of the country where the landraces were collected, three letters of the landrace short acronym, and a number (when multiple contigs from the same landrace were obtained).

### 3.2. Genome Annotation

The sequences comparisons between the various isolates identified two clearly separated groups (Figure 1). The nearly complete genome of a representative isolate was selected for each group (respectively, CG-Nmb for the MEaV-1 and MG-Men-9 for MEaV-2) and was analyzed in detail to determine the genome organization. As shown in Figure 1, this analysis predicted in both cases seven open reading frames encoding for polypeptides with a mass ranging from 6 to 258 kDa (Table 5). The genomic organization, the size, and the number of ORFs identified revealed a similarity with the subgroup II ampeloviruses. The recognized members of this subgroup are grapevine leafroll-associated virus 4 (GLRaV-4), pineapple mealybug wilt-associated virus 1 and 3 (PMWaV-1 and -3),

<footnote n="1">¹ Isolates were named considering two letters of the code of the country where the landraces were collected—three letters of the landrace short acronym—and a number (when multiple contigs from the same landrace were obtained)</footnote>
air potato Ampelovirus 1 (AiPoV1), and plum bark necrosis and stem pitting-associated virus (PBNSPaV) [15].

The ORFs *1a* and *1b* encode the replication-associated proteins (Figure 1). As expected, a conserved domain search identified the following two replication-associated domains in the *1a* protein: a methyltransferase motif (*MTR*; pfam01660) in the N-terminal part, and a viral helicase (*HEL*; superfamily 1, pfam01443) in the C-proximal part. The RNA methyltransferase domain was found in a wide range of ssRNA viruses and is known to be involved in mRNA capping.

The *1b* ORF encodes the RNA-dependent RNA polymerase (*RdRp* 2 superfamily, pfam00978) expressed through a +1 ribosomal frameshift, as is seen in other members of the *Closteroviridae* family. Together with the *MTR* domain located at the *1a* protein N-terminus, the *RdRp* forms the defining unique feature of the alpha-like ensemble of viruses [21].

Downstream of the polymerase, a small ORF2 encodes a hypothetical protein with a predicted molecular mass varying from 5 to 7 kDa depending on the isolate, and which is lacking known functional domains. The numbers and the sizes of the hypothetical proteins occurring downstream of the polymerase in the family *Closteroviridae* are highly variable. For example, Actinidia virus 1 (AcV-1) has three hypothetical proteins of predicted molecular masses of 13.6, 25.4 and 5.7 kDa [22]. However, in subgroup II ampeloviruses, only one hypothetical protein between 5 and 6 kDa is reported, sometimes overlapping with the 5'-proximal region of ORF3. The small size and the high proportion of the hydrophobic amino acids of the proteins encoded by the MG-Men-9 and CG-Nmb isolates (36/48 aa and 42/67 aa, respectively) could indicate that they are similar to the transmembrane proteins often present in the *Closteroviridae* members, where they have been proposed to function as cell-to-cell movement proteins [23].

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic representation of the genomic organization of representative isolates CG-Nmb (MEaV-1) and MG-Men-9 (MEaV-2) (top) and structure of the defective variants (dRNA) identified (bottom). The RNA genome is drawn as a black line and the predicted open reading frames (ORFs) are represented by colored rectangles. Annotations and ORF numbers are given inside and above rectangles, respectively. Abbreviations: *MTR*—methyltransferase, *HEL*—helicase, *RdRp*—RNA-dependent RNA polymerase, *HSP70h*—heat shock protein 70 homolog, *HSP90h*—heat shock protein 90 homolog, *CP*—coat protein, *CPd*—minor coat protein.</figcaption></figure>

<table><caption>**Table 5.** Genome length, position of ORFs and size of the encoded proteins for the representative isolates of the two species CG-Nmb (MEaV-1) and MG-Men-9 (MEaV-2).</caption><thead><tr><th rowspan="2">Isolates</th><th rowspan="2">Genome Length (Nt)</th><th colspan="7">Length of Open Reading Frames (Orfs) (Nt) and Molecular Mass of Encoded Proteins (KDa)</th></tr><tr><th>1a</th><th>1b</th><th>P5-P7</th><th>HSP70h</th><th>HSP90h</th><th>CP</th><th>ORF6</th></tr></thead><tbody><tr><td rowspan="2">CG-Nmb (MEaV-1)</td><td rowspan="2">13,616 nt</td><td>6789 nt</td><td>1503 nt</td><td>204 nt</td><td>1707 nt</td><td>1686 nt</td><td>774 nt</td><td>717 nt</td></tr><tr><td>254.6 kDa</td><td>56.6 kDa</td><td>7.6 kDa</td><td>62.2 kDa</td><td>63.4 kDa</td><td>27.7 kDa</td><td>27.1 kDa</td></tr><tr><td rowspan="2">MG-Men-9 (MEaV-2)</td><td rowspan="2">13,752 nt</td><td>6975 nt</td><td>1593 nt</td><td>147 nt</td><td>1608 nt</td><td>1698 nt</td><td>816 nt</td><td>639 nt</td></tr><tr><td>258.8 kDa</td><td>60.4 kDa</td><td>5.3 kDa</td><td>58.4 kDa</td><td>64.3 kDa</td><td>29.3 kDa</td><td>24 kDa</td></tr></tbody></table>
ORF3 encodes a heat shock protein 70 homolog (*HSP70h*; cd10170) from the NBD-sugar-kinase-HSP70-actin superfamily. The functions postulated for *HSP70h* are as follows: the mediation of cell-to-cell movement through the plasmodesmata, involvement in the assembly of multisubunit complexes for genome replication and/or subgenomic RNAs synthesis, and the assembly of viral particles [9]. The size of this ORF is variable, being 62.2 kDa for MEaV-1 and 58.49 kDa for MEaV-2, as shown in the Supplementary Material (Table S1).

ORF4 codes for a heat shock protein 90 homolog (*HSP90h*, pfam03225), which is found to partially overlap with the 3'-proximal region of ORF3 and the 5'-proximal region of ORF5 as is the case for subgroup II ampeloviruses.

The coat protein (Closter coat superfamily, pfam01785) is encoded by ORF5.

It has been shown that downstream of the coat protein, some members of the family *Closteroviridae* have a variable number of ORFs encoding accessory proteins, some of which may have a functional conserved domain [24], while others show no or only very limited identity to any other proteins. The translated amino acid sequences of the protein encoded by the ORF6 have shown similarity with the p24 protein of the pineapple mealybug wilt-associated virus-1 (25.1%) and the minor coat protein (*CPd*) of the grapevine leafroll-associated virus-5 (20.8%) from the TrEMBL database of UniProt. This suggests that, although they share <20% identity with their respective, the proteins coded by the ORF6 could be considered as minor coat proteins. This ORF6 overlaps partially with the 3'-proximal region of ORF5.

All of the isolates sequenced share the same genome organization, although the contigs obtained for the isolates MG-Mia-403 and MG-Mian-2-2 were incomplete, missing ORF5 and ORF6.

The reads belonging to four defective RNA (D-RNAs) molecules were identified (Figure 1). Three of them were detected in the long Java accession (RE-Ljv) originating from the Réunion island (two for MEaV-1 (RE-Ljv-D-RNA-1 and RE-Ljv-D-RNA-2, accession numbers MW306827 and MW306828, respectively), one for MEaV-2 (RE-Ljv-D-RNA-3, accession number MW306829) and the remaining one in the Menatana accession originating from Madagascar (MG-Men-D-RNA-4 belonging to MEaV-2, accession number MW306830)). The deleted zones cover different genome regions and consist of most of the ORF5 and the ORF6 for the RE-Ljv-D-RNA-1, most of ORF3, ORF4, ORF5 and the start of ORF6 for the RE-Ljv-D-RNA-2, and most of the ORF3 and ORF4 for the RE-Ljv-D-RNA-3. In the case of the MG-Men-D-RNA-4, the deletion is at the beginning of the 1b ORF and does not alter the coding frame. The presence of molecules bearing these deletions in the original plant material was validated by the sequencing of amplicons spanning each of the individual deletions (see below). The proportion of deleted RNAs, as compared to the non-deleted genomic RNA, was estimated by comparing the read counts for the deletion borders and for the corresponding regions of the undeleted genome. This proportion was found to be variable, corresponding to about 4%, 37%, 27% and 3.8% for RE-Ljv-D-RNA-1, RE-Ljv-D-RNA-2, RE-Ljv-D-RNA-3 and MG-Men-D-RNA-4, respectively.

### 3.3. Phylogenetic Analysis of *Cassava Isolates with Members of the Family Closteroviridae*

A multiple alignment of the *HSP70h* protein of the thirteen isolates, and of the corresponding protein of known members of the family *Closteroviridae*, was used to perform pairwise amino acid comparisons and to generate a phylogenetic tree that allowed to first address the taxonomic position of the cassava isolates in the family *Closteroviridae*.

The resulting tree (Figure 2) placed them with strong bootstraps support in a clade with members of the genus *Ampelovirus*, specifically in the subgroup II clade adjacent to GLRaV-4 (NC_016416) and PMWaV-1 and 3 (NC_010178 and DQ_399259), confirming the taxonomical relatedness of these cassava isolates to the subgroup II of the *Ampelovirus* genus.

The molecular identity at the amino acid level for the *RdRp*, *HSP70h* and *CP* of the cassava isolates with members of the family *Closteroviridae* is systematically lower than 55% (Table 6), well below the 75% molecular species demarcation criterion for these
taxonomically relevant genes in this family [9]. This shows that the isolates from cassava presented in this study are novel agents that do not belong to a known species in the family *Closteroviridae*.

<figure><img src="image_5.png" /></figure>

**Figure 2.** Phylogenetic analysis of the aligned amino acid sequences of the *HP70h* (ORF3) of the thirteen isolates from cassava and of selected members of the family *Closteroviridae* (see the Supplementary Material Table S2 for detailed information on these viruses). Green arrows indicate isolate sequences obtained in this study. Maximum likelihood phylogenetic tree was reconstructed in Mega X using Poisson model with uniform distribution for amino acid sequence alignments. Bootstrap values are indicated at the main branch nodes. The bar represents the number of amino acid substitution per site.
<table><caption>**Table 6.** Amino acid sequence identity of the RdRp, HSP70h and CP of MEaV-2 (isolate MG-Men-9) with those of other selected *Closteroviridae* members.</caption><thead><tr><th rowspan="2">Genus</th><th rowspan="2">Representative Members</th><th colspan="3">Proteins</th></tr><tr><th>RdRp</th><th>HSP70H</th><th>CP</th></tr></thead><tbody><tr><td><i>Ampelovirus</i></td><td>MEaV-1 (isolate CG-Nmb)</td><td>66.5%</td><td>69.9%</td><td>80.1%</td></tr><tr><td><i>Ampelovirus</i></td><td>Grapevine<br>leafroll-associated virus 4<br>(NC_016416)</td><td>41.7%</td><td>51.6%</td><td>47.6%</td></tr><tr><td><i>Closterovirus</i></td><td>Citrus tristeza virus<br>(NC_001661)</td><td>27.2%</td><td>25%</td><td>12.5%</td></tr><tr><td><i>Crinivirus</i></td><td>Potato yellow vein virus<br>(NC_006062 and<br>NC_006063)</td><td>25.7%</td><td>26.5%</td><td>17.2%</td></tr><tr><td><i>Velarivirus</i></td><td>Grapevine<br>leafroll-associated virus 7<br>(NC_016436)</td><td>27.5%</td><td>27.7%</td><td>13.6%</td></tr></tbody></table>

### 3.4. Diversity of the New Isolates

To investigate the phylogenetic relationships between the various reconstructed *Closteroviridae* genomic sequences from cassava, multiple alignments of the near-complete genomic sequences and of the amino acid sequences of the three taxonomically relevant proteins were used to reconstruct phylogenetic trees (Figure 3a–d).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Phylogenetic trees reconstructed using the amino acid sequences of the following three taxonomically relevant proteins for the family *Closteroviridae*: (a) RdRp; (b) HSP70h; (c) CP; and (d) the whole genome nucleotide sequences. Maximum likelihood phylogenetic trees were reconstructed in Mega X using the GTR+GI model for nucleotide sequences alignments and a Poisson model with uniform distribution for amino acid sequence alignments. Only bootstrap values of more than 70% are mentioned on nodes. The scales provide branch distance for the given number of substitutions per site.</figcaption></figure>
From these analyses, it is clear that the cassava isolates form two well-separated phylogenetic groups (MEaV-1 and MEaV-2) that show 41.4 ± 0.3% average nucleotide divergence based on the whole genome and, respectively, 32.3 ± 1.7% (*RdRp*), 30.8 ± 2.0% (*HSP70h*) and 19 ± 2.2% (*CP*) average amino acid divergence (Table 6 and see the Supplementary Material Tables S3–S6 for additional information).

The first group (MEaV-1) clusters two isolates identified in two cassava landraces from Madagascar (MG-Men-241 and MG-Mia-403), one isolate identified in an accession from Reunion (RE-Ljv), two isolates identified in accessions from Mayotte (MY-6mb-4 and MY-Ren-5) and two isolates from DR Congo (CG-Nmb and CG-Kah). The following two isolates from this group are divergent from the others and form a distinct subgroup: MY-Ren-5 and MG-Mia-403. Indeed, the MEaV-1 average pairwise divergence is found to be 17.1 ± 0.2% based on the whole genome, and, respectively, 7.1 ± 0.7% (*RdRp*), 6.0 ± 0.6% (*HSP70h*) and 7.4 ± 1.0% (*CP*) for the various proteins (Table 6). However, if excluding the divergent isolates MG-Mia-403 and MY-Ren-5, these values fall down to 7.9 ± 0.2% for the whole genome and to 0.8 ± 0.4% (*RdRp*), 2.3 ± 0.4% (*HSP70h*) and 3.2 ± 0.6% (*CP*) for the various proteins, respectively. The average pairwise distance between these two divergent isolates, MG-Mia-403 and MY-Ren-5, is only 2.3% for the whole genome (Table 6).

The second group (MEaV-2) is composed of four isolates identified in two landraces from Madagascar (MG-Men-9, MG-Mia-2-2, MG-Mia-10 and MG-Mia-362) and two isolates identified in two landraces from Mayotte (MY-6mb-6 and MY-Ren-8). The isolate MG-Mian-2-2 forms a sub-cluster of its own. The MEaV-2 average pairwise divergence is 12.3 ± 0.1% for the whole genome, and, respectively, 7.3 ± 0.7% (*RdRp*) and 2.8 ± 0.4% (*HSP70h*) for the various genes (Table 7) (the average *CP* divergence was not computed because the *CP* gene is not covered by the MG-Mian-2-2 contig). If excluding the divergent MG-Mian-2-2 isolate, the average pairwise divergence for the whole genome falls to 3.8 ± 0.1%, and to 1.1 ± 0.3% (*RdRp*), 1.0 ± 0.3% (*HSP70h*) and 3.1 ± 0.8% (*CP*) for the various proteins, respectively (Table 7).

<table><caption>**Table 7.** Intergroup and intragroup average pairwise divergence and standard deviation calculated for the three taxonomically relevant proteins (*RdRp*, *HSP70h* and *CP*) and for the nearly complete genomes.</caption><thead><tr><th></th><th></th><th>Within MEaV-2 with Mian 2-2</th><th>Within MEaV-2 without Mian2-2</th><th>Within MEaV-1 with Divergents</th><th>Within MEaV-1 without Divergents</th><th>Within MEaV-1 Divergents</th><th>Between MEaV-1 and MEaV-2</th></tr></thead><tbody><tr><td>RdRp</td><td>aa divergences</td><td>7.3 +/- 0.7%</td><td>1.1 +/- 0.3%</td><td>7.1 +/- 0.7%</td><td>0.8 +/- 0.4%</td><td>0.8%</td><td>32.3 +/- 1.7%</td></tr><tr><td>HSP70h</td><td>aa divergences</td><td>2.8 +/- 0.4%</td><td>1.0 +/- 0.3%</td><td>6.0 +/- 0.6%</td><td>2.3 +/- 0.4%</td><td>0.7%</td><td>30.8 +/- 2.0%</td></tr><tr><td>CP</td><td>aa divergences</td><td>na</td><td>3.1 +/- 0.8%</td><td>7.4 +/- 1.0%</td><td>3.2 +/- 0.6%</td><td>na</td><td>19 +/- 2.2%</td></tr><tr><td>Genome</td><td>nt divergences</td><td>12.3 +/- 0.1%</td><td>3.8 +/- 0.1%</td><td>17.1 +/- 0.2%</td><td>7.9 +/- 0.2%</td><td>2.3%</td><td>41.4 +/- 0.3%</td></tr></tbody></table>

MEaV-1 divergents: MY-Ren-5 and MG-Mia-403. na: information not available as the contig/scaffold does not extend to the region concerned.

Comparisons across the trees based on the *RdRp*, *HSP70h* and *CP* proteins reveal variations in the branching topology that point to recombination events.

Considering the taxonomic differentiation between these two isolate groups (MEaV-1 and MEaV-2), largely supported by the phylogenetic trees, and the average divergence rates well over the 25% species demarcation criterion [14,15] for two of the three proteins (Table 6), it can be considered that the analyzed isolates form two distinct species, tentatively named "Manihot esculenta-associated ampelovirus 1 and 2" (MEaV-1 and MEaV-2), respectively.

MEaV-1 is detected from the four regions included in the study, while MEaV-2 is observed only from Madagascar and Mayotte (Figure 4).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Geographic repartition of the two novel cassava ampelovirus species (MEaV-1 and MEaV-2) according to the country of origin of positive samples.</figcaption></figure>

**3.5. Recombination Analysis**

Intraspecific and interspecific recombination events were investigated among the cassava isolates, and the recombination events supported that more than four of the nine algorithms integrated in RDP4 were considered as possible events [19].

Only one such recombination event was reliably detected using this criterion (Figure 5). This recombination event involves only the MEaV-2 isolates (two isolates MG-Mia-362 and MG-Mia-10 from Madagascar, and one isolate MY-6mb-6 from Mayotte) and is detected by seven out of the nine algorithms, with a strong probability (p-value = 2.75 × 10⁻⁹¹). The predicted recombination breakpoints are identified within ORF1a (nucleotide position 906 of the isolate MG-Mia-362) and within ORF4 (nucleotide position 8997 of the isolate MG-Mia-362).

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Recombination event detected among isolates from MeaV-2 by the RDP4 program.</figcaption></figure>
### 3.6. Validation of Defective RNAs

The presence of defective molecules was validated by the sequencing of RT-PCR amplicons spanning each of the identified deletions (Figure 6). The alignment of the four defective clones to their HTS-defective references, as shown in the Supplementary Material (Figures S7–S10), confirmed the presence of deletion zones at their predicted location. On the reference parental genomes (RE-Ljv and MG-Men-9), the deleted zones were found to be, respectively, 568 nt, 3761 nt, 648 nt and 444 nt for the D-RNA1, 2, 3 and 4.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Gel red-stained 2% agarose gel showing RT-PCR products from: (a) D-RNA1 expected amplicon size 380 bp; (b) D-RNA2 and D-RNA3, expected amplicon size 380 bp and 400 bp, respectively; and (c) D-RNA4, expected amplicon size 500 bp. M1: GeneRuler 1 kb plus ladder (Thermo Scientific, SM1343). M2: fast DNA ladder 10 kb (NEB N3238S).</figcaption></figure>

These findings confirm that the D-RNAs reported here are not artifacts from high-throughput sequencing analysis, but continuous sequences that exist in the analyzed cassava samples.

## 4. Discussion

We have reported here the genome characterization from field-grown cassava plants from the DR Congo, Madagascar, Mayotte and Reunion of thirteen new *Closteroviridae* isolates, belonging to two potentially new ampelovirus species.

The analysis of the sequences of the corresponding contigs (ranging between 10,417 and 13,752 nucleotides in length) revealed seven open reading frames. In subgroup II, the length of the untranslated regions (UTRs) varies between 218 and 353 nt at the 5' end, and between 125 and 132 nt at the 3' end [25–30]. The UTR sizes reported here for the various genomes vary between 210 and 253 nt for the 5'end, and between 157 and 166 nt for the 3' end. The comparison of these values, and the unambiguous assignation of MEaV-1 and MEaV-2 to the subgroup II of ampelovirus, show that the genomes reported here are unlikely to miss much more than a hundred nucleotides at their 5' end and probably even less at the 3' end, since some of the reported 3' UTRs are already longer than the longest ones reported to date in the subgroup. This is confirmed when comparing the complete genome sizes; GLRaV-4 (isolate Man086, KJ810572) and PMWaV-1 (isolate HN, KJ872494) are the two subgroup II ampeloviruses with, respectively, the longest and the shortest genome (13,858 nt and 13,069 nt, respectively) [25,27]. The longest contigs for the representative isolates of MEaV-1 and MEaV-2 are, respectively, 13,616 nt (CG-Nmb) and 13,752 nt (MG-Men-9) long, suggesting that near-complete genomes, comprising the totality of the coding sequences and missing limited terminal non-coding nucleotide sequences, have been obtained, although the extremities of the genomes were not determined by RACE.

Based on the divergence of their *CPs*, which was determined to be below the species threshold (19 ± 2.2%), MEaV-1 and MEaV-2 could be considered as a single species.
However, the official formulation of the taxonomic criterion for species demarcation in the family *Closteroviridae* and in ampelovirus genus is as follows: "Amino acid sequence of relevant gene products (polymerase, *CP* and *HSP70h*) differing by more than 25%" [14,15]. Therefore, there is an ambiguity on the fact of whether this criterion should be met for at least one of the three proteins or by all three simultaneously. Such a situation has previously been reported for Rehmannia virus 1 (ReV1) [31], for which the *RdRp* shows only 11% divergence with that of tobacco virus 1, while the *HSP70h* and the *CP* show 26% and 38% divergence, respectively. ReV1 was nevertheless accepted as a valid new *Closterovirus* species, suggesting that, as no recombination event was identified between them, MEaV-1 and MEaV-2 could similarly be considered as two distinct *Ampelovirus* species. This notion is further supported by the observation that all of the proteins diverge by more than 25%, between MEaV-1 and MEaV-2, except for the *CP* (and for a few comparisons involving one specific isolate, for the *HSP90h*).

Recombination is one of the mechanisms by which viruses evolve. Several studies have reported recombination events in *Closteroviridae* members [24,32,33]. In this study, we report a recombination event involving MEaV-2 isolates. No evidence of recombination between MEaV-1 and MEaV-2 was detected.

In most cases, the analyzed cassava plants contained complex mixed infections involving either several isolates from a single MEaV species, or isolates belonging to the two species. Such complex infection patterns are probably due to the vegetative propagation practices used for cassava cultivation. In two plants with these complex mixed infections, defective molecules were identified and confirmed by RT-PCR, cloning and sequencing. Defective RNAs belong to the category of virus-associated molecules that are not required for normal virus propagation, but can sometimes affect the accumulation of the helper virus and symptoms expression [34–36]. The presence of defective RNAs is reported in three genera of the *Closteroviridae* family [28,37–40], and is therefore not unexpected here. The considerable advances in virus characterization by HTS-based approaches are now revealing that, in addition to genomic and subgenomic RNAs, plants infected with viruses from the family *Closteroviridae* may contain several different subviral defective RNAs whose role is not known [35]. Further studies are needed to assess their impact on epidemiology and pathogenicity.

The sampling campaign conducted in the DR Congo collected stems/cuttings from cassava plants showing typical symptoms of CBSD on their leaves and stems. These symptoms consisted of a yellow blotchy pattern on mature lower leaves and brown-black marks ('streaks') on green stem portions. The collected cuttings were planted in a glasshouse and, after approximately six months, two of the planted cuttings that have previously expressed foliar symptoms in a farmer's field remained asymptomatic in the glasshouse and were used for the present study. The HTS results have confirmed the presence of CBSV in one of the samples, but not in the second one. Furthermore, the typical symptoms of CBSD and/or CMD that were observed on the samples from Mayotte are to be connected, in both cases, to the identification of the corresponding causal virus(es). No virus-like symptoms were observed on the samples from Reunion or from Madagascar, the Reunion sample being even selected in this study as a healthy control plant. Additional experiments of artificial inoculation and larger epidemiological studies are needed to clarify the symptomatology of MEaV-1 and MEaV-2, and to estimate the synergistic interaction between coinfecting viruses such as CBSVs and/or CMGs. The symptoms mentioned in this study consist only of foliar symptoms, while no observation was made on below-ground organs.

There is a need to develop a diagnostic test in order to be able to evaluate the distribution and prevalence of these new viral agents in other regions of the world, and to evaluate their impact on the yield. This evaluation of the phytosanitary risk should further be completed through wider surveys on symptomatic and asymptomatic plants in both agricultural and natural ecosystems to gain insight into the genetic variability of these new viruses, and their biological significance and impact [41].
**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/v13061030/s1; Figure S1: multiple sequence alignment of D-RNA1 clone (referred to in the figure as D-RNA1_Sanger) with the defective reference from HTS (MEaV-1_RE-LJV-D-RNA-1). The term "DELETIONZONE" has been inserted into the HTS-defective reference to identify the zone where the deletion is located. This zone consists of 568 nucleotide sequences and is indicated by the rectangle. Figure S2: multiple sequence alignment of D-RNA2 clone (referred to in the figure as D-RNA2_Sanger) with the defective reference from HTS (MEaV-1_RE-LJV-D-RNA-2). The term "DELETIONZONE" has been inserted into the HTS-defective reference to identify the zone where the deletion is located. This zone consists of 3761 nucleotide sequences and is indicated by the rectangle. Figure S3: multiple sequence alignment of D-RNA3 clones (two clones were sequenced and are referred to in the figure as D-RNA3a-Sanger and D-RNA3b-Sanger) with the defective reference from HTS (MEaV-2_RE-LJV-D-RNA3). The term "DELETIONZONE" has been inserted into the reference HTS-defective to identify the zone where the deletion has been identified. This zone consists of 648 nucleotide sequences and is indicated by the rectangle. Figure S4: multiple sequence alignment of D-RNA4 clone (referred to in the figure as D-RNA4-Sanger) with the defective reference from HTS (MEaV-2_MG-Mena-D-RNA-4). The term "DELETIONZONE" has been inserted into the reference HTS-defective to identify the zone where the deletion is located. This zone consists of 444 nucleotide sequences and is indicated by the rectangle. Table S1: length of ORF and molecular mass of 3 encoded proteins (RdRp, HSP70h and CP) of the 15 isolates of the cassava ampelovirus reported in this study. Table S2: accession numbers and abbreviations for selected Closterovirus used in the phylogenetic reconstruction of the HSP70h nucleotide sequences. Table S3: estimates of evolutionary divergence between amino acid (aa) sequences for the HSP70h protein among cassava ampeloviruses isolates reported in this study. Table S4: estimates of evolutionary divergence between amino acid (aa) sequences for the HSP70h protein among cassava ampeloviruses isolates reported in this study. Table S5: estimates of evolutionary divergence between amino acid (aa) sequences for the coat protein among cassava ampeloviruses isolates reported in this study. Table S6: estimates of evolutionary divergence between nucleotide sequences for the whole genomes of cassava ampelovirus isolates reported in this study.

**Author Contributions:** Conceptualization: S.M., T.C., Y.K. and F.M.; methodology: S.M., T.C., Y.K., E.B., S.S.; validation: S.M., T.C., A.G.B., J.-M.L., D.F.; formal analysis: Y.K., T.C., A.G.B., J.-M.L., D.F., S.S., C.F.; investigation: Y.K., E.B., A.M. resources: S.M., C.B., H.V., J.-M.L., E.B., A.M., S.R.; writing—original draft preparation: Y.K., S.M., T.C.; writing—review and editing: all authors; supervision: S.M., H.V., C.B., T.C., D.F., J.-M.L.; funding: S.M., C.B., H.V., T.C., D.F., J.-M.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** Y.K. and work in DR Congo were supported by the iCARE project funded by ARES (granted to Sébastien Massart, Claude Bragard and Hervé Vanderschuren). This study was co-funded in part by the European Union (ERDF, INTERREGV), the Conseil Régional de la Réunion and CIRAD, and conducted in part on the Plant Protection Platform (3P, IBISA). SS was supported by LEAP-Agri (A Long-term EU–Africa research and innovation Partnership on food and nutrition security and sustainable Agriculture), CASSANDRA FNRS-LEAP-288 project.

**Data Availability Statement:** The datasets of genome sequences generated and analyzed during this study are available in Genbank repository under the following accession numbers: MT773584, MT773585, MT773586, MT773587, MT773588, MT773589, MT773590, MT773591, MT773592, MT773593, MT773594, MT773595, MT773596, MW306827, MW306828, MW306829, MW306830.

**Acknowledgments:** The Genotoul Platform (INRAE, Toulouse, France) is gratefully acknowledged for Illumina Myseq sequencing. The authors thank Marc Fuchs (Cornell University, USA) for a valuable exchange on *Closteroviridae* taxonomy and species demarcation criteria.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Food and Agriculture Organization of the United Nations. FAOSTAT Statistical Database; FAO: Rome, Italy, 2017.

2. Burns, A.; Gleadow, R.; Cliff, J.; Zacarias, A.; Cavagnaro, T. Cassava: The Drought, War and Famine Crop in a Changing World. Sustainability **2010**, *2*, 3572–3607. [CrossRef]

3. Lobell, D.B.; Burke, M.B.; Tebaldi, C.; Mastrandrea, M.D.; Falcon, W.P.; Naylor, R.L. Prioritizing Climate Change Adaptation Needs for Food Security in 2030. *Science* **2008**, *319*, 607–610. [CrossRef] [PubMed]
4. Jarvis, A.; Ramirez-Villegas, J.; Campo, B.V.H.; Navarro-Racines, C. Is Cassava the Answer to African Climate Change Adaptation? *Trop. Plant. Biol.* **2012**, *5*, 9–29. [CrossRef]

5. Fermont, A.; Van Asten, P.; Tittonell, P.; Van Wijk, M.; Giller, K. Closing the cassava yield gap: An analysis from smallholder farms in East Africa. *Field Crop. Res.* **2009**, *112*, 24–36. [CrossRef]

6. Rey, C.; Vanderschuren, H.V. Cassava Mosaic and Brown Streak Diseases: Current Perspectives and Beyond. *Annu. Rev. Virol.* **2017**, *4*, 429–452. [CrossRef] [PubMed]

7. Owor, B.; Legg, J.P.; Okao-Okuja, G.; Obonyo, R.; Ogenga-Latigo, M.W. The effect of cassava mosaic geminiviruses on symptom severity, growth and root yield of a cassava mosaic virus disease-susceptible cultivar in Uganda. *Ann. Appl. Biol.* **2004**, *145*, 331–337. [CrossRef]

8. Dolja, V.V.; Karasev, A.V.; Koonin, E.V. Molecular biology and evolution of Closterovirus: Sophisticated Build-up of Large RNA Genomes. *Annu. Rev. Phytopathol.* **1994**, *32*, 261–285. [CrossRef]

9. Martelli, G.P.; Agranovsky, A.A.; Bar-joseph, M. Family Closteroviridae. In *Ninth Report of the International Committee on Taxonomy of Viruses*; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier: New York, NY, USA, 2012; pp. 987–1001.

10. Karasev, A.V. Genetic diversity and evolution of Closteroviruses. *Annu. Rev. Phytopathol.* **2000**, *38*, 293–324. [CrossRef]

11. Blouin, A.G.; Pearson, M.N.; Chavan, R.R.; Woo, E.N.Y.; Lebas, B.S.M.; Veerakone, S.; Ratti, C.; Biccheri, R.; MacDiarmid, R.M.; Cohen, D. Viruses of Kiwifruit (Actinidia species). *J. Plant. Pathol.* **2013**, *95*, 221–235.

12. Dey, K.K.; Sugikawa, J.; Kerr, C.; Melzer, M.J. Air potato (*Dioscorea bulbifera*) plants displaying virus-like symptoms are co-infected with a novel potyvirus and a novel ampelovirus. *Virus Genes* **2018**, *55*, 117–121. [CrossRef]

13. Verdin, E.; Wipf-Scheibel, C.; Gognalons, P.; Aller, F.; Jacquemond, M.; Tepfer, M. Sequencing viral siRNAs to identify previously undescribed viruses and viroids in a panel of ornamental plant samples structured as a matrix of pools. *Virus Res.* **2017**, *241*, 19–28. [CrossRef] [PubMed]

14. Fuchs, M.; Bar-Joseph, M.; Candresse, T.; Maree, H.J.; Martelli, G.P.; Melzer, M.J.; Menzel, W.; Minafra, A.; Sabanadzovic, S. ICTV Report Consortium ICTV Virus Taxonomy Profile: Closteroviridae. *J. Gen. Virol.* **2020**, *101*, 364–365. [CrossRef] [PubMed]

15. Candresse, T.; Fuchs, M. Closteroviridae. eLS **2020**, 1–10. [CrossRef]

16. Fuchs, M.; Bar-Joseph, M.; Candresse, T. ICTV Virus Taxonomy Profile: Closteroviridae. *J. Gen. Virol.* **2020**, *2*, [CrossRef]

17. Marais, A.; Faure, C.; Bergey, B.; Candresse, T. Viral Double-Stranded RNAs (dsRNAs) from Plants: Alternative Nucleic Acid Substrates for High-Throughput Sequencing. In *Methods in Molecular Biology*; Springer Science and Business Media LLC: Clifton, NJ, USA, 2018; Volume 1746, pp. 45–53.

18. Palanga, E.; Filloux, D.; Martin, D.P.; Fernandez, E.; Gargani, D.; Ferdinand, R.; Zabré, J.; Bouda, Z.; Neya, J.B.; Sawadogo, M.; et al. Metagenomic-Based Screening and Molecular Characterization of Cowpea-Infecting Viruses in Burkina Faso. *PLoS ONE* **2016**, *11*, e0165188. [CrossRef]

19. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. *Virus Evol.* **2015**, *1*, vev003. [CrossRef]

20. Scussel, S.; Candresse, T.; Marais, A.; Claverie, S.; Hoareau, M.; Azali, H.A.; Verdin, E.; Tepfer, M.; Filloux, D.; Fernandez, E.; et al. High-throughput sequencing of complete genomes of ipomoviruses associated with an epidemic of cassava brown streak disease in the Comoros Archipelago. *Arch. Virol.* **2019**, *164*, 2193–2196. [CrossRef]

21. Koonin, E.V.; Dolja, V.V.; Morris, T.J. Evolution and Taxonomy of Positive-Strand RNA Viruses: Implications of Comparative Analysis of Amino Acid Sequences. *Crit. Rev. Biochem. Mol. Biol.* **1993**, *28*, 375–430. [CrossRef]

22. Blouin, A.G.; Biccheri, R.; Khalifa, M.E.; Pearson, M.N.; Pollini, C.P.; Hamiaux, C.; Cohen, D.; Ratti, C. Characterization of Actinidia virus 1, a new member of the family Closteroviridae encoding a thaumatin-like protein. *Arch. Virol.* **2017**, *163*, 229–234. [CrossRef]

23. Peremyslov, V.V.; Pan, Y.-W.; Dolja, V.V. Movement Protein of a Closterovirus Is a Type III Integral Transmembrane Protein Localized to the Endoplasmic Reticulum. *J. Virol.* **2004**, *78*, 3704–3709. [CrossRef]

24. Koloniuk, I.; Thekke-Veetil, T.; Reynard, J.-S.; Pleško, I.M.; Přibylová, J.; Brodard, J.; Kellenberger, I.; Sarkisova, T.; Špak, J.; Lamovšek, J.; et al. Molecular Characterization of Divergent Closterovirus Isolates Infecting Ribes Species. *Viruses* **2018**, *10*, 369. [CrossRef]

25. Adiputra, J.; Jarugula, S.; Naidu, R.A. Intra-species recombination among strains of the ampelovirus Grapevine leafroll-associated virus 4. *Virol. J.* **2019**, *16*, 1–13. [CrossRef]

26. Dey, K.K.; Green, J.C.; Melzer, M.; Borth, W.; Hu, J.S. Mealybug Wilt of Pineapple and Associated Viruses. *Horticulturae* **2018**, *4*, 52. [CrossRef]

27. Yu, N.; Luo, Z.; Fan, H.; Zhang, Z.; Li, X.; Wang, J.; Liu, Z.; He, F. Complete genomic sequence of a Pineapple mealybug wilt-associated virus-1 from Hainan Island, China. *Eur. J. Plant. Pathol.* **2015**, *141*, 611–615. [CrossRef]

28. Melzer, M.J.; Sether, D.M.; Karasev, A.V.; Borth, W.; Hu, J.S. Complete nucleotide sequence and genome organization of pineapple mealybug wilt-associated virus-1. *Arch. Virol.* **2008**, *153*, 707–714. [CrossRef] [PubMed]

29. Green, J.C.; Rwahnih, M.A.; Olmedo-Velarde, A.; Melzer, M.J.; Hamim, I.; Borth, W.B.; Brower, T.M.; Wall, M.; Hu, J.S. Further genomic characterization of pineapple mealybug wilt-associated viruses using high-throughput sequencing. *Trop. Plant. Pathol.* **2020**, *45*, 64–72. [CrossRef]

30. Sether, D.M.; Melzer, M.J.; Borth, W.B.; Hu, J.S. Genome organization and phylogenetic relationship of Pineapple mealybug wilt associated virus-3 with family Closteroviridae members. *Virus Genes* **2009**, *38*, 414–420. [CrossRef] [PubMed]
31. Kwon, S.-J.; Jin, M.; Cho, I.-S.; Yoon, J.-Y.; Choi, G.-S. Identification of rehmannia virus 1, a novel putative member of the genus Closterovirus, from *Rehmannia glutinosa*. *Arch. Virol.* **2018**, *163*, 3383–3388. [CrossRef]

32. Martín, S.; Sambade, A.; Rubio, L.; Vives, M.C.; Moya, P.; Guerri, J.; Elena, S.F.; Moreno, P. Contribution of recombination and selection to molecular evolution of Citrus tristeza virus. *J. Gen. Virol.* **2009**, *90*, 1527–1538. [CrossRef]

33. Rubio, L.; Guerri, J.; Moreno, P. Genetic variability and evolutionary dynamics of viruses of the family Closteroviridae. *Front. Microbiol.* **2013**, *4*, 151. [CrossRef]

34. Li, X.H.; Heaton, L.A.; Morris, T.J.; Simon, A.E. Turnip crinkle virus defective interfering RNAs intensify viral symptoms and are generated de novo. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 9173–9177. [CrossRef] [PubMed]

35. Bar-Joseph, M.; Mawassi, M. The defective RNAs of Closteroviridae. *Front. Microbiol.* **2013**, *4*, 132. [CrossRef]

36. White, K.A.; Morris, T.J. RNA determinants of junction site selection in RNA virus recombinants and defective interfering RNAs. *RNA* **1995**, *1*, 1029–1040.

37. Rubio, L.; Yeh, H.-H.; Tian, T.; Falk, B.W. A Heterogeneous Population of Defective RNAs Is Associated with Lettuce infectious yellows virus. *Virology* **2000**, *271*, 205–212. [CrossRef]

38. Eliasco, E.; Livieratos, I.C.; Müller, G.; Guzman, M.; Salazar, L.F.; Coutts, R.H.A. Sequences of defective RNAs associated with potato yellow vein virus. *Arch. Virol.* **2005**, *151*, 201–204. [CrossRef]

39. Menzel, W.; Goetz, R.; Lesemann, D.E.; Vetten, H.J. Molecular characterization of a closterovirus from carrot and its identification as a German isolate of Carrot yellow leaf virus. *Arch. Virol.* **2009**, *154*, 1343–1347. [CrossRef] [PubMed]

40. Mongkolsiriwattana, C.; Chen, A.Y.; Ng, J.C. Replication of Lettuce chlorosis virus (LCV), a crinivirus in the family Closteroviridae, is accompanied by the production of LCV RNA 1-derived novel RNAs. *Virology* **2011**, *420*, 89–97. [CrossRef]

41. Massart, S.; Candresse, T.; Gil, J.; Lacomme, C.; Predajna, L.; Ravnikar, M.; Reynard, J.-S.; Rumbou, A.; Saldarelli, P.; Škorić, D.; et al. A Framework for the Evaluation of Biosecurity, Commercial, Regulatory, and Scientific Impacts of Plant Viruses and Viroids Identified by NGS Technologies. *Front. Microbiol.* **2017**, *8*, 45. [CrossRef] [PubMed]