Article

# Characterization of Six Ampeloviruses Infecting Pineapple in Reunion Island Using a Combination of High-Throughput Sequencing Approaches

Delphine Massé<sup>1,2,*</sup>, Thierry Candresse<sup>3,□</sup>, Denis Filloux<sup>4,5</sup>, Sébastien Massart<sup>6</sup>, Nathalie Cassam<sup>1</sup>, Bruno Hostachy<sup>1</sup>, Armelle Marais<sup>3,□</sup>, Emmanuel Fernandez<sup>4,5</sup>, Philippe Roumagnac<sup>4,5,□</sup>, Eric Verdin<sup>7</sup>, Pierre-Yves Teycheney<sup>8,□</sup>, Jean-Michel Lett<sup>8</sup> and Pierre Lefeuvre<sup>8,□</sup>

1. ANSES—LSV RAPT, F-97410 St. Pierre, La Réunion, France; nathalie.cassam@anses.fr (N.C.); bruno.hostachy@anses.fr (B.H.)

2. UMR PVBMT, Université de La Réunion, F-97410 St. Pierre, La Réunion, France

3. INRAe, UMR 1332 Biologie du Fruit et Pathologie, Université Bordeaux, CS20032, F-33882 Villenave d'Ornon, France; thierry.candresse@inrae.fr (T.C.); armelle.marais-colombel@inrae.fr (A.M.)

4. CIRAD, UMR PHIM, F-34090 Montpellier, France; denisfilloux@cirad.fr (D.F.); emmanuel.fernandez@cirad.fr (E.F.); philippe.roumagnac@cirad.fr (P.R.)

5. PHIM Plant Health Institute, Université Montpellier, CIRAD, INRAE, Institut Agro, IRD, F-34090 Montpellier, France

6. Plant Pathology Laboratory, Gembloux Agro-Bio Tech, University of Liège, 5030 Gembloux, Belgium; sebastien.massart@uliege.be

7. INRAe, UR407 Unité de Pathologie Végétale, CS 60094, F-84140 Montfavet, France; eric.verdin@inrae.fr

8. CIRAD, UMR PVBMT, F-97410 St. Pierre, La Réunion, France; pierre-yves.teycheney@cirad.fr (P.-Y.T.); jean-michel.lett@cirad.fr (J.-M.L.); pierre.lefeuvre@cirad.fr (P.L.)

* Correspondence: delphine.masse@anses.fr

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

check for
updates

**Citation:** Massé, D.; Candresse, T.; Filloux, D.; Massart, S.; Cassam, N.; Hostachy, B.; Marais, A.; Fernandez, E.; Roumagnac, P.; Verdin, E.; et al. Characterization of Six Ampeloviruses Infecting Pineapple in Reunion Island Using a Combination of High-Throughput Sequencing Approaches. *Viruses* **2024**, *16*, 1146. https://doi.org/10.3390/v16071146

Academic Editor: Svetlana
Y. Folimonova

Received: 19 June 2024

Revised: 9 July 2024

Accepted: 9 July 2024

Published: 16 July 2024

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

**Copyright:** © 2024 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:** The cultivation of pineapple (*Ananas comosus*) is threatened worldwide by mealybug wilt disease of pineapple (MWP), whose etiology is not yet fully elucidated. In this study, we characterized pineapple mealybug wilt-associated ampeloviruses (PMWaVs, family *Closteroviridae*) from a diseased pineapple plant collected from Reunion Island, using a high-throughput sequencing approach combining Illumina short reads and Nanopore long reads. Reads co-assembly resulted in complete or near-complete genomes for six distinct ampeloviruses, including the first complete genome of pineapple mealybug wilt-associated virus 5 (PMWaV5) and that of a new species tentatively named pineapple mealybug wilt-associated virus 7 (PMWaV7). Short reads data provided high genome coverage and sequencing depths for all six viral genomes, contrary to long reads data. The 5' and 3' ends of the genome for most of the six ampeloviruses could be recovered from long reads, providing an alternative to RACE-PCRs. Phylogenetic analyses did not unveil any geographic structuring of the diversity of PMWaV1, PMWaV2 and PMWaV3 isolates, supporting the current hypothesis that PMWaVs were mainly spread by human activity and vegetative propagation.

**Keywords:** *Ananas comosus*; pineapple mealybug wilt-associated viruses; ampelovirus; high-throughput sequencing; short and long reads; molecular diversity

## 1. Introduction

Pineapple (*Ananas comosus* L. [Merr.]) is a tropical plant in the family *Bromeliaceae*. It is grown extensively throughout the tropics and subtropics for its edible fruits. In 2021, worldwide pineapple production was estimated at 28.6 million metric tons, making it the third largest tropical fruit production [1]. Among the many pests and diseases that hamper the cultivation of pineapple, mealybug wilt of pineapple (MWP) is considered the most important and complex [2]. Typical symptoms of MWP include reddening and/or yellowing of the leaves, downward curling of margin tips, leaf-tip dieback, wilting and
root decay [3–5]. Since its first description in Hawaii in 1910 [6], MWP has been described
in most pineapple growing areas and reported to cause yield reductions of up to 55% [7,8].

The involvement of a virus in the etiology of MWP was first supported by immuno-
logical evidence [9]. Long, flexuous, rod-shaped virus particles and high molecular weight
double stranded RNAs (i.e., nucleic acid forms originating mostly from viruses) were
later observed in infected plants [10,11], and associated viruses were assigned to genus
*Ampelovirus* [2,12] in the family *Closteroviridae* [13]. Ampeloviruses have single stranded
positive-sense RNA genomes ranging in size between 13.0 and 18.5 kb, with 7 to 12 open
reading frames (ORFs) [14]. They are distributed in two subgroups (I and II) depending on
genome organisation and phylogenetic relationships [15,16]. Thirteen species are currently
assigned to the *Ampelovirus* genus, including *Ampelovirus unananas*, *Ampelovirus duananas*
and *Ampelovirus triananas*, whose exemplar isolates are pineapple mealybug wilt-associated
virus 1 (PMWaV1), PMWaV2 and PMWaV3, respectively [15]. Additional ampeloviruses
were described in pineapple but not yet included in the taxonomy by the International Com-
mittee on Taxonomy of Viruses (ICTV): PMWaV4, which was later reassigned as a strain
of PMWaV1 [17]; PMWaV5 [7,17], for which only a partial genome sequence is available;
and more recently PMWaV6, for which a complete genome sequence was obtained using
high-throughput sequencing (HTS) [18]. Ampeloviruses are transmitted by mealybugs
*(Pseudococcidae)* and scale insects (*Coccidae*), in a semi-persistent mode [14,19]. There is
no evidence of seed or mechanical transmission. Long distance dissemination primarily
occurs through the exchange of infected propagation materials and germplasm [20].

The etiology of MWP has not yet been fully elucidated. Biotic and abiotic factors,
such as environmental factors, the presence of mealybugs and co-infections with distinct
ampeloviruses, and/or with viruses from other genera may be involved. In Hawaii, only
pineapple plants infected by PMWaV2 and infested with mealybugs were reported to dis-
play MWP symptoms [21]. However, another study conducted in Australia suggested that
the presence of PMWaV1 and PMWaV3 together with other viruses such as badnaviruses,
was involved in the etiology of MWP [7]. Additional viruses infecting pineapple were
recently characterized [22–25], further increasing the number of virus combinations pos-
sibly involved in the etiology of MWP, and challenging serological or PCR-based viral
diagnostics. HTS-based approaches, which have enabled the discovery and detection of
new viruses and viral strains in a large range of plants [26–28], have been successfully
applied to the characterization of several ampeloviruses infecting pineapple [18,29] and
have a definite potential for comprehensive diagnosis of pineapple-associated viruses.

The current diversity of pineapple viruses and the still-unresolved etiology of MWP
argue for a better characterization of the pineapple virome. Hence, we compared and
assessed two different HTS strategies for the characterization of ampelovirus genomes
using total RNA extracted from the leaves of a pineapple plant showing MWP symptoms.
The first strategy is based on cDNA synthesis from ribodepleted total RNAs and Illumina
short reads sequencing, whereas the second is based on cDNA synthesis from poly-A tailed
RNAs and Nanopore long read sequencing. Using the combination of these two methods,
we obtained near complete genome sequences of PMWaV2, PMWaV3, PMWaV5, and
PMWaV6, and complete genome sequences of PMWaV1 and a putative new ampelovirus
tentatively named pineapple mealybug wilt-associated virus 7 (PMWaV7). Phylogenetic
analyses of PMWaV1, PMWaV2 and PMWaV3 suggested that recurrent exchanges of
infected germplasm promoted the worldwide dissemination of these viruses.

**2. Materials and Methods**

2.1. *Plant Material*

A plant of cultivar "Queen Victoria" showing reddening and leaf tip dieback symptoms was collected in 2016 from Cirad's experimental station at Bassin Plat in Saint Pierre (Reunion Island, France). Two full leaves from this plant (referred to as sample 16-1) were cut into ~0.5 cm² pieces and stored at -80 °C until further use.
## 2.2. RNA Extraction

Total RNAs were extracted from the 16-1 leaf sample using the RNeasy Plant Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions, including an optional DNAse treatment (DNAse I, Thermo Fisher Scientific, Waltham, MA, USA). The suitability of extracted RNAs quality for HTS sequencing was assessed by measuring the A260/280 (~2.0) and A260/230 (2.0 to 2.2) ratios using a NanoVue spectrophotometer (Cytiva, Washington, DC, USA) and the RNA Integrity Number (RIN, >6.0) value using an Agilent TapeStation and an RNA 6000 Pico chip (Agilent Technologies, Santa Clara, CA, USA). Extracted RNAs were quantified using a Qubit Fluorometer 2.0 (Qubit RNA Broad Range Assay Kit, Thermo Fisher Scientific, Waltham, MA, USA). Total RNA samples were stored at -80 °C until further use.

## 2.3. High-Throughput Sequencing

Two HTS approaches were used in the study. In the first approach (hereafter referred to as "Illumina short reads approach"), two Illumina sequencing runs were performed. For both runs, total RNAs were first ribodepleted prior to cDNA synthesis and sequencing. The first run was performed according to Marais et al. [30]. Briefly, ribosomal RNAs were removed from the total RNAs using a Ribominus Plant Kit (Invitrogen, Fisher Scientific, Illkirch, France). A cDNA library was generated using the Illumina TruSeq Stranded total RNA library prep kit (Illumina, San Diego, CA, USA) and sequenced on an Illumina NextSeq500 sequencer (2 × 150 bp) (GIGA-Genomics Facility, Université de Liège, Liège, Belgium). For the second run, libraries preparation and sequencing were performed by Genewiz (now Azenta Life Sciences, South Plainfield, NJ, USA) using the Ribozero plant rRNA depletion kit and the TruSeq stranded total RNA library Prep kit (Illumina, San Diego, CA, USA), followed with a paired-ends 2 × 150 bp run on an Illumina NovaSeq sequencer.

In the second approach, hereafter referred to as "Nanopore long reads approach", nanopore sequencing was performed using the MinION device (Mk1B, Oxford Nanopore Technologies, Oxford, UK). For this, a cDNA library was generated from the total RNAs using poly(dT) primers, with the cDNA-PCR Barcoding kit SQK-PCB109 (Oxford Nanopore Technologies, Oxford, UK) according to the manufacturer's instructions, except that Sera-Mag Select Size Selection beads (Cytiva, Washington, DC, USA) were used for the nucleic acid purification step. Two distinct sequencing runs were performed. In the first one, the RNA sample was treated in quadruplicate using distinct barcodes. A cDNA library with equimolar concentrations of products from each replicate was produced. In the second run, a single cDNA library was prepared from the RNA sample. In both runs, libraries were loaded on FLO-MIN 106D R9.4.1 flowcells and sequencing was performed using a Mk1B MinION device and monitored using the MinKNOW v21.11.8 software (ONT) for 16 and 28 h, respectively.

## 2.4. Data Analysis and Assembly of Viral Genomes

Sequencing data were uploaded on the Galaxy instance of the Migale bioinformatics server (https://migale.inrae.fr) and analysed using a homemade workflow described below. Demultiplexing and quality control of short reads were performed using Trimmomatic v0.38.1 [31]. After adapters and tags removal using default settings, reads were quality trimmed using a minimum quality score of 25 and a minimum read length of 80 nt. For long reads, highly accurate base calling was performed using the Guppy v6.0.1 software (ONT), and reads that passed the default quality threshold (mean read PHRED score of 7) were processed with PoreChop v0.2.4 [32] for adapters removal and demultiplexing when needed. Reads shorter than 80 nt were removed using NanoFilt v0.1.0 [33].

Hybrid de novo assembly of short and long reads was performed using SPAdes v3.13 [34] on the SouthGreen bioinformatic platform [35]. Then, viral contigs were identified using BLASTn and BLASTx searches against the reference viral sequences database retrieved from GenBank in April 2022 (v211). When possible, viral contigs were merged using the EGAssembler [36] before being polished using Pilon v1.23 [37] after mapping
back short reads using the BWA-MEM algorithm (BWA v0.7017.1) [38]. Each set of trimmed and demultiplexed reads was individually mapped against polished viral contig references using BWA-MEM for short reads and Minimap2 v2.17 [39] for long reads, respectively. After discarding non-primary alignments, mapping statistics were obtained using Samtools v1.9 [40]. Subsamplings of the mapped reads were performed in order to estimate the coverage of viral contigs as a function of total reads number. To this end, we evaluated the number of reads mapped per position for sets of decreasing sequencing efforts with 100 random replicates per set. Sequencing depth (the number of times a position was covered with a read) and breadth of coverage (the proportion of the genome that is covered with reads) were computed for each subsample.

## 2.5. Detection of PMWaVs by RT-PCR

RT-PCR diagnosis was performed in order to confirm the presence of the viruses detected by HTS. The detection of PMWaV1, PMWaV2 and PMWaV3 was performed by a new multiplex RT-PCR using virus-specific primers [2,7] (Supplementary Table S1). Additional simplex RT-PCRs were carried out to detect PMWaV5, PMWaV6 and PMWaV7 using virus-specific primers (Supplementary Table S1) designed from available complete viral genome sequences.

Reactions were performed using the One-Step RT-PCR kit (QIAGEN, Hilden, Germany) in a Veriti 96-Well Fast Thermal Cycler (Applied Biosystems™, Thermo Fisher Scientific, Waltham, MA, USA). For the multiplex detection of PMWaV1, PMWaV2 and PMWaV3, the reaction mix contained 2.5 µL of 5X RT-PCR Buffer, 0.5 µL of 10 mM dNTPs, 1.25 µL of a primer mix containing PMWaV1 and PMWaV3 primers at a concentration of 10 mM, 2 µL of a primer mix containing the PMWaV2 specific primers at a concentration of 10 mM, 2 µL of RT-Taq DNA polymerase mix, 2 µL of total plant RNA extracts and nuclease-free water for a total reaction volume of 25 µL. For the simplex RT-PCR detection of PMWaV5, PMWaV6 or PMWaV7, the reaction mix contained 5 µL of 5X RT-PCR buffer, 1 µL of 10 mM dNTPs, 1.5 µL of a primer mix containing each specific primer at a concentration of 10 mM, 1 µL of RT-Taq DNA polymerase mix, 2 µL of total plant RNA extracts and nuclease-free water for a total reaction volume of 25 µL. RT-PCR conditions were a reverse-transcription step of 30 min at 50 °C followed by an initial denaturation step of 15 min at 95 °C and 35 cycles of 30 s at 94 °C, 30 s at 58 °C (for PMWaV1, PMWaV2 and PMWaV3) or 60 °C (for PMWaV5, PMWaV6, and PMWaV7), and 1 min at 72 °C, then a final extension of 10 min at 72 °C. Amplification products were analyzed by electrophoresis on 2.5% agarose gels in a Tris-acetate-EDTA (TAE) buffer and visualized under UV light following staining with ethidium bromide.

## 2.6. RACE PCR

Rapid amplification of cDNA 3'-ends (3'-RACE) and 5'-ends (5'-RACE) experiments were performed for viruses for which genomic 5' and 3' ends were missing from the assemblies of HTS reads. The SMARTer RACE 5'/3' Kit (Takara, Beijing, China) was used following the manufacturer's instructions. cDNAs were amplified by PCR using the universal primer A mix provided by the manufacturer and a virus-specific primer designed from the assembled genome sequence (Supplementary Table S2). The resulting PCR amplicons were analyzed by electrophoresis as described above. Bands corresponding to amplification products of the expected sizes were excised from the gels under UV light, purified using the Nucleospin PCR and Gel Purification Kit (Macherey-Nagel, Düren, Germany) and ligated into the CloneJet PCR cloning vector (Thermo Scientific, Waltham, MA, USA). Inserts of the recombinant vectors were sequenced by standard Sanger sequencing (Macrogen Inc., Seoul, Republic of Korea) using a primer walking approach.

## 2.7. *Search for Recombination and Phylogenetic Analyses*

Sequence alignments were performed using MAFFT v7.450 in Geneious Prime soft v2021.1 [41]. Inter-species recombination signals were searched using the RDP [42],
GENECONV [43], BOOTSCAN [44], MAXCHI [45], CHIMERA [46], SISCAN [47] and 3SEQ [48] methods implemented in the RDP5 program [49], using default settings. Recombination events were considered significant when detected (p-value < 0.05) by at least three different recombination detection methods.

For phylogenetic analyses, subsets of nucleotide sequences encoding viral RNA-dependent RNA polymerase (RdRp), heat shock protein 70 homolog (HSP70h) and coat protein (CP) sequences were obtained from the complete genome sequences, and conceptually translated into amino acid sequences. Homologous sequences of beet yellows virus (genus *Closterovirus*) were used as outgroups. Sequences alignments were performed using MAFFT v7.450 [41]. Maximum likelihood phylogenetic reconstruction was performed using FastTree v2.1.11 [50]. The LG substitution model with gamma-distributed rate among sites was used and the Shimodaira-Hasegawa-like test (SH-like) for branch support (1000 replicates) was performed.

### 3. Results

#### 3.1. Identification of Known and Novel Ampeloviruses

HTS performed on sample 16-1 generated 136.9 M paired reads (representing ~20.1 Gb) and 4.9 M reads (representing ~3.3 Gb) from the Illumina short reads and Nanopore long reads approaches, respectively. After trimming, the size of long reads ranged from 80 to 8053 nt with a mean size of 615 nt, and that of short reads ranged from 80 to 150 nt with a mean size of 122 nt. A total of 97,021 contigs (ranging in size between 80 and 56,804 nt) was obtained from de novo assembly performed on both short and long reads using the hybrid assembly approach. After similarity searches against viral databases, 183 of these contigs were identified as viral sequences. After aligning and merging these 183 viral contigs, a second similarity search was performed against viral databases using BLASTn and revealed the presence of six supercontigs (referred to as Contig-A to Contig-F) ranging in size between ~13 and 18 kbp, with similarities to viruses from the genus *Ampelovirus* (Table 1), especially PMWaV1, PMWaV2, PMWaV3 and PMWaV6 for which similarities ranged between 79% and 99% nucleotide identity (E-values of 0) (Table 1).

<table><caption>**Table 1.** Characteristics of the contigs with similarities to ampeloviruses.</caption><thead><tr><th>Supercontigs ID (Length bp)</th><th>Number of Contigs</th><th>Contigs Length (bp)</th><th>Number of Reads</th><th>Cov ≥ 10<sup>1</sup> (%)</th><th>Cov ≥ 100<sup>1</sup> (%)</th><th>Mean Depth</th><th>Most Similar Virus on NCBI (BLASTn)</th><th>Accession Numbers</th><th>Identity (%)<sup>2</sup></th><th>Query Cover (%)</th><th>Virus Species</th></tr></thead><tbody><tr><td>Contig-A (13,074)</td><td>67</td><td>86–13,073</td><td>42,491</td><td>99.9</td><td>98.3</td><td>415</td><td>Pineapple mealybug wilt-associated virus 1</td><td>OP860292</td><td>86.3</td><td>99.3</td><td>PMWaV1</td></tr><tr><td>Contig-B (16,199)</td><td>38</td><td>81–5816</td><td>200,516</td><td>100.0</td><td>99.9</td><td>1683</td><td>Pineapple mealybug wilt-associated virus 2</td><td>OP860299</td><td>99.4</td><td>99.2</td><td>PMWaV2</td></tr><tr><td>Contig-C (13,229)</td><td>59</td><td>80–5040</td><td>28,779</td><td>99.8</td><td>98.3</td><td>285</td><td>Pineapple mealybug wilt-associated virus 3</td><td>MN539274</td><td>96.3</td><td>99.4</td><td>PMWaV3</td></tr><tr><td>Contig-D (12,971)</td><td>8</td><td>159–6146</td><td>19,534</td><td>99.8</td><td>90.6</td><td>197</td><td>Pineapple mealybug wilt-associated virus 5</td><td>EF467920</td><td>84.1</td><td>13.1</td><td>PMWaV5</td></tr><tr><td>Contig-E (17,440)</td><td>4</td><td>466–7117</td><td>7304</td><td>97.5</td><td>1.1</td><td>41</td><td>Pineapple mealybug wilt-associated virus 6</td><td>OP86029</td><td>99.1</td><td>99.1</td><td>PMWaV6</td></tr><tr><td>Contig-F (18,388)</td><td>7</td><td>389–18,092</td><td>23,269</td><td>100.0</td><td>85.1</td><td>157</td><td>Grapevine leafroll-associated virus 3</td><td>KY073324</td><td>79.3</td><td>0.5</td><td>PMWaV7</td></tr></tbody></table>

Other contigs with similarities with viruses from genera *Vitivirus* (family *Betaflexiviridae* [25]), *Sadwavirus* (family *Secoviridae* [22–24]) and *Badnavirus* (family *Caulimoviridae* [21]) genera, were also obtained (Supplementary Table S3).

<footnote n="1">¹ Percent genome coverage at 10X and 100X. ² E value was always of 0 except for PMWaV7 with 8 × 10⁻⁷.</footnote>
The presence of a strand-switching primer (SSP) sequence and poly-A tail within some of the long reads allowed the identification of the putative 5'-ends for five out of the six assembled genomes (except for Contig-D) and the putative 3'-ends of Contig-A and Contig-F. Using RACE PCR and custom primers (Supplementary Table S2) designed from the sequences of the contigs, the 5'-ends of Contig-D, Contig-E and Contig-F and the 3'-end of Contig-F were confirmed. Our attempts to obtain the sequences of 3'-ends of Contig-B, Contig-C, Contig-D and Contig-E remained unsuccessful. The presence of the six identified ampeloviruses in sample 16-1 was confirmed using RT-PCR assays performed with virus-specific primers (Supplementary Table S1).

## 3.2. *Species Identification and Genome Organization*

The sizes of the six assembled ampelovirus genomes, including the 5' and 3' ends, ranged between 12,971 and 18,388 nt (Table 1 and Figure 1). These genomes shared a common organisation typical of either ampelovirus subgroup I (with seven ORFs) or subgroup II (with 10 or 11 ORFs) (Figure 1 and Supplementary Table S4).

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

**Figure 1.** Genomic organization of the six ampelovirus isolates whose genomes were sequenced in this work. ORFs are represented as coloured boxes. Size scale (in nucleotides) is provided at the bottom of the figure. The acronyms of ampeloviruses from subgroup I and subgroup II are coloured in blue and orange, respectively. Abbreviations are as follows: L-Pro, leader papain-like protease; Mtr, methyltransferase; Hel, helicase; RdRp, RNA-dependent RNA polymerase; HSP70h, heat-shock protein 70 homolog; CP, coat protein; CPm, minor coat protein; other ORFs are named according to the weight in kDa of the protein they encode, preceded with the letter "p" for protein.

The RdRp, HSP70h and CP amino acid sequences conceptually translated from the six assembled ampelovirus genomes were used in phylogenetic analyses to assign these genomes to existing or new ampelovirus species, following the current ICTV criteria for species demarcation in genus *Ampelovirus* of less than 75% amino acid (aa) identity [14]. RdRps, HSP70hs, and CPs encoded by Contigs-A, -B, -C and -E displayed 86.4 to 99.2% identity with homologous proteins of previously characterized isolates of PMWaV1, PMWaV2, PMWaV3, and PMWaV6, respectively (Table 2), providing evidence that Contigs-A, -B, -C and -E can be safely assigned to isolates of PMWaV1, PMWaV2, PMWaV3 and PMWaV6, respectively, named hereafter PMWaV1-RUN, PMWaV2-RUN, PMWaV3-RUN and PMWaV6-RUN, respectively. The genomes of the four PMWaV6 isolates from Hawaii for which
genome sequences are publicly available, and that of PMWaV6-RUN varied in size due to insertion/deletion, notably between the ORFs encoding the RdRp and p6. The genome of PMWaV6-RUN (17,355 nt) displayed a 470 nt deletion similar to that of isolate 6/S1-1 from Hawaii. None of the other three isolates of PMWaV6 described in Hawaii (S1, S2-1 and S3-1), whose sizes range from 17,225 to 18,213 nt, displayed this deletion.

<table><caption>**Table 2.** Pairwise nucleotide and amino acid sequence identities in the RNA-dependant RNA polymerase (RdRp); heat shock protein 70 homolog (HSP70h); coat protein (CP) of PMWaV1, PMWaV2, PMWaV3, PMWaV5, PMWaV6 and PMWaV7 isolates whose genome were sequenced in this work; and the most closely related ampeloviruses.</caption><thead><tr><th>Protein</th><th>Supercontig</th><th>Virus Species</th><th>Most Closely Related Virus (BLASTp)</th><th>Nucleotides Identities (%)</th><th>Amino Acid Identities (%) <sup>1</sup></th></tr></thead><tbody><tr><td rowspan="6">RdRp</td><td>Contig-A</td><td>PMWaV1</td><td>NC_010178—Pineapple mealybug wilt-associated virus 1</td><td>90.6</td><td>90.7</td></tr><tr><td>Contig-B</td><td>PMWaV2</td><td>NC_043105—Pineapple mealybug wilt-associated virus 2</td><td>98.9</td><td>98.4</td></tr><tr><td>Contig-C</td><td>PMWaV3</td><td>NC_043406—Pineapple mealybug wilt-associated virus 3</td><td>97.6</td><td>98.6</td></tr><tr><td>Contig-D</td><td>PMWaV5</td><td>EF467922—Pineapple mealybug wilt-associated virus 5</td><td>85.6</td><td>85.7</td></tr><tr><td>Contig-E</td><td>PMWaV6</td><td>MW269512—Pineapple mealybug wilt-associated virus 6</td><td>99.5</td><td>98.8</td></tr><tr><td>Contig-F</td><td>PMWaV7</td><td>NC_004667—Grapevine leafroll-associated virus 3</td><td>55.2</td><td>43.0</td></tr><tr><td rowspan="6">HSP70h</td><td>Contig-A</td><td>PMWaV1</td><td>NC_010178—Pineapple mealybug wilt-associated virus 1</td><td>88.9</td><td>86.4</td></tr><tr><td>Contig-B</td><td>PMWaV2</td><td>NC_043105—Pineapple mealybug wilt-associated virus 2</td><td>97.5</td><td>95.9</td></tr><tr><td>Contig-C</td><td>PMWaV3</td><td>NC_043406—Pineapple mealybug wilt-associated virus 3</td><td>97.2</td><td>94.8</td></tr><tr><td>Contig-D</td><td>PMWaV5</td><td>EF467920—Pineapple mealybug wilt-associated virus 5</td><td>84.7</td><td>91.8</td></tr><tr><td>Contig-E</td><td>PMWaV6</td><td>MW269512—Pineapple mealybug wilt-associated virus 6</td><td>99.3</td><td>99.1</td></tr><tr><td>Contig-F</td><td>PMWaV7</td><td>NC_004667—Grapevine leafroll-associated virus 3</td><td>55.2</td><td>41.2</td></tr><tr><td rowspan="6">CP</td><td>Contig-A</td><td>PMWaV1</td><td>NC_010178—Pineapple mealybug wilt-associated virus 1</td><td>91.1</td><td>94.9</td></tr><tr><td>Contig-B</td><td>PMWaV2</td><td>NC_043105—Pineapple mealybug wilt-associated virus 2</td><td>99.5</td><td>99.2</td></tr><tr><td>Contig-C</td><td>PMWaV3</td><td>NC_043406—Pineapple mealybug wilt-associated virus 3</td><td>97.2</td><td>97.3</td></tr><tr><td>Contig-D</td><td>PMWaV5</td><td>NC_043406—Pineapple mealybug wilt-associated virus 3</td><td>66.6</td><td>60.6</td></tr><tr><td>Contig-E</td><td>PMWaV6</td><td>MW269512—Pineapple mealybug wilt-associated virus 6</td><td>99.4</td><td>98.6</td></tr><tr><td>Contig-F</td><td>PMWaV7</td><td>NC_022072—Blackberry vein banding-associated virus</td><td>46.9</td><td>34.3</td></tr></tbody></table>

The RdRp and HSP70h conceptually translated from the complete genome sequence corresponding to Contig-D shared 85.7% and 91.8% aa identity, respectively, with their counterparts encoded by the only publicly available and partial PMWaV5 genome sequence (EF467920) [6], providing a strong indication that Contig-D corresponds to the complete genome of a Reunion isolate of PMWaV5 tentatively named PMWaV5-RUN (OQ850040), making it the first complete genome of a PMWaV5 isolate. Sequence comparisons of the CP sequence of this isolate with that of the partially sequenced genome could not be carried out, since this partial sequence (EF467920) does not encompass the ORF encoding the CP. The genome of PMWaV5-RUN displays the typical organisation of subgroup II ampeloviruses, with seven putative ORFs (Supplementary Table S4). ORF1a encodes a putative polyprotein of 226.8 kDa with the domains of a papain-like protease (L-Pro), a methyltransferase (Mtr) and a helicase (Hel) involved in replication [51]. ORF1b encodes a putative 52.3 kDa protein with the conserved motifs of an RNA-dependent RNA polymerase (RdRp). ORF2 encodes a putative movement protein of 5.7 kDa [52]. ORF3, ORF4, ORF5 and ORF6 encode a putative HSP70h (58.4 kDa), p61 (61.2 kDa), CP (30 kDa) and coat protein minor (CPm; 24.3 kDa), respectively.

The RdRp and HSP70h conceptually translated from the complete genome sequence corresponding to Contig-F (18,388 nt) were most closely related to those of grapevine leafroll-associated virus 3 (GRLaV3, NC_004667), with 43.0% and 41.2% aa identities, respectively, whereas the CP was most closely related to that of the blackberry vein banding-associated virus (BVBaV, NC_022072), with 34.3% aa identity (Table 2). Therefore, following the criterion for the demarcation of ampelovirus species, the complete genome sequence assembled from Contig-F belongs to the new subgroup I ampelovirus (Figure 2), for which, we propose the name pineapple mealybug wilt-associated virus 7 (PMWaV7) according to the current species name nomenclature. ORF1a encodes a putative polyprotein of 219 kDa

<footnote n="1">¹ Values below the current ICTV species discrimination criteria (<75% aa identity) are shown in bold.</footnote>
with the domains L-Pro, Mtr and Hel; ORF1b encodes a putative RdRp of 53.8 kDa; ORFs 2 and 3 encode putative proteins of 3.3 kDa (p3) and 4.6 kDa (p5), respectively; ORF4, 5, 6 and 7 encodes a putative HSP70h (60.2 kDa), p57 (56.8 kDa), CP (34.3 kDa) and CPm (52.8 kDa), respectively (Supplementary Table S4). ORF8 and ORF9 encode putative proteins of 16.9 kDa (p20) and 26.5 kDa (p27), which are potentially involved in systemic transport and RNA silencing suppression, respectively [51]. No alkylation B (AlkB) domain was found in the replication-related proteins of the Reunion isolates of PMWaV1, PMWaV2, PMWaV3, PMWaV5, PMWaV6 and PMWaV7. No inter-species recombination was detected between the isolates of these six PMWaVs.

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

**Figure 2.** Maximum-likelihood phylogenetic trees showing the placement of isolates of pineapple mealybug wilt-associated viruses (PMWaV1 to 7) within genus *Ampelovirus*, using the amino acid sequences of the RdRp (A), HSP70h (B) and CP (C). Sequences obtained during this study are indicated in bold. Values at nodes represent bootstrap supports of the branches to their left. The scale bars represents the number of substitutions per site. Beet yellows virus (BYV, *Closterovirus*) was used as an outgroup and solely the *Ampelovirus* clade was presented for convenience. *Ampelovirus* subgroups I and II are shown in blue and yellow, respectively. Partial sequences are indicated with asterisks. PMWaV4 was recently classified as a strain of PMWaV1 and is indicated with a hashtag. For virus acronyms, see Supplementary Table S5.

### 3.3. Phylogenetic Analyses

Phylogenetic analyses using RdRp, HSP70h, or CP aa sequences showed the placement of PMWaV2-RUN, PMWaV6-RUN and PMWaV7-RUN in *Ampelovirus* subgroup I and that of PMWaV1-RUN, PMWaV3-RUN and PMWaV5-RUN in subgroup II (Figure 2).

We took advantage of the availability in public databases of PMWaV1, PMWaV2 and PMWaV3 HSP70h sequences originating from samples collected in various locations covering most pineapple growing regions to carry out a phylogeographic analysis of the diversity of these viruses. We used 48, 32 and 21 publicly available HSP70h aa sequences for PMWaV1, PMWaV2 and PMWaV3, respectively (Figure 3). Mean sequence identities (aa) were 95.5%, 95.9% and 92.8% for the PMWaV1, PMWaV2 and PMWaV3 isolates, respectively.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Maximum-likelihood phylogenetic trees built from the comparison of HSP70h amino acid sequences of pineapple mealybug wilt-associated virus 1 (A), 2 (B) and 3 (C) isolates from Reunion with sequences available from GenBank. Sequences are colored according to their geographical origin and the color legend is indicated on the map on top of the figure. Isolates from Reunion are indicated in brown and in bold. HSP70h sequences from complete genome sequences are indicated with two black asterisks. Bootstrap values equal to or greater than 70% are indicated by solid black circles. PMWaV3 (YP_0099666005), PMWaV6 (QZB90243) and PMWaV1 (YP_001642339) isolates from Hawaii were used to root the PMWaV1, 2 and 3 trees, respectively (not presented on the figures). The scale bars show the number of substitutions per site.</figcaption></figure>

There was no clear geographic structuring of the diversity, regardless of the considered *Ampelovirus* species. Sequences originating from given locations, such as PMWaV1 sequences originating from Hawaii, were sometimes scattered throughout the phylogenies
and/or were more closely related to sequences originating from distant locations than to sequences originating from nearby locations. For example, the isolate PMWaV1-RUN was more closely related to isolates from Hawaii and Australia than to an isolate from nearby Mauritius (Figure 3A). Likewise, the isolate PMWaV2-RUN described in this work was grouped with highly similar sequences (minimum percentage aa identity of 99.51%) originating from Cuba, Mauritius, Thailand, Taiwan, Brazil and Hawaii, but was more distantly related to previously described isolates from Reunion [53] (Figure 3B). It must be noted however, that none of the branch of the PMWaV2 tree showed strong bootstrap support. Lastly, the isolate PMWaV3-RUN was grouped with isolates from various locations such as Australia, Cuba, Taiwan and Hawaii (Figure 3C).

### 3.4. Comparison between the Illumina Short Reads and Nanopore Long Reads Approaches

All the reads obtained using the two HTS approaches were mapped back to the six assembled ampelovirus genomes and the coverage statistics were calculated (Table 3). Between 7235 and 198,147 short reads and 17 and 2374 long reads were mapped, depending on the virus considered. For Illumina short reads, all the genomes were entirely covered, and 97.1% to 99.8% of their length was covered at a 10× sequencing depth. PMWaV2-RUN displayed the highest mean sequencing depth (N = 1606). Due to the much lower number of reads obtained with the Nanopore long reads approach, sequencing depths were drastically lower compared to the Illumina short reads approach (Table 3). Some genomes presented with relatively large mean sequencing depths, such as those of PMWaV2-RUN (N = 77), and PMWaV1-RUN (N = 53.3). However, the genomes of PMWaV5-RUN and PMWaV6-RUN were poorly covered, with mean sequencing depths of 0.4 and 1.5, respectively.

<table><caption>**Table 3.** Reads mapping statistics for Illumina short and Nanopore long reads approaches (after quality control).</caption><thead><tr><th>Virus Species</th><th>Approach</th><th>Number of Reads</th><th>Cov ≥ 1<sup>1</sup> (%)</th><th>Cov ≥ 10<sup>1</sup> (%)</th><th>Mean Depth</th><th>Number Viral Bases/ Million Sequenced Bases</th></tr></thead><tbody><tr><td rowspan="2">PMWaV1</td><td>Illumina short reads</td><td>41,417</td><td>100</td><td>99.8</td><td>361.0</td><td>310.6</td></tr><tr><td>Nanopore long reads</td><td>1074</td><td>100</td><td>82.1</td><td>53.3</td><td>234.7</td></tr><tr><td rowspan="2">PMWaV2</td><td>Illumina short reads</td><td>198,147</td><td>100</td><td>100</td><td>1606</td><td>1710.1</td></tr><tr><td>Nanopore long reads</td><td>2374</td><td>99.7</td><td>72.0</td><td>77.0</td><td>420</td></tr><tr><td rowspan="2">PMWaV3</td><td>Illumina short reads</td><td>28,561</td><td>100</td><td>99.7</td><td>275.0</td><td>238.9</td></tr><tr><td>Nanopore long reads</td><td>218</td><td>94.8</td><td>35.5</td><td>10.1</td><td>45.3</td></tr><tr><td rowspan="2">PMWaV5</td><td>Illumina short reads</td><td>19,517</td><td>99.9</td><td>99.8</td><td>197.0</td><td>168.1</td></tr><tr><td>Nanopore long reads</td><td>17</td><td>18.8</td><td>0.0</td><td>0.4</td><td>1.8</td></tr><tr><td rowspan="2">PMWaV6</td><td>Illumina short reads</td><td>7235</td><td>99.9</td><td>97.1</td><td>39.4</td><td>45.2</td></tr><tr><td>Nanopore long reads</td><td>69</td><td>53.9</td><td>3.7</td><td>1.5</td><td>8.7</td></tr><tr><td rowspan="2">PMWaV7</td><td>Illumina short reads</td><td>22,937</td><td>99.9</td><td>99.7</td><td>148.8</td><td>179.9</td></tr><tr><td>Nanopore long reads</td><td>332</td><td>94.0</td><td>30.4</td><td>8.3</td><td>51.7</td></tr></tbody></table>

We further analyzed the sequencing depths of Illumina short reads and Nanopore long reads along the genomes of PMWaV1, 2, 3, 5, 6 and 7 (Figure 4).

The sequencing depths for Illumina short reads was relatively stable along genome sequences (ranging from 39 to 1605), except for PMWaV1 and, to a lesser extent, PMWaV7 (Figure 4). On the contrary, large variations in coverage were observed for Nanopore long reads along the genomes of all six studied viruses (ranging from 0.4 to 77), with relatively higher depths for the 5' and 3' ends of the genomes, except for PMWaV3 and PMWaV5 (Figure 4D). Indeed, seven of the twelve genome extremities sequences were obtained (5' end essentially) using Nanopore long reads sequencing (from 1 to 108 reads covering the 5' and 3' ends), whereas three of these 12 genome extremities sequences were obtained using Illumina short reads (from 5 to 254 reads covering the 5' and 3' ends). Very few reads were obtained for PMWaV5, and none could be mapped to the extremities of its genome using either Nanopore long or Illumina short reads sequencing.

<footnote n="1">¹ Percent genome coverage at 1X and 10X.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Mean sequencing depths for Illumina short reads (red curves) and Nanopore long reads (blue curves) along the genomes of PMWaV1 (A), PMWaV2 (B), PMWaV3 (C), PMWaV5 (D), PMWaV6 (E) and PMWaV7 (F). The x-axis represents the genomic position (in bases, b) from the 5' end to 3' end and the y-axis represents the fold coverage (Log10 scale). Colored curves were obtained after smoothing (window size of 20) of the raw depth (light grey). Mean depth values are indicated on the right of each curve.</figcaption></figure>

Regarding the proportion of viral bases in relation to the total number of sequenced bases, important differences were observed between viral genomes (Table 3). For Illumina short reads, this proportion ranged between 34.4 bases per million of sequenced bases for
PMWaV6 and 1292.0 for PMWaV2. For Nanopore long reads, it ranged between 1.6 for PMWaV5 and 373.5 for PMWaV2. If these proportions were considered proxies of the viral loads, PMWaV2 would be the virus with the highest viral load in the sampled plant, and PMWaV5 and PMWaV6 the viruses with the lowest viral loads.

In order to evaluate the correlation between the sequencing effort (i.e., the amount of raw sequences obtained per sample) and genome coverage, we calculated the sequencing depth (i.e., the number of times a nucleotide is read during sequencing) and breadth (i.e., the proportion of nucleotide positions in the consensus sequence relative to the length of the complete genome sequence) for sets of subsampled reads (Figure 5A,B). We obtained the distribution of the expected genome coverage (y-axis on Figure 5A,B) for every virus and for increasing sequencing efforts (x-axis on Figure 5A,B). A total of 15.2 Gb of sequenced data for the Illumina short reads, corresponding to the overall number of bases obtained through this approach after quality control, revealed that all ampelovirus genomes sequences exhibited an almost complete 10× depth (a minimum of 97% for PMWaV6), but a decrease in coverage breadth was observed when reducing the sequencing effort (Figure 5). Importantly, the slope of the curve (representing the rate at which bases are covered when increasing the number of reads) was mostly similar for all six viruses (Figure 5A). For a sequencing effort of 1 Gb, corresponding to a common Illumina run in which 96 samples would have been multiplexed, 10× sequencing depth would range between 29.3% and 98.6% for all viral genomes except PMWaV6. Although not all genome regions would be properly covered, at least five of the viral genomes would have been detected (with an estimation of ~1200 to 13,000 minimum reads for all of the six species). Only PMWaV6 may remain less detected, with only ~475 reads expected at 1 Gb. For this virus, the interpolation curve showed that the breadth of coverage at a 10× sequencing depth reached 50% of the genome for a 4.7 Gb sequencing effort (Figure 5A).

<figure><img src="image_7.png" /><figcaption>**Figure 5.** In interpolation of sequence coverage (10×, y-axis) for each PMWaV for the Illumina short (A) and Nanopore long (B) reads approaches according to the number of sequenced bases (x-axis in Log10 scale).</figcaption></figure>

For the Nanopore long reads approach, a 3 Gb sequencing effort corresponding to the total number of sequenced bases obtained after quality control when combining the two runs, did not allow the complete sequencing of any viral genome and resulted in a 10× sequencing depth for 0 to 82% of the nucleotide positions depending on the virus considered (Figure 5B). Very few reads were obtained for PMWaV5 (N = 17), for which
no position was covered with a 10× depth. Contrary to the results obtained using the Illumina short reads, the slopes of the curves differed slightly between viruses (Figure 5B), suggesting that some genomes were preferentially sequenced over others. If the sequencing effort was reduced to 1 Gb, the proportion of the viral genomes sequenced at a 10× depth would be drastically reduced: no position of the genomes of PMWaV5 and PMWaV6 would be covered with a 10× coverage and it would be below 10% for that of PMWaV3 and PMWaV7. PMWaV1 and PMWaV2 would display 32.9% to 47.1% genome breadths with a 10× depth, respectively. The estimated number of reads that would be obtained for each of the six viruses targeted by this study would range from ~5 to 800.

**4. Discussion**

In this study, we recovered near complete genome sequences of isolates from six ampeloviruses infecting a wilt-diseased pineapple plant in Reunion using two different HTS approaches. One sequence obtained in this study represents the first complete genome sequence of PMWaV5, and another is the complete genome sequence of a new ampelovirus, for which we propose the name pineapple wilt-associated virus 7 (PMWaV7). The remaining four viruses (PMWaV1, PMWaV2, PMWaV3 and PMWaV6) were described previously [24], and three of them (PMWaV1, PMWaV2 and PMWaV3) were previously reported in Reunion from wilt-diseased pineapple plants [53]. Our work provides the first evidence that PMWaV5, PMWaV6 and PMWaV7 isolates are also present in Reunion, adding to the complexity of the etiology of MWP in Reunion. Mixed viral infections have been repeatedly reported in pineapple. Indeed, contigs with similarities to viruses from the *Vitivirus* (family *Betaflexiviridae*, [25]), *Sadwavirus* (family *Secoviridae*, [22–24]) and *Badnavirus* (family *Caulimoviridae*, [21]) genera were also identified from the 16-1 pineapple sample used in this work, providing evidence that this plant was infected by at least nine viruses from four distinct families. Vegetative propagated crops such as pineapple are known to accumulate viruses [17,24,25,54] because they do not undergo sexual reproduction, which acts as natural sanitation against viruses since the majority of plant viruses are not seed-transmitted [55]. Co-infections by large numbers of viruses may favour synergistic effects, which in the case of pineapple, could influence the severity of MWP symptoms.

Pineapple was disseminated throughout the tropics and subtropics from its South American centre of origin (Paraguay, southern Brazil and northern Argentina) in less than 600 years through navigation routes [56–58], likely resulting in the spread of pineapple pests and diseases worldwide. Further intensification of plant material exchanges may also have contributed to the spread of these pests and pathogens. Indeed, the phylogeographic analyses reported in this work did not provide evidence for a geographical structuring of the diversity of PMWaV1, PMWaV2 and PMWaV3 isolates, supporting the hypothesis that PMWVs were mainly disseminated through recurrent exchanges of infected germplasm. The analysis of a larger number of samples collected worldwide, including ancient samples conserved in herbaria throughout the world, could help refine this scenario and better understand the migration routes of PMWVs and other pineapple viruses. A similar approach proved successful for unravelling the origin of viruses of grasses [59], grapevine [60] and cassava [61].

Our previous studies on PMWVs [18,24,62,63] highlighted the lack of knowledge about the diversity of viruses infecting pineapple and the potential of HTS to help fill this gap. In this study, we used a combination of HTS Illumina short and Nanopore long reads to tackle this issue. Whereas the first approach requires access to Illumina sequencing devices usually affordable for medium to large laboratories, the second uses the Oxford Nanopore MinION sequencing device, which is easy to implement even in small research facilities with average levels of equipment. A probable limitation of the Nanopore long reads sequencing approach described in this work lies in the use of a cDNA-PCR barcoding kit designed for the sequencing of polyadenylated RNAs. Although this approach is suitable for mRNA sequencing, it is not suitable for not-polyadenylated RNA, such as those of the ampelovirus genomes. The sequencing was nevertheless effective because
of the priming of poly-T primer in A-rich regions of ampelovirus genomes, although
probably not as much as with polyadenylated RNA templates, resulting in the observed
reduced efficiency of the Nanopore long reads approach. The use of random priming
for cDNA synthesis could be a more efficient option for ampelovirus sequencing and
viruses' discovery in general. Nevertheless, the combination of the Illumina short reads
and Nanopore long reads approaches enabled the assembly of complete or near complete
genomes of PMWaV1, PMWaV2, PMWaV3, PMWaV5, PMWaV6 and PMWaV7, including
the sequences of most 5' ends and some of the 3' ends of these viral genomes without
having to resort to time-consuming RACE PCRs. Therefore, the approach described in this
paper paves the way for the development of HTS-based viral indexing and sequencing
of complete viral genomes, for which technical and management guidelines have been
produced to cover the process of implementing HTS technologies in a research or diagnostic
laboratory (selection, development, verification and validation) to detect plant pathogens
and pests [64].

**5. Conclusions**

Our results underscore the existence of important viral diversity among ampeloviruses
infecting pineapple in Reunion. Additional work is now required to assess the biology of
these viruses, their role in the etiology of MWP, implement better risk assessment and to
design appropriate disease management methods for the control of MWP [65].

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v16071146/s1. Table S1: Sequence of the primers used in RT-PCR assays for the detection of PMWaVs. Table S2: Identification of the 5' and 3' ends of the genome sequences of PMWaV species using RACE PCR and Nanopore long reads approaches. Table S3: Read mapping statistics for Illumina short and Nanopore long reads approaches for contigs other than *Ampeloviruses*. Table S4: Genomic organization of PMWaV sequences described on pineapple in Reunion. Table S5: List of sequences used for pairwise identity comparisons and phylogenetic analysis. (Ref. [66] is cited in Supplementary Materials Table S1).

**Author Contributions:** D.M., J.-M.L. and P.L. designed the experiments; D.M., N.C., E.F. and S.M. performed the experiments; D.M., P.L., D.F. and T.C. analyzed the data; D.M., P.L., P.-Y.T. and J.-M.L. wrote the initial version of the manuscript. All authors have revised the initial draft. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was co-funded by the European Union (ERDF, INTERREG VI), the Con-seil Régional de La Réunion, CIRAD and ANSES, and conducted on the Plant Protection Platform (3P, IBISA).

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

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** New sequences described in this paper were deposited on GenBank under accession numbers OQ850037-OQ850042.

**Acknowledgments:** We are grateful to the INRAE MIGALE bioinformatics facility (MIGALE, INRAE, 2020. Migale bioinformatics Facility, doi:10.15454/1.5572390655343293E12) for providing resources through its Galaxy server. This study was supported by the European Cost Action FA1407.

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

**References**

1. FAOSTAT. Food and Agriculture Organization of the United Nations. In *FAOSTAT Statistics Database*; FAO: Rome, Italy, 2023.

2. Sether, D.M.; Melzer, M.J.; Busto, J.; Zee, F.; Hu, J.S. Diversity and mealybug transmissibility of ampeloviruses in pineapple. *Plant Dis.* **2005**, *89*, 450–456. [CrossRef] [PubMed]

3. Carter, W. The influence of plant nutrition on susceptibility of pineapple plants to mealybug wilt. *Phytopathology* **1945**, *35*, 316–323.

4. Carter, W. (Ed.) The systemic phytotoxemias: Mealybug wilt of pineapple. In *Insects in Relation to Plant Disease*; Interscience: New York, NY, USA, 1962; pp. 238–265.

5. Carter, W. Mealybug wilt of pineapple; a reappraisal. *Ann. N. Y. Acad. Sci.* **1963**, *105*, 741–764. [CrossRef]
6. Larsen, L.D. Diseases of pineapple: Hawaii Sugar Planters' Association. *Pathol. Physiol. Ser. Exp. Stn. Bull.* **1910**, *10*, 1-72.

7. Gambley, C.F.; Steele, V.; Geering, A.D.W.; Thomas, J.E. The genetic diversity of ampeloviruses in Australian pineapples and their association with mealybug wilt disease. *Australas. Plant Pathol.* **2008**, *37*, 95–105. [CrossRef]

8. Sether, D.M.; Hu, J.S. Yield Impact and Spread of Pineapple mealybug wilt associated virus-2 and Mealybug Wilt of Pineapple in Hawaii. *Plant Dis.* **2002**, *86*, 867–874. [CrossRef] [PubMed]

9. Ito, K. Additional immunological evidence supporting the virus nature of mealybug wilt. *Pineapple Res. Inst. News* **1962**, *10*, 158–162.

10. Gunasinghe, U.; German, T. Association of virus-particles with mealybug-wilt of pineapple. *Phytopathology* **1986**, *76*, 1073.

11. Gunasinghe, U.; German, T. Further characterization of virus associated with mealybug-wilt of pineapple. *Phytopathology* **1987**, *77*, 1776.

12. Melzer, M.J.; Karasev, A.V.; Sether, D.M.; Hu, J.S. Nucleotide sequence, genome organization and phylogenetic analysis of pineapple mealybug wilt-associated virus-2. *J. Gen. Virol.* **2001**, *82*, 1–7. [CrossRef]

13. Gunasinghe, U.; German, T. Purification and partial characterization of a virus from pineapple. *Phytopathology* **1989**, *79*, 1337–1341. [CrossRef]

14. Martelli, G.P.; Agranovsky, A.A.; Bar-Joseph, M.B.D.; Boscia, D.; Candresse, T.; Coutts, R.H.A.; Dolja, V.V.; Falk, B.W.; Gonsalves, D.; Jelkmann, W.; et al. ICTV Study Group on closteroviruses and allied viruses. *The family Closteroviridae revised. Arch. Virol.* **2002**, *147*, 2039–2044. [CrossRef] [PubMed]

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

16. Maree, H.J.; Almeida, R.P.P.; Bester, R.; Chooi, K.M.; Cohen, D.; Dolja, V.V.; Fuchs, M.F.; Golino, D.A.; Jooste, A.E.C.; Martelli, G.P.; et al. Grapevine leafroll-associated virus 3. *Front. Microbiol.* **2013**, *4*, 82. [CrossRef] [PubMed]

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

18. Larrea-Sarmiento, A.; Olmedo-Velarde, A.; Wang, X.; Borth, W.; Matsumoto, T.K.; Suzuki, J.Y.; Wall, M.M.; Melzer, M.; Hu, J. A novel ampelovirus associated with mealybug wilt of pineapple (*Ananas comosus*). *Virus Genes* **2021**, *57*, 464–468. [CrossRef]

19. Herrbach, E.; Le Maguet, J.; Hommay, G. CHAPTER 11: Virus Transmission by Mealybugs and Soft Scales (Hemiptera: Coccoidea). In *Vector-Mediated Transmission of Plant Pathogens*; American Phytopathological Society (APS): Saint Paul, MN, USA, 2016; pp. 147–161. [CrossRef]

20. Sether, D.M.; Borth, W.B.; Melzer, M.J.; Hu, J. Spatial and Temporal Incidences of Pineapple mealybug wilt-associated viruses in Pineapple Planting Blocks. *Plant Dis.* **2010**, *94*, 196–200. [CrossRef]

21. Sether, D.M.; Hu, J.S. Closterovirus infection and mealybug exposure are necessary for the development of mealybug wilt of pineapple disease. *Phytopathology* **2002**, *92*, 928–935. [CrossRef]

22. Larrea-Sarmiento, A.; Geering, A.D.W.; Olmedo-Velarde, A.; Wang, X.; Borth, W.; Matsumoto, T.K.; Suzuki, J.Y.; Wall, M.M.; Melzer, M.; Moyle, R.; et al. Genome sequence of pineapple secovirus B, a second sadwavirus reported infecting Ananas comosus. *Arch. Virol.* **2022**, *167*, 2801–2804. [CrossRef]

23. Larrea-Sarmiento, A.; Olmedo-Velarde, A.; Green, J.C.; Al Rwahnih, M.; Wang, X.; Li, Y.H.; Wu, W.; Zhang, J.; Matsumoto, T.K.; Suzuki, J.Y.; et al. Identification and complete genomic sequence of a novel sadwavirus discovered in pineapple (*Ananas comosus*). *Arch. Virol.* **2020**, *165*, 1245–1248. [CrossRef]

24. Larrea-Sarmiento, A.E.; Olmedo-Velarde, A.; Wang, X.; Borth, W.; Domingo, R.; Matsumoto, T.K.; Suzuki, J.Y.; Wall, M.M.; Melzer, M.J.; Hu, J. Genetic Diversity of Viral Populations Associated with Ananas Germplasm and Improvement of Virus Diagnostic Protocols. *Pathogens* **2022**, *11*, 1470. [CrossRef]

25. Massé, D.; Filloux, D.; Candresse, T.; Massart, S.; Marais, A.; Verdin, E.; Cassam, N.; Fernandez, E.; Roumagnac, P.; Teycheney, P.Y.; et al. Identification of a novel vitivirus from pineapple in Reunion Island. *Arch. Virol.* **2022**, *167*, 2355–2357. [CrossRef]

26. Adams, I.P.; Glover, R.H.; Monger, W.A.; Mumford, R.; Jackeviciene, E.; Navalinskiene, M.; Samuitiene, M.; Boonham, N. Next-generation sequencing and metagenomic analysis: A universal diagnostic tool in plant virology. *Mol. Plant Pathol.* **2009**, *10*, 537–545. [CrossRef]

27. Roossinck, M.J. Deep Sequencing for Discovery and Evolutionary Analysis of Plant Viruses. *Virus Res.* **2017**, *239*, 82–86. [CrossRef]

28. Wu, Q.; Ding, S.-W.; Zhang, Y.; Zhu, S. Identification of Viruses and Viroids by Next-Generation Sequencing and Homology-Dependent and Homology-Independent Algorithms. *Annu. Rev. Phytopathol.* **2015**, *53*, 425–444. [CrossRef]

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

30. Marais, A.; Murolo, S.; Faure, C.; Brans, Y.; Larue, C.; Maclot, F.; Massart, S.; Chiumenti, M.; Minafra, A.; Romanazzi, G.; et al. Sixty Years from the First Disease Description, a Novel Badnavirus Associated with Chestnut Mosaic Disease. *Phytopathology* **2021**, *111*, 1051–1058. [CrossRef]

31. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics* **2014**, *30*, 2114–2120. [CrossRef]

32. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Completing bacterial genome assemblies with multiplex MinION sequencing. *Microb. Genom.* **2017**, *3*, e000132. [CrossRef]
33. De Coster, W.; D'Hert, S.; Schultz, D.T.; Cruts, M.; Van Broeckhoven, C. NanoPack: Visualizing and processing long-read sequencing data. *Bioinformatics* **2018**, *34*, 2666-2669. [CrossRef]

34. 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]

35. Rouard, M.D.A.; Droc, G.; Sarah, G.; Dufayard, J.F.; Guignon, V.; Hamelin, C.; Homa, F.; de Lamotte, F.; Larmande, P.; Lariviere, D.; et al. The South Green Bioinformatics platform. In Proceedings of the Plant and Animal Genome, PAG XXII, San Diego, CA, USA, 11-15 January 2014.

36. Masoudi-Nejad, A.; Tonomura, K.; Kawashima, S.; Moriya, Y.; Suzuki, M.; Itoh, M.; Kanehisa, M.; Endo, T.; Goto, S. EGassembler: Online bioinformatics service for large-scale processing, clustering and assembling ESTs and genomic DNA fragments. *Nucleic Acids Res.* **2006**, *34*, W459-W462. [CrossRef]

37. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. *PLoS ONE* **2014**, *9*, e112963. [CrossRef]

38. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **2009**, *25*, 1754-1760. [CrossRef]

39. Li, H. Minimap2: Pairwise alignment for nucleotide sequences. *Bioinformatics* **2018**, *34*, 3094-3100. [CrossRef] [PubMed]

40. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Subgroup, G.P.D.P. The Sequence Alignment/Map format and SAMtools. *Bioinformatics* **2009**, *25*, 2078-2079. [CrossRef]

41. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. *Mol. Biol. Evol.* **2013**, *30*, 772-780. [CrossRef]

42. Martin, D.; Rybicki, E. RDP: Detection of recombination amongst aligned sequences. *Bioinformatics* **2000**, *16*, 562-563. [CrossRef]

43. Padidam, M.; Sawyer, S.; Fauquet, C.M. Possible Emergence of New Geminiviruses by Frequent Recombination. *Virology* **1999**, *265*, 218-225. [CrossRef] [PubMed]

44. Mika, O.; Salminen, J.K.C.; Burke, D.S.; Mccutchan, F.E. Identification of Breakpoints in Intergenotypic Recombinants of HIV Type 1 by Bootscanning. *AIDS Res. Hum. Retroviruses* **1995**, *11*, 1423-1425. [CrossRef]

45. Smith, J.M. Analyzing the mosaic structure of genes. *J. Mol. Evol.* **1992**, *34*, 126-129. [CrossRef] [PubMed]

46. Posada, D.; Crandall, K.A. Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 13757-13762. [CrossRef] [PubMed]

47. Gibbs, M.J.; Armstrong, J.S.; Gibbs, A.J. Sister-Scanning: A Monte Carlo procedure for assessing signals in recombinant sequences. *Bioinformatics* **2000**, *16*, 573-582. [CrossRef] [PubMed]

48. Boni, M.F.; Posada, D.; Feldman, M.W. An Exact Nonparametric Method for Inferring Mosaic Structure in Sequence Triplets. *Genetics* **2007**, *176*, 1035-1047. [CrossRef] [PubMed]

49. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. *Virus Evol.* **2021**, *7*, veeaa087. [CrossRef] [PubMed]

50. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree: Computing large minimum evolution trees with profiles instead of a distance matrix. *Mol. Biol. Evol.* **2009**, *26*, 1641-1650. [CrossRef] [PubMed]

51. Dolja, V.V.; Kreuze, J.F.; Valkonen, J.P.T. Comparative and functional genomics of closteroviruses. *Virus Res.* **2006**, *117*, 38-51. [CrossRef] [PubMed]

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

53. Massé, D.; Cassam, N.; Hostachy, B.; Iskra-Caruana, M.-L.; Darnaudery, M.; Lefeuvre, P.; Lett, J.-M. First Report of Three Pineapple Mealybug Wilt-Associated Viruses in Queen Victoria Pineapples in Reunion Island. *Plant Dis.* **2021**, *105*, 715. [CrossRef]

54. Moreno, I.; Rodríguez-Arévalo, K.A.; Tarazona-Velásquez, R.; Kondo, T. Occurrence and distribution of pineapple mealybug wilt-associated viruses (PMWaVs) in MD2 pineapple fields in the Valle del Cauca Department, Colombia. *Trop. Plant Pathol.* **2023**, *48*, 217-225. [CrossRef]

55. Bradamante, G.; Mittelsten Scheid, O.; Incarbone, M. Under siege: Virus control in plant meristems and progeny. *Plant Cell* **2021**, *33*, 2523-2537. [CrossRef] [PubMed]

56. Paull, R.E.; Duarte, O. *Tropical Fruits Volume I. Crop Production Science, Horticulture*; CABI Publishing: Wallingford, UK, 2011.

57. Zanella, C.M.; Janke, A.; Palma-Silva, C.; Kaltchuk-Santos, E.; Pinheiro, F.G.; Paggi, G.M.; Soares, L.E.S.; Goetze, M.; Büttow, M.V.; Bered, F. Genetics, evolution and conservation of Bromeliaceae. *Genet. Mol. Biol.* **2012**, *35*, 1020-1026. [CrossRef] [PubMed]

58. Joy, P.; Anjana, R. *Evolution of Pineapple*; Pineapple Research Station: Muvattupuzha, India, 2015; pp. 670-686.

59. Malmstrom, C.M.; Shu, R.; Linton, E.W.; Newton, L.A.; Cook, M.A. Barley yellow dwarf viruses (BYDVs) preserved in herbarium specimens illuminate historical disease ecology of invasive and native grasses. *J. Ecol.* **2007**, *95*, 1153-1166. [CrossRef]

60. Rwahnih, M.A.; Rowhani, A.; Golino, D.A.; Islas, C.M.; Preece, J.E.; Sudarshana, M.R. Detection and genetic diversity of Grapevine red blotch-associated virus isolates in table grape accessions in the National Clonal Germplasm Repository in California. *Can. J. Plant Pathol.* **2015**, *37*, 130-135. [CrossRef]
61. Rieux, A.; Campos, P.; Duvermy, A.; Scussel, S.; Martin, D.; Gaudeul, M.; Lefeuvre, P.; Becker, N.; Lett, J.-M. Contribution of historical herbarium small RNAs to the reconstruction of a cassava mosaic geminivirus evolutionary history. *Sci. Rep.* **2021**, *11*, 21280. [CrossRef] [PubMed]

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

63. Sether, D.M.; Melzer, M.J.; Borth, W.B.; Hu, J.S. Genome organization and phylogenetic relationship of Pineapple mealybug wilt associated virus-3 with family Closteroviridae members. *Virus Genes* **2009**, *38*, 414–420. [CrossRef] [PubMed]

64. Massart, S.; Adams, I.; Al Rwahnih, M.; Baeyen, S.; Bilodeau, G.J.; Blouin, A.G.; Boonham, N.; Candresse, T.; Chandellier, A.; De Jonghe, K. Guidelines for the reliable use of high throughput sequencing technologies to detect plant pathogens and pests. *Peer Community J.* **2022**, *2*, e62. [CrossRef]

65. Fontdevila Pareta, N.; Khalili, M.; Maachi, A.; Rivarez, M.P.S.; Rollin, J.; Salavert, F.; Temple, C.; Aranda, M.A.; Boonham, N.; Botermans, M.; et al. Managing the deluge of newly discovered plant viruses and viroids: An optimized scientific and regulatory framework for their characterization and risk analysis. *Front. Microbiol.* **2023**, *14*, 1181562. [CrossRef]

66. Gambley, C.F.; Geering, A.D.W.; Thomas, J.E. Development of an immunomagnetic capture-reverse transcriptase-PCR assay for three pineapple ampeloviruses. *J. Virol. Methods* **2009**, *155*, 187–192. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.