# Molecular diversity of yam virus Y and identification of banana mild mosaic virus isolates infecting yam (Dioscorea spp.)

Mame Boucar Diouf, Sébastien Guyader, Marie-Michelle Nopoly, Olyvia Gaspard, Denis Filloux, Thierry Candresse, Armelle Marais, Marie Umber,

Pierre Yves Teycheney

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

To cite this version:

Mame Boucar Diouf, Sébastien Guyader, Marie-Michelle Nopoly, Olyvia Gaspard, Denis Filloux, et al. Molecular diversity of yam virus Y and identification of banana mild mosaic virus isolates infecting yam (Dioscorea spp.). Archives of Virology, 2023, 168 (7), pp.180. 10.1007/s00705-023-05809-3. hal-04132918

HAL Id: hal-04132918

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

Submitted on 22 Jun 2023

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Molecular diversity of yam virus Y and identification of banana mild mosaic virus isolates infecting yam (Dioscorea spp.)

Mame Boucar Diouf<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> · Sébastien Guyader<fnref n="1" /> · Marie-Michelle Nopoly<fnref n="1" /> · Olyvia Gaspard<fnref n="1" /> · Denis Filloux<fnref n="4" />,<fnref n="5" /> · Thierry Candresse<fnref n="6" /> · Armelle Marais<fnref n="6" /> · Pierre-Yves Teycheney<fnref n="7" />,<fnref n="8" /> · Marie Umber<fnref n="1" />

Received: 17 February 2023 / Accepted: 2 May 2023
© The Author(s) 2023

## Abstract

Two members of the family *Betaflexiviridae* associated with yam (*Dioscorea* spp.) have been described so far: yam latent virus (YLV) and yam virus Y (YVY). However, their geographical distribution and molecular diversity remain poorly documented. Using a nested RT-PCR assay, we detected YVY in *D. alata*, *D. bulbifera*, *D. cayenensis*, *D. rotundata*, and *D. trifida* in Guadeloupe, and in *D. rotundata* in Côte d'Ivoire, thus extending the known host range of this virus and geographical distribution. Using amplicon sequencing, we determined that the molecular diversity of YVY in the yam samples analyzed in this work ranged between 0.0 and 29.1% and that this diversity is partially geographically structured. We also identified three isolates of banana mild mosaic virus (BanMMV) infecting *D. alata* in Guadeloupe, providing the first evidence for BanMMV infection in yam.

Yam (*Dioscorea* spp.) is widely cultivated in the intertropical belt for its nutritional and pharmaceutical properties [1–3] and economic value [4]. Twenty-five viruses infecting yam are currently recognized by the International Committee on Taxonomy of Viruses (ICTV) [5], including one member of family *Betaflexiviridae*, yam latent virus (YLV) [6], which is assigned to genus *Carlavirus*. An additional

virus, yam virus Y (YVY) [7], has all the hallmarks of viruses of a member of the family *Betaflexiviridae*, but has not been formally recognized as a member of this family. YLV has been reported only in *Dioscorea opposita* in China, whereas YVY has been reported in *D. rotundata* in germplasm collections in Nigeria and Ghana, where it was found mostly in mixed infections with yam mosaic virus (YMV), a potyvirus [7]. The symptomatology associated with YLV or YVY infection has not been elucidated. Overall, the molecular diversity, host range, geographical distribution, impact on production, and transmission mode of YLV and YVY remain largely unknown.

In order to fill these knowledge gaps, we undertook a comprehensive prevalence survey of yam-infecting *Betaflexiviridae* members using three sampling sources (Supplementary Table S1). The first source, referred to as ‘Guadeloupe’, included 896 leaf samples of the species *D. alata*, *D. bulbifera*, *D. cayenensis*, *D. esculenta*, *D. rotundata*, and *D. trifida* collected between November 2019 and December 2020 from 18 yam fields scattered across Guadeloupe, as described by Diouf *et al.* [8]. Considering that the relationship between leaf symptoms and viral infection in yam is poorly documented, symptoms, or lack thereof, were not taken into consideration during sample collection. A second source, referred to as ‘BRC-TP’, included 406 leaf samples from the same species, collected between 2014 and 2019 from the *in vitro* yam collection maintained by the Biological Resources

Handling Editor: Ioannis E. Tzanetakis.

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

Marie Umber
marie.umber@inrae.fr

<footnote n="1">1 INRAE, UR ASTRO, 97170 Petit-Bourg, Guadeloupe, France</footnote>

<footnote n="2">2 CIRAD, UMR AGAP Institut, 97130 Capesterre Belle Eau, France</footnote>

<footnote n="3">3 UMR AGAP Institut, Univ Montpellier, CIRAD, INRAE, Institut Agro, 97130 Capesterre Belle Eau, France</footnote>

<footnote n="4">4 CIRAD, UMR PHIM, 34090 Montpellier, France</footnote>

<footnote n="5">5 PHIM Plant Health Institute, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, 34090 Montpellier, France</footnote>

<footnote n="6">6 INRAE, Univ. Bordeaux, UMR BFP, CS20032, 33882 Villenave d'Ornon Cedex, France</footnote>

<footnote n="7">7 CIRAD, UMR PVBMT, 97410 Saint-Pierre, La Réunion, France</footnote>

<footnote n="8">8 UMR PVBMT, Université de la Réunion, 97410 Saint-Pierre, La Réunion, France</footnote>
Center for Tropical Plants (BRC-TP) in Guadeloupe [9]. A third and last source, referred to as 'CNRA', included 38 leaf samples of *D. rotundata* collected in 2018 from the yam germplasm collection maintained under open field conditions by the Food Crops Research Station of the Centre National de Recherche Agronomique (CNRA) in Bouaké (Côte d'Ivoire) [10]. Total nucleic acids (TNAs) were purified from each leaf sample using extraction procedure 2, described by Foissac *et al.* [11], and were used as template for cDNA synthesis as described by Umber *et al.* [9]. This was followed by a nested PCR using the five inosine-containing degenerate primers (PDO) designed by Foissac *et al.* [11] (Supplementary Table S2). These primers allow the detection of viruses of the genera *Bamivirus*, *Capillovirus*, *Foveavirus*, and *Trichovirus* (family *Betaflexiviridae*), resulting in the amplification of a 362-bp fragment containing the conserved motifs II, III, IV, and V of the RNA-directed RNA polymerases (RdRps) of positive-stranded RNA viruses [12]. The results of the prevalence survey showed that 24.6% (220/896) of the samples from the 'Guadeloupe' sampling tested positive (Table 1). Interestingly, the highest prevalence in this sampling was in plants of the African species *D. cayenensis* and *D. rotundata* (96.2% and 67.8%, respectively). The prevalence was similarly high (78.9%; 30/38) in the 'CNRA' sampling, which consisted exclusively of samples from the species *D. rotundata*, but it was significantly lower (14.3%; 58/406) in the 'BRC-TP' sampling.

The molecular diversity of yam-infecting members of the family *Betaflexiviridae* was also investigated. For this, 31 amplification products (18, 11, and 2 from the 'Guadeloupe', 'BRC-TP', and 'CNRA' samplings, respectively; Supplementary Table S1) were cloned in pGEM-T Easy Vector (Promega, Charbonnières, France) according to the manufacturer's protocol. One to six clones per amplification product were randomly chosen and sequenced (Genewiz, Leipzig, Germany), resulting in a total of 80 nucleotide sequences of 310 nt length after primer sequences were removed: 50, 22 and 8 sequences originated from the 'Guadeloupe', 'BRC-TP', and 'CNRA' samplings, respectively (GenBank accession numbers: OP495428-OP495507). This dataset was completed with nucleotide sequences of reference isolates for each of the 15 recognized genera of the family *Betaflexiviridae*, banana virus X (BVX-Som; AY710267; unassigned *Betaflexiviridae*) and YLV-SG1 (KJ789130), sequences available for YVY isolates YVY-Mak and YVY-Dan (MK782910 and MK782911), identified in Nigeria by high-throughput sequencing (HTS) [7], and for dioscorea virus A (DioVA-PE5.1; LC467961), a putative member of family *Betaflexiviridae* reported in yam in Brazil (Hayashi *et al.*, unpublished). In addition, sequences from a YVY isolate (IgC10YVYDr) and a BanMMV isolate (IgB10BanMMVDa) identified in a *D. rotundata* plant from Benin and in a *D. alata* plant collected in Guadeloupe, respectively, were

<table><caption>Table 1 Prevalence of *Betaflexiviridae* members in yam samples from the three samplings analyzed in this work</caption><thead><tr><th rowspan="2">Sampling<br/>Yam species</th><th colspan="3">Guadeloupe</th><th colspan="3">BRC-TP</th><th colspan="3">CNRA</th></tr><tr><th>Number of<br/>indexed sam-<br/>ples</th><th>PDO primers</th><th>YVY-specific<br/>primers</th><th>Ban-<br/>MMV-<br/>specific<br/>primers</th><th>Number of<br/>indexed sam-<br/>ples</th><th>PDO prim-<br/>ers</th><th>YVY-<br/>specific<br/>primers</th><th>Ban-<br/>MMV-<br/>specific<br/>primers</th><th>Number of<br/>indexed sam-<br/>ples</th><th>PDO prim-<br/>ers</th><th>YVY-<br/>specific<br/>primers</th><th>Ban-<br/>MMV-<br/>specific<br/>primers</th></tr><tr><th></th><th>N</th><th>%</th><th>N</th><th>%</th><th>N</th><th>%</th><th>N</th><th>%</th><th>N</th><th>%</th><th>N</th><th>%</th><th>N</th><th>%</th></tr></thead><tbody><tr><td><i>D. alata</i></td><td>720</td><td>135</td><td>18.8</td><td>129</td><td>17.9</td><td>3</td><td>0.4</td><td>141</td><td>36</td><td>25.5</td><td>35</td><td>24.8</td><td>1</td><td>0.7</td><td>-</td></tr><tr><td><i>D. bulbifera</i></td><td>16</td><td>1</td><td>6.3</td><td>1</td><td>6.3</td><td>0</td><td>0.0</td><td>4</td><td>1</td><td>25.0</td><td>1</td><td>25.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td><i>D. cayenensis</i></td><td>26</td><td>25</td><td>96.2</td><td>25</td><td>96.2</td><td>0</td><td>0.0</td><td>49</td><td>10</td><td>20.4</td><td>7</td><td>14.3</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td><i>D. esculenta</i></td><td>31</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>8</td><td>1</td><td>12.5</td><td>1</td><td>12.5</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td><i>D. rotundata</i></td><td>87</td><td>59</td><td>67.8</td><td>59</td><td>67.8</td><td>0</td><td>0.0</td><td>27</td><td>6</td><td>22.2</td><td>5</td><td>18.5</td><td>0</td><td>0.0</td><td>38</td></tr><tr><td><i>D. trifida</i></td><td>16</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>173</td><td>4</td><td>2.3</td><td>4</td><td>2.3</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Other species</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>4</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Total</td><td>896</td><td>220</td><td>24.6</td><td>214</td><td>23.9</td><td>3</td><td>0.3</td><td>406</td><td>58</td><td>14.3</td><td>53</td><td>13.1</td><td>1</td><td>0.2</td><td>38</td></tr></tbody></table>

N, number of positive samples; %, percentage of positive samples relative to the number of indexed samples
also included. These sequences were obtained from contigs assembled *de novo* from yam RNASeq transcriptomic Illumina reads generated by the ARCAD project and referenced under numbers SRX1967863 and SRX1967874, respectively, in the SRA GenBank database [13]. A nucleotide sequence alignment was performed using MAFFT v7.490 [14] and checked visually to ensure that none of the introduced gaps would disrupt conceptually translated protein sequences. A maximum-likelihood phylogenetic tree constructed using IQ-TREE v2.2.0 [15], with the best nucleotide substitution model (TIM3+F+I+G4) chosen by the program based on the Bayesian information criterion, and branch support was estimated with 10,000 ultrafast bootstrap replicates [16] and 1,000 SH-aLRT replicates. The sequence identity matrix was obtained using Clustal Omega v1.2.4 [17].

The phylogenetic analysis confirmed that all 80 sequences determined in this work belonged to viruses of the family *Betaflexiviridae* (Fig. 1). Most of the sequences (77/80) formed a clade distinct from other genera in the family *Betaflexiviridae*. This clade also included YVY-Dan, YVY-Mak, and DioVA-PE5.1 and had the CRB603YVYDc sequence in a basal position. This sequence appeared particularly divergent relative to the other sequences in this clade, sharing with them, on average, only 74.8% nucleotide sequence identity (range, 71.9-77.7%) and 82.3% amino acid sequence identity (79.6-83.5%), while the other sequences in the clade shared, on average, 87.7% nucleotide sequence identity (range 79.0-100.0%) and 96.0% amino sequence identity (89.3-100%). Further studies are therefore required to decide on the taxonomic status of the CRB603YVYDc isolate. Analysis of the IgC10YVYDr sequence, which corresponds to a nearly complete viral genome, showed that it shared 96.0% and 94.9% amino acid sequence identity with both YVY-Mak and YVY-Dan isolates in the analyzed complete RdRp and coat protein (CP) region, respectively. Considering that these figures are well above the species demarcation threshold established by the ICTV for the family *Betaflexiviridae* (80% amino acid sequence identity in the RdRp or CP; [18]), IgC10YVYDr has all the hallmarks of a YVY isolate. The other 76 sequences in this clade, which were determined in this study and cover part of the RdRp coding region, shared, on average, 85.6% (range, 82.9-87.7%) and 86.9% (range, 84.2-96.1%) nucleotide sequence identity with those of YVY-Mak and YVY-Dan isolates, respectively (96.9% [94.2-99.0%] and 97.3% [94.2-100%] amino acid sequence identity), and 84.3% (range, 81.3-86.8%) nucleotide sequence identity with that of IgC10YVYDr (92.7% [89.3-95.1%] amino acid identity). Since these identity scores are all in the same range, the 76 sequences generated in this study should most likely also be considered isolates of YVY. The same can be concluded for isolate DioVA-PE5.1, which shared 85.2% and 87.4% nucleotide sequence identity with YVY-Mak and YVY-Dan, respectively (97.1%

amino acid sequence identity in the encoded protein for both isolates). The remaining three sequences out of the 80 determined in this study (Kin8BanMMVDa-1, Kin8Ban-MMVDa-2, and Kin2BanMMVDa), which were obtained from *D. alata* plants of the 'Guadeloupe' sampling, shared 79.4-80.0% nucleotide sequence identity and 87.4% amino acid sequence identity with the reference sequence of banana mild mosaic virus (BanMMV; AF314662) and 98.4-99.0% nucleotide sequence identity and 100% amino acid sequence identity with IgB10BanMMVDa. Sequence comparisons of the RdRp and CP proteins from the nearly complete genome sequence of IgB10BanMMVDa and the BanMMV reference isolate showed that IgB10BanMMVDa unambiguously corresponds to a BanMMV isolate (86.3% and 80.2% amino acid sequence identity to AF314662 in the complete RdRp and CP, respectively). These figures strongly support the notion that the three partial sequences determined in this study belong to BanMMV isolates and provide evidence that BanMMV infects *D. alata* in Guadeloupe.

The YVY isolates characterized in this work displayed moderately high rates of nucleotide sequence diversity, with an average sequence dissimilarity of 12.9% (0-22.3%), approaching those previously reported for other members of the family *Betaflexiviridae* [11, 19, 20], but with a rather low average amino acid sequence dissimilarity of 3.9% (0-10.7%). Statistical analysis was performed on the 80 YVY sequences used in the phylogenetic analysis described above, including the 76 YVY sequences generated in this work but excluding the CRB603YVYDc sequence, to investigate a possible geographical and/or historical structuring of the molecular diversity of YVY. Four groups of sequences were defined: 'Africa' (eight sequences from the 'CNRA' sampling, one from Benin, and two from Nigeria), 'BRC-TP' (21 sequences), 'Guadeloupe' (47 sequences), and 'S. America' (one sequence from Brazil). An analysis of molecular variance (AMOVA) was carried out on single nucleotide polymorphisms (SNPs), using R v4.2.1 [21] in order to assess sequence differentiation among these four groups. This analysis revealed a between-group variance of 7.97 (significantly greater than zero; *P* = 0.001), showing that the groups of isolates have diverged from each other. This result was confirmed by a discriminant analysis of principal components (DAPC; Fig. 2), which looks for SNP patterns that maximize between-group distance while minimizing within-group distance. The first discriminant axis appeared to separate the 'Africa' group from a metagroup composed of sequences from sources located in the New World ('BRC-TP', 'Guadeloupe', and 'S. America' groups), and the second axis separated the latter three groups, thereby showing that isolates in these groups form genetically differentiated pools. Sequence comparisons revealed an average between-group nucleotide dissimilarity of 14.1-15.2%, while the average within-group dissimilarity was lower (10.2% and 10.5%
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>
►Fig. 1 Maximum-likelihood phylogenetic tree showing the relationships between betaflexivirid sequences obtained in this study (in bold) and those of previously characterized isolates. Phylogenetic analysis was performed based on the amplified sequences (310 nt) corresponding to nucleotide positions 4612–4921 of the genome of YVY-Mak and coding for part of the RdRp. Potato virus X (PVX) was used to root the tree. Keys are provided for the coloured sequence names according to their taxonomic position in the family *Betaflexiviridae* (subfamilies recognized by the ICTV), and to the coloured dots referring for the location of the samples from which the sequences originated. Significant branch support (ultrafast bootstrap values above 95% of 10,000 replicates and SH-aLRT above 85% of 1,000 replicates) is indicated by an orange rectangle across the supported branches. Da. *Dioscorea alata*; Db, *D. bulbifera*; Dc, *D. cayenensis*; Dr, *D. rotundata*; Dt, *D. trifida*

for ‘Africa’ and ‘Guadeloupe’, respectively), except for the
‘BRC-TP’ group, whose average within-group nucleotide
dissimilarity (14.1%) was in the same range as the average
between-group dissimilarity.

The prevalence of YVY and BanMMV in the three
samplings was investigated. For this, two species-specific
primer pairs were designed based on sequence comparisons
between all available sequences of YVY and BanMMV yam
isolates, in order to capture the molecular diversity of both
viruses, and used for virus indexing. Primers YVY-F and
YVY-R (Supplementary Table S2) were used in nested RT-
PCR experiments, using the PCR products generated from
the first PCR step of the PDO-PCR assay as a template.
This YVY-specific RT-PCR assay produced amplicons of
275 bp corresponding to nucleotide positions 4640-4914 in
the YVY-Mak genome. Conditions for the nested PCR were
an initial denaturation step of 2 min at 95° C, followed by
30 cycles of 30 s at 95°C, 30 s at 54°C, and 30 s at 72°C,
and a final elongation of 5 min at 72°C. The primers Yam-
BanF and YamBanR (Supplementary Table S2) were used
in RT-PCR experiments, producing amplicons of 226 bp
corresponding to nucleotide positions 4552-4777 of the Ban-
MMV reference isolate (AF314662). PCR conditions were
the same as above except that the annealing temperature
was increased to 57°C. Indexing results showed that most
of the samples used in this study that were positive with
broad-range PDO primers also tested positive when using
YVY-specific primers (Table 1; Supplementary Table S1).
This proportion was 97.3% for samples of the ‘Guadeloupe’
sampling, 91.4% for those of the ‘BRC-TP’ sampling, and
96.7% for those of the ‘CNRA’ sampling (Table 1), show-
ing the efficiency of this specific detection tool. In contrast,
the prevalence of BanMMV was very low, with only 0.3%
and 0.2% of the samples from the ‘Guadeloupe’ and ‘BRC-
PT’ samplings, respectively, testing positive. All 38 samples
from the ‘CNRA’ sampling were negative for BanMMV,
strengthening the observation that only *D. alata* seems to
be prone to infection by this virus.

Our results extend the known natural host range of
YVY to several additional yam species, including *D.
alata*, *D. bulbifera*, *D. cayenensis*, and *D. trifida* in addi-
tion to *D. rotundata*, which had been reported previously
[7], and its known geographical distribution to Benin,
Côte d'Ivoire and Guadeloupe, where it is widespread.
Whether YVY was introduced into Guadeloupe through
infected yam imported from Africa during the slave trade
[22], as hypothesized for YMV and yam mild mosaic virus
(YMMV) [23, 24], remains to be investigated and requires
extensive phylogeographical studies with large samplings
in Africa. The ‘BRC-TP’ collection includes yam acces-
sions originating from various parts of the world (the
Caribbean, South America, Africa, and the South Pacific;
Supplementary Table S1) that were introduced before
accurate viral diagnostic tools were available, making it
impossible to know whether accessions in this collection
were infected by YVY prior to their introduction into Gua-
deloupe or became contaminated while being maintained
under field conditions in Guadeloupe prior to their intro-
duction into the ‘BRC-TP’ *in vitro* collection.

The YVY isolates characterized in this work display levels of molecular diversity similar to those reported for other yam-infecting viruses, such as YMV, YMMV, and Dioscorea mosaic associated virus (DMaV) [23, 25, 26]. However, diversity among YVY isolates could be partitioned into four distinct genetic pools by the DAPC analysis (Fig. 2), reflecting a genetic differentiation also seen in the phylogenetic tree in Figure 1 in the form of origin-specific subclades of isolates. The molecular diversity among infected plants was moderate, whereas within-plant molecular diversity was low, contrary to the situation encountered for other yam-infecting viruses such as DMaV [26]. Some YVY sequences with 100% identity were amplified from distinct samples collected from distinct species (Fig. 1), suggesting plant-to-plant transmission of this virus. Whether such transmission occurs mechanically through the use of infected tools for fragmenting yam tubers before planting or involves biological vectors remains to be elucidated. YVY might also be transmitted vertically through vegetative propagation and the use of infected tubers as planting material, as demonstrated for other yam viruses [8]. Since most YVY-infected plants reported in this work were coinfect ed with other viruses such as YMMV and YaV1 [8, 10], we could not associate YVY with specific symptoms.

The use of primers YVY-F and YVY-R allowed the
detection of YVY in 94.2% (290/308) of the samples that
tested positive using degenerate PDO primers (Supple-
mentary Table S1), as well as in six samples that tested
negative with the PDO primers. On the other hand, 5.8%
(18/308) of the samples tested positive when using the
PDO primers and negative when using the YVY-specific
primers. Four of the latter samples were found to be
<figure><img src="image_6.png" /><figcaption>**Fig. 2** Ordination plot of the discriminant analysis of principal components performed on 76 YVY sequences obtained in this work and four YVY sequences obtained from GenBank. Sequences were assigned to four groups ('Africa', 'BRC-TP', 'Guadeloupe', and 'S. America') according to their sampling location and are coloured</figcaption></figure>

accordingly on the plot. Ellipses delineate the confidence interval around the group centroids. A between-group percent sequence similarity matrix is shown in the bottom right corner, with the standard deviation in parentheses and within-group dissimilarities shown along the matrix diagonal. NA, not applicable

infected by BanMMV, suggesting that the remaining 14 could have been infected by another virus of the family *Betaflexiviridae*, including potentially uncharacterized YVY variants that escaped detection with the YVY-specific primers described in this work.

Here, we report the first occurrence of BanMMV isolates infecting yam. BanMMV was detected in three *D. alata* samples collected from the same plot, in one *D. alata* hybrid collected in Guadeloupe and used for RNASeq analysis and in one *D. alata* accession from the 'BRC-TP' sampling, originating in Vanuatu (CRB767; Supplementary Table S1). The latter was introduced into Guadeloupe in 2003 and maintained under field conditions until it was introduced *in vitro* in 2014. Therefore, it may have been infected by BanMMV either before or after its introduction into Guadeloupe. There is experimental evidence of plant-to-plant transmission of BanMMV in banana in Guadeloupe [19], although the actual mode of transmission (mechanical and/or vector-borne) of this virus has not been elucidated. A possible scenario might involve vector-borne transmission of BanMMV from banana to yam, considering that BanMMV is widespread in banana in Guadeloupe and that both crops are often present together on

farm plots, a situation that could favour host jumps from banana to yam.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s00705-023-05809-3.

**Acknowledgements** The authors wish to thank Chantal Faure (UMR BFP, INRAE and University of Bordeaux, Villenave d'Ornon, France) for help with nucleic acid extraction, Rose-Marie Gomez (UR ASTRO, INRAE, Guadeloupe, France) for technical help, and Yacouba Bakayoko and Amani Michel Kouakou for providing samples from the CNRA germplasm collection.

**Funding** This work was funded by the Agence Nationale de la Recherche (ANR) and the Guadeloupe Region through the NetBiome program under project SafePGR (ANR-11-EBIM-0002-06) and by the European Commission and the Guadeloupe Region through the Malin project (2018-FED-1088) under the EU's Seventh Framework Programme.

**Data availability statement** The data presented in this study are available in Supplementary Tables S1–S2. The full dataset is openly available at Recherche Data Gouv, V1, https://doi.org/10.57745/GNUXVA, UNF:6:uRiYq5/k7MnOABwsIvhK9A==. The nucleotide sequences reported in this work have been deposited in the GenBank database under accession numbers OP495428—OP495507.
## Declarations

**Conflict of interest** The authors declare that they have no conflict of interest.

**Ethical approval** This article does not contain any studies with human participants or animals performed by any of the authors.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Ayensu ES, Coursey DG (1972) Guinea Yams: the botany, ethno-botany, use and possible future of yams in West Africa. Econ Bot 26:301–318

2. Padhan B, Panda D (2020) Potential of neglected and underutilized yams (*Dioscorea* spp.) for improving nutritional security and health benefits. Front Pharmacol 11:496. https://doi.org/10.3389/fphar.2020.00496

3. Lebot V, Lawac F, Legendre L (2023) The greater yam (Dioscorea alata L.): a review of its phytochemical content and potential for processed products and biofortification. J Food Compos Anal 115:104987. https://doi.org/10.1016/j.jfca.2022.104987

4. Amusa TA, Isiwu EC, Oketoobo EA (2018) Economics of yam (*Dioscorea rotundata*) production among small holder farmers in Abia state, Nigeria. https://aksujaeerd.com/viewpdf/articles/publications/d/5530.pdf. Accessed 3 Nov 2022

5. Diouf MB, Festus R, Silva G, Guyader S, Umber M, Seal S, Teycheney P-Y (2022) Viruses of yams (Dioscorea spp.): current gaps in knowledge and future research directions to improve disease management. Viruses 14:1884. https://doi.org/10.3390/v14091884

6. Lebas BSM, Ochoa-Corona FM, Elliott DR, Tang Z, Alexander JR (2005) Partial characterization of a Carla-like virus infecting yam (*Dioscorea* spp.) from China. Plant Dis 89:912–912. https://doi.org/10.1094/PD-89-0912A

7. Silva G, Bömer M, Rathnayake AI, Sewe SO, Visendi P, Oyekanmi JO, Quain MD, Akomeah B, Kumar PL, Seal S (2019) Molecular characterization of a new virus species identified in yam (Dioscorea spp.) by high-throughput sequencing. Plants 8:167. https://doi.org/10.3390/plants8060167

8. Diouf MB, Guyader S, Gaspard O, Francius E, Teycheney P-Y, Umber M (2022) Epidemiology of yam viruses in Guadeloupe: role of cropping practices and seed-tuber supply. Viruses 14:2366. https://doi.org/10.3390/v14112366

9. Umber M, Filloux D, Gélabale S, Gomez RM, Marais A, Gallet S, Gamiette F, Pavis C, Teycheney P-Y (2020) Molecular viral diagnosis and sanitation of yam genetic resources: implications for safe yam germplasm exchange. Viruses 12:1101. https://doi.org/10.3390/v12101101

10. Bakayoko Y, Kouakou AM, Kouassi AB, Gomez RM, Dibi KEB, Essis BS, N'Zué B, Adebola P, N'Guetta AS-P, Umber

M (2021) Detection and diversity of viruses infecting African yam (*Dioscorea rotundata*) in a collection and F1 progenies in Côte d'Ivoire shed light to plant-to-plant viral transmission. Plant Pathol 70:1486–1495. https://doi.org/10.1111/ppa.13393

11. Foissac X, Svanella-Dumas L, Gentit P, Dulucq M-J, Marais A, Candresse T (2005) Polyvalent degenerate oligonucleotides reverse transcription-polymerase chain reaction: a polyvalent detection and characterization tool for trichoviruses, capilloviruses, and foveaviruses. Phytopathology 95:617–625. https://doi.org/10.1094/PHYTO-95-0617

12. Koonin EVY (1991) The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J Gen Virol 72:2197–2206. https://doi.org/10.1099/0022-1317-72-9-2197

13. Sarah G, Homa F, Pointet S, Contreras S, Sabot F, Nabholz B, Santoni S, Sauné L, Ardisson M, Chantret N et al (2017) A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives. Mol Ecol Resour 17:565–580. https://doi.org/10.1111/1755-0998.12587

14. Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol 30:772–780. https://doi.org/10.1093/molbev/mst010

15. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R (2020) IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol 37:1530–1534. https://doi.org/10.1093/molbev/msaa015

16. Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS (2018) UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol 35:518–522. https://doi.org/10.1093/molbev/msx281

17. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins D (2011) Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7:539. https://doi.org/10.1038/msb.2011.75

18. *Betaflexiviridae*—Positive sense RNA viruses—positive sense RNA viruses (2011) - ICTV. https://talk.ictvonline.org/ictv-reports/ictv_9th_report/positive-sense-rna-viruses-2011/w/posrna_viruses/241/betaflexiviridae. Accessed 26 May 2021

19. Teycheney P-Y, Laboureau N, Iskra-Caruana M-L, Candresse T (2005) High genetic variability and evidence for plant-to-plant transfer of Banana mild mosaic virus. J Gen Virol 86:3179–3187. https://doi.org/10.1099/vir.0.81197-0

20. Marais A, Svanella-Dumas L, Barone M, Gentit P, Faure C, Charlot G, Ragozzino A, Candresse T (2012) Development of a polyvalent RT-PCR detection assay covering the genetic diversity of Cherry capillovirus A. Plant Pathol 61:195–204. https://doi.org/10.1111/j.1365-3059.2011.02488.x

21. R Core Team (2022) R: a language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria

22. Sharif BM, Burgarella C, Cormier F, Mournet P, Causse S, Van Nguyen K, Kaoh J, Rajaonah MT, Lakshan SR, Waki J et al (2020) Genome-wide genotyping elucidates the geographical diversification and dispersal of the polyploid and clonally propagated yam (*Dioscorea alata*). Ann Bot 126:1029–1038. https://doi.org/10.1093/aob/mcaa122

23. Bousalem M, Dallot S, Fuji S, Natsuaki KT (2003) Origin, worldwide dispersion, bio-geographical diversification, radiation and recombination: an evolutionary history of Yam mild mosaic virus (YMMV). Infect Genet Evol 3:189–206. https://doi.org/10.1016/S1567-1348(03)00085-6

24. Mendoza AR, Margaria P, Nagata T, Winter S, Blawid R (2022) Characterization of yam mosaic viruses from Brazil reveals a new phylogenetic group and possible incursion from the African
continent. Virus Genes 58:294–307. https://doi.org/10.1007/s11262-022-01903-x

25. Bousalem M, Douzery EJP, Fargette D (2000) High genetic diversity, distant phylogenetic relationships and intraspecies recombination events among natural populations of Yam mosaic virus: a contribution to understanding potyvirus evolution. J Gen Virol 81:243–255. https://doi.org/10.1099/0022-1317-81-1-243

26. Umber M, Filloux D, Svanella-Dumas L, Bonheur L, Acina-Mambole I, Gomez R-M, Faure C, Anzala F, Pavis C, Roumagnac P, Marais A, Theil S, Candresse T, Teycheney P-Y

(2022) Host range and molecular variability of the sadwavirus dioscorea mosaic associated virus. Arch Virol 167:917–922. https://doi.org/10.1007/s00705-022-05379-w

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.