# Pineapple mycobiome related to fruitlet core rot occurrence and the influence of fungal species dispersion patterns

Manon Vignassa, Jean-Christophe Meile, Frédéric Chiroleu, Christian Soria, Charlène Leneveu-Jenvrin, Sabine Schorr-Galindo, Marc Chillet

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

## To cite this version:

Manon Vignassa, Jean-Christophe Meile, Frédéric Chiroleu, Christian Soria, Charlène Leneveu-Jenvrin, et al.. Pineapple mycobiome related to fruitlet core rot occurrence and the influence of fungal species dispersion patterns. Journal of Fungi, 2021, 7 (3), pp.175. 10.3390/jof7030175. hal-05177322

HAL Id: hal-05177322

https://hal.science/hal-05177322v1

Submitted on 22 Jul 2025

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Pineapple Mycobiome Related to Fruitlet Core Rot Occurrence and the Influence of Fungal Species Dispersion Patterns

Manon Vignassa<sup>1,2,*</sup>, Jean-Christophe Meile<sup>1,2</sup>, Frédéric Chiroleu<sup>3</sup>, Christian Soria<sup>1,2</sup>, Charlène Leneveu-Jenvrin<sup>4</sup>, Sabine Schorr-Galindo<sup>2</sup> and Marc Chillet<sup>1,2</sup>

**Abstract:** Fruitlet Core Rot (FCR) is a fungal disease that negatively impacts the quality of pineapple, in particular the 'Queen Victoria' cultivar. The main FCR causal agent has been identified as *Fusarium ananatum*. This study focused on the correlation between FCR disease occurrence, fungal diversity, and environmental factors. FCR incidence and fungal species repartition patterns were spatially contextualized with specific surrounding parameters of the experimental plots. The myco-biome composition of healthy and diseased fruitlets was compared in order to search for potential fungal markers. A total of 240 pineapple fruits were sampled, and 344 fungal isolates were identified as belonging to 49 species among 17 genera. FCR symptom distribution revealed a significant gradient that correlated to that of the most abundant fungal species. The association of wind direction and the position of proximal cultivated crops sharing pathogens constituted an elevated risk of FCR incidence. Five highly represented species were assayed by Koch's postulates, and their pathogenicity was confirmed. These novel pathogens belonging to *Fusarium fujikuroi* and *Talaromyces purpureo-genus* species complexes were identified, unravelling the complexity of the FCR pathosystem and the difficulty of apprehending the pathogenesis over the last several decades. This study revealed that FCR is an airborne disease characterized by a multi-partite pathosystem.

**Keywords:** *Fusarium*; pathogenic fungi; pineapple disease; *Talaromyces*

## 1. Introduction

*Ananas comosus* has become a highly appreciated fruit and is considered the best tropical fruit in terms of world exportation volumes, with 3.2 million tons in 2019 [1]. Pineapple production is affected by numerous pests and diseases, including Fruitlet Core Rot (FCR). The infection leads to brown discoloration of the flesh, impacting fruit acceptability by the consumers. FCR occurrence in the field of tropical and subtropical areas impacts pineapple post-harvest quality and both local and export markets. Current strategies focus on the creation of new cultivars trying to combine desired sensorial qualities and an optimal level of resistance to pathogens, as with the 'MD-2' cultivar [2]. Nevertheless, the 'Queen Victoria' cultivar still remains a favorite, especially for its flavor, in the Indian Ocean, as well as in Europe and some parts of Asia, despite its high susceptibility to a broad spectrum of pathogens [3]. Disease incidence and severity have increased in the

**Citation:** Vignassa, M.; Meile, J.-C.; Chiroleu, F.; Soria, C.; Leneveu-Jenvrin, C.; Schorr-Galindo, S.; Chillet, M. Pineapple Mycobiome Related to Fruitlet Core Rot Occurrence and the Influence of Fungal Species Dispersion Patterns. *J. Fungi* 2021, 7, 175. https://doi.org/10.3390/jof7030175

Academic Editor: Lei Cai

Received: 28 January 2021

Accepted: 24 February 2021

Published: 28 February 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 (http://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD), UMR Qua-lisud, F-97410 Saint-Pierre, F-97410 Réunion, France; meile@cirad.fr (J.-C.M.); christian.soria@cirad.fr (C.S.); marc.chillet@cirad.fr (M.C.)</footnote>

<footnote n="2">² Unité Mixte de Recherche Qualisud, Université de Montpellier, Avignon Université, Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Institut Agro, Institut de Recherche pour le Développement, Université de La Réunion, Montpellier, France; sabine.galindo@umontpellier.fr</footnote>

<footnote n="3">³ Centre de coopération internationale en recherche agronomique pour le développement CIRAD, UMR PVBMT, 97410 Saint-Pierre, F-97410 La Réunion, France; frederic.chiroleu@cirad.fr</footnote>

<footnote n="4">⁴ Unité Mixte de Recherche Qualisud, Université de La Réunion, Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Université de Montpellier, Institut Agro, Avignon Université, Sainte-Clotilde, France; charlene.leneveujenvrin@yahoo.fr</footnote>

<footnote>* Correspondence: manon.vignassa@cirad.fr; Tel.: +262-2-62-49-28-62</footnote>
producing regions, resulting in significant fruit quality depreciation and thus economic concerns [4]. The FCR pathosystem was originally considered as the relation between *Anas comosus* var. *comosus* and the fungus *Fusarium verticillioides* (syn. *Fusarium moniliforme*) associated with *Penicillium funiculosum* [5–7]. However, the evolution of molecular tools led to the identification of *Fusarium ananatum* and *Talaromyces funiculosus* as causal agents of FCR in China, South Africa, Mauritius, and Reunion Island [4,8–11]. Contamination by soil-borne pathogens (such as *P. funiculosum*) occurs at pre-flowering stages, and the fungus remains latent in a blossom cup. During fruit development, causal agents spread into septal nectaries, but propagation is restricted by the lignification of young plant tissues. The accumulation of sugar following fruit ripening finally enables the fungus to counter-balance the plant defenses and then colonize host tissues [12]. These events trigger the browning of the flesh, whose expansion is limited to the fruitlet, without reaching the core or the skin of the pineapple fruit. Those symptoms are defined as ‘black spots’, thus constituting a typical marker of FCR-infected fruitlets. Nevertheless, the disease incidence could not be directly evaluated during pineapple cultivation or after harvesting due to the absence of external symptoms. The complexity of the FCR pathosystem relies on the presence of both diseased and healthy fruitlets within the same pineapple fruit, thus suggesting that pathogenesis could be driven by a specific microbiota. Thus far, the critical biotic parameters that mediate pineapple immunity or susceptibility remain unknown.

Plant diseases are closely related to numerous biotic and abiotic factors. Fungal species can be found in association with various plant species leading to a large-scale distribution and high levels of diversity according to their habitat [13–15]. Considering the wide range of pineapple crop systems and abiotic environments, fungi identified as belonging to *Trichoderma* spp. have been isolated in higher relative abundances in leaf and leaf sheath tissues when compared to roots and rhizospheric areas. By contrast, these plant sections are characterized by a microbiota, at least 70% of which is composed of bacterial species [16]. Nonetheless, there has been no study about fungal taxa richness and diversity associated with pineapple. Among abiotic factors, previous work has shown a link between FCR incidence and specific climatic conditions, which allowed for the elucidation of a predictive model based on a pluviothermic index. Predictions showed that, in subtropical regions (such as Reunion Island), the combination of a high altitude, high rainfall and cool temperatures (16 °C to 27 °C) led to an elevated risk of harvesting fruits with FCR symptoms [17,18]. In spite of this, the impact of environmental factors on the fruit mycobiome and pathogen dispersion patterns remains a ‘black box’ in the fate of this plant–pathogen interaction. Further studies on inoculum availability have reported the prevalence of disease severity related to the spatial distribution pattern on mango and banana crops [19,20]. Current ecological concerns have induced the development of alternative strategies for managing the spreading of fungal pathogens. To this end, multi-scale assays are conducted in order to fit with environmental changes and their impact on microbial epidemics [21,22]. The few studies dealing with short-distance dispersal patterns have mostly focused on vineyard [23], mango [19], and apple orchards [24]. Those data demonstrated the impact of soil and row management on conidia dispersion dynamics. One of the relevant aspects of estimating the infection pressure is the relation between the pineapple plot and its surrounding environment, notably the concomitant crops. Indeed, fungal isolates belonging to *Fusarium* genus are able to infect a broad spectrum of hosts such as fruits, small grain cereals, and horticultural plants [25–27].

Based on these assumptions, the main goal of the present study was to characterize the fungal flora of healthy and naturally infected pineapple fruitlets. From the data obtained, fungal species distribution was linked to the proximal environment of the pineapple plots. We also investigated the potential co-occurrence between fungal species as potential key factors for disease development. Thanks to a context-specific approach, we report here for the first time the contribution of five new pathogenic fungal species belonging to *Fusarium* and *Talaromyces* genera in FCR incidence in pineapple cv. 'Queen Victoria'.
## 2. Materials and Methods

### 2.1. Plant Material and Agricultural Practices

The experiment was conducted on the 'Queen Victoria' pineapple cultivar in accordance with the conventional cultural practices locally applied [28]. The plots were located in the southwest of Reunion Island at the Experimental Research Station of CIRAD Bassin-Plat, Saint-Pierre, located at an altitude of 150 meters above sea level (21°19'21''S, 55°29'26''E). The surface was divided into 2 plots with perpendicular row directions. Thereby, row orientations were North-West/South-East (N-W/S-E) and North-East/South-West (N-E/S-W) for Plot 1 and Plot 2, respectively (Figure 1). Each of the 2 experimental plots were subdivided into 3 sections (2 rows per section) along an N-E/S-W axis corresponding to rows adjacent to the mango orchard (1), central rows (2), and rows proximal to jackfruit trees (3). Each section was subdivided into 4 quadrats of 4 m long (n = 24). Meteorological data from flower induction treatment (FIT) in April 2018 to harvest in November 2018 were considered [29]. The wind data collection was performed by the Météo-France weather station located in Pierrefonds [30], Saint-Pierre, approximatively 6 km away from the plots (Table 1). A total of 240 fruits were harvested at the C1 stage, corresponding to 1/4 yellow fruit according to the shell's color. The ripening was completed at 19 °C in a cold room until the C4 stage, when pineapples became entirely yellow.

### 2.2. FCR Symptom Occurrence and Fruitlets Sampling

Ten fruits per quadrat were collected. Each fruit was cut, and the number of black spots (infected fruitlets) was counted and recorded. FCR occurrence per fruit was then calculated for each section of each of the 2 plots.

Fruitlets sampling was performed on 96 fruits randomly chosen (4 fruits per quadrat). For each pineapple fruit, 1 healthy-looking fruitlet (HF) and 1 naturally infected fruitlet (IF) were sampled with sterile equipment for a total of 192 fruitlets.

### 2.3. Isolation of Cultivatable Fungal Flora from Fruitlet and Soil Samples

#### 2.3.1. Fruitlet Tissues

The pineapple fruitlets collected (Section 2.2.) were deposited on Sabouraud glucose agar Petri dish (Biokar diagnostic, Solabia, Allonne, France) supplemented with 100 mg L⁻¹ chloramphenicol at 27 °C in the dark. Fungal colony growth and isolation were processed daily following similar conditions until 5 days of incubation. Pure strains isolated were observed for morphological (colony aspect and color) and microscopic (macroconidia and microconidia) characterizations (Zeiss, Germany, model Axiostar Plus, magnification x100). These observations enabled the taxonomic assignation for each isolate at the genus level. Consequently, a conidia solution was prepared by adding 2 mL of SPW (Saline Peptone Water, Condalab, Torrejón de Ardoz, Madrid, Spain) to the Petri dish and gently scratching the mycelium surface with a sterile spreader. From the conidia solution, 1 mL was transferred into sterile cryovials with the addition of an equal volume of 40% glycerol (Sigma-Aldrich, Darmstadt, Germany) and stored at -80 °C. Another volume of 500 µL was recovered for the molecular identification procedure and stored at -20 °C.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Environmental context of pineapple plots located in Reunion Island. Plot 1: North-West/South-East-directed rows; Plot 2: South-West/North-East-directed rows. (1) Plot sections proximal to the mango orchard, (2) central section of the plots, and (3) plot sections near the jackfruit trees. Quadrats are shown as yellow rectangles. White asterisks indicate coring points.</figcaption></figure>

<table><caption>**Table 1.** Meteorological factors recorded from FIT (Flower Induction Treatment) in April 2018 to harvest in November 2018. Source of satellite image: https://www.geoportail.gouv.fr/.</caption><thead><tr><th></th><th>Min Temp (°C)</th><th>Average Temp (°C)</th><th>Max Temp (°C)</th><th>Average Rainfall (mm)</th><th>Number of day(s) with Rain</th><th>Wind Direction</th><th>Max Gust Speed (km h<sup>-1</sup>)</th></tr></thead><tbody><tr><td>April</td><td>20.12</td><td>23.19</td><td>27.85</td><td>15.78 *</td><td>10</td><td>south</td><td>116*</td></tr><tr><td>May</td><td>18.42</td><td>22.20</td><td>27.48</td><td>0.66</td><td>5</td><td>east-south-east</td><td>85</td></tr><tr><td>June</td><td>17.10</td><td>20.68</td><td>25.72</td><td>2.47</td><td>8</td><td>south-east</td><td>73</td></tr><tr><td>July</td><td>15.53</td><td>19.50</td><td>24.64</td><td>4.46</td><td>12</td><td>south-east</td><td>80</td></tr><tr><td>August</td><td>17.07</td><td>21.10</td><td>26.60</td><td>0.00</td><td>0</td><td>south-east</td><td>75</td></tr><tr><td>September</td><td>16.64</td><td>21.15</td><td>27.60</td><td>0.15</td><td>1</td><td>south-east</td><td>80</td></tr><tr><td>October</td><td>17.65</td><td>21.75</td><td>27.79</td><td>3.85</td><td>8</td><td>east-south-east</td><td>88</td></tr><tr><td>November</td><td>18.98</td><td>23.49</td><td>28.66</td><td>0.03</td><td>1</td><td>south-east</td><td>69</td></tr></tbody></table>

### 2.3.2. Soil

Soil sampling was conducted at 3 points 10 m apart on the 2 opposite junctions between the pineapple field and the surrounding crops: The jackfruit trees (North-East) and the mango orchard (South-West) (Figure 1). For each coring point (white asterisks), 1 g of a 5 cm soil layer was suspended in a sterile tube containing 9 mL of sterile physiological saline water with 9 g of NaCl supplemented with 1% Tween® 80 (Merck, Darmstadt, Germany). Soil suspensions were then placed into a shaking incubator for 30 min at 25 °C and 200 rpm. Two technical repetitions per coring point were performed. For microbial assays, serial decimal dilutions were conducted to 10⁻⁶ in physiological saline water supplemented with 1% Tween® 80. Fungi were isolated on PDA (Potato Dextrose Agar, BD Difco™, Le Pont-de-Claix, France) supplemented with cycloheximide (10 mg L⁻¹, Merck, Darmstadt, Germany) and on Sabouraud glucose agar supplemented with Chloramphenicol 100 mg L⁻¹, after incubation for 5 days at 27 °C and 4 days at 30 °C, respectively. Fungal

<footnote>(*) The tropical cyclone Fakir passed near the eastern coast of Reunion Island from April 24th to April 26th 2018.</footnote>
colonies were removed with a sterile scalpel and transferred into 1.5 mL sterile tubes with 500 µL of SPW. Isolates were finally stored at -20 °C prior to molecular procedures.

## 2.4. Molecular characterization of cultivatable fungal flora from fruitlets and soil samples

According to the morphological analysis, genus-specific DNA primers were used to infer phylogenetic relationships between closely related genotypes of strains isolated from fruitlets.

PCR amplification of fungal DNA was performed on conidia solutions with no DNA extraction procedure. For *Fusarium* sp. genus, the EF1α (TEF-1α) reference gene region was amplified by PCR. For *Aspergillus* sp., *Penicillium* sp., and *Talaromyces* sp. genus, the β-tubulin gene region was targeted for PCR amplification. Afterward, when no amplification could be obtained with those specific primers, an amplification of the Internal Transcribed Spacer (ITS) region was conducted with the ITS1F/ITS4 primer pair. The list of the primers used in this study is given in Table 2. The identification of soil fungal isolates was performed by PCR amplification of the Internal Transcribed Spacer (ITS) region.

PCR amplification reactions were performed in a final volume of 50 µL containing 0.3 µM of each primer, all the deoxyribonucleotide triphosphates (dNTPs) at 0.2 mM, 2 mM MgCl₂, 10 µL of 5x GoTaq Flexi buffer, 1.25 U GoTaq® G2 Flexi DNA polymerase (Promega, Charbonnières-les-Bains, France), and 2 µL of conidia solution. Conditions for the PCR amplification of the TEF-1α region were an initial denaturation at 95 °C for 3 min, 33 cycles at 95 °C for 30 s, 55 °C for 30 s, 72 °C for 1 min, and a final extension stage of 5 min at 72 °C. Similarly, PCR amplification reactions for the β-tubulin region were established as follows: An initial denaturation at 95 °C for 2 min, 30 cycles at 95 °C for 30 s, 62 °C for 30 s, and 72 °C for 30 s, and a final extension at 72 °C for 5 min. PCR amplification reactions for the ITS region were carried out as follows: An initial denaturation at 95 °C for 2 min, 40 cycles at 95 °C for 15 s, 53 °C for 30 s, and 72 °C for 45 s, and a final extension at 72 °C for 5 min. The PCR reactions were performed in a thermocycler (Veriti, Applied Biosystems, United Kingdom). PCR products were analyzed with the Qiaxcel® Advanced System (Qiagen, Hilden, Germany) using size markers 250 bp - 4 kb.

The PCR products were sent to Macrogen (Amsterdam, The Netherlands) for purification and sequencing. The DNA sequences obtained were aligned with SnapGene v5.0 software, and identification was performed using a BLASTn similarity search. Sequences having a percentage of identity of at least 98% and 95% for fruit and soil isolates, respectively, and those with the lowest E-values were considered as belonging to the same species.

<table><caption>**Table 2.** PCR primers sequences (forward and reverse) used for fungal strains identification.</caption><thead><tr><th>Sequence ID of Primers Pair</th><th>Target Locus</th><th>Sequence Forward (5' 3')</th><th>Sequence Reverse (5' 3')</th><th>Product Length (bp)</th><th>References</th></tr></thead><tbody><tr><td>ef1 / ef2</td><td>Translation Elongation Factor-1α</td><td>ATGGGTAAGGAAGAACAGAC</td><td>GGAAGTACCAGTGATCATGTT</td><td>380-680</td><td>[31-33]</td></tr><tr><td>Bt2a / Bt2b</td><td>β-tubulin</td><td>GGTAACCAAATCGGTGCTGCTTTTC</td><td>ACCCTCAGTGTAGTGACCCTTGGC</td><td>250-500</td><td>[34]</td></tr><tr><td>ITS1F / ITS2</td><td rowspan="2">Internal Transcribed Spacer 1</td><td rowspan="2">CTTGGTCATTTAGAGGAAGTAA</td><td>GCTGCGTTCTTCATCGATGC</td><td rowspan="2">145-695</td><td rowspan="2"></td></tr><tr><td>ITS1F / ITS4</td><td>TCCTCCGCTTATTGATATGC</td></tr><tr><td>GC-ITS1F / ITS4</td><td>Entire Internal Transcribed Spacer</td><td>CGCCCGCCGCGCGCGCGCGCGCGGGGGGACGGGGGGCTTGGTCATTTAGAGGAAGTAA</td><td>TCCTCCGCTTATTGATATGC</td><td>600-800</td><td>[35,36]</td></tr></tbody></table>
## 2.5. Koch's *Postulates*

### 2.5.1. Controlled Inoculations and Plant Material

Among fungal strains identified on infected fruitlets, 5 highly represented species were selected for pathogenicity assays. For each of the 5 species, 2–3 isolates were randomly tested according to the following correspondence: *Fusarium oxysporum* strains BP369 and BP460; *Fusarium proliferatum* strains BP114, BP429, and BP436; *Fusarium sacchari* strains BP138 and BP575; *Talaromyces stollii* strains BP054, BP185, and BP462; *Talaromyces amestolkiae* strains BP002, BP257, and BP605. First, strains of interest were grown on PDA for 7 days at 27 °C in the dark. For each isolate, an inoculum solution was prepared with sterile water as described in Section 2.3.1, and a final concentration was normalized at 10⁵ conidia per mL. Inoculations were performed on cv. 'Queen Victoria' pineapple fruits were harvested at the C1 stage and grown according to locally recommended conventional agricultural practices. For each of the 13 tested strains, 3 fruits were inoculated. For each fruit, 3 fruitlets were inoculated by injecting 25 µL of inoculum in a blossom cup from the upper, median, or basal fruits parts. Inoculations were also conducted with *Fusarium ananatum* strain BP383 as a positive control and with sterile water as a negative control (H₂O), and a 2nd class of negative control fruits (Mock) was not inoculated for comparison to the initial fungal load of fruitlets. Fruits were then incubated for 7 days at 19 °C in a cold room. Thus, 16 different conditions (13 tested strains + 1 positive and 2 negative controls) were assayed on a total of 144 fruitlets from 48 pineapple fruits.

### 2.5.2. Fruit Sampling

After incubation, inoculated fruitlets were sampled. The 3 infected or control fruitlets of each fruit were pooled and placed in a sealed sterile lab blender bag and mixed with 10 mL of SPW in a stomacher for 2 min at maximal speed. Prior to storage at -80 °C, 5 mL of lysate were recovered and supplemented with an equal volume of 40% glycerol.

### 2.5.3. DNA Extraction

For each fruit sample and each pure fungal strain used for inoculation, 2 mL of lysate were collected in 2 mL sterile tubes and centrifuged at 14,000 g for 2 min. DNA extractions were performed on the biomass pellets with the FastDNA SPIN kit and the FastPrep-24 Instrument (MP Biomedicals®, Ilkirch, France) using Lysing Matrix A and Lysis Buffer CLS-Y in accordance with the manufacturer's instructions.

### 2.5.4. PCR-Denaturing Gradient Gel Electrophoresis (DGGE)

PCR amplifications were performed using the GC-ITS1F and ITS2 DNA primers (Table 2). A 40 bp GC-clamp was added to the 5' end of the ITS1F primer in order to ensure that the DNA fragment remains partially double-stranded and that the region screened was in the lowest melting domain [37]. This method was limited by the fragment length, which should not exceed 500 bp. In spite of this, a large fragment length resulted in a poor gel resolution complicating the recovery of DNA bands. For these reasons, PCR-DGGE was performed on the ITS1 region (~200 bp), conferring a good resolution but discrimination of fungal strains only at genus level. PCR reaction was performed in a final volume of 50 µL containing 0.6 µM of each primer, all the deoxyribonucleotide triphosphates (dNTPs) at 200 µM, 2 mM of MgCl₂, 10 µL of 5× GoTaq Flexi buffer, 1.25 U of GoTaq® G2 Flexi DNA polymerase, and 1 µL of extracted DNA. PCR amplification reactions were carried out as follows: An initial denaturation at 95 °C for 2 min, 40 cycles at 95 °C for 15 s, 60 °C for 30 s, and 72 °C for 20 s, and a final extension at 72 °C for 5 min. The PCR reactions were performed in a thermocycler (Veriti, Applied Biosystems, United Kingdom). PCR products were then analyzed as previously described.

The PCR products were separated by DGGE using a Cleaver Scientific system (Cleaver Scientific, Rugby, Warwickshire, United Kingdom). Briefly, 30 µL of PCR amplicons were loaded onto 8% (*w/v*) polyacrylamide gels (acrylamide: N,N-methylene bis-
acrylamide, 37.5:1, Sigma-Aldrich, Darmstadt, Germany) in a 1× TAE buffer (40 mM Tris-HCl pH 7.4, 20 mM sodium acetate, 1.0 mM Na₂-EDTA). Electrophoresis was performed at 60 °C using a denaturing gradient ranging from 20% to 60% (100% denaturant corresponding to 7 M urea and 40% *v/v* formamide, Carlo Erba Reagents, Val-de-Reuil, France). The gels were run at 20 V for 10 min and then at 75 V for 16 h. After electrophoresis, the gels were stained for 1 h with ethidium bromide solution (50 µg mL⁻¹ in 1× TAE), rinsed for 1 h in distilled water, and then visualized with a UV transilluminator (Bio-Rad®, Marnes-la-Coquette, France). Bands of interest were excised from the DGGE gels with a sterile scalpel as previously described [38]. Briefly, DNA of each band was then eluted in 100 µL of TE buffer (10 mM TrisHCl; 1 mM EDTA; pH 7.4, Sigma-Aldrich, Darmstadt, Germany) at 4 °C overnight. DNA was precipitated by adding 10 µL of sodium acetate (3 M, pH 5), 1 µL of glycogen (Molecular Grade, Roche Diagnostics, Meylan, France), and 300 µL of 100% ethanol and centrifuged at 15,000 g for 30 min at 4 °C. The supernatant was discarded, DNA pellets were washed with 500 µL of 70% ethanol, and, after 5 min of centrifugation, the DNA pellets were air-dried for 1 h. Finally, the DNA was re-suspended in 20 µL of TE buffer (10 mM TrisHCl; 1 mM EDTA; pH 7.4) and stored at -20 °C. PCR amplification reaction was performed on purified DNA with the ITS1F/ITS2 primer set. Reactions were carried out as follows: An initial denaturation at 95 °C for 2 min, 40 cycles at 95 °C for 15 s, 53 °C for 30 s, and 72 °C for 45 s, and a final extension at 72 °C for 5 min. PCR products were analyzed and sequenced as previously described in Section 2.4.

## 2.6. Computational Analysis

### 2.6.1. Phylogenetic Analysis

The determination of phylogenetic relationships between fungal strains was conducted by considering each genomic region. Multiple alignments of the nucleotide sequences were performed with the MEGAX computer program and the MUSCLE automatic alignments method [39,40]. Subsequently, the jModelTest v2.1.7 program was used to select the nucleotide substitution model for each sequence set [41]. The maximum-Likelihood analysis was inferred with the PhyML v3.0 program implemented with 'ape' R package according to the best fitting model [42,43]. Branch supports were tested with 100 bootstrap replications. Trees were then visualized and edited with the the FigTree v1.4.4 program.

### 2.6.2. DGGE Band Pattern

Electrophoretic profiles were analyzed with CLIQS 1D Pro software (TotalLab, Newcastle-Upon-Tyne, United Kingdom). In the DGGE gel, for each sample (lane), the DNA band presence and its relative intensity were recorded. The Unweighted Pair-Group Method using Arithmetic Average (UPGMA) and Pearson coefficient correlation were used to build the dendrogram.

### 2.6.3. Statistical Analysis

Data were analyzed with the R statistical language v4.0.3 (R Core Team, 2020). The effect of row direction and proximal crops on the mean number of infected-fruitlets per fruit was tested with a Deviance test on a generalized linear model with Poisson distribution, and a Tukey's multiple comparison test. Principal Component Analysis (PCA) was computed with 'FactoMineR' on the centered-fungal species abundances with no scaling [44]. The 6 combinations between proximal crop and plot were considered as individuals. Results of PCA were then visualized with the 'factoextra' package [45]. Correlation between variables and each principal component of the factorial plane was tested at *p* <0.05. Hierarchical clustering of species abundances in healthy and infected tissues was performed according to the Raup-Crick dissimilarity method from the 'vegan' package [46]. Species co-occurrence was computed with the 'cooccur' package based on a probabilistic
model of calculations between expected and observed pairwise frequencies [47,48]. Correlation between species co-occurrence and fruitlet phenotype was assessed by a Mantel test on distance matrixes of each dataset.

## 3. Results

### 3.1. Fungal Flora of Healthy and Naturally Infected Fruitlets

The exploration of the mycobiome resulted in the identification of 170 isolates from healthy-looking fruitlets (HF), and 174 isolates from infected samples (IF) exhibiting black spots. In all, the 344 isolates belonged to 17 genera and represented a total of 49 different species. The relative abundance of the isolated fungal species identified five highly represented genera that corresponded to *Fusarium*, *Talaromyces*, *Aspergillus*, *Phialemoniopsis*, and *Trichoderma*, together supporting 90% (153 isolates) and 94% (164 isolates) of the cultivatable fungal flora of HF and IF, respectively (Table 3). The remaining isolates belonged to 12 genera (*Pestalotiopsis*, *Lasiodiplodia*, *Penicillium*, *Clonostachys*, *Diaporthe*, *Epicoccum*, *Bionectria*, *Cosmospora*, *Curvularia*, *Davidiella*, *Glomus*, and *Rhizopus*) with low relative abundances ranging from 0.58% to 2.35%. Although these were sparse isolates, *Lasiodiplodia*, *Bionectria*, *Cosmospora*, *Davidiella*, and *Glomus* were exclusively recorded in HF. Contrariwise, isolates corresponding to *Curvularia* and *Rhizopus* could only be observed in IF. The contextualization of predominant species among studied conditions revealed a decrease in relative abundances for *F. proliferatum* and *T. amestolkiae*, from 16.5% and 7.6% in HF to 12.6% and 5.2% in IF samples, respectively. By contrast, the other highly represented species were characterized by an increase in their relative abundances in the IF samples. This trend was notably observed for *T. stollii* and *P. curvata*, which extended from 4.7% and 3.5% in HF to 9.2% and 6.9% in IF, respectively (Figure 2). Particularly, only one *Trichoderma* species (*T. paraviridescens*) was isolated, showing a low relative abundance (0.6%) in HF, while this genus was supported by five species in IF (*T. asperellum*, *T. erinaceum*, *T. harzianum*, *T. paraviridescens*, and *T. trixiae*) contributing to 4.6% of the identified fungal flora. Although *F. ananatum* was the most frequently described FCR-associated pathogen [10,11,49], its occurrence and relative abundances were similar in HF and IF, with 11.2% and 10.3%, respectively. Interestingly, *F. proliferatum*, *F. ananatum*, *F. fujikuroi*, *F. circinatum*, *F. oxysporum*, *F. sacchari*, and *F. verticillioides* all belonged to the *Fusarium fujikuroi* species complex (FFSC) [49], which represented 43.5% and 37.9% of HF and IF mycobiomes, respectively. Similarly, the *Talaromyces purpureogenus* species complex was also highly represented following the identification of *T. stollii*, *T. amestolkiae*, and *T. purpureogenus*, together gathering 15.3% of isolates in HF and 17.2% in IF. Thus, those two species complexes together contributed to 58.8% and 55.2% of the cultivatable fungal flora of HF and IF, respectively.

<table><caption>**Table 3.** Identification of fungal species and their respective occurrence in healthy and naturally infected fruitlets.</caption><thead><tr><th rowspan="2">Genus</th><th rowspan="2">Species</th><th colspan="2">Occurrence</th><th rowspan="2">Top Match GenBank<br>Accession Number *<br>(Frequency in data)</th></tr><tr><th>Healthy Fruitlet (n<br>= 96)</th><th>Infected Fruitlet<br>(n = 96)</th></tr></thead><tbody><tr><td><strong>Fusarium</strong></td><td><em>proliferatum</em></td><td>28</td><td>22</td><td>KF993985.1 (12)</td></tr><tr><td>Fusarium</td><td>ananatum</td><td>19</td><td>18</td><td>MT010996.1 (36)</td></tr><tr><td>Talaromyces</td><td>stollii</td><td>8</td><td>16</td><td>JX315634.1 (24)</td></tr><tr><td>Fusarium</td><td>verticillioides</td><td>12</td><td>12</td><td>MT594370.1 (10)</td></tr><tr><td>Talaromyces</td><td>amestolkiae</td><td>13</td><td>9</td><td>KJ413360.1 (13)</td></tr><tr><td>Fusarium</td><td>equiseti</td><td>10</td><td>11</td><td>MN589630.1 (4)</td></tr><tr><td>Phialemoniopsis</td><td>curvata</td><td>6</td><td>12</td><td>AB278180.1 (15)</td></tr><tr><td>Fusarium</td><td>oxysporum</td><td>9</td><td>9</td><td>CP053267.1 (4)</td></tr><tr><td>Fusarium</td><td>sp.</td><td>10</td><td>5</td><td>JF740861.1 (5)</td></tr><tr><td>Fusarium</td><td>solani</td><td>5</td><td>8</td><td>MT594367.1 (3) / MK968891.1 (3)</td></tr><tr><td>Talaromyces</td><td>purpureogenus</td><td>5</td><td>5</td><td>KJ528885.1 (3) / MF476006.1 (3)</td></tr><tr><td>Aspergillus</td><td>flavus</td><td>2</td><td>6</td><td>MN955851.1 (3)</td></tr></tbody></table>
<table><tr><td>Fusarium</td><td>incarnatum</td><td>3</td><td>5</td><td>MK752398.1 (2) / MN882829.1 (2)</td></tr><tr><td>Fusarium</td><td>fujikuroi</td><td>6</td><td>1</td><td>MF281281.2 (4)</td></tr><tr><td>Aspergillus</td><td>niger</td><td>3</td><td>4</td><td>KY357318.1 (2) / MN788116.1 (2)</td></tr><tr><td>Fusarium</td><td>chlamydosporum</td><td>3</td><td>1</td><td>KJ125830.1 (3)</td></tr><tr><td>Talaromyces</td><td>funiculosus</td><td>2</td><td>2</td><td>AB893941.1 (3)</td></tr><tr><td>Penicillium</td><td>sp.</td><td>3</td><td>1</td><td>EU330619.1 (2)</td></tr><tr><td>Pestalotiopsis</td><td>trachicarpicola</td><td>2</td><td>2</td><td>MN295594.1 (4)</td></tr><tr><td>Fusarium</td><td>sacchari</td><td>0</td><td>3</td><td>MN193868.1 (2)</td></tr><tr><td>Trichoderma</td><td>erinaceum</td><td>0</td><td>2</td><td>MK109820.1 (2)</td></tr><tr><td>Fusarium</td><td>napiforme</td><td>1</td><td>1</td><td>MH862670.1 (2)</td></tr><tr><td>Trichoderma</td><td>paraviridescens</td><td>1</td><td>1</td><td>MF782827.1 (1) / MK418756.1 (1)</td></tr><tr><td>Epicoccum</td><td>sorghinum</td><td>1</td><td>1</td><td>MF782827.1 (1) / MK418756.1 (1)</td></tr><tr><td>Trichoderma</td><td>sp.</td><td>0</td><td>2</td><td>KX449479.1 (1) / MK870964.1 (1)</td></tr><tr><td>Aspergillus</td><td>awamori</td><td>2</td><td>0</td><td>KY416558.1 (2)</td></tr><tr><td>Fusarium</td><td>falciforme</td><td>2</td><td>0</td><td>MT251175.1 (2)</td></tr><tr><td>Lasiodiplodia</td><td>mahajangana</td><td>2</td><td>0</td><td>MH057188.1 (2)</td></tr><tr><td>Trichoderma</td><td>asperellum</td><td>0</td><td>1</td><td>KX538815.1 (1)</td></tr><tr><td>Fusarium</td><td>circinatum</td><td>0</td><td>1</td><td>MK334369.1 (1)</td></tr><tr><td>Fusarium</td><td>cortaderiae</td><td>0</td><td>1</td><td>AH012626.2 (1)</td></tr><tr><td>Aspergillus</td><td>fumigatus</td><td>0</td><td>1</td><td>MH844690.1 (1)</td></tr><tr><td>Fusarium</td><td>graminearum</td><td>0</td><td>1</td><td>MK460853.1 (1)</td></tr><tr><td>Trichoderma</td><td>harzianum</td><td>0</td><td>1</td><td>JN116710.1 (1)</td></tr><tr><td>Diaporthe</td><td>kongii</td><td>0</td><td>1</td><td>KR024740.1 (1)</td></tr><tr><td>Curvularia</td><td>lunata</td><td>0</td><td>1</td><td>MN971669.1 (1)</td></tr><tr><td>Diaporthe</td><td>masirevicii</td><td>0</td><td>1</td><td>MF668289.1 (1)</td></tr><tr><td>Aspergillus</td><td>novoparasiticus</td><td>0</td><td>1</td><td>MH279415.1 (1)</td></tr><tr><td>Aspergillus</td><td>oryzae</td><td>0</td><td>1</td><td>MN648727.1 (1)</td></tr><tr><td>Rhizopus</td><td>stolonifer</td><td>0</td><td>1</td><td>MF461025.1 (1)</td></tr><tr><td>Trichoderma</td><td>trixiae</td><td>0</td><td>1</td><td>MN889512.1 (1)</td></tr><tr><td>Pestalotiopsis</td><td>vismiae</td><td>0</td><td>1</td><td>KP747694.1 (1)</td></tr><tr><td>Clonostachys</td><td>wenpingii</td><td>0</td><td>1</td><td>NR_119651.1 (1)</td></tr><tr><td>Cosmospora</td><td>butyri</td><td>1</td><td>0</td><td>KU204560.1 (1)</td></tr><tr><td>Lasiodiplodia</td><td>citricola</td><td>1</td><td>0</td><td>KU530119.1 (1)</td></tr><tr><td>Glomus</td><td>clarum</td><td>1</td><td>0</td><td>AY035654.1 (1)</td></tr><tr><td>Fusarium</td><td>dlamini</td><td>1</td><td>0</td><td>MN173109.1 (1)</td></tr><tr><td>Fusarium</td><td>ficicrescens</td><td>1</td><td>0</td><td>KP662895.1 (1)</td></tr><tr><td>Bionectria</td><td>ochroleuca</td><td>1</td><td>0</td><td>EU552110.1 (1)</td></tr><tr><td>Clonostachys</td><td>rosea</td><td>1</td><td>0</td><td>MH047188.1 (1)</td></tr><tr><td>Clonostachys</td><td>sp.</td><td>1</td><td>0</td><td>MH681594.1 (1)</td></tr><tr><td>Davidiella</td><td>sp.</td><td>1</td><td>0</td><td>KX621979.1 (1)</td></tr><tr><td>Diaporthe</td><td>sp.</td><td>1</td><td>0</td><td>MH220834.1 (1)</td></tr><tr><td>Lasiodiplodia</td><td>theobromae</td><td>1</td><td>0</td><td>KR260829.1 (1)</td></tr><tr><td>Aspergillus</td><td>welwitschiae</td><td>1</td><td>0</td><td>MH374611.1 (1)</td></tr></table>

<footnote>(*) The listed accession numbers correspond to the most frequent top match for each species.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Relative abundances of fungal isolates sampled from healthy and infected pineapple fruitlets at genus and species taxonomic levels.</figcaption></figure>

## 3.2. *Correlation between the Diversity of Fruitlet Mycobiomes and FCR Incidence*

To determine the prevalence of FCR over the experimental plots, the number of IF per fruit was reported on a total of 240 pineapple fruits. The deviance test evidenced significant interaction between plot orientation and adjacent crops (p < 0.0001) on the average of IF per pineapple fruit. Considering Plot 1, the mean IF counts per pineapple fruit (mean ± sd) were 6.65 ± 4.36 for the row section proximal to the mango orchard, 4.92 ± 4.16 for rows localized in the center of the plot, and 1.57 ± 1.92 for rows near the jackfruit trees. Similar patterns were observed over rows of Plot 2, however, showing the highest infestation level with a mean of 8.55 ± 4.80 IF per pineapple fruit grown near the mango orchard, 5.82 ± 4.32 IF for pineapples in the plot center, and 4.75 ± 2.81 IF for fruits proximal to the jackfruit trees (Figure 3). Among Plot 1, the main significant differences in the average IF number were observed between mango orchard and jackfruit trees sections (p = 3.03 × 10⁻⁷) and between center and jackfruit trees sections (p = 5.88×10⁻⁵). Considering Plot 2, the average IF number per fruit were significantly different between mango orchard and jackfruit trees sections (p = 0.0002) and between center and mango orchard sections (p = 0.013). Moreover, the comparison data in jackfruit sections also evidenced a significant rise of the average IF number per fruit in Plot 2 compared to Plot 1 (p = 0.0001).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Scatter and box plots of numbers of FCR-infected fruitlets (IF) per pineapple fruit in the function of the adjacent crop according to plot orientation (Plot 1: Red, Plot 2: Blue). Black dots with bars show mean numbers of IF per pineapple with asymptotic confidence intervals at 95%. Letters indicate significant average differences between plot x sections (n = 40) according to Tukey's multiple comparison test (p ≤ 0.013).</figcaption></figure>

Considering the distribution profile of FCR symptoms, the distribution of the identi-
fied fungal species over the plots was investigated to search for a correlation. The first
factorial axis (Dim1) of the PCA computed on the abundance of the 344 identified isolates
according to plot sections explained 55.6% of the variance (Figure 4). Jackfruit trees and
mango orchard sections of Plot 1 contributed to Dim1 with 48.9% and 22.8%, respectively.
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Principal component analysis (PCA) of combinations between proximal crop and plot direction. Only the 15 most contributory fungal species are shown (contributions ranging from 23.60% to 1.08%). Plot 1: North-West/South-East-directed rows; Plot 2: South-West/North-East-directed rows; mango orchard or jackfruit trees as adjacent crops, respectively; center: No adjacent crops. The size of each dot is proportional to the average number of infected fruitlets per fruit for each section of plots.</figcaption></figure>

Symmetrically, Dim1 was supported by *F. proliferatum*, *F. ananatum*, *T. stolii*, *F. equiseti*, and *F. oxysporum* with respective contributions of 30.4%, 23.7%, 9.02%, 8.4%, and 7.68%. Correlation between variables and the first component evidenced that *F. ananatum* (R = 0.99) and *F. proliferatum* (R = 0.96) were positively correlated to PC1 at *p* = 0.0002 and *p* = 0.002, respectively. Thus, the mango orchard section of Plot 1 was associated with elevated abundances of *F. proliferatum*, *F. ananatum*, *T. stolii*, *F. equiseti*, and *F. oxysporum* by opposition to the jackfruit trees section of the same plot. Additionally, 16.8% of the variance was explained by the Dim2. The construction of PCA factorial plane established that the center section of Plot 1 (57.71%) and mango orchard proximal rows of Plot 2 (18.1%) were the main contributive individuals of Dim2 in association with *F. incarnatum* (15.0%), *A. flavus* (12.2%), *F. solani* (9.0%), and *F. verticillioides* (7.5%). Moreover, *A. flavus* (R = 0.95), *F. incarnatum* (R = 0.84), *F. dlamini* (R = 0.83), *Diaporthe* sp. (R = 0.83), *D. masirevicii* (R = 0.83), *D. kongii* (R = 0.83), *C. butyri* (R = 0.83), *C. wenpingii* (R = 0.83), *A. fumigatus* (R = 0.83), and *F. verticillioides* (R = 0.82) were the main variables positively correlated to Dim2 at *p* ≤ 0.04. This demonstrates that center and mango orchard sections of Plot 1 were characterized by different structuration in fungal communities even if their average number of IF per fruit were not significantly different. In addition, the jackfruit trees section of Plot 1 that presented the lowest FCR incidence also demonstrated the poorest fungal diversity and abundances. Considering that the most contributive fungal species were also among the most abundant in fruitlet samples, this suggested that fungal species repartition was closely linked with FCR occurrence inside the plot sections. Consequently, relations between host tissues and fungal species prevalence were tested through the establishment of hierarchical clustering based on isolate counts in the 24 quadrats (Figure 5). Two distinct groups could be observed: Group 1 highlighted co-occurrence depicting similar profiles of highly represented fungal species, corresponding to identification in at least 25% of each sample set.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Heat map representation of fungal species frequencies in healthy and infected fruitlets of each quadrat. Group 1 corresponds to species identified in at least 25% of healthy and infected fruitlet sets, and Group 2 shows species identified in less than 25% of healthy and infected fruitlet sets.</figcaption></figure>

HF (n = 24) seems mainly characterized by *F. ananatum*, *F. proliferatum*, *F. equiseti*, *T. stollii*, *T. amestolkiae*, *F. verticillioides*, *F. oxysporum*, and *P. curvata*, together gathering 61.8% of the fungal flora. Interestingly, IF clustering (n = 24) shared the same species in addition to *A. flavus* and *F. solani*, which all contributed to 70.7% of the overall mycobiome diversity. In both types of fruitlet sets, a second group (Group 2) contained low abundance species identified in a limited number of quadrats and were thus considered as non-prevalent species for the rest of the analyses. Fungal communities were thus explored according to a combinatorial approach. Pairwise effect sizes were computed and revealed positive co-occurrence patterns in both HF and IF (only significant species combinations were displayed). Both matrixes showed 19 and 22 positive pair profiles for HF and IF fungal
species datasets, respectively (Figure 6). Positive associations evidenced species pairs with significantly greater and large frequency than predicted by the model. By contrast, random associations were defined by observed values nearly equal to expectations and with no significant difference.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Matrixes of fungal species co-occurrences determined in healthy (a) and infected fruitlets (b) based on a pairwise approach. Significant positive associations are indicated in blue with the corresponding number of observations. Random associations are indicated in grey. The correlation between matrixes was assayed by the Mantel test (*** significance at *p* ≤0.001).</figcaption></figure>

F. proliferatum was found in a positive combination with all individuals belonging to Group 1. Despite being the most prevalent species in the dataset, F. proliferatum had the highest number of pairwise associations in both tissues by positively co-occurring with 10 species (69 observations) in HF and with 8 species in IF (47 observations). Alternatively, co-occurrence patterns of several species were contrasted by presenting fewer combinations in IF samples than in HF samples. This trend was reported for F. equiseti and T. amestolkiae, which were represented 28 and 25 times in associations with 5 and 4 species in HF but only combined 9 and 5 times with 2 and 1 species in IF, respectively. On the contrary, 5 species displayed more positive combinations in the IF set than in the HF set. As an example, P. curvata was described to be positively associated only with F. proliferatum in HF (5 observations), while a combination with 6 species (31 observations) was highlighted in IF samples. This profile was also observed, to a lesser extent, for T. stollii (HF: 1 association observed 6 times; IF: 5 associations observed 31 times), F. solani (HF: 1 association observed 5 times; IF: 4 associations observed 19 times), and F. oxysporum (HF: 1 association observed 6 times; IF: 5 associations observed 23 times). Association patterns suggested that the fungal flora associated with HF and IF was structured. However, the Mantel test processed on distance matrixes demonstrated a significant relationship between the mycobiomes of the HF and IF samples (*p* = 0.0001). This revealed that species associations did not depend on the fruitlet phenotype (healthy or infected) when considered through a pairwise approach.

Among plants pathogenic species, strains may exhibit a contrasting level of virulence, partly due to a genetic polymorphism inherited from co-evolution. TEF-1α and ITS gene sequence datasets represented, respectively, 101 and 98 Fusarium DNA sequences, with at least 98% of identity according to the reference nucleotide sequence from NCBI (Table S1). These loci were submitted to phylogenetic analyses following the Maximum Likelihood (ML) approach, with a bootstrap of 100 replications. Over the TEF-1α region, 6 clades related to the 6 different species represented in this sequence dataset were discriminated without the distinction of strains isolated from HF or IF samples, thus describing a high level of homology (Figure 7).
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Phylogenetic tree of *Fusarium* species computed by Maximum Likelihood (ML) analysis (TIM2ef+G as the best fit model) and based on a TEF-1α sequences dataset. Only ML bootstrap branches that support values greater than 70% are shown. Strains isolated from infected (IF) and healthy fruitlets (HF) are indicated in red and black, respectively. The tree is rooted with *F. graminearum* accession KU885947.1.</figcaption></figure>
The same method was performed over the 61 *Talaromyces* strains. All of the 48 isolates characterized under the β-tubulin region corresponded to *T. amestolkiae* and *T. stollii*. As for the *Fusarium* strains, 2 clades were distinguished in concordance with the represented species (Figure 8). Strains were closely related inside clades, regardless of the phenotype (infected or healthy) of the fruitlets. This indicated that species virulence and pathogenicity were not linked to a polymorphism of TEF-1 α, ITS, or β-tubulin regions. Moreover, sequence analysis on TEF-1α and β-tubulin loci enabled a higher taxonomic resolution when compared to ITS (Figure S1; Figure S2).
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Phylogenetic tree of *Talaromyces* species computed by Maximum Likelihood (ML) analysis (TPM2+G as the best fit model) and based on a β-tubulin sequences dataset. Only ML bootstrap branches that support values greater than 70% are shown. Strains isolated from infected (IF) and healthy fruitlets (HF) are indicated in red and black, respectively. The tree is rooted with *T. purpureogenus* accession MN233617.1.</figcaption></figure>
### 3.3. Potential Sources of Inoculum and Dispersion Patterns

The distribution of fungal species over the experimental plots suggested the presence of an important inoculum source. Fungal flora of the soil at the interface between pineapple and jackfruit trees or the mango orchard was studied in order to elucidate its origin. Along the two junctions, 26 (mango orchard) and 33 (jackfruit trees) isolates were identified as belonging to 27 genera and 35 species. Most of the identified species belonged to the Trichocomaceae (32.2%) and Nectriaceae (10.2%) families (Figure 9A). The Trichocomaceae family was represented by *Talaromyces* (teleomorph of *Penicillium*), *Penicillium*, and *Aspergillus* species, which were detected in lower abundances in the soil related to jackfruit trees (19.2%) than in the mango orchard soil (45.6%), respectively.

*Talaromyces* was strictly represented by *T. purpureogenus*, while 8 *Penicillium* species were isolated: *P. janthinellum*, *P. sumatraense*, *P. citrinum*, *P. guanacastense*, *P. multicolor*, *P. ochrochloron*, *P. sclerotigenum*, and *P. spinulosum*. This family was also supported by 4 *Aspergillus* species: *A. sclerotiorum*, *A. subramanianii*, *A. flavipes*, and *A. tardicrescens*. Moreover, the Nectriaceae family was only represented by *Fusarium* species. The fungal isolates that belonged to *Fusarium* genus were detected at higher frequencies close to the jackfruit trees (15.4%) when compared to the mango orchard interface (6.1%) (Figure 9B, Table S2). Considering the diversity, two distinct species were identified in each junction: *F. equiseti*/*F. polyphialidicum* at the jackfruit trees and *F. oxysporum*/*F. solani* at the mango orchard interface. Remarkably, among all fungal soil isolates, no *F. ananatum* was identified. Subsequent isolates belonged to the *Cladosporium*, *Plectosphaerella*, *Trichoderma*, *Diaporthe*, *Mortierella*, *Xepicula*, *Bartalinia*, *Roussoella*, *Acrostalagmus*, *Cystofilobasidiales*, *Pestalotiopsis*, *Phoma*, *Pleosporales*, *Pyrenochaeta*, *Heterocephalum*, *Edenia*, *Purpureocillium*, *Aureobasidium*, *Fennellia*, *Lecanicillium*, *Robillarda*, *Peniophora*, and *Leptosphaeria* genera and were found with relative abundances ranging from 3% to 11.5%. Several contrasted profiles were noticed with the strict characterization of *Lecanicillium*, *Heterocephalum*, *Aureobasidium*, *Edenia*, *Fennellia*, *Leptosphaeria*, *Peniophora*, *Purpureocillium*, *Robillarda*, and *Talaromyces* genera in samples from the soil interface with the mango orchard. In contrast, *Diaporthe*, *Acrostalagmus*, *Cystofilobasidiales*, *Pestalotiopsis*, *Phoma*, *Plectosphaerella*, *Pleosporales*, *Pyrenochaeta*, *Roussoella*, and *Xepicula* genera were exclusively detected in the jackfruit tree junctions.

To evaluate the influence of abiotic factors on the fungal species spreading, meteorological parameters were obtained for the seven-month period from the flower induction treatment to harvest. Average temperatures ranged from 19.5 °C to 23.49 °C, with a low average rainfall below 5 mm. Data also demonstrated that the major wind direction was south-east during anthesis and fruit development stages, with a gust speed between 69 and 88 km h⁻¹ (Table 1). The dynamics of these factors typically described the setting of the austral winter.
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Relative abundances of fungal isolates identified from soil sampled on either side of the pineapple plots (jackfruit trees or mango orchard) at (a) family, (b) genus, and species taxonomic levels.</figcaption></figure>
### 3.4. *Koch's Postulates for the Determination of New FCR Pathogens*

In IF samples, species belonging to the *Fusarium fujikuroi* species complex were represented with high frequencies and concomitantly with two *Talaromyces* species related to the *Talaromyces purpureogenus* species complex. Koch's postulates were verified for all of the following species: *F. proliferatum*, *F. sacchari*, *F. oxysporum*, *T. stollii*, and *T. amestolkiae*, in order to test for their pathogenicity.

After 7 days of incubation on PDA Petri dishes, all species presented fast growth rates with a mean diameter (av. ± sd) reaching 7.25 cm, 7.77 ± 0.69 cm, 7.75 ± 0.70 cm, and 7.90 ± 0.19 cm for *F. ananatum*, *F. oxysporum*, *F. proliferatum*, and *F. sacchari*, respectively (Figure 10 A1-H1). *T. stollii* and *T. amestolkiae* showed a slightly slower growth rate with respective mean diameters of 6.66 ± 0.08 cm and 6.68 ± 0.90 cm (Figure 10 I1-N1). *Fusarium* colonies were circular with a thick and aerial white mycelium, except for strain BP575 (*F. sacchari*), which showed a short aerial mycelium. Species can be distinguished by their coloration when observed from a bottom range from orange (*F. sacchari*), to violet (*F. oxysporum*), to ochre/brown (*F. proliferatum*).

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Morphological observations of fungal species isolated from pineapple-infected fruitlets and assayed in Koch's postulates. (1) Colonies were incubated on potato dextrose agar (PDA) for 7 days at 27 °C in the dark. Each colony is shown from top (left) to bottom (right). (2) Conidia. *Fusarium sacchari* strains (a) BP138 and (b) BP575; *Fusarium ananatum* strain (c) BP383; *Fusarium oxysporum* strains (d) BP369 and (e) BP460; *Fusarium proliferatum* strains (f) BP114, (g) BP429, and (h) BP436; *Talaromyces stollii* strains (i) BP054, (j) BP185, and (k) BP462; *Talaromyces amestolkiae* strains (l) BP002, (m) BP257, and (n) BP605. Scale bar: 20 µm.</figcaption></figure>

*Talaromyces* species showed circular to oval colonies with short and rough mycelium characterized by a green coloration in the oldest zone and with a white ring on the edges. The isolates identified as *T. amestolkiae* could be distinguished by a red pigmentation visible at the bottom of the Petri dishes. This was in addition to the secretion of a red exudate in the center of the colonies (Figure 10 L1-N1). Microscopic observations of *Fusarium* species isolates showed straight (Figure 10 A2-C2) and curved fusiform macroconidia (Figure 10 D2, E2), while microconidia were ovoidal (Figure 10 D2, E2) or ellipsoidal (Figure 10 A2-C2, F2-H2). *Talaromyces amestolkiae* and *Talaromyces stollii* exhibited aggregated conidia with an ellipsoidal shape (Figure 10 I2-N2).
Subsequently to sporulation on PDA, controlled inoculations were performed by injecting a normalized conidia solution into a blossom cup of fruitlets. After 7 days of incubation, fruitlets from low, median, and upper sections of pineapple fruit were extracted and cut into two equal parts. Samples inoculated with sterile water as well as the non-inoculated fruitlets (Mock) did not show any visible FCR symptom. All 5 tested species (13 isolates) led to the development of a vitreous aspect or black spot (Figure 11). Following *F. oxysporum* inoculations, black spot development was observed in both pineapple fruit sections.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** FCR symptoms observed 7 days post-inoculation (dpi) following controlled inoculations of five putative pathogenic fungal species in the blossom cup of fruitlets from the low (L), median (M), and upper (U) parts of 'Queen Victoria' pineapple fruit. *F. ananatum* strain BP383 was used as positive control; H₂O: Sterile water-inoculated; mock: No inoculation.</figcaption></figure>

Similar observations were recorded for *F. sacchari*. Interestingly, *F. proliferatum* and *T. stollii* caused strong dark spot development in the fruitlet flesh. Through *T. stollii*-inoculated fruitlets, symptoms came out with a red pigmentation that was visible in samples from the median section only. Inoculations with 3 different isolates of *T. amestolkiae* led to slight necrosis associated with a vitreous appearance. Comparable symptoms were also reported after *F. ananatum* conidia injection, considered in this experiment as a positive control [10].

To fulfill Koch's postulates, the potential causal agents were recovered by a molecular approach (PCR-DGGE) that allowed for a culture-independent characterization of microbial communities associated with a semi-quantitative determination of band patterns. To this end, *Fusarium*- and *Talaromyces*-inoculated fruitlets were distinguished, but the identification of species could not be achieved. First, a reference gel was performed on pure fungal strains using migration conditions that were similar to that of the inoculated samples (Figure S3). According to the UPGMA construction method, fungal patterns were structured depending on the genus of the tested strain (Figure 12). Two groups were then identified. One group contained *Fusarium*-inoculated samples and presented the
‘*Fusarium* marker’ according to the reference profiles. The second group was supported by *Talaromyces*-inoculated fruitlets with matching markers and control samples.

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

**Figure 12.** Dendrogram of PCR-DGGE profiles constructed by the UPGMA method. Red rectangles show *Fusarium* markers in *Fusarium*-inoculated fruitlets. Green rectangles show *Talaromyces* markers of *Talaromyces*-inoculated fruitlets. Red arrows indicate naturally present *Fusarium* sp. in control fruitlets. Green arrows indicate naturally present *Talaromyces* sp. in control fruitlets. BP138 and BP575: *Fusarium sacchari*; BP383: *Fusarium ananatum*; BP369 and BP460: *Fusarium oxysporum*; BP114, BP429, and BP436: *Fusarium proliferatum*; BP054, BP185, and BP462: *Talaromyces stollii*; BP002, BP257, and BP605: *Talaromyces amestolkiae*. H₂O: sterile water-inoculated samples; mock: No inoculation. Letters label DNA bands, which were sequenced.

The DNA bands recovered from the gel were sequenced, and the identification results confirmed that the targeted bands corresponded to the *Fusarium* sp. and *Talaromyces* sp. that were used for inoculations (Table S3). Another difficulty relied on the presence of an important natural core mycobiome partly involving *Fusarium* and *Talaromyces* species. The fungal diversity of fruitlets could be appreciated in negative control samples that naturally presented fungal banding patterns corresponding to *Fusarium* (red asterisks) and *Talaromyces* (green asterisks) migration profiles. Fungi-inoculated samples, however, showed an intense band that was not detectable in symptom-free fruitlets that corresponded to the inoculated fungi.

**4. Discussion**

By analyzing pineapple cultivatable fungal flora, the composition of the mycobiome associated with FCR could be characterized. The data generated in the present study showed that an interaction between biotic and abiotic factors may impact disease spreading over a production area. Moreover, the determination of species contribution to FCR disease led to the discovery of new causal agents.

Our results demonstrate that structuration of the fungal flora of fruitlets was mainly composed of two species complexes that were recently characterized as *Fusarium fujikuroi* and *Talaromyces purpureogenus* [49,50]. This type of association was previously described in FCR through the isolation of *F. verticillioides* and *T. funiculosus* recovered from FCR symptoms [8,51], even though *T. funiculosus* had never been defined as prevalent in Réunion Island before. Previous sampling over eight production areas in the island showed a predominance of *F. ananatum* (72%), *F. oxysporum* (6.67%), *F. proliferatum* (0.67%), and *T. stollii* (20.67%) [52]. Conjointly, conidia with a fusiform shape were observed in a blossom cup of symptomless pineapples related to the 'Queen Victoria' cultivar, while ellipsoidal
conidia were observed in the 'MD-2' defined as the FCR-resistant cultivar of reference [53]. However, the identification of species such as *Disporotrichum dimorphosporum*, *Galactomyces candidus*, and *Clavispora lusitaniae* in freshly cut pineapple (cv. 'Queen Victoria') revealed a higher level of fungal diversity in the fruit flesh [54].

By focusing on a specific plot environment, 49 fungal species were identified through 344 isolates. The infection level was reported over the three plot sections defined for the two experimental plots. Results showed that rows proximal to the mango orchard presented an elevated number of FCR symptoms, and this number was significantly lower in the rows proximal to the jackfruit trees. This suggested the existence of a disease gradient ranging from the mango orchard (high FCR incidence zone) to the jackfruit trees (low FCR incidence zone). Several soil-borne *Fusarium* species can cross long distances through water, especially irrigation water, under the form of chlamydospores. As an example, *F. oxysporum* causing *Fusarium* wilt of banana spread through the soil conductance, leading to pathogen invasion [55]. Nevertheless, the composition of the fungal flora of soil samples from both junctions was not correlated with the disease gradient observed in fruits over row sections. In fact, FCR pathogens such as *Fusarium* species were surprisingly mainly identified over the jackfruit tree interface, and members of the *Talaromyces* genus were poorly represented across soil samples. Mango orchard junctions only demonstrated the presence of *Talaromyces purpureogenus* (one of the four components of the *T. purpureogenus* species complex) as a potential new contributor to FCR. Interestingly, *F. proliferatum* and *F. sacchari* were previously identified as pathogens causing mango malformation disease and leading to an anarchic development of inflorescence [56,57]. Even though this pathology had never been reported in the various mango cultivars of this orchard, it is proposed that mango trees could represent a prevalent habitat for some FCR pathogens. In fact, across a landscape, fungi can switch from an endophyte lifestyle on a specific host to a pathogenic behavior in another host [58,59]. By considering the interaction between crop proximity (adjacent cultures were only 2 m away) and meteorological data, it was suggested that wind could help the dispersion of conidia from a mango tree canopy to pineapple rows, and it was hypothesized that FCR was an airborne disease. This profile also emphasized the significance of row orientation in relation to the main wind direction. The FCR gradient was slightly more pronounced in Plot 2, and the mean level of infestation per fruit was significantly higher. This suggested that the inoculum widespread from the mango tree canopy to the pineapple plot was carried by the wind further than for the plot having rows perpendicular to the major wind direction. In fact, the N-W/S-E-directed rows (Plot 1) showed a strong FCR gradient with a significant difference between plot sections. Thus, plants proximal to the mango orchard appeared as a 'buffer zone' for the shared pathogens and limited the wind dispersion of conidia over the plot. These observations are consistent with Parnell et al., who demonstrated that the spatial configuration of host in a landscape is the keystone determinant of plant-disease epidemics [60].

Among FCR pathogens, *F. ananatum* has been described the most. *F. ananatum* isolates from pineapples of Costa Rica and Ecuador showed a capacity to produce fumonisins B1 (FB1), fumonisins B2 (FB2), fumonisins B3 (FB3), and beauvericin (BEA). However, the production levels observed under in vitro conditions were only in low concentrations, thus constituting a limited risk for food contamination [61]. Following FCR sampling in Reunion Island, *F. proliferatum* and *F. oxysporum* have been described as mycotoxin producers of FB1, FB2, and BEA, while no toxigenic profile was detected for *T. stollii*. Significant amounts of FB1 and BEA have also been detected in naturally infected fruitlets in comparison to healthy ones [52]. It was thus suggested that these species could also be involved in FCR pathogenesis with a significant impact on food safety.

The characterization of the mycobiome related to both healthy and infected fruitlets resulted in the determination of significant pairwise co-occurrence between *Fusarium* and *Talaromyces* species. Although *F. ananatum* was described as the first causal agent of FCR in Reunion Island [11,52], our data showed a prevalence of *F. proliferatum* (50 strains), which seemed to define a hub over fungal interactions in pineapple. The prevalence of
this species was also demonstrated following the sampling of diseased leaves and pine-
apple fruits in Malaysia [62]. Interestingly, this study identified *F. proliferatum* as the main
causal agent of leaf spot disease and pineapple fruit rot, characterized by the development
of brown necrosis in the fruit flesh and skin.

Across related abundance profiles and distribution patterns, *F. proliferatum*, *F. sacchari*, *F. oxysporum*, *T. stollii*, and *T. amestolkiae* were investigated for their pathogenicity. Following the injection of conidia in a blossom cup, all of these species showed a capacity to induce symptoms related to FCR. Several symptoms observed after seven days of incubation were only slightly pronounced, particularly for *F. oxysporum* and *T. amestolkiae*. This could be explained by the time variation for fungal development and symptom appearance, which may fluctuate according to the studied species. Surprisingly, the evolution of fungal flora of minimally processed pineapples also revealed the implication of *T. amestolkiae* in fruit spoilage over cold storage conditions [54]. Taken together, these elements illustrate that specific structural arrangement of microbial communities may lead to pathogenesis and/or to fruit quality depreciation [63]. We, therefore, propose the FCR pathosystem to be considered as multi-partite and hypothesize that other species related to *Fusarium fujikuroi* and *Talaromyces purpureogenus* species complexes could be involved. Thus, the presence of multiple pathogens in a production area may increase the epidemiological risks. This partially explained the difficulties faced over the past decades to elucidate the infectious process leading to pineapple FCR susceptibility.

Conventional verification of Koch's postulates had to be adapted for pineapple fruits
due to the constraints of dealing with non-microbial-free tissues for in vivo inoculation.
As exposed with DGGE migration profiles, the initial (natural) microbial loads in fruitlets
contain a high level of fungal diversity as numerous DNA bands that could be detected
in symptomless fruitlets. As a comparison, 5-9 bands could be visualized in DGGE lanes
of negative controls (H₂O and Mock), while 1-5 species per fruitlets were isolated by a
classical microbial procedure. The recovery of the pathogens from inoculated fruitlets was
achieved following the establishment of reference migration profiles for each tested spe-
cies. Over the 13 strains, strong DNA signals corresponded to a highly abundant inocu-
lated species. This showed that the tested strains had been able to extend in the blossom
cup of all biological replicates, except for one fruit inoculated with the *F. sacchari* (strain
BP138). The inoculation of the tested species also modified the basal fungal communities.
Numerous bands with various intensity were detected and contributed to the discrimina-
tion between *Fusarium*- and *Talaromyces*-inoculated mycobiomes. Interestingly, the fungal
flora profiles of negative control samples were related to *Talaromyces*-inoculated samples
rather than to *Fusarium*-inoculated fruitlets. This suggested a pathogen-specific evolution
of the fruitlet mycobiome that is much more noticeable following inoculations with
*Fusarium* species. The over-representation of these species may thus involve potentially
synergistic or antagonistic microbial interactions. As an example, *F. verticillioides*, as a part
of the endophytic fungal flora of maize, was described for its capacity to lower corn smut
disease severity by limiting the biomass development of the pathogen *Ustilago maydis* [64].
In pineapple, the determination of positive pairwise co-occurrences led to the assumption
that mycobiomes of healthy and naturally infected fruitlets were significantly correlated
and that species pair combinations thus did not diverge from fruitlet phenotypes. More-
over, phylogenetic trees evidenced strong relationships among strains recovered from
both HF and IF, proving that pathogenesis did not correlate with a fungal genetic poly-
morphism. Nevertheless, the interactions inside fungal communities are complex and
may imply numerous species with various communication strategies, as observed in nu-
merous studies [65-68]. It is important to consider that the present approach for the de-
tected fungal species did not enable a determination of either the form (micro/macro-
conidia and mycelium) or their respective concentration in the fruitlets. The colonization
of host tissues could be achieved by prevalent species forming an extensive mycelial net-
work and could potentially modulate the development of specialized structures of the
remaining fungal species [69]. The establishment of cooperative or competitive fungal in-
teractions partly relies on chemical recognition [66,68]. Interestingly, *T. stollii* and *T. amestolkiae* were defined as red-pigment-producing fungi that correspond to azaphilone extrolites [70]. This corroborates the observations performed following the in vitro culture of *T. amestolkiae* and the in vivo development of *T. stollii* [50]. These secondary metabolites are related to antiviral, antimicrobial, antifungal, and other biological activities [71]. This illustrates the necessity to determine how chemical communication between fungi may contribute to FCR establishment. In addition, the present study focused on fungal flora without considering the role of yeast and bacteria, which may also influence the microbial dialog during pathogenesis. In tomato, Wei et al. showed that endophytic microbial com-munities sharing an ecological niche with the *Ralstonia solanacearum* pathogen exhibited a low disease incidence [72]. As frequently exemplified, *F. oxysporum* has been described for its capacities to secrete and accumulate fusaric acid in order to inhibit the production by *Pseudomonas fluorescens* of 2,4-diacetylphloroglucinol acting as an antibiotic [73]. Biologi-cal association between microbial species represented a source for the determination of biocontrol agents. In pineapple, the antagonistic activity of *Trichoderma asperellum* against *F. guttiforme* was already described for the control of fusariosis [74].

In addition, both relative abundance and spatial distribution of fungal species are
also fluctuating in accordance with abiotic factors as precipitation, relative humidity, tem-
perature, and wind speed [75]. In the present study, pineapples were cultivated during
austral winter, which is considered the prevalent season for FCR incidence in Reunion
Island [18]. The SIMPIÑA model was developed to predict the quality of 'Queen Victoria'
following total soluble sugar estimation, agricultural practices, and production area [76].
Our data highlighted that the understanding of plant disease outbreaks requires the eval-
uation of both biotic and abiotic environments. *Fusarium* species belonging to the FFSC,
such as *F. sacchari*, *F. verticillioides*, and *F. proliferatum*, were reported as causal agents of
several tropical crops, including sugarcane and mango (as previously demonstrated),
which constitute preponderant cultures in Reunion Island [77-79]. This supports the po-
tential cross-contribution between shared pathogens and abiotic factors on fungal epi-
demics in tropical and subtropical regions. The evaluation of influence parameters and
their interactions may lead to the definition of risk factors for FCR occurrence. These var-
iables could then be incremented into the VICTORIA database to adapt disease manage-
ment strategies to various plot contexts [80,81].

**5. Conclusions**

Our study provided a new understanding of the FCR pathosystem by demonstrating the implications of *F. proliferatum*, *F. oxysporum*, *F. sacchari*, *T. stollii*, and *T. amestolkiae* on pathogenesis. We demonstrated that the mycobiome composition of fruitlets is influenced by adjacent crops that share common pathogens and may cause elevated levels of FCR incidence. It appears essential to consider the plot environment, especially row orientation, in relation to the major wind direction. Finally, similar fungal populations were described in healthy and infected fruitlets, suggesting that pineapple susceptibility may be inherent to chemical communication between pathogens that modulate strategies of host tissue colonization. Determining in vitro and in vivo pathogens interacting-behaviors would be necessary to understand the virulence factors of FCR-pathogens and subsequently investigate the molecular crosstalk of host-multi-pathogens interaction.
**Supplementary Materials:** The following are available online at www.mdpi.com/2309-608X/7/3/175/s1, Figure S1: Phylogenetic tree of *Fusarium* species computed by Maximum Likelihood (ML) analysis (JC+I+G as the best fit model) and based on the ITS sequences dataset. Only ML bootstrap branches that support values greater than 70% are shown. Strains isolated from infected and healthy fruitlets are indicated in red and black, respectively. The tree is rooted with *F. sporotrichoides* accession MT218410.1, Figure S2: Phylogenetic tree of *Talaromyces* species computed by Maximum Likelihood (ML) analysis (K80+G as the best fit model) and based on the ITS sequences dataset. Only ML bootstrap branches that support values greater than 70% are shown. Strains isolated from infected and healthy fruitlets are indicated in red and black, respectively. The tree is rooted with *T. cecidicola* accession MT182957.1, Figure S3: Reference DGGE profiles over ITS-1 region of DNA extracted from pure strains evaluated in Koch's postulates. (Mix) indicates the combined migration of *F. ananatum* strain BP383, *F. sacchari* strain BP575, *F. oxysporum* strain BP460, *F. proliferatum* strain BP429, *T. stollii* strain BP054, and *T. amestolkiae* strain BP002, Table S1: Identification of fungal species isolated from healthy and naturally infected pineapple fruitlets, Table S2: Identification of fungal species isolated from soil coring points of cultures proximal to pineapple plots. Table S3: Identification of fungal species from DNA extracted from DGGE gels corresponding to *Fusarium-* and *Talaromyces*-inoculated pineapple fruitlets.

**Author Contributions:** Conceptualization, M.V., J.-C.M., and S.S.-G.; methodology, M.V. and C.L.-J.; formal analysis, M.V. and F.C.; investigation, M.V. and C.S.; data curation, M.V.; writing—original draft preparation, M.V.; writing—review and editing, M.V., J.-C.M., C.L.-J., C.S., S.S.-G., F.C., and M.C.; visualization, M.V. and F.C.; supervision, J.-C.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** Funding for this research, including M.V.'s PhD grant, was provided by the European Regional Development Fund (FEDER Réunion, PO-COSAQ).

**Acknowledgments:** We acknowledge Pierre Lefeuvre (CIRAD, UMR PVBMT, Reunion Island) for his useful advice and help in phylogenetics. We also kindly thank Bernard Abufera (CIRAD, UPR GECO, Reunion Island) and Georget Tullus (CIRAD, UPR GECO, Reunion Island) for pineapple field maintenance.

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

## References

1. FAO. *Major Tropical Fruits—Preliminary Market Results 2019*; Food and Agriculture Organization of the United Nations: Rome, Italy, 2020; p. 24.

2. Bartholomew, D.P.; Hawkins, R.A.; Lopez, J.A. *Hawaii Pineapple: The Rise and Fall of an Industry*. *HortScience* **2012**, *47*, 1390–1398, doi:10.21273/HORTSCI.47.10.1390.

3. Nassr, M.S.; Abu-Naser, S.S. *Knowledge Based System for Diagnosing Pineapple Diseases*. *Int. J. Acad. Pedagog. Res.* **2018**, *2*, 12–17.

4. Gu; Zhan, R.L.; Zhan, L.B.; Gong, D.Q.; Jia, Z.W. *First Report of Fusarium Ananatum Causing Pineapple Fruitlet Core Rot in China*. *Plant Dis.* **2015**, *99*, 1653, doi:10.1094/PDIS-03-15-0279-PDN.

5. Rohrbach, K.G.; Pfeiffer, J.B. *Susceptibility of Pineapple Cultivars to Fruit Diseases Incited by Penicillium Funiculosum and Fusarium Moniliforme*. *Phytopathology* **1976**, *66*, 1386–1390, doi:10.1094/Phyto-66-1386.

6. Rohrbach, K.G.; Johnson, M.W. *Pests, Diseases and Weeds. In The Pineapple: Botany, Production and Uses*; Department of Tropical Plant and Soil Science, CTAHR, University of Hawaii at Manoa: Honolulu, HI, USA, 2003; pp. 203–251, doi:10.1079/9780851995038.0000.

7. Seifert, K.A.; Aoki, T.; Baayen, R.P.; Brayford, D.; Burgess, L.W.; Chulze, S.; Gams, W.; Geiser, D.; de Gruyter, J.; Leslie, J.F.; et al. *The Name Fusarium Moniliforme Should No Longer Be Used*. *Mycol. Res.* **2003**, *107*, 643–644, doi:10.1017/S095375620323820X.

8. Lim, T.-K. *Role of Penicillium Funiculosum* Strains in the Development of Pineapple Fruit Diseases. *Phytopathology* **1980**, *70*, 663, doi:10.1094/Phyto-70-663.

9. Petty, G.J.; Tustin, H.A.; Dicks, H.M. *Control of black spot disease / fruitlet core rot in queen pineapple with integrated mealybug, pineapple fruit mite and fungus control programmes*. *Acta Hortic.* **2006**, 143–149, doi:10.17660/ActaHortic.2006.702.17.

10. Jacobs, A.; Van Wyk, P.S.; Marasas, W.F.O.; Wingfield, B.D.; Wingfield, M.J.; Coutinho, T.A. *Fusarium Ananatum Sp. Nov.* in the Gibberella Fujikuroi Species Complex from Pineapples with Fruit Rot in South Africa. *Fungal Biol.* **2010**, *114*, 515–527, doi:10.1016/j.funbio.2010.03.013.

11. Barral, B.; Chillet, M.; Minier, J.; Léchaudel, M.; Schorr-Galindo, S. *Evaluating the Response to Fusarium Ananatum Inoculation and Antifungal Activity of Phenolic Acids in Pineapple*. *Fungal Biol.* **2017**, *121*, 1045–1053, doi:10.1016/j.funbio.2017.09.002.

12. Mourichon, X. *Pineapple Fruitlet Core Rot (Black Spot) and Leathery Pocket: Review and Prospects*. *Acta Hortic.* **1997**, 501–508.
13. Maltz, M.R.; Treseder, K.K.; McGuire, K.L. Links between Plant and Fungal Diversity in Habitat Fragments of Coastal Shrubland. *PLoS ONE* **2017**, *12*, e0184991, doi:10.1371/journal.pone.0184991.

14. Hugoni, M.; Luis, P.; Guyonnet, J.; el Zahar Haichar, F. Plant Host Habitat and Root Exudates Shape Fungal Diversity. *Mycorrhiza* **2018**, *28*, 451–463, doi:10.1007/s00572-018-0857-5.

15. Hawksworth, D.L.; Lücking, R. Fungal Diversity Revisited: 2.2 to 3.8 Million Species. In *The Fungal Kingdom*; Heitman, J., Howlett, B.J., Crous, P.W., Stukenbrock, E.H., James, T.Y., Gow, N.A.R., Eds.; ASM Press: Washington, DC, USA, 2017; pp. 79–95, ISBN 978-1-68367-082-7.

16. dos Santos Souza, C.R.; de Oliveira Barbosa, A.C.; Fortes Ferreira, C.; Vidigal Duarte Souza, F.; de Souza Rocha, L.; de Souza, E.H.; de Oliveira, S.A.S. Diversity of Microorganisms Associated to Ananas Spp. from Natural Environment, Cultivated and Ex Situ Conservation Areas. *Sci. Hortic.* **2019**, *243*, 544–551, doi:10.1016/j.scienta.2018.09.015.

17. Rohrbach, K.G.; Taniguchi, G. Effects of Temperature, Moisture, and Stage of Inflorescence Development on Infection of Pineapple by Penicillium Funiculosum and Fusarium Moniliforme Var. Subglutinans. *Phytopathology* **1984**, *74*, 995–1000, doi:10.1094/Phyto-74-995.

18. Fournier, P.; Benneveau, A.; Hardy, C.; Chillet, M.; Léchaudel, M. A Predictive Model Based on a Pluviothermic Index for Leathery Pocket and Fruitlet Core Rot of Pineapple Cv. 'Queen.' *Eur. J. Plant Pathol.* **2015**, *142*, 449–460, doi:10.1007/s10658-015-0625-8.

19. Gamliel-Atinsky, E.; Sztejnberg, A.; Maymon, M.; Shtienberg, D.; Freeman, S. Inoculum Availability and Conidial Dispersal Patterns of *Fusarium Mangiferae*, the Causal Agent of Mango Malformation Disease. *Phytopathology* **2009**, *99*, 160–166, doi:10.1094/PHYTO-99-2-0160.

20. Rieux, A.; Soubeyrand, S.; Bonnot, F.; Klein, E.K.; Ngando, J.E.; Mehl, A.; Ravigne, V.; Carlier, J.; de Lapeyre de Bellaire, L. Long-Distance Wind-Dispersal of Spores in a Fungal Plant Pathogen: Estimation of Anisotropic Dispersal Kernels from an Extensive Field Experiment. *PLoS ONE* **2014**, *9*, e103225, doi:10.1371/journal.pone.0103225.

21. Aylor, D.E. SPREAD OF PLANT DISEASE ON A CONTINENTAL SCALE: ROLE OF AERIAL DISPERSAL OF PATHOGENS. *Ecology* **2003**, *84*, 1989–1997, doi:10.1890/01-0619.

22. Rieux, A.; Lenormand, T.; Carlier, J.; de Lapeyre de Bellaire, L.; Ravigné, V. Using Neutral Cline Decay to Estimate Contemporary Dispersal: A Generic Tool and Its Application to a Major Crop Pathogen. *Ecol. Lett.* **2013**, *16*, 721–730, doi:10.1111/ele.12090.

23. Fredrikson, L.; Skinkis, P.A.; Peachey, E. Cover Crop and Floor Management Affect Weed Coverage and Density in an Establishing Oregon Vineyard. *HortTechnology* **2011**, *21*, 208–216, doi:10.21273/HORTTECH.21.2.208.

24. Charest, J.; Dewdney, M.; Paulitz, T.; Philion, V.; Carisse, O. Spatial Distribution of *Venturia Inaequalis* Airborne Ascospores in Orchards. *Phytopathology* **2002**, *92*, 769–779, doi:10.1094/PHYTO.2002.92.7.769.

25. Osborne, L.E.; Stein, J.M. Epidemiology of Fusarium Head Blight on Small-Grain Cereals. *Int. J. Food Microbiol.* **2007**, *119*, 103–108, doi:10.1016/j.ijfoodmicro.2007.07.032.

26. Dita, M.; Barquero, M.; Heck, D.; Mizubuti, E.S.G.; Staver, C.P. Fusarium Wilt of Banana: Current Knowledge on Epidemiology and Research Needs Toward Sustainable Disease Management. *Front. Plant Sci.* **2018**, *9*, 1468, doi:10.3389/fpls.2018.01468.

27. van der Does, H.C.; Constantin, M.E.; Houterman, P.M.; Takken, F.L.W.; Cornelissen, B.J.C.; Haring, M.A.; van den Burg, H.A.; Rep, M. Fusarium Oxysporum Colonizes the Stem of Resistant Tomato Plants, the Extent Varying with the R-Gene Present. *Eur. J. Plant Pathol.* **2019**, *154*, 55–65, doi:10.1007/s10658-018-1596-3.

28. Fournier, P. La Culture de l'ananas Victoria à La Réunion Pour l'exportation. In *Recueil de Bonnes Pratiques*; CIRAD: Reunion Island, France, 2011;

29. Smart Information System (SmartIS)-METEOR-CIRAD-UR AÏDA Available online: https://smartis.re/METEOR (accessed on 25 November 2020).

30. Météo France Available online: https://donneespubliques.meteofrance.fr/?fond=produit&id_produit=129&id_rubrique=52 (accessed on 30 September 2020).

31. O'Donnell, K.; Cigelnik, E.; Nirenberg, H.I. Molecular Systematics and Phylogeography of the *Gibberella Fujikuroi* Species Complex. *Mycologia* **1998**, *90*, 465–493, doi:10.1080/00275514.1998.12026933.

32. O'Donnell, K.; Kistler, H.C.; Cigelnik, E.; Ploetz, R.C. Multiple Evolutionary Origins of the Fungus Causing Panama Disease of Banana: Concordant Evidence from Nuclear and Mitochondrial Gene Genealogies. *Proc. Natl. Acad. Sci. USA* **1998**, *95*, 2044–2049, doi:10.1073/pnas.95.5.2044.

33. Liu, X.; Xing, M.; Kong, C.; Fang, Z.; Yang, L.; Zhang, Y.; Wang, Y.; Ling, J.; Yang, Y.; Lv, H. Genetic Diversity, Virulence, Race Profiling, and Comparative Genomic Analysis of the Fusarium Oxysporum f. Sp. Conglutinans Strains Infecting Cabbages in China. *Front. Microbiol.* **2019**, *10*, 1373, doi:10.3389/fmicb.2019.01373.

34. Glass, N.L.; Donaldson, G.C. Development of Primer Sets Designed for Use with the PCR to Amplify Conserved Genes from Filamentous Ascomycetes. *Appl. Environ. Microbiol.* **1995**, *61*, 1323–1330, doi:10.1128/AEM.61.4.1323-1330.1995.

35. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification And Direct Sequencing Of Fungal Ribosomal Rna Genes For Phylogenetics. In *PCR Protocols*; Academic Press: Cambridge, MA, USA, 1990; pp. 315–322 ISBN 978-0-12-372180-8.

36. Gardes, M.; Bruns, T.D. ITS Primers with Enhanced Specificity for Basidiomycetes—Application to the Identification of Mycorrhizae and Rusts. *Mol. Ecol.* **1993**, *2*, 113–118, doi:10.1111/j.1365-294X.1993.tb00005.x.
37. Sheffield, V.C.; Cox, D.R.; Lerman, L.S.; Myers, R.M. Attachment of a 40-Base-Pair G + C-Rich Sequence (GC-Clamp) to Genomic DNA Fragments by the Polymerase Chain Reaction Results in Improved Detection of Single-Base Changes. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 232–236, doi:10.1073/pnas.86.1.232.

38. Hamdouche, Y.; Meile, J.C.; Nganou, D.N.; Durand, N.; Teyssier, C.; Montet, D. Discrimination of Post-Harvest Coffee Processing Methods by Microbial Ecology Analyses. *Food Control* **2016**, *65*, 112–120, doi:10.1016/j.foodcont.2016.01.022.

39. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549, doi:10.1093/molbev/msy096.

40. Edgar, R.C. MUSCLE: Multiple Sequence Alignment with High Accuracy and High Throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797, doi:10.1093/nar/gkh340.

41. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. JModelTest 2: More Models, New Heuristics and Parallel Computing. *Nat. Methods* **2012**, *9*, 772–772, doi:10.1038/nmeth.2109.

42. Guindon, S.; Dufayard, J.-F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. *Syst. Biol.* **2010**, *59*, 307–321, doi:10.1093/sysbio/syq010.

43. Paradis, E.; Schliep, K. Ape 5.0: An Environment for Modern Phylogenetics and Evolutionary Analyses in R. *Bioinformatics* **2019**, *35*, 526–528, doi:10.1093/bioinformatics/bty633.

44. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. *J. Stat. Softw.* **2008**, *25*, doi:10.18637/jss.v025.i01.

45. Kassambara, A.; Mundt, F. *Factoextra*; 2020;

46. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O'Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. *Vegan: Community Ecology Package*; 2020.

47. Veech, J.A. A Probabilistic Model for Analysing Species Co-Occurrence: Probabilistic Model. *Glob. Ecol. Biogeogr.* **2013**, *22*, 252–260, doi:10.1111/j.1466-8238.2012.00789.x.

48. Griffith, D.M.; Veech, J.A.; Marsh, C.J. Cooccur: Probabilistic Species Co-Occurrence Analysis in R. *J. Stat. Softw.* **2016**, *69*, doi:10.18637/jss.v069.c02.

49. Barral, B.; Chillet, M.; Doizy, A.; Grassi, M.; Ragot, L.; Léchaudel, M.; Durand, N.; Rose, L.J.; Viljoen, A.; Schorr-Galindo, S. Diversity and Toxigenicity of Fungi That Cause Pineapple Fruitlet Core Rot. *Toxins* **2020**, *12*, 339, doi:10.3390/toxins12050339.

50. Wigmann, É.F.; Behr, J.; Vogel, R.F.; Niessen, L. MALDI-TOF MS Fingerprinting for Identification and Differentiation of Species within the Fusarium Fujikuroi Species Complex. *Appl. Microbiol. Biotechnol.* **2019**, *103*, 5323–5337, doi:10.1007/s00253-019-09794-z.

51. Yilmaz, N.; Houbraken, J.; Hoekstra, E.S.; Frisvad, J.C.; Visagie, C.M.; Samson, R.A. Delimitation and Characterisation of Talaromyces Purpurogenus and Related Species. *Persoonia Mol. Phylogeny Evol. Fungi* **2012**, *29*, 39–54, doi:10.3767/003158512X659500.

52. Oxenham, B. Etiology of Fruitlet Core Rot of Pineapple in Queensland. *Qld. J. Agric. Sci.* **1962**, *19*, 27–31.

53. Barral, B.; Chillet, M.; Léchaudel, M.; Larlaud, M.; Verdeil, J.-L.; Conéjéro, G.; Schorr-Galindo, S. An Imaging Approach to Identify Mechanisms of Resistance to Pineapple Fruitlet Core Rot. *Front. Plant Sci.* **2019**, *10*, 1065, doi:10.3389/fpls.2019.01065.

54. Leneveu-Jenvrin, C.; Quentin, B.; Assemat, S.; Hoarau, M.; Meile, J.-C.; Remize, F. Changes of Quality of Minimally-Processed Pineapple (Ananas Comosus, Var. 'Queen Victoria') during Cold Storage: Fungi in the Leading Role. *Microorganisms* **2020**, *8*, 185, doi:10.3390/microorganisms8020185.

55. Ou, Y.; Penton, C.R.; Geisen, S.; Shen, Z.; Sun, Y.; Lv, N.; Wang, B.; Ruan, Y.; Xiong, W.; Li, R.; et al. Deciphering Underlying Drivers of Disease Suppressiveness Against Pathogenic Fusarium Oxysporum. *Front. Microbiol.* **2019**, *10*, 2535, doi:10.3389/fmicb.2019.02535.

56. Lima, C.S.; Pfenning, L.H.; Costa, S.S.; Campos, M.A.; Leslie, J.F. A New *Fusarium* Lineage within the *Gibberella Fujikuroi* Species Complex Is the Main Causal Agent of Mango Malformation Disease in Brazil. *Plant Pathol.* **2009**, *58*, 33–42, doi:10.1111/j.1365-3059.2008.01946.x.

57. Zhan, R.-L.; Yang, S.-J.; Ho, H.-H.; Liu, F.; Zhao, Y.-L.; Chang, J.-M.; He, Y.-B. Mango Malformation Disease in South China Caused by Fusarium Proliferatum: Mango Abnormal Growth. *J. Phytopathol.* **2010**, *158*, 721–725, doi:10.1111/j.1439-0434.2010.01688.x.

58. García-Guzmán, G.; Heil, M. Life Histories of Hosts and Pathogens Predict Patterns in Tropical Fungal Plant Diseases. *New Phytol.* **2014**, *201*, 1106–1120, doi:10.1111/nph.12562.

59. Chen, L.; Zhou, J.; Zeng, T.; Miao, Y.; Mei, L.; Yao, G.; Fang, K.; Dong, X.; Sha, T.; Yang, M.; et al. Quantifying the Sharing of Foliar Fungal Pathogens by the Invasive Plant *Ageratina Adenophora* and Its Neighbours. *New Phytol.* **2020**, *227*, 1493–1504, doi:10.1111/nph.16624.

60. Parnell, S.; Gottwald, T.R.; Gilligan, C.A.; Cunniffe, N.J.; van den Bosch, F. The Effect of Landscape Pattern on the Optimal Eradication Zone of an Invading Epidemic. *Phytopathology* **2010**, *100*, 638–644, doi:10.1094/PHYTO-100-7-0638.

61. Stępień, Ł.; Koczyk, G.; Waśkiewicz, A. Diversity of Fusarium Species and Mycotoxins Contaminating Pineapple. *J. Appl. Genet.* **2013**, *54*, 367–380, doi:10.1007/s13353-013-0146-0.

62. Ibrahim, N.F.; Mohd, M.H.; Mohamed Nor, N.M.I.; Zakaria, L. Characterization of *Fusarium* Spp. Associated with Pineapple Fruit Rot and Leaf Spot in Peninsular Malaysia. *J. Phytopathol.* **2017**, *165*, 718–726, doi:10.1111/jph.12611.
63. Leneveu-Jenvrin, C.; Quentin, B.; Assemat, S.; Remize, F. Maintaining Physicochemical, Microbiological, and Sensory Quality of Pineapple Juice (Ananas Comosus, Var. 'Queen Victoria') through Mild Heat Treatment. *Processes* **2020**, *8*, 1186, doi:10.3390/pr8091186.

64. Rodriguez Estrada, A.E.; Hegeman, A.; Corby Kistler, H.; May, G. In Vitro Interactions between Fusarium Verticillioides and Ustilago Maydis through Real-Time PCR and Metabolic Profiling. *Fungal Genet. Biol.* **2011**, *48*, 874–885, doi:10.1016/j.fgb.2011.06.006.

65. Hynes, J.; Müller, C.T.; Jones, T.H.; Boddy, L. Changes in Volatile Production During the Course of Fungal Mycelial Interactions Between Hypholoma Fasciculare and Resinicium Bicolor. *J. Chem. Ecol.* **2006**, *33*, 43–57, doi:10.1007/s10886-006-9209-6.

66. Combès, A.; Ndoye, I.; Bance, C.; Bruzaud, J.; Djediat, C.; Dupont, J.; Nay, B.; Prado, S. Chemical Communication between the Endophytic Fungus Paraconiothyrium Variabile and the Phytopathogen Fusarium Oxysporum. *PLoS ONE* **2012**, *7*, e47313, doi:10.1371/journal.pone.0047313.

67. Marmann, A.; Aly, A.; Lin, W.; Wang, B.; Proksch, P. Co-Cultivation—A Powerful Emerging Tool for Enhancing the Chemical Diversity of Microorganisms. *Mar. Drugs* **2014**, *12*, 1043–1065, doi:10.3390/md12021043.

68. Barriuso, J.; Hogan, D.A.; Keshavarz, T.; Martínez, M.J. Role of Quorum Sensing and Chemical Communication in Fungal Biotechnology and Pathogenesis. *FEMS Microbiol. Rev.* **2018**, *42*, 627–638, doi:10.1093/femsre/fuy022.

69. Lu, Z.; Tombolini, R.; Woo, S.; Zeilinger, S.; Lorito, M.; Jansson, J.K. In Vivo Study of Trichoderma-Pathogen-Plant Interactions, Using Constitutive and Inducible Green Fluorescent Protein Reporter Systems. *Appl. Environ. Microbiol.* **2004**, *70*, 3073–3081, doi:10.1128/AEM.70.5.3073-3081.2004.

70. Lagashetti, A.C.; Dufossé, L.; Singh, S.K.; Singh, P.N. Fungal Pigments and Their Prospects in Different Industries. *Microorganisms* **2019**, *7*, 604, doi:10.3390/microorganisms7120604.

71. Osmanova, N.; Schultze, W.; Ayoub, N. Azaphilones: A Class of Fungal Metabolites with Diverse Biological Activities. *Phytochem. Rev.* **2010**, *9*, 315–342, doi:10.1007/s11101-010-9171-3.

72. Wei, Z.; Yang, T.; Friman, V.-P.; Xu, Y.; Shen, Q.; Jousset, A. Trophic Network Architecture of Root-Associated Bacterial Communities Determines Pathogen Invasion and Plant Health. *Nat. Commun.* **2015**, *6*, 8413, doi:10.1038/ncomms9413.

73. Schouten, A.; van den Berg, G.; Edel-Hermann, V.; Steinberg, C.; Gautheron, N.; Alabouvette, C.; de Vos, C.H. (Ric); Lemanceau, P.; Raaijmakers, J.M. Defense Responses of *Fusarium Oxysporum* to 2,4-Diacetylphloroglucinol, a Broad-Spectrum Antibiotic Produced by *Pseudomonas Fluorescens*. *Mol. Plant-Microbe Interact.* **2004**, *17*, 1201–1211, doi:10.1094/MPMI.2004.17.11.1201.

74. Santos, C.; Ventura, J.A.; Costa, H.; Fernandes, P.M.B.; Lima, N. MALDI-TOF MS to Identify the Pineapple Pathogen Fusarium Guttiforme and Its Antagonist Trichoderma Asperellum on Decayed Pineapple. *Trop. Plant Pathol.* **2015**, *40*, 227–232, doi:10.1007/s40858-015-0027-7.

75. Crandall, S.G.; Gilbert, G.S. Meteorological Factors Associated with Abundance of Airborne Fungal Spores over Natural Vegetation. *Atmos. Environ.* **2017**, *162*, 87–99, doi:10.1016/j.atmosenv.2017.05.018.

76. Dorey, E.; Fournier, P.; Léchaudel, M.; Tixier, P. Modeling Sugar Content of Pineapple under Agro-Climatic Conditions on Reunion Island. *Eur. J. Agron.* **2016**, *73*, 64–72, doi:10.1016/j.eja.2015.10.010.

77. Zheng, Q.; Ploetz, R. Genetic Diversity in the Mango Malformation Pathogen and Development of a PCR Assay: *PCR Assay for Mango Malformation Pathogen*. *Plant Pathol.* **2002**, *51*, 208–216, doi:10.1046/j.1365-3059.2002.00677.x.

78. Lin, Z.; Xu, S.; Que, Y.; Wang, J.; Comstock, J.C.; Wei, J.; McCord, P.H.; Chen, B.; Chen, R.; Zhang, M. Species-Specific Detection and Identification of Fusarium Species Complex, the Causal Agent of Sugarcane Pokkah Boeng in China. *PLoS ONE* **2014**, *9*, e104195, doi:10.1371/journal.pone.0104195.

79. Filières Végétales—DAAF de La Réunion Available online: https://daaf.reunion.agriculture.gouv.fr/Filieres-Vectetales (accessed on 11 December 2020).

80. Smart Information System (SmartIS)-VICTORIA-CIRAD-UR AÏDA Available online: https://smartis.re/VICTORIA (accessed on 23 October 2020).

81. Dorey, E.; Fournier, P.; Lechaudel, M.; Tixier, P. Validity of the Pineapple Crop Model SIMPIÑA across the Climatic Gradient in Réunion Island. *Eur. J. Agron.* **2015**, *62*, 1–12, doi:10.1016/j.eja.2014.09.004.