# Comparison of the Virome of Quarantined Sugarcane Varieties and the Virome of Grasses Growing near the Quarantine Station

Jean H Daugrois, Denis Filloux, Charlotte Julian, Lisa Claude, Romain

Ferdinand, Emmanuel Fernandez, Hugo Fontes, Philippe C Rott, Philippe

Roumagnac

► To cite this version:

Jean H Daugrois, Denis Filloux, Charlotte Julian, Lisa Claude, Romain Ferdinand, et al.. Comparison of the Virome of Quarantined Sugarcane Varieties and the Virome of Grasses Growing near the Quarantine Station. Viruses, 2021, 13 (5), pp.922. 10.3390/v13050922. hal-03247843

HAL Id: hal-03247843

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

Submitted on 3 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

# Comparison of the Virome of Quarantined Sugarcane Varieties and the Virome of Grasses Growing near the Quarantine Station

Jean H. Daugrois <fnref n="1" />,<fnref n="2" />, Denis Filloux <fnref n="1" />,<fnref n="2" />, Charlotte Julian <fnref n="1" />,<fnref n="2" />, Lisa Claude <fnref n="1" />,<fnref n="2" />, Romain Ferdinand <fnref n="1" />,<fnref n="2" />, Emmanuel Fernandez <fnref n="1" />,<fnref n="2" />, Hugo Fontes <fnref n="3" />,<fnref n="4" /> ID, Philippe C. Rott <fnref n="1" />,<fnref n="2" /> ID and Philippe Roumagnac <fnref n="1" />,<fnref n="2" />,* ID

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

check for
updates

**Citation:** Daugrois, J.H.; Filloux, D.; Julian, C.; Claude, L.; Ferdinand, R.; Fernandez, E.; Fontes, H.; Rott, P.C.; Roumagnac, P. Comparison of the Virome of Quarantined Sugarcane Varieties and the Virome of Grasses Growing near the Quarantine Station. *Viruses* **2021**, *13*, 922. https://doi.org/10.3390/v13050922

Academic Editor: Henryk Czosnek

Received: 13 April 2021

Accepted: 10 May 2021

Published: 16 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:** Visacane is a sugarcane quarantine station located in the South of France, far away from sugarcane growing areas. Visacane imports up to 100 sugarcane varieties per year, using safe control and confinement measures of plants and their wastes to prevent any risk of pathogen spread outside of the facilities. Viruses hosted by the imported material are either known or unknown to cause disease in cultivated sugarcane. Poaceae viruses occurring in plants surrounding the quarantine glasshouse are currently unknown. These viruses could be considered as a source of new sugarcane infections and potentially cause new sugarcane diseases in cases of confinement barrier failure. The aim of this study was to compare the plant virome inside and outside of the quarantine station to identify potential confinement failures and risks of cross infections. Leaves from quarantined sugarcane varieties and from wild Poaceae growing near the quarantine were collected and processed by a metagenomics approach based on virion-associated nucleic acids extraction and library preparation for Illumina sequencing. While viruses belonging to the same virus genus or family were identified in the sugarcane quarantine and its surroundings, no virus species was detected in both environments. Based on the data obtained in this study, no virus movement between quarantined sugarcane and nearby grassland has occurred so far, and the confinement procedures of Visacane appear to be properly implemented.

**Keywords:** viral metagenomics; plant quarantine; sugarcane

## 1. Introduction

Sugarcane (*Saccharum* interspecific hybrids) belongs to the Poaceae family and is cultivated in tropical and subtropical areas for the production of sugar and ethanol. To be competitive on the global market and for diversification of products, sugarcane farmers need to use the best-performing varieties and breeders need genetic resources to create them. Consequently, moving plants across international borders for exchange between breeders or directly to producers is an essential step of this process. Plant movement involves biosafety to protect local production by avoiding pathogen introduction and spread. Safe and accurate quarantine procedures need to be implemented, including accurate disease testing [1].

To control disease introduction through sugarcane germplasm movement in the outermost French European regions and for partners in Africa, a sugarcane quarantine station was established in 1971 by the Institut de Recherche en Agronomie Tropicale in

<footnote n="1">¹ CIRAD, UMR PHIM, 34090 Montpellier, France; jean-heinrich.daugrois@cirad.fr (J.H.D.); denis.filloux@cirad.fr (D.F.); charlotte.julian@cirad.fr (C.J.); claude.lisa83@hotmail.fr (L.C.); romain.ferdinand@cirad.fr (R.F.); emmanuel.fernandez@cirad.fr (E.F.); philippe.rott@cirad.fr (P.C.R.)</footnote>

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

<footnote n="3">³ Tour du Valat, Research Institute for the Conservation of Mediterranean Wetlands, 13200 Arles, France; fontes@tourduvalat.org</footnote>

<footnote n="4">⁴ Institut Méditerranéen de Biodiversité et Ecologie, UMR CNRS-IRD, Avignon Université, Aix-Marseille Université, IUT d'Avignon, 337 chemin des Meinajariés, Site Agroparc BP 61207, 84911 Avignon, France</footnote>

<footnote>* Correspondence: philippe.roumagnac@cirad.fr</footnote>
Nogent sur Marne, France [2]. This quarantine station moved to Montpellier in 1978 and became Cirad's sugarcane quarantine in 1984. In 2010, the sugarcane quarantine service was named Visacane (http://visacane.cirad.fr/, accessed on 31 March 2021). Nowadays, Visacane works with research organizations and sugarcane companies from tropical and subtropical locations around the world [3]. It is a go-through sugarcane quarantine located in the South of France, far away from sugarcane-growing areas. The major objective of Visacane is to supply sugarcane germplasm free from potentially damaging pathogens [3]. Viruses are the most detected pathogens in sugarcane quarantine [3,4]. Imported planting material can host viruses [4–6] that are either known for damaging sugarcane crops (the overall majority of sugarcane mastreviruses, poaceviruses, poleroviruses, and potyviruses) or whose effect on sugarcane growth is unknown (ampeloviruses, badnaviruses, and at least one umbravirus and one mastrevirus) [4,6–9]. On the other hand, plant viruses occurring in Poaceae growing around a quarantine glasshouse are usually unknown. If virus-infected, these plants are potential sources of new infections in quarantine in case of passing of viruses through the confinement barriers.

Because Visacane is located outside a sugarcane growing area, the risk of sugarcane infection by surrounding local viruses should be nonexistent or minimal. Nevertheless, the purpose of quarantining is also to avoid the release or escape from quarantine of viruses that could potentially threaten the local vegetation. Infrastructures and processes are rarely 100% safe and a flow may exist between a quarantine glasshouse and its surrounding environment. Importation of infected plants in a new area is a risk for disease emergence and newly introduced viruses may invade communities of native plants [10]. Additionally, virus species infecting several plant species of a same plant family is a common feature, especially in the Poaceae family. For example, maize, sorghum, sugarcane, and several other grasses can host sugarcane mosaic virus [11]. Maize yellow mosaic virus, a polerovirus initially described in maize, was also recently found in sugarcane in Nigeria and China [12,13]. Sorghum mosaic virus, a potyvirus from sorghum, is the causal agent of sugarcane mosaic in the USA [14] and in China [15]. Grassland or wild plants can also serve as reservoirs for known pathogenic viruses such as wheat streak mosaic virus in the Czech Republic [16] or for new emerging plant viruses as illustrated by the epidemic of maize streak virus in Africa [17]. Wild Poaceae from different grasslands host viruses from the genera *Mastrevirus* (*Geminiviridae* family), *Tritimovirus* (*Potyviridae* family) and *Luteovirus* (*Luteoviridae* family) [18–20]. Among these, viruses from the genera *Polerovirus* (*Luteoviridae*), *Potyvirus* (*Potyviridae*) and *Mastrevirus* (*Geminiviridae*) cause diseases in sugarcane [8].

Consequently, to investigate the potential flow of plant viruses between a quarantine and its environment, deciphering the plant virome of both entities is necessary. At the present time, the best methodology for virome identification in a specific environment is viral metagenomics, and different approaches are available [21,22]. Virion-associated nucleic acid (VANA) metagenomics has already been used to identify known and new viruses in sugarcane [5,7,23] and in wild Poaceae [24–26]. The objective of this study was to use this methodology for the characterization and comparison of the virome of sugarcane varieties located at the quarantine glasshouse of Visacane and the virome of wild Poaceae present in the surrounding environment of this sugarcane quarantine station.

## 2. Material and Methods

### 2.1. *Quarantined Sugarcane Sampling*

The youngest fully developed leaf (also called the top visible dewlap leaf or TVD leaf) was sampled from 20 sugarcane plants representing 19 diseased varieties currently maintained in the quarantine glasshouse of Visacane at Cirad in Montpellier (Table 1). Five grams of tissue from the bottom part of each leaf blade were cut in small pieces (approximately 2 × 2 mm) with sterilized instruments and stored in plastic bags at -80 °C until further use.
<table><caption>**Table 1.** List of plant samples processed by virion-associated nucleic acid metagenomics and number of reads obtained for each species/location.</caption><thead><tr><th colspan="3">Quarantined Sugarcane</th><th colspan="2">Wild Poaceae from Flat Stone Land</th><th colspan="2">Wild Poaceae from Sloped Land</th><th colspan="2">Wild Poaceae from Grassland</th></tr><tr><th>Variety</th><th>Number of Reads</th><th>Species</th><th>Number of Reads</th><th>Species</th><th>Number of Reads</th><th>Species</th><th>Number of Reads</th></tr></thead><tbody><tr><td>KN88260</td><td>169,914</td><td>Aegilops triuncialis (3) <sup>a</sup></td><td>1,228,780</td><td>Aegilops triuncialis (3)</td><td>455,952</td><td>Avena sterilis (2)</td><td>376,016</td></tr><tr><td>R579</td><td>337,700</td><td>Aegilops geniculata (3)</td><td>1,051,536</td><td>Aegilops geniculata (3)</td><td>425,570</td><td>Bothriochloa barbinodis (4)</td><td>1,190,788</td></tr><tr><td>KN88104</td><td>364,446</td><td>Anisantha madritensis (6)</td><td>1,661,880</td><td>Avena sterilis (2)</td><td>249,804</td><td>Brachypodium phoenicoides (5)</td><td>885,284</td></tr><tr><td>Q112</td><td>189,526</td><td>Avena sterilis (2)</td><td>800,428</td><td>Bothriochloa barbinodis (3)</td><td>981,718</td><td>Bromopsis erecta (3)</td><td>1,028,474</td></tr><tr><td>B46364-USA51/1</td><td>301,906</td><td>Bothriochloa ischaemum (6)</td><td>931,920</td><td>Trachynia distachya (2)</td><td>374,270</td><td>Cynodon dactylon (2)</td><td>447,196</td></tr><tr><td>KN88147</td><td>483,278</td><td>Trachynia distachya (3)</td><td>1,088,638</td><td>Brachypodium phoenicoides (2)</td><td>859,170</td><td>Festuca sp. (5)</td><td>728,278</td></tr><tr><td>BJ79038</td><td>244,194</td><td>Bromopsis erecta (2)</td><td>1,102,196</td><td>Bromopsis erecta (3)</td><td>152,906</td><td>Gastridium ventricosum (6)</td><td>1,065,580</td></tr><tr><td>KN8924</td><td>625,654</td><td>Bromus lanceolatus (3)</td><td>516,490</td><td>Catapodium rigidum (2)</td><td>681,122</td><td>Helictochloa bromoides (1)</td><td>150,384</td></tr><tr><td>KN8843</td><td>391,108</td><td>Catapodium rigidum (2)</td><td>922,976</td><td>Cynodon dactylon (2)</td><td>539,866</td><td></td><td></td></tr><tr><td>GT9</td><td>596,222</td><td>Cynodon dactylon (2)</td><td>52,6150</td><td>Phalaris minor (6)</td><td>1,161,018</td><td></td><td></td></tr><tr><td>B46364-PAK155</td><td>386,752</td><td>Elytrigia intermedia (6)</td><td>1,270,630</td><td></td><td></td><td></td><td></td></tr><tr><td>X</td><td>269,742</td><td>Hainardia cylindrical (6)</td><td>998,856</td><td></td><td></td><td></td><td></td></tr><tr><td>USDA</td><td>598,072</td><td>Holcus lanatus (2)</td><td>386,482</td><td></td><td></td><td></td><td></td></tr><tr><td>LF653661</td><td>189,262</td><td>Hordeum murinum subsp. Leporinum (6)</td><td>806,796</td><td></td><td></td><td></td><td></td></tr><tr><td>NA021668</td><td>236,036</td><td>Lolium rigidum (6)</td><td>1,353,330</td><td></td><td></td><td></td><td></td></tr><tr><td>SP701284</td><td>267,242</td><td>Melica ciliata (6)</td><td>520,800</td><td></td><td></td><td></td><td></td></tr><tr><td>CR9821</td><td>186,898</td><td>Piptatherum miliaceum (1)</td><td>284,228</td><td></td><td></td><td></td><td></td></tr><tr><td>Q140</td><td>218,214</td><td>Rostraria cristata (3)</td><td>602,560</td><td></td><td></td><td></td><td></td></tr><tr><td>FG087484</td><td>128,726</td><td>Vulpia ciliata (6)</td><td>570,698</td><td></td><td></td><td></td><td></td></tr><tr><td>FR95433</td><td>125,826</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="3">Total (samples) reads</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>(20)</td><td>6,310,178</td><td>(74)</td><td>12,561,594</td><td>(28)</td><td>4,388,287</td><td>(28)</td><td>5,872,000</td></tr></tbody></table>

<footnote><sup>a</sup> The number in parentheses represents the number of samples collected per plant species.</footnote>
## 2.2. Wild Poaceae Sampling

Wild Poaceae were collected from three areas surrounding the quarantine glasshouse, including a flat stone land adjacent to the quarantine glasshouse, a sloped land slightly further from the quarantine glasshouse, and a grassland at least 25 m distant from the glasshouse (Figure 1). One hundred and thirty plants representing 25 wild Poaceae species were collected in these three areas (Table 1). A local botanical expert (co-author Hugo Fontes) confirmed the identity of these species, named in accordance with the Euro+Med PlantBase (https://www.emplantbase.org/home.html, accessed on 31 March 2021). Five grams of leaf and stem tissue from each wild Poaceae were cut in small pieces and stored as described above for sugarcane.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Aerial view of Cirad's campus showing the sugarcane quarantine (Visacane) and the sampling locations used for analysis of the virome of Poaceae in the environment of the sugarcane quarantine facilities. 1, Research laboratories; 2, plant pathology glasshouse; 3, sugarcane quarantine glasshouse; 4 and 5, flat stone land (framed in red) adjacent to the quarantine glasshouse; 6, sloped land (framed in blue) near the quarantine glasshouse, and 7, grassland (framed in green) near the quarantine glasshouse. This image was retrieved from Google Earth.</figcaption></figure>

## 2.3. *Virion-Associated Nucleic Acid-Based Viral Metagenomics*

Each of the 150 collected samples was processed using the virion-associated nucleic acid (VANA)-based viral metagenomics approach [27]. Briefly, 1 g of frozen leaf tissue was ground in Hanks' buffered salt solution (HBSS) (1:10) with four ceramic beads (MP Biomedicals, USA) using a tissue homogenizer (MP biomedicals, USA). The homogenized plant extracts were centrifuged and supernatants were filtered through a 0.45 µm filter and centrifuged at 148,000× g for 2.5 h at 4 °C to concentrate viral particles. Nonencapsidated nucleic acids were eliminated by DNase and RNase incubation at 37 °C for 1.5 h. Total RNA and DNA was then extracted using the NucleoSpin kit (Macherey Nagel). Reverse
transcription was performed with the SuperScript III reverse transcriptase (Invitrogen), cDNAs were purified with the QIAquick PCR Purification Kit (Qiagen) and complementary strands were synthesized using the Klenow DNA polymerase I. Double-stranded DNA was amplified by random PCR amplification. Samples were barcoded during reverse transcription and PCR steps were performed using homemade 26-nt Dodeca Linkers and PCR multiplex identifier primers. PCR products were purified using NucleoSpin gel and PCR clean-up (Macherey Nagel). Wild Poaceae and sugarcane samples were sequenced in two independent Illumina HiSeq runs in order to avoid index-hopping contamination [28]. Bioinformatics analyses were performed as described previously [27]. Briefly, demultiplexing was performed with the agrep command-line tool to assign reads to the samples from which they originated [29]. Adaptors were removed and the reads were filtered for quality (q30 quality and read length > 45 nt) using Cutadapt 1.9 [30]. The cleaned reads were assembled de novo into contigs using SPAdes 3.6.2 [31]. Putative virus reads and plant virus reads obtained using either BLASTn or BLASTx [32] from sugarcane and wild Poaceae with e-values < 0.001 were retained. Finally, only putative plant viruses with more than 10 reads per sample were taken into consideration for further analysis. Cleaned reads have been deposited in the sequence read archive of GenBank (accession number PRJNA721112).

## 2.4. Partial Genomic Characterization of a Novel Sugarcane Umbravirus

RT-PCR reactions were carried out to extend the genomic sequence of a novel umbravirus isolated from sugarcane variety BJ790038, for which four short contigs (131–325 nt) were recovered using the VANA-based metagenomics approach. PCR reactions were performed using the Qiagen OneStep RT-PCR Kit. Ten primers were designed (Table S1) based on the four contig sequences. The 25 µL RT-PCR reaction mix consisted of 1 µL of eluted RNA, 14.5 µL of RNase-free water, 1 µL of RNase inhibitor (RNase-Out, Invitrogen), 5 µL of RT-PCR buffer (5X), 0.5 µL of dNTPmix (10 mM), 1 µL of each primer (10 µM), and 1 µL of the RT-PCR enzyme mix. The RT-PCR program was as follows with the extension time (Ext) for each primer pair listed in Table S1: 50 °C for 30 min, 95 °C for 15 min, 35 cycles at 94 °C for 1 min, annealing temperature 55 °C for 1 min, and 72 °C for Ext with a final 72 °C extension for 10 min. PCR products were analyzed by electrophoresis using a 1.2% agarose gel in TAE buffer stained with ethidium bromide and visualized under UV light. Amplification products were sequenced using the Sanger method (Genewiz, Leipzig, Germany).

## 2.5. Phylogenetic Analyses

Contigs produced by assembly of Illumina reads were utilized as queries to perform BLASTn and BLASTx searches [32] using the NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 31 March 2021). Forty-eight contigs assigned to plant viruses obtained from both compartments (within and outside the quarantine) were selected for further comparisons. These 48 sequences were subsequently aligned to reference genomes using ClustalW with default settings [33]. Visual inspection of alignment quality and BLAST graphic summaries was done and one contig presenting one chimeric end was trimmed. The resulting alignments were used to infer neighbor-joining phylogenetic trees using the MegaX software [34]. The Jukes-Cantor nucleotide substitution model with 1000 bootstrap replicates for branch support was applied. The 48 contigs obtained from quarantined sugarcane varieties and wild Poaceae samples, assigned to the *Caulimoviridae*, *Closteroviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae*, and *Tombusviridae* families, that were used for the phylogenetic analyses described above are listed in the Supplementary Material.
### 3. Results

#### 3.1. Analysis of the Virome of Quarantined Sugarcane and of Wild Poaceae Growing Outside the Sugarcane Quarantine Glasshouse

Total numbers of Illumina reads obtained for each of the 20 samples of quarantined sugarcane varieties varied from 125,826 (FR95433) to 625,654 (KN8924), with an average number of 315,508 reads per sample (Table 1). Retained virus reads represented 24% of total reads whereas the 1483,466 plant virus reads represented 97.6% of the virus reads (Table 2), with an average of 74,173 plant virus reads per quarantined sugarcane sample. The number of reads per putative plant virus family ranged from 1339 (Tombusviridae) to 833,570 (Geminiviridae) (Table 2).

<table><caption>**Table 2.** Viruses identified in wild Poaceae and quarantined sugarcane varieties by virion-associated nucleic acid metagenomics.</caption><thead><tr><th rowspan="2">Putative Virus</th><th rowspan="2">Virus Host *</th><th colspan="2">Wild Poaceae</th><th colspan="2">Quarantined Sugarcane</th></tr><tr><th>Number of Samples</th><th>Number of Reads</th><th>Number of Samples</th><th>Number of Reads</th></tr></thead><tbody><tr><td>Geminiviridae</td><td>Plant</td><td>2</td><td>48,512</td><td>8</td><td>833,570</td></tr><tr><td>Potyviridae</td><td>Plant</td><td>2</td><td>15,029</td><td>3</td><td>414,737</td></tr><tr><td>Tombusviridae</td><td>Plant</td><td>13</td><td>3578</td><td>2</td><td>1339</td></tr><tr><td>Luteoviridae</td><td>Plant</td><td>2</td><td>270</td><td>2</td><td>83,338</td></tr><tr><td>Closteroviridae</td><td>Plant</td><td>1</td><td>368</td><td>7</td><td>10,260</td></tr><tr><td>Alphasatellitidae</td><td>Plant</td><td></td><td></td><td>3</td><td>140,222</td></tr><tr><td>Amalgaviridae</td><td>Plant</td><td>4</td><td>657</td><td></td><td></td></tr><tr><td>Bromoviridae</td><td>Plant</td><td>1</td><td>25</td><td></td><td></td></tr><tr><td>Aspiviridae</td><td>Plant</td><td>1</td><td>19</td><td></td><td></td></tr><tr><td>Sobemovirus</td><td>Plant</td><td>1</td><td>20</td><td></td><td></td></tr><tr><td>Trifolium-ass. DNA virus</td><td>Plant</td><td>3</td><td>587</td><td></td><td></td></tr><tr><td>Caulimoviridae</td><td>Plant/insect</td><td>2</td><td>44</td><td>9</td><td>37,192</td></tr><tr><td>Partitiviridae</td><td>Plant/fungus</td><td>89</td><td>89,391</td><td>5</td><td>210</td></tr><tr><td>Endornaviridae</td><td>Plant/fungus/oomycete</td><td>2</td><td>66</td><td></td><td></td></tr><tr><td>Rhabdoviridae</td><td>Plant/vertebrate/invertebrate</td><td>4</td><td>493</td><td></td><td></td></tr><tr><td>Reoviridae</td><td>Plant/fungus/vertebrate/invertebrate</td><td>1</td><td>413</td><td></td><td></td></tr><tr><td>Chrysoviridae</td><td>Fungus</td><td>29</td><td>4197</td><td></td><td></td></tr><tr><td>Narnaviridae</td><td>Fungus</td><td>1</td><td>15</td><td></td><td></td></tr><tr><td>Hypoviridae</td><td>Fungus</td><td>6</td><td>560</td><td></td><td></td></tr><tr><td>Gammaflexiviridae</td><td>Fungus</td><td>2</td><td>477</td><td></td><td></td></tr><tr><td>Metaviridae</td><td>Fungus</td><td>2</td><td>125</td><td></td><td></td></tr><tr><td>Botybirnavirus</td><td>Fungus</td><td>21</td><td>2554</td><td></td><td></td></tr><tr><td>Tymovirales (unclassified)</td><td>Fungus</td><td>2</td><td>75</td><td></td><td></td></tr><tr><td>Unclassified fungal viruses</td><td>Fungus</td><td>31</td><td>35,671</td><td></td><td></td></tr><tr><td>Totiviridae</td><td>Fungus/protozoan</td><td>40</td><td>60,090</td><td></td><td></td></tr><tr><td>Genomoviridae</td><td>Fungus/human/mammal/bird</td><td>30</td><td>6678</td><td>1</td><td>15</td></tr><tr><td>Microviridae</td><td>Bacteria/spiroplasma</td><td>26</td><td>35,410</td><td></td><td></td></tr><tr><td>Myoviridae</td><td>Bacteria/archaea</td><td>5</td><td>1419</td><td>1</td><td>13</td></tr><tr><td>Podoviridae</td><td>Bacteria/archaea</td><td>5</td><td>6673</td><td></td><td></td></tr><tr><td>Siphoviridae</td><td>Bacteria/archaea</td><td>5</td><td>152,598</td><td></td><td></td></tr><tr><td>Unclassified bacterial viruses</td><td>Bacteria</td><td>3</td><td>108</td><td></td><td></td></tr><tr><td>Iflaviridae</td><td>Insect</td><td>1</td><td>220</td><td></td><td></td></tr><tr><td>Bidnaviridae</td><td>Insect</td><td>1</td><td>32</td><td></td><td></td></tr><tr><td>Parvoviridae</td><td>Insect/vertebrate</td><td>10</td><td>3271</td><td></td><td></td></tr><tr><td>Baculoviridae</td><td>Insect/decapod</td><td>1</td><td>19</td><td>1</td><td>18</td></tr><tr><td>Peribunyaviridae</td><td>Insect/rodent</td><td>1</td><td>16</td><td></td><td></td></tr><tr><td>Unclassified arthropod viruses</td><td>Arthropod</td><td>18</td><td>20,812</td><td></td><td></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th rowspan="2">Putative Virus</th><th rowspan="2">Virus Host *</th><th colspan="2">Wild Poaceae</th><th colspan="2">Quarantined Sugarcane</th></tr><tr><th>Number of Samples</th><th>Number of Reads</th><th>Number of Samples</th><th>Number of Reads</th></tr></thead><tbody><tr><td>Retroviridae</td><td>Vertebrate</td><td>2</td><td>70</td><td></td><td></td></tr><tr><td>Circoviridae</td><td>Bird/mammal</td><td>26</td><td>8208</td><td></td><td></td></tr><tr><td>Iridoviridae</td><td>Amphibia/fish/invertebrate</td><td>1</td><td>23</td><td></td><td></td></tr><tr><td>Phycodnaviridae</td><td>Alga</td><td>15</td><td>450</td><td>1</td><td>19</td></tr><tr><td>Mimiviridae</td><td>Amoebae</td><td>10</td><td>277</td><td></td><td></td></tr><tr><td>Nimaviridae</td><td>Crustacean</td><td>1</td><td>15</td><td></td><td></td></tr><tr><td>Arenaviridae</td><td>Rodent</td><td>1</td><td>19</td><td></td><td></td></tr><tr><td>unclassified</td><td>Vertebrate/invertebrate/crustacean</td><td>11</td><td>3915</td><td></td><td></td></tr><tr><td>unclassified</td><td>environmental</td><td>8</td><td>96,793</td><td></td><td></td></tr><tr><td>unclassified</td><td>unknown</td><td>6</td><td>33,824</td><td></td><td></td></tr></tbody><tfoot><tr><td colspan="2">Total reads</td><td></td><td>634,566</td><td></td><td>1,520,933</td></tr></tfoot></table>

Total numbers of Illumina reads obtained for each of the 130 samples of wild Poaceae collected outside the sugarcane quarantine varied from 40 (*Helictochloa bromoides* from grassland) to 414,878 (*Bromopsis erecta* from flat stone land), with an average number of 175,553 reads per sample. Retained virus reads represented 2.8% of the 22,821,881 reads obtained for the three sampled locations (Table 1), whereas the 69,065 plant virus reads represented 10.9% of the virus reads, with an average of 531 plant virus reads per wild Poaceae sample. The overall majority of viral sequences were assigned to viruses infecting bacteria, fungi, and arthropoda (Table 2).

Contigs from quarantined sugarcane varieties were assigned to seven virus families, i.e., *Alphasatellitidae*, *Caulimoviridae*, *Closteroviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae*, and *Tombusviridae* (Table 2). With the exception of 11 contigs, the contigs obtained from the quarantined sugarcane varieties were all assigned to sugarcane viruses already known to infect the sugarcane plants maintained in the quarantine glasshouse of Visacane at Cirad in Montpellier. The 11 contigs that were not assigned to a known sugarcane virus were produced from two quarantined sugarcane varieties and all shared highest identity with umbraviruses (see below). Plant virus contigs obtained from wild Poaceae samples were distributed into 11 plant virus families (*Amalgaviridae*, *Aspiviridae*, *Bromoviridae*, *Caulimoviridae*, *Closteroviridae*, *Endornaviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae*, *Retroviridae*, and *Tombusviridae*), one unclassified genus (*Sobemovirus*) and one unclassified plant-associated virus, Trifolium-associated circular DNA virus (Table 2). Contigs assigned to *Closteroviridae*, *Caulimoviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae* and *Tombusviridae* families were each produced from wild Poaceae samples and from quarantined sugarcane samples. Specifically, contigs assigned to the *Badnavirus* (*Caulimoviridae*), *Mastrevirus* (*Geminiviridae*) and *Umbravirus* (*Tombusviridae*) genera were present in both quarantined sugarcane samples and wild Poaceae samples.

## 3.2. Phylogenetic Relationships of Plant Viral Sequences Assigned to the Closteroviridae, Luteoviridae, and Potyviridae Families

While the 101 *Closteroviridae* contigs obtained from seven quarantined sugarcane varieties (Q112, B46364-USA51/1, B46364-PAK155, BJ790038, X, LF653661, and NA021668) were assigned to the *Ampelovirus* genus, one *Closteroviridae* contig obtained from a grassland sample of *Brachypodium phoenicoides* was assigned to the *Closterovirus* genus using BLASTx searches (Figure 2A and Table S2). The highest identity score for this 535-nt-long contig obtained from *B. phoenicoides* was obtained with raspberry leaf mottle virus (accession number QOS14265, highest percent identity = 36%, e-value = $3 \times 10^{-10}$), suggesting that it represented a novel closterovirus that needs further characterization. The phylogenetic analyses also showed that the sugarcane *Ampelovirus* contigs grouped together but were also distributed into several subgroups, thus indicating that sugarcane-associated am-

<footnote>* as reported by ViralZone, SIB Swiss Institute of Bioinformatics.</footnote>
peloviruses are diverse and may represent several divergent variants of the same species, or even several species (Figure 2A).

A

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

B

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

**Figure 2.** (A) Neighbor-joining phylogenetic tree of *Closteroviridae* contigs obtained from one wild Poaceae species (bold characters colored in blue), sugarcane varieties (bold characters colored in red), and representative closteroviruses and ampeloviruses genome sequences (colored in black). (B) Neighbor-joining phylogenetic tree of *Luteoviridae* contigs obtained from two wild Poaceae species (bold characters colored in blue), two plants from one sugarcane variety (bold characters colored in red), and representative *Luteoviridae* genome sequences (colored in black). Scale bars represent the number of substitutions per site; branch labels correspond to the consensus support (%). The 12 contigs produced from quarantined sugarcane varieties and wild Poaceae samples, assigned to the *Closteroviridae* and *Luteoviridae* families, are listed in the Supplementary Material.
Two contigs from an *Aegilops triuncialis* plant, one contig from a *Bothriochloa barbinodis* plant, and two large contigs from two sugarcane varieties (B46364-USA51/1 and B46364-PAK155) were assigned to the *Luteoviridae* family using BLAST searches (Table S3). The large contig from each sugarcane variety shared 99% identity with sugarcane yellow leaf virus (SCYLV, accession number KY052166) and represented 97.7% of the complete genome sequence of this SCYLV isolate. The two overlapping contigs from *A. triuncialis* and *B. barbinodis* shared 71% identity among each other, suggesting that they could be variants of the same virus. Identities of the three wild Poaceae contigs with SCYLV (accession number KY052166) ranged from 39 to 43% and from 70.8 to 74.9% with unclassified luteoviruses isolated from insects (Fagle virus, accession number MK440658 and Norway luteo-like virus 1, accession number MF141065). The phylogenetic analysis allowed us to confirm that the *Luteoviridae* contigs obtained from the two wild grasses were highly divergent from plant luteoviruses (Figure 2B).

Complete genome sequences of viral isolates belonging to the *Potyviridae* family were obtained from both quarantined sugarcane varieties and wild Poaceae samples. Three complete genome sequences, each produced from a different sugarcane variety, were assigned to the genus *Potyvirus*. On the other hand, two complete genome sequences obtained from two wild Poaceae samples were assigned to the genus *Tritimovirus* (Table S4). These results confirmed the perennial conservation of two potyviruses (sorghum mosaic virus (SrMV) isolate USA51/1 and sugarcane mosaic virus (SCMV) isolate PAK155) that were inoculated in 2002 in Montpellier on two plants of sugarcane variety B46364. The complete genome sequence obtained from sugarcane variety GT9 was also phylogenetically related to sorghum mosaic virus (Figure 3A), and the genome nucleotide sequences of the two virus isolates were 98% identical. The two complete genome sequences recovered from a sloped land sample of *Brachypodium phoenicoides* and a grassland sample of *Gastridium ventricosum* shared 99% nucleotide identity. This suggested that both sequences belonged to the same *Tritimovirus* species. Furthermore, these two Poaceae tritimovirus sequences shared 67.6% nucleotide identity with oat necrotic mottle virus (accession number AY377938), which is below the species demarcation (<76% nucleotide identity) of the *Potyviridae* family [35]. Consequently, these sequences could represent a novel species of the *Tritimovirus* genus.

### 3.3. Phylogenetic Relationships of Plant Viral Sequences Assigned to the Mastrevirus, Badnavirus, and Umbravirus Genera

Mastrevirus contigs were obtained from eight sugarcane varieties and from two *Anisantha madritensis* plant samples collected from the flat stone land (Table S5). One complete genome sequence of sugarcane streak Egypt virus (SSEV) and five complete genome sequences of sugarcane white streak virus (SWSV) were produced from one sugarcane variety (USDA) and from five sugarcane varieties (USDA, R579, X, KN88147, and KN8924), respectively. The five SWSV complete genome sequences shared 91–99% nucleotide identities among each other. Additionally, three partial SWSV genome sequences were obtained from three other sugarcane plant samples (KN88260, KN88104, and KN8843). The complete genome nucleotide sequences that were recovered from the two plants of *A. madritensis* were 99.9% identical. These two sequences had 75% identity with sorghum arundinaceum-associated virus (MK546381), a mastrevirus recently reported from the Réunion Island. They also shared 44.3–46.2% identity with the SSEV and SWSV genome sequences retrieved from the quarantined sugarcane plant samples. In a phylogenetic tree constructed with entire genome sequences of known mastreviruses and those obtained in this study, the two sequences from the wild Poaceae plants were also located at a unique branch (Figure 3B). According to the current species demarcation (<78% nucleotide identity) for the genus *Mastrevirus* [36], the *A. madritensis*-derived mastrevirus is therefore likely to be a novel species.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** (A) Neighbor-joining phylogenetic tree of *Potyviridae* contigs obtained from wild Poaceae (bold characters colored in blue), sugarcane varieties (bold characters colored in red), and representative potyviruses and tritimovirus genome sequences (colored in black). (B) Neighbor-joining phylogenetic tree of *Geminiviridae* contigs obtained from wild Poaceae species (bold characters colored in blue), sugarcane varieties (bold characters colored in red), and representative mastrevirus genome sequences (colored in black). Scale bars represent the number of substitutions per site; branch labels correspond to the consensus support (%). The 16 contigs obtained from quarantined sugarcane varieties and wild Poaceae samples, assigned to the *Potyviridae* and *Geminiviridae* families, are listed in the Supplementary Material.</figcaption></figure>
Eighty-five contigs assigned to the *Badnavirus* genus (*Caulimoviridae*) were obtained from nine sugarcane samples (B46364 USA51/1, B46364 PAK155, BJ790038, KN8843, KN88104, LF653661, Q112, R579, and X), from which two contigs, each covering the full badnavirus genome, were produced for sugarcane variety BJ790038 and variety X (Table S6). For the wild Poaceae, only one sample of *B. phoenicoides* collected from grassland contained four badnavirus contigs (Table S6) that shared 74.5–77% nucleotide identity with canna yellow leaf mottle virus (CaYMV, accession number KX255725) and 30 to 72% nucleotide identity with the badnavirus contigs obtained from quarantined sugarcane varieties. Based on phylogenetic analyses, the two longest contigs of *B. phoenicoides* clustered with CaYMV and formed a group that was apart from any of the badnaviruses currently known to infect sugarcane (Figure 4A and Figure S1) [37]. The contigs produced from sugarcane varieties R579, B46364 USA51/1, B46364 PAK155, and BJ79038 were related to sugarcane baciliform Guadeloupe A virus (SCBGAV, accession number NC_038382) (Figure 4A). The contigs obtained from the variety X clustered with sugarcane bacilliform Guadeloupe D virus (SCBGDV, accession number NC_013455) (Figure 4A and Figure S1) and contigs obtained from varieties KN88104, KN8843, and LF653661 grouped with banana streak CA virus (accession number HQ593111) (Figure 4A and Figure S1). The contigs from sugarcane variety Q112 matched with two different badnavirus genomes, i.e., SCBGAV and sugarcane bacilliform virus isolate Iscam (accession number JN377534) (data not shown).

Contigs assigned by BLAST searches to the *Tombusviridae* family, including members of the *Umbravirus* genus and unclassified umbraviruses, were obtained from five wild Poaceae species and two sugarcane varieties (BJ79038 and KN8924). These wild Poaceae samples included *Trachynia distachya*, *Bothriochloa barbinodis*, *B. phoenicoides*, *Phalaris minor*, and *Festuca* sp. plants from both the sloped land and the grassland areas (Table S7). *Umbravirus* contigs recovered from wild Poaceae plants shared 45–72% nucleotide identities with a sugarcane umbra-like virus genome (accession number MN868593). Four contigs obtained from two *P. minor* plants, one *Festuca* sp. plant, and one *B. phoenicoides* plant had 66–69% nucleotide identity with strawberry-associated virus A (accession number MK211274). The four contigs also clustered in a neighbor-joining phylogenetic tree, thus suggesting that a single novel umbravirus infected these wild Poaceae plants (Figure 4B). Furthermore, while the contig recovered from the other *Festuca* sp. plant shared 71% nucleotide identity with Patrinia mild mottle virus (accession number MH922775), contigs from *T. distachya* and *B. barbinodis* were 100% identical and had 71–73% nucleotide identity with Ethiopia-maize-associated virus (accession number MF415880). Finally, two partial genome sequences of umbraviruses (1931 nt and 2039 nt, Table S7) were recovered from sugarcane varieties (BJ79038 and KN8924), respectively, using VANA-based reads and RT-PCR assays. These two partial genomes were 71% identical. The 2039-nt-long partial genome from the sugarcane variety KN8924 shared 97% nucleotide identity with Ethiopia-maize-associated virus (accession number MF415880), suggesting that this still-unclassified umbra-like virus, initially isolated from maize, was also infecting sugarcane. The 1931-nt-long partial genome from sugarcane variety BJ79038 shared 72% nucleotide identity with Ethiopia-maize-associated virus (accession number MF415880) and 71% with sugarcane umbra-like virus (accession number MN868593).
<figure><img src="image_8.png" /><figcaption>**Figure 4.** (A) Neighbor-joining phylogenetic trees of *Caulimoviridae* contigs obtained from one wild Poaceae species (bold characters colored in blue), sugarcane varieties (bold characters colored in red), and representative badnaviruses genome sequences (colored in black). (B) Neighbor-joining phylogenetic trees of *Tombusviridae* contigs obtained from wild Poaceae species (bold characters colored in blue), sugarcane varieties (bold characters colored in red), and representative classified and unclassified umbravirus genome sequences (colored in black). Scale bars represent the number of substitutions per site; branch labels correspond to the consensus support (%). The 15 contigs produced from quarantined sugarcane varieties and wild Poaceae samples, assigned to the *Caulimoviridae* and *Tombusviridae* families are listed in the Supplementary Material.</figcaption></figure>
## 4. Discussion

The plant virome of the wild Poaceae collected near the quarantine glasshouse included 11 virus families, one unclassified genus (*Sobemovirus*), and one unclassified plant-virus (trifolium-associated circular DNA virus). In contrast, only seven families formed the plant virome of the quarantined sugarcane varieties. Most virus sequence reads (85%) obtained from wild Poaceae were attributed to non-plant viruses (bacterial viruses: *Microviridae* viruses, bacteriophages, and caudovirales; bird and mammal viruses: *Circoviridae* viruses and *Genomoviridae* viruses; fungal viruses: viruses belonging to the families *Chrysoviridae*, *Partitiviridae*, *Totiviridae*, and *Genomoviridae*; vertebrate and insect viruses: *Partoviridae* viruses). Fungal virus reads related to *Chrysoviridae*, *Partitiviridae*, and *Totiviridae* were previously identified in wild plants, including Poaceae [26,38]. The frequency of this type of reads was very low (1%) in quarantined sugarcane samples for which more than 99% of the virus reads were associated to plant viruses, suggesting that plants growing in quarantine are protected from infections by environmental microorganisms. However, six of the seven virus families found in quarantined sugarcane were also represented in wild Poaceae species growing outside the restricted area. These families included *Caulimoviridae*, *Closteroviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae*, and *Tombusviridae*. The sequences obtained herein for the *Closteroviridae*, *Luteoviridae*, and *Potyviridae* families belonged to different virus genera when comparing virus isolates from quarantined plants and from plants growing in the outside quarantine environment. Nevertheless, even if the virus genera between quarantined and non-quarantined plants are different, the risk of cross infection cannot be eliminated. For the contigs assigned to the three other families, the genera were common between quarantined sugarcane and wild Poaceae. These common genera were identified as *Badnavirus* (*Caulimoviridae* family), *Mastrevirus* (*Geminiviridae* family), and *Umbravirus* (*Tombusviridae* family).

Major plant species infected by badnaviruses have a tropical or sub-tropical origin (aglaonema, alpinia, banana, bougainvillea, cacao, canna, citrus, codonopsis, dracaena, jujube, kalanchoe, pagoda, pineapple, piper, shefflera, stilbocarpa, sugarcane, sweet potato, taro, wisteria, yacon, yam, and yucca) [39]. Badnaviruses were also reported in Europe in rubus [40], grapevine [41], birch [42], and fig [43], but not in wild Poaceae. Badnaviruses are frequently detected in sugarcane growing in tropical and sub-tropical locations of all continents with relative high genetic variability and numerous species groups [44]. The badnavirus sequences found in *B. phoenicoides* were different but close to the sequence identified in sugarcane variety X that belonged to badnavirus group 1 subgroup D [37], thus indicating potential virus transfer between the sugarcane quarantine in Montpellier and surrounding areas. However, the probability that sugarcane is the source of the badnavirus identified in *B. phoenicoides* near Visacane is very low because the same badnavirus was also found in *B. phoenicoides* samples collected at Villeveyrac (12 reads) and Pellissanne (12 reads), two locations that are, respectively, 30 and 120 km distant from the sugarcane quarantine in Montpellier (Roumagnac, unpublished results). However, further investigations are needed to identify the potential risk for sugarcane to be infected by the putative *B. phoenicoides* badnavirus discovered in this study.

Mastreviruses occur in cultivated and wild Poaceae in several continents [18,45–47]. For instance, wild Poaceae host various mastreviruses in Nigeria and are considered as reservoirs for maize streak virus [47]. In Australia, three mastreviruses were identified in wild Poaceae but not in cultivated Poaceae [46]. Up to now, only three mastreviruses have been reported in Europe. Wheat and barley dwarf viruses were initially reported in only a few European countries [48], but were found in additional European countries a few years later [49,50]. A survey undertaken in Germany revealed occurrence of at least three mastreviruses in cultivated Poaceae in this country: barley dwarf virus, oat dwarf virus, and wheat dwarf virus [51]. To our knowledge, the mastrevirus identified in the wild Poaceae *A. madritensis* in our study is the first report of this virus. It could be the cause of new diseases of cultivated Poaceae in Europe. This virus is different from currently known mastreviruses of sugarcane and, although genetically distant from sugarcane streak Egypt
virus (SSEV) and sugarcane white streak virus (SWSV), it might be able to cause disease in sugarcane in cases of accidental spread to healthy sugarcane plants.

Umbraviruses and unclassified umbra-like viruses were the final viruses shared by quarantined sugarcane and wild Poaceae. Umbra-like viruses, which are similar to umbraviruses but do not contain all of the umbravirus genomic features [9], have been previously found in sugarcane in Florida [7] and South Africa [9]. Umbra-like virus sequences have, so far, not been reported in other Poaceae with the exception of Ethiopia-maize-associated virus that was initially considered as an unclassified virus [52]. In this study, sequences related to sugarcane and maize umbra-like viruses were clearly identified in sugarcane varieties KN8924 and BJ79038, as well as in wild Poaceae species (*B. barbinodis* and *T. distachya*). Specifically, Ethiopia-maize-associated virus was detected in the quarantine glasshouse in the sugarcane variety KN8924 that originated from Sudan. Consequently, umbra-like virus diversity does not appear related to plant species and this type of virus could spread among different plant hosts, such as between plants in quarantine and in surrounding area compartments. Several other umbravirus isolates found in wild Poaceae in this study differed genetically from isolates close to the sugarcane and maize umbra-like viruses. Isolates of these putative novel umbraviruses were distributed in two phylogenetic subgroups suggesting that two unclassified umbraviruses, not related to any currently known sugarcane umbravirus, were circulating among several wild Poaceae species growing near the sugarcane quarantine greenhouse in the South of France. These two viruses may represent two novel species of the genus *Umbravirus* while the International Committee on Taxonomy of Viruses recommends a species demarcation threshold of <70% nt sequence identity for umbraviruses [53].

**5. Conclusions**

The area surrounding the sugarcane quarantine of Cirad in the South of France hosts a great diversity of viruses in wild Poaceae, which appears greater than the one occurring in quarantined sugarcane varieties. As these two environments share viruses from the same family and even from the same genus, the risk of cross contamination seems higher than expected. However, we showed in this study that no identical or highly similar virus sequences were present in quarantined material and in plants growing in its surrounding environment. This suggested that plant virus movement between quarantined sugarcane and plants of nearby locations has not, so far, occurred. While a large majority of plant viruses are transmitted by insect vectors [54], these results imply that the confinement procedures of Visacane (including insect-proof structures and security portals) appear to be properly implemented.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/v13050922/s1, Figure S1: Neighbor-joining phylogenetic trees of Caulimoviridae contigs obtained from sugarcane varieties (bold characters colored in red) and representative badnaviruses genome sequences (colored in black). Scale bars represent the number of substitutions per site; Branch label corresponds to the consensus support (%). Table S1: Primer pairs used for amplification of the partial genome sequence of the umbra-like sequence isolated from sugarcane variety BJ790038, Table S1: Wild Poaceae sample collected from grassland and sugarcane samples from the quarantine glasshouse infected by *Closteroviridae* viruses. Table S3: Wild Poaceae samples collected from the sloped land and the grassland and sugarcane samples from the quarantine glasshouse infected by *Luteoviridae* viruses. Table S4: Wild Poaceae samples collected from the sloped land, the grassland and sugarcane samples from the quarantine glasshouse infected by *Potyviridae* viruses. Table S5: Wild Poaceae samples collected from the flat stone land and sugarcane samples from the quarantine glasshouse infected by mastreviruses (*Geminiviridae* family). Table S6: Selected contig sequences of wild Poaceae samples collected from the grassland and sugarcane samples from the quarantine glasshouse potentially infected by badnaviruses (*Caulimoviridae* family). Table S7: Wild Poaceae samples collected from the sloped land, the grassland and from sugarcane samples located in the quarantine glasshouse, containing *Umbravirus* or umbra-like virus contigs. Contigs in bold were used for phylogenetic analyses. Material S1: Sequences of the 48 contigs obtained from quarantined sugar-
cane varieties and wild Poaceae samples, assigned to the *Caulimoviridae*, *Closteroviridae*, *Geminiviridae*, *Luteoviridae*, *Potyviridae*, and *Tombusviridae* families, that were used for the phylogenetic analyses.

**Author Contributions:** J.H.D., D.F., and P.R. conceived and designed the experiments; J.H.D., D.F., C.J., L.C., R.F., E.F., H.F., and P.R. collected the sugarcane and wild Poaceae samples, C.J., L.C., H.F., and E.F. performed the experiments; J.H.D., D.F., H.F., and P.R. analyzed the data; J.H.D., P.C.R., and P.R. wrote the paper which was edited and approved by all authors. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Agropolis Foundation (E-Space flagship program, grant number 1504-004).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Cleaned Illumina reads determined as part of this study have been deposited in the sequence read archive of GenBank (accession number PRJNA721112).

**Acknowledgments:** We are grateful to the CIRAD Southgreen bioinformatics platform and INRAE MIGALE bioinformatics facilities for providing computing resources.

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

## References

1. Martin, R.R.; Constable, F.; Tzanetakis, I.E. Quarantine regulations and the impact of modern detection methods. *Annu. Rev. Phytopathol.* **2016**, *54*, 189-205. [CrossRef]

2. Rott, P.; Bousquet, J.-F.; Muller, E.; Chatenet, M. La quarantaine de canne à sucre du CIRAD à Montpellier. *Agric. Dév.* **1997**, *13*, 22-28.

3. Guinet-Brial, I.; Girard, J.-C.; Roumagnac, A.; Daugrois, J.-H.; Fernandez, E.; Rott, P. Visacane: An innovative quarantine tool for the exchange of pest and disease-free sugarcane germplasm. In *International Society of Sugar Cane Technologists, Proceedings of the XXVIIIth Congress, São Paulo, Brazil, 24-27 June 2013*; Sociedade dos Técnicos Açucareiros e Alcooleiros do Brasil (STAB); São Paulo, Brazil; The XXVIIIth ISSCT Organising Committee: São Paulo, Brazil, 2013; pp. 890-896.

4. Fernandez, E.; Ferdinand, R.; Filloux, D.; Guinet, I.; Julian, C.; Rott, P.; Roumagnac, P.; Daugrois, J.-H. Improvements in virus detection at CIRAD's sugarcane quarantine using both viral metagenomics and PCR-Based approaches. *Sugar Tech* **2020**, *22*, 498-503. [CrossRef]

5. Candresse, T.; Filloux, D.; Muhire, B.; Julian, C.; Galzi, S.; Fort, G.; Bernardo, P.; Daugrois, J.-H.; Fernandez, E.; Martin, D.P.; et al. Appearances can be deceptive: Revealing a hidden viral infection with deep sequencing in a plant quarantine context. *PLoS ONE* **2014**, *9*, e102945. [CrossRef]

6. Roumagnac, P.; Cirad, F.; Mollov, D.; Daugrois, J.-H.; Filloux, D.; Usda-Ars, U. Viral Metagenomics and Sugarcane Pathogens. In *Understanding the Behaviour and Improving the Welfare of Dairy Cattle*; Burleigh Dodds Science Publishing: Cambridge, UK, 2018; Volume 2, pp. 183-200.

7. Filloux, D.; Fernandez, E.; Comstock, J.C.; Mollov, D.; Roumagnac, P.; Rott, P. Viral metagenomic-based screening of sugarcane from florida reveals occurrence of six sugarcane-infecting viruses and high prevalence of sugarcane yellow leaf virus. *Plant Dis.* **2018**, *102*, 2317-2323. [CrossRef]

8. Rott, P.; Bailey, R.A.; Comstock, J.C.; Croft, B.J.; Saumtally, A.S. *A Guide to Sugarcane Diseases*; Cirad: Montpellier France, 2000.

9. Tahir, M.N.; Bolus, S.; Grinstead, S.C.; McFarlane, S.A.; Mollov, D. A new virus of the family Tombusviridae infecting sugarcane. *Arch. Virol.* **2021**, *166*, 961-965. [CrossRef] [PubMed]

10. Jones, R.A.C.; Coutts, B.A. Spread of introduced viruses to new plants in natural ecosystems and the threat this poses to plant biodiversity. *Mol. Plant Pathol.* **2015**, *16*, 541-545. [CrossRef] [PubMed]

11. Redinbaugh, M.G.; Zambrano, J.L. Control of Virus Diseases in Maize. *Adv. Clin. Chem.* **2014**, *90*, 391-429. [CrossRef]

12. Sun, S.-R.; Chen, J.-S.; Yang, J.; Huang, X.-C.; Huang, M.-T.; Gao, S.-J. First Report of Maize yellow mosaic virus Infecting Sugarcane in China. *Plant Dis.* **2019**, *103*, 2482. [CrossRef]

13. Yahaya, A.; Dangora, D.B.; Alabi, O.J.; Zongoma, A.M.; Kumar, P.L. Detection and diversity of maize yellow mosaic virus infecting maize in Nigeria. *VirusDisease* **2019**, *30*, 538-544. [CrossRef]

14. Grisham, M.P. Strains of sorghum mosaic virus causing sugarcane mosaic in Louisiana. *Plant Dis.* **1994**, *78*, 729-732. [CrossRef]

15. Chen, J.; Chen, J. Sugarcane mosaic disease in Zhejiang Province was caused by sorghum mosaic virus and sugarcane mosaic virus. *Bing Du Xue Bao* **2002**, *18*, 362-366.

16. Chalupniková, J.; Kundu, J.K.; Singh, K.; Bartaková, P.; Beoni, E. Wheat streak mosaic virus: Incidence in field crops, potential reservoir within grass species and uptake in winter wheat cultivars. *J. Integr. Agric.* **2017**, *16*, 523-531. [CrossRef]

17. Fargette, D.; Konaté, G.; Fauquet, C.; Muller, E.; Peterschmitt, M.; Thresh, J.M. Molecular Ecology and Emergence of Tropical Plant Viruses. *Annu. Rev. Phytopathol.* **2006**, *44*, 235-260. [CrossRef] [PubMed]
18. Claverie, S.; Ouattara, A.; Hoareau, M.; Filloux, D.; Varsani, A.; Roumagnac, P.; Martin, D.P.; Lett, J.-M.; Lefeuvre, P. Exploring the diversity of Poaceae-infecting mastreviruses on Reunion Island using a viral metagenomics-based approach. *Sci. Rep.* **2019**, *9*, 1–11. [CrossRef]

19. Ellis, M.H.; Rebetzke, G.J.; Moore, C.S.; Hyles, J.E.; Kelman, W.M. Detection of Wheat streak mosaic virus in four pasture grass species in Australia. *Plant Pathol.* **2004**, *53*, 239. [CrossRef]

20. Malmstrom, C.M.; McCullough, A.J.; Johnson, H.A.; Newton, L.A.; Borer, E.T. Invasive annual grasses indirectly increase virus incidence in California native perennial bunchgrasses. *Oecologia* **2005**, *145*, 153–164. [CrossRef]

21. Kumar, A.; Murthy, S.; Kapoor, A. Evolution of selective-sequencing approaches for virus discovery and virome analysis. *Virus Res.* **2017**, *239*, 172–179. [CrossRef]

22. Roossinck, M.J.; Martin, D.P.; Roumagnac, P. Plant Virus metagenomics: Advances in virus discovery. *Phytopathology* **2015**, *105*, 716–727. [CrossRef]

23. Boukari, W.; Alcalá-Briseño, R.I.; Kraberger, S.; Fernandez, E.; Filloux, D.; Daugrois, J.-H.; Comstock, J.C.; Lett, J.-M.; Martin, D.P.; Varsani, A.; et al. Occurrence of a novel mastrevirus in sugarcane germplasm collections in Florida, Guadeloupe and Réunion. *Virol. J.* **2017**, *14*, 146. [CrossRef]

24. François, S.; Bernardo, P.; Filloux, D.; Roumagnac, P.; Yaverkovski, N.; Froissart, R.; Ogliastro, M. A novel itera-like densovirus isolated by viral metagenomics from the Sea Barley Hordeum marinum. *Genome Announc.* **2014**, *2*, 2. [CrossRef] [PubMed]

25. Kraberger, S.; Farkas, K.; Bernardo, P.; Booker, C.; Argüello-Astorga, G.R.; Mesléard, F.; Martin, D.P.; Roumagnac, P.; Varsani, A. Identification of novel Bromus- and Trifolium-associated circular DNA viruses. *Arch. Virol.* **2015**, *160*, 1303–1311. [CrossRef]

26. Bernardo, P.; Charles-Dominique, T.; Barakat, M.; Ortet, P.; Fernandez, E.; Filloux, D.; Hartnady, P.A.; Rebelo, T.; Cousins, S.R.; Mesleard, F.; et al. Geometagenomics illuminates the impact of agriculture on the distribution and prevalence of plant viruses at the ecosystem scale. *ISME J.* **2018**, *12*, 173–184. [CrossRef] [PubMed]

27. François, S.; Filloux, D.; Fernandez, E.; Ogliastro, M.; Roumagnac, P. Viral metagenomics approaches for high-resolution screening of multiplexed arthropod and plant viral communities. *Breast Cancer* **2018**, *1746*, 77–95. [CrossRef]

28. van der Valk, T.; Vezzi, F.; Ormestad, M.; Dalén, L.; Guschanski, K. Estimating the rate of index hopping on the Illumina HiSeq X platform. *bioRxiv* **2018**, *179028*. [CrossRef]

29. Wu, S.; Manber, U. A fast approximate pattern-matching tool. In Proceedings of the Usenix Winter 1992 Technical Conference, San Francisco, CA, USA, 20–24 January 1992; pp. 153–162.

30. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet* **2011**, *17*, 10. [CrossRef]

31. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. *J. Comput. Biol.* **2012**, *19*, 455–477. [CrossRef] [PubMed]

32. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. *J. Mol. Biol.* **1990**, *215*, 403–410. [CrossRef]

33. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; Mcgettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. *Bioinformatics* **2007**, *23*, 2947–2948. [CrossRef]

34. Stecher, G.; Tamura, K.; Kumar, S. Molecular Evolutionary Genetics Analysis (MEGA) for macOS. *Mol. Biol. Evol.* **2020**, *37*, 1237–1239. [CrossRef]

35. Wylie, S.J.; Adams, M.; Chalam, C.; Kreuze, J.; López-Moya, J.J.; Ohshima, K.; Praveen, S.; Rabenstein, F.; Stenger, D.; Wang, A.; et al. ICTV Virus Taxonomy Profile: Potyviridae. *J. Gen. Virol.* **2017**, *98*, 352–354. [CrossRef] [PubMed]

36. Muhire, B.; Martin, D.P.; Brown, J.K.; Navas-Castillo, J.; Moriones, E.; Zerbini, F.M.; Rivera-Bustamante, R.; Malathi, V.G.; Briddon, R.W.; Varsani, A. A genome-wide pairwise-identity-based proposal for the classification of viruses in the genus Mastrevirus (family Geminiviridae). *Arch. Virol.* **2013**, *158*, 1411–1424. [CrossRef]

37. Muller, E.; Dupuy, V.; Blondin, L.; Bauffe, F.; Daugrois, J.-H.; Nathalie, L.; Iskra-Caruana, M.-L. High molecular variability of sugarcane bacilliform viruses in Guadeloupe implying the existence of at least three new species. *Virus Res.* **2011**, *160*, 414–419. [CrossRef]

38. Roossinck, M.J.; Saha, P.; Wiley, G.B.; Quan, J.; White, J.D.; Lai, H.; Chavarría, F.; Shen, G.; Roe, B.A. Ecogenomics: Using massively parallel pyrosequencing to understand virus ecology. *Mol. Ecol.* **2010**, *19*, 81–88. [CrossRef]

39. Sastry, K.S.; Mandal, B.; Hammond, J.; Scott, S.W.; Briddon, R.W. *Encyclopedia of Plant Viruses and Viroids*; Springer: Berlin/Heidelberg, Germany, 2019.

40. Diaz-Lara, A.; Mosier, N.J.; Keller, K.E.; Martin, R.R. A variant of Rubus yellow net virus with altered genomic organization. *Virus Genes* **2014**, *50*, 104–110. [CrossRef] [PubMed]

41. Basso, M.F.; Fajardo, T.V.M.; Saldarelli, P. Grapevine virus diseases: economic impact and current advances in viral prospection and management. *Rev. Bras. Frutic.* **2017**, *39*, e-411. [CrossRef]

42. Rumbou, A.; Candresse, T.; Marais, A.; Theil, S.; Langer, J.; Jalkanen, R.; Büttner, C. A novel badnavirus discovered from *Betula* sp. affected by birch leaf-roll disease. *PLoS ONE* **2018**, *13*, e0193888. [CrossRef]

43. Minafra, A.; Chiumenti, M.; Elbeaino, T.; Digiaro, M.; Bottalico, G.; Pantaleo, V.; Martelli, G. Occurrence of fig Badnavirus 1 in fig trees from different countries and in symptomless seedlings. *J. Plant Pathol.* **2012**, *94*, 4–105.

44. Wu, X.; Alabi, O.J.; Damaj, M.B.; Sun, S.; Mirkov, T.E.; Fu, H.; Chen, R.; Gao, S. Prevalence and RT/RNase H Geneology of Sugarcane Bacilliform Virus Isolates from China. *J. Phytopathol.* **2016**, *164*, 595–607. [CrossRef]
45. Agindotan, B.O.; Domier, L.L.; Bradley, C.A. Detection and characterization of the first North American mastrevirus in switchgrass. *Arch. Virol.* **2015**, *160*, 1313–1317. [CrossRef]

46. Hadfield, J.; Martin, D.P.; Stainton, D.; Kraberger, S.; Owor, B.E.; Shepherd, D.N.; Lakay, F.; Markham, P.G.; Greber, R.S.; Briddon, R.W.; et al. Bromus catharticus striate mosaic virus: A new mastrevirus infecting Bromus catharticus from Australia. *Arch. Virol.* **2010**, *156*, 335–341. [CrossRef]

47. Oluwafemi, S.; Kraberger, S.; Shepherd, D.N.; Martin, D.P.; Varsani, A. A high degree of African streak virus diversity within Nigerian maize fields includes a new mastrevirus from Axonopus compressus. *Arch. Virol.* **2014**, *159*, 2765–2770. [CrossRef] [PubMed]

48. Manurung, B.; Witsack, W.; Mehner, S.; Grüntzig, M.; Fuchs, E. The epidemiology of Wheat dwarf virus in relation to occurrence of the leafhopper Psammotettix alienus in Middle-Germany. *Virus Res.* **2004**, *100*, 109–113. [CrossRef]

49. Achon, M.A.; Serrano, L.; Ratti, C.; Rubies-Autonell, C. First Detection of Wheat dwarf virus in Barley in Spain Associated with an Outbreak of Barley Yellow Dwarf. *Plant Dis.* **2006**, *90*, 970. [CrossRef] [PubMed]

50. Lemmetty, A.; Huusela-Veistola, E. First report of wheat dwarf virus in winter wheat in finland. *Plant Dis.* **2005**, *89*, 912. [CrossRef]

51. Schubert, J.; Habekuß, A.; Kazmaier, K.; Jeske, H. Surveying cereal-infecting geminiviruses in Germany—Diagnostics and direct sequencing using rolling circle amplification. *Virus Res.* **2007**, *127*, 61–70. [CrossRef]

52. Adams, I.P.; Braidwood, L.A.; Stomeo, F.; Phiri, N.; Uwumukiza, B.; Feyissa, B.; Mahuku, G.; Wangi, A.; Smith, J.; Mumford, R.; et al. Characterising maize viruses associated with maize lethal necrosis symptoms in sub Saharan Africa. *bioRxiv* **2017**, 161489.

53. King, A.M.Q.; Adams, M.J.; Carstens, E.B.; Lefkowitz, E.J. *Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses*; Elsevier: Amsterdam, The Netherlands, 2012.

54. Dietzgen, R.G.; Mann, K.S.; Johnson, K.N. Plant virus–insect vector interactions: Current and potential future research directions. *Viruses* **2016**, *8*, 303. [CrossRef] [PubMed]