1. Introduction
Banana is a staple food for most tropical and subtropical regions worldwide, considered one of the most relevant agricultural activities in developing countries with economic and social impact especially in Asia, Africa, Latin America and the Caribbean, Ásia, África [
1]. There are more than 50 banana subgroups of which the “Cavendish” subgroup accounts for approximately 90% of the international fruit trade, followed by plantains with 21% of the market share [
2,
3,
4]. In Brazil, approximately 60% of the banana production is based on the “Prata” subgroup, with cultivars ‘Prata-Anã’ and ‘Pacovan’ being the most representative. In the South and Southeastern Regions of the Country, the cultivars of the “Cavendish” subgroup have a well-established consumer market, whereas the cultivation of plantains is mostly appreciated in the North and Northeast Regions [
5].
One of the main limitations on current banana production is
Fusarium wilt, considered the most destructive disease of the crop, and it is caused by the fungus
Fusarium oxysporum f. sp.
cubense (Foc). In addition to the “
forma specialis” classification where each isolate has a unique and normally related host range, Foc isolates can be also classified according to its race, meaning, its ability to infect particular bananas cultivars. Thus, Foc is classified into three physiological races according to its virulence to the following cultivars: ‘race 1’: isolates virulent to the cultivars ‘Gros Michel’ and ‘Maçã’ (syn: Silk); race 2: isolates virulent to ‘Bluggoe’ and other cooking bananas (plantains); race 4: isolates virulent to cultivars of the “Cavendish” subgroup [
3,
6,
7,
8]. Races 1 and 2 are present in Brazil, with greater relevance for ‘race 1’. ‘Race 3’, is not a pathogen of
Musa spp., and has ambiguity in its taxonomy since this race causes symptoms in heliconias, a wild relative of bananas [
9].
Foc ‘race 4’ is subdivided into subtropical (ST4) and tropical (TR4) following the usual nomenclature adopted in recent research and official database of nomenclature of the fungus. ST4 causes production losses in subtropical regions of Australia, South Africa and the Canary Islands, commonly associated with plantations exposed to abiotic stresses [
3]. TR4 is found in many countries in Asia, such as India, China, Israel, Jordan, Laos, Lebanon, Malaysia, Myanmar, Pakistan, Philippines, Sulawesi, Sumatra, Vietnam, Turkey, Mayotte, Taiwan and Indonesia [
3,
8,
10,
11,
12,
13,
14,
15], in Africa in Mozambique and Oceania in Australia, and Queensland [
16].
TR4 has not yet been reported in Brazil and is considered a quarantine pest by the Brazilian Ministry of Agriculture and Livestock. However, the disease was officially identified in Latin America in August 2019 in Colombia, in April 2021 in Peru and in early 2023 in Venezuela [
17,
18,
19,
20].
For Foc disease, the classification into physiological races does not correlate with the genetic variability and also with existing lineages and phylogenetic groups, but is extremely important for management purposes [
15]. In an attempt to obtain a reliable molecular diagnosis capable of distinguishing several races of the pathogen, a family of effector genes known as “SIX” (secreted in xylem), which are reported as virulence factors, previously used successfully to differentiate the three races of
F. oxysporum f. sp.
lycopersici (Fol) x tomato (
Solanum lycopersicum) [
21], were also used to distinguish Foc isolates in bananas [
22,
23,
24,
25].
In these studies, the SIX7 and SIX8 homologues were detected only in the Foc race 4 isolates (TR4 and ST4), while Foc-SIX8a was present in all TR4 isolates, Foc-SIX8b was present in all ST4 isolates, but was not detected in TR4 isolates [
22,
23]. Likewise, a set of assays was able to distinguish relevant Foc races and VCGs, in which sets of primers specifically amplified regions of the SIX6 gene in the Foc race 1, SIX1 gene in TR4, SIX8 in the subtropical race 4, SIX9 / SIX10 in the Foc VCG 0121 and SIX13 in the Foc VCG 0122 [
25]. To investigate the evolution of SIX genes, whole genome sequencing data generated for 23 Foc genetic lineages were queried for 14 known SIX genes (SIX1 – SIX14), where SIX gene profile variation, including the presence of SIX-specific homologues, correlated with race structure and evidence of horizontal transmission of SIX genes in Foc, was identified [
24].
In Brazil, the genetic variability of Foc was studied by means of a representative survey of Foc isolates from the main banana production regions and their levels of variability and aggressiveness, estimated; in addition, the first detailed map of the distribution of Foc was elaborated demonstrating the diversity among isolates, especially as to their distribution and ability to infect different banana cultivars [
26]. In our work, the 214 isolates analyzed by SSR markers were separated into 52 distinct haplotypes and distributed in different banana-producing regions of Brazil, but no Foc population structure based on host or geographic regions, was found. Recently, Batista et al. [
27] quantified the genetic variability present in Foc isolates from the Brazilian population to investigate the presence of TR4 and to determine the genetic structure based on the SIX gene profile, VCG tests, IGS PCR-RFLP and SSR markers. In their study, the authors showed that Brazil is, in fact, free of Foc TR4 isolates. In addition, this study identified a high diversity of VCG compared to some Asian countries and a trend of association between the profile of SIX, clades and VCG genes.
In the present study, new collections of Foc isolates were carried out in “Cavendish” and “Prata-Anã” production areas, especially in the South and Southeast Regions of Brazil. The isolates collected in the current study together with isolates belonging to the collection of Embrapa Cassava and were investigated for the presence of putative SIX pathogenicity effectors using PCR analysis. In addition, among the isolates evaluated, two highly aggressive pathotypes of Foc with virulence to cultivar ‘Grande Naine’, were selected, which will be used as a new standard for the selection of banana genotypes resistant to Fusarium wilt by the Embrapa banana genetic breeding program.
2. Materials and Methods
2.1. Foc isolates
The regions of São Paulo and Santa Catarina were selected because they are large banana producers in Brazil. After empirical research at the site, an increase in the incidence of Fusarium wilt was found in cultivars of the ‘Cavendish’ subgroup, previously unaffected by the disease. Thus, the cultivars ‘Cavendish’ and ‘Prata-Anã’ that showed symptoms of Fusarium wilt in the field were selected for sampling.
The description of sampling and conservation of the isolates followed the same protocols performed by Costa et al. [
26]. At least three samples were collected per field in four plantations sampled by municipality in each state (
Table S1). For indirect isolation, fragments of stems infected with Foc (pseudoestemics) were disinfected and then deposited in potato dextrose agar (PDA) culture medium. The isolations were incubated at 25º C for 8 days. At the end, the isolates obtained were cataloged and preserved in sterile distilled water and also on filter paper with silica gel and became part of Embrapa’s collection of Foc isolates.
The isolate Foc 0801, characterized as Race 1, previously used as a standard of aggressiveness by Embrapa’s banana genetic improvement program for the selection of genotypes resistant to Foc [
28,
29], and the isolate Foc 218A, characterized as Subtropical race 4 (ST4) [
30] were selected for this study (
Table S1).
2.2. Bioassay
The virulence/aggressiveness trials were conducted at the Plant Phytopathology Laboratory and under greenhouse conditions at Embrapa Mandioca and Fruticultura, located in Cruz das Almas, Bahia, Brazil. For the inoculation trials two cultivars with known behavior to Foc “race 1” were used, ‘Grand Naine’ (Cavendish), (resistant) and ‘Prata-Anã’ (Prata), (moderately susceptible). The plantlets were obtained by in vitro culture and acclimated on commercial substrate and coconut fiber (5:1) for 60 days. After this period, they were transferred to polyethylene bags with dimensions of 20 x 30 cm until reaching 50 cm in height for the bioassay.
For the inoculation of plantlets of Grand Naine and Prata-Anã, the 52 isolates of Foc were individually grown in substrate sand/cornmeal. To produce the substrate, washed sand and cornmeal (5:1 ratio) were used in addition of 150 mL of distilled water for each 300 g of substrate. The mixture was added to plastic bags and autoclaved at 120 °C for 20 minutes. Foc mycelia discs were transferred to the substrate and incubated in a growth chamber at 25 °C for 30 days. Serial dilutions were performed to quantify the colonies on the substrate. The dilutions 104 and 105 were plated in BDA culture medium and after 48 hours the colony forming units were counted. Subsequently, the inoculum was adjusted to 106 colony-forming units per gram of substrate (UFC g-1). For the inoculation, 10 grams of the inoculum were distributed in four holes around the plantlets at a depth of 10 cm. The control treatment was performed with 10 g of the sand/cornmeal per orifice, without the presence of Foc [
31].
A trial of completely randomized design with factorial scheme (52 x 2) was used for the greenhouse experiment, totaling 104 treatments, being 52 isolates of Foc and two banana genotypes (‘Grand Naine’ and ‘Prata-Anã’). For each isolate, 5 replicates were used. The incubation period was considered as the time elapsed between the inoculation and the appearance of the first external symptoms, at which time the disease severity was evaluated at intervals of seven days up to 85 days after inoculation (DAI).
For the assessment of the external symptoms, the following scale was used according to Dita et al. [
32]: (0) no external symptoms; (1) initial yellowing on old leaves; (2) yellowing of old leaves with initial discoloration on young leaves; (3) all leaves with intense yellowing and (4) dead plant. At 85 DAI, the plants were removed from the substrate and evaluated for internal symptoms by cross-section of the rhizome, following the scale proposed by Dita et al. [
32]: (0) no symptoms; (1) initial rhizome necrosis; (2) mild necrosis on the rhizome; (3) intense necrosis symptoms on the rhizome; (4) rhizome with most internal tissues showing necrosis; and (5) totally necrotic rhizome.
The disease index (DI) for external symptoms (EDI), and internal (IDI), was calculated according to the formula proposed by McKinney [
33], being DI (%) = 100.∑[(f.v)/(n.x)], whereas, f is the number of plants with the same score; v is the observed score; n is the number of plants evaluated; and x is the maximum score from the scale.
2.3. Root bleaching and coloring
At 85 days after inoculation, after assessing the internal symptoms of the disease, the root fragments were collected and washed and taken immediately to the laboratory to begin the whitening process. The root clarification technique was used to analyze ‘Grand Naine’ plants inoculated with isolates 218A and 0801. Root clarification and staining techniques were used according to Phillips and Haymann [
34] and Brundrett et al. [
35]. Clarification was performed by immersing the root system cuttings of approximately 2 cm length in a 10% KOH solution for 48 hours. Afterwards, the KOH was discarded, washed with water to remove all solution, and the samples neutralized in 1% HCl for 5 minutes. The staining was conducted with trypan blue in 0.05% solution of lactoglycerol (2:1:1 lactic acid, glycerol, water) at room temperature for 1 hour. After this period, the dye was discarded and the roots stored in lactoglycerol solution until observation under light microscope.
2.4. DNA extraction
At The Foc isolates used in this study were grown in BDA medium for seven days. After growth, the mycelium was collected, dried at room temperature and macerated in liquid nitrogen into a fine powder using a mortar and pestle. The DNA was extracted according to the methodology described by Doyle and Doyle [
36] with modifications by addition of 700 μL of CTAB extraction buffer (1% CTAB, 0.7 mM NaCl, 20 mM Tris-HCl, 10 mM EDTA, 1% 2-mercaptoethanol) previously heated at 60-65 ° C for 45 min.
The samples were then incubated in a water bath (60-65 °C) for one hour, vortexed every 10 minutes. Subsequently, 500 μL of the chloroform: alcohol isoamylic solution (24: 1) was added, centrifuged at 10,000 rpm for sedimentation for 10 minutes, and the supernatant transferred to new microtubes containing 500 μL isopropanol and 200 μL 3M sodium acetate. The DNA was pelleted by centrifugation on 10,000 rpm for 12 minutes, followed by two washes with 500 µL of ethanol (70%). The pellets were dried at room temperature for 30 minutes and resuspended in 100 µL of TE buffer (Tris-HCl 10 mM, 1 mMEDTA, pH 8 0) with RNase A (Invitrogen) and incubated for at 37 °C for 30 min. The amount of gDNA was checked by 1% agarose gel electrophoresis, using ethidium bromide and UV-light.
The DNA of isolate 701 of Tropical Race 4 (TR4) was kindly granted by Plant Pathology Department at the Universidade Federal de Viçosa (UFV), Viçosa-MG, Brazil, and was used as a positive control in this study.
2.5. PCR analysis
The 52 Foc isolates from this study were used for PCR analysis using the SIX1-F/R (AJ608702), SIX2 (GQ268949), SIX3 (AM234063), SIX4 (HQ260602), SIX5 (FJ767863), SIX6 (FJ755835), SIX7 (FJ755836), SIX8 (FJ755837) primer set according to SIX gene data published in GenBank [
22] (
Table 1). Additionally, the new set of primers for the SIX8 homologs, being SIX8a-F/R (KF548063) and SIX8b-F/R (KF58064) was used [
23]. A set of primers which amplify a region of the intergenic spacer (IGS) region of the nuclear ribosomal operon generated by Dita et al. [
37] to specifically detect Tropical Race 4 (TR4) was also used. In addition, primers Foc - 1 / Foc - 2 (EU379562) developed from a Random Amplified Polymorphic (RAPD) DNA marker reported for specific detection of the breed Foc 4 has been tested [
38]. To evaluate the quality of the Foc DNA, the translation elongation factor 1a (TEF-1a) gene was also used with the primers EF1 and EF2 [
39] (
Table 1).
PCR parameters were as follows: denaturation at 94 °C for 10 min, followed by 35 cycles of denaturation at 94 °C for 45 s; annealing at 60 °C for 45 s; and extension at 72 °C for 1 min, with a final extension at 72 °C for 10 min in a thermocycler (Applied Biosystems). The final volume of the reactions was 25 μL consisting of 1 x GoTaq DNA Polymerase (Promega, Corporation, Madison, WI), 10 nmoles dNTPs (Promega), 4 μM of each primer and 10 ng of DNA. The resulting PCR products were visualized on 1% agarose using ethidium bromide and UV-light. For PCR with the TR4-specific primer the following program was used: 95°C for 2 min and 35 cycles of 95°C for 30 s, 60°C for 30 s and 72°C for 1 min, followed by an additional extension time for 10 min at 72°C [
37].
2.5. Data analysis
The disease index data were submitted to analysis of variance at the 5% significance level and for the factors that presented significant differences, the interactions unfolded by the Skott-Knott test at 5% of significance. The heat maps were developed using the ‘devtools’ and ‘gplots’ packages in the R software [
40].
A derivative analysis of PCA was used for the comparison of genotype vs genotype (G x G) interaction, and performed to explain the variation of the isolates in the cultivars ‘Prata-Anã’ and ‘Grand-Naine’, using the ‘Ggebiplot’ package in the R software. The isolates were compared based on the band profile (bp – base pair) for the primers related to the genes. For the construction of the dendrogram a binary matrix was established for the presence (1) and absence (0) of the gene.
The correlation among the different parameters such as internal and external disease indexes (IDI and EDI) and the presence of band patterns associated with presence of different SIX genes was evaluated based on the Multiple Factorial Analysis (MFA), for the mixed analysis of qualitative and quantitative data obtained using the ‘FactoMiner ’ package implemented in R software.
4. Discussion
According to a study of a population of Foc isolates in Brazil carried out by Costa et al. [
26], the population structure of Foc was characterized by high genetic variability. Although a sexual phase is not known in the reproductive structure of the
Fusarium oxysporum species, other evolutionary mechanisms may be related, leading to the emergence of more aggressive variants capable of infecting resistant cultivars [
26]. In addition to the existence of physiological races and vegetative compatibility groups, the occurrence of cryptic recombination and the presence of mobile chromosomes that allow selective pathogenicity, have been proven [
42]. In the present study, some isolates of Foc were virulent when cultivating ‘Grand Naine’, belonging to the Cavendish subgroup and a significant number were extremely aggressive in the cultivation of ‘Prata-Anã’, while isolate 218A was more aggressive in both cultivars, which reinforces the idea that this is a new Foc pathotype in Brazil [
26,
28,
30,
43].
The phylogenetic basis which explains the variations in the genetic structure of Foc regarding the pathogenic specificity has been associated with horizontal gene transfer [
24,
39,
44]. The data suggests that chromosome 14 of lineage-specific genomic regions (LS) harbor genes responsible for pathogenicity that are dispensable and can be transferred laterally between pathogens, including several polyphyletic clonal lines within the species complex [
42], although, to date, there is no published evidence that chromosome 14 exists in Foc. Recent studies indicate that many genes are required for general pathogenicity in relation to plants and are widely distributed throughout the
F. oxysporum species complex, but the specific pathogenicity of the Foc host can be attributed to fungal effectors such as SIX genes [
25,
45,
46]. Furthermore, a study has stated that the presence of the SIX1a allele is linked to the virulence of Foc ‘TR4’ in bananas from the Cavendish subgroup based on mutagenesis analyzes [
47]. Other studies reported that the effector gene SIX8 is required for virulence of Foc ‘TR4’ to Cavendish banana [
48].
In this study, the SIX1 gene was present in all tested isolates, including avirulent isolates for both cultivars (
Figure S1). In a similar study Czislowski et al. [
24] identified the diversity of the same gene and found nine types of sequences in total, even though their role in pathogenicity is still unknown. Meldrum et al. al. [
22] and Fraser-Smith et al. [
23] assessing the presence of putative genes in pathogenic Foc isolates showed that the SIX7 and SIX8 genes were present only in race 4 isolates and absent in race 1 and 2; similarly, the homologous gene SIX8 differ - (TR4) from subtropical race 4 (ST4) so that the SIX8a gene was present only in TR4 isolates and the SIX8b gene only in ST4 isolates. Furthermore, a comparison of SIX8a with SIX8b revealed 95% identity. However, in this study, new perspectives on the role of SIX genes in Foc virulence, were noticed.
PCR amplification with the set of primers SIX7, SIX8, SIX8a and SIX8b showed a band for most of the evaluated isolates, but there was no presence of amplification of bands for the set of specific primers of TR4. Based on PCR reactions, we confirmed the absence of TR4 strains in our collection from different regions of Brazil. In addition, the results confront the proposition that ‘Foc-SIX8’ is a suitable candidate for the molecular differentiation of isolates of race ‘4’ among races ‘1’ and ‘2’, as well as the differentiation of the tropical race 4 (TR4) of subtropical race 4 (ST4) as stated by Meldrum et al. [
22], Fraser-Smith et al. [
23] and An et al. [
48]. Thus, it is assumed that copies of the SIX8 gene may perform different virulence functions in these isolates, since unlike other SIX genes, which are found as single copies, there are nine copies of SIX8 and four copies of SIX8b in the genome of Fol 4287 [
22,
42]. In general, one must take into account the non-specificity of the primers used to evaluate these genes.
The presence of association between virulence and/or aggressiveness parameters, and the specific presence of some of the SIX genes, were noticed in this study, since the presence of the SIX8b allele was strongly related to the behavior of avirulence in the cultivar ‘Grand Naine’ (Cavendish subgroup) (
Figure 6). It therefore may be inferred that this gene possibly performs functions related to specificity in relation to the host.
The understanding of the host-pathogen interaction is essential to explain the aspects of variation in the virulence and aggressiveness of the isolates in the cultivars ‘Prata-Anã’ and ’Grand Naine’ in this study. The presence of both the virulent isolate 218A and the avirulent isolate 0801, were thus present within the tissues, but the structures produced by an isolate in which the cultivar is resistant were predominantly chlamydospores, which are structures of resistance, indicating that barriers have been imposed to prevent the advancement of hyphae and colonization of the tissue and therefore, the pathogen needed to produce resistance structures. The formation of chlamydospores is correlated with the survival capacity of the fungus as a result of plant defense. This intense production of chlamydospores occurs as defense of the pathogen in relation to the plant, since these structures can guarantee the survival of the fungus for many years [
3,
8,
49].
In another study, the infection process of two commercial Australian cotton varieties differed in susceptibility to
Fusarium oxysporum f. sp.
vasinfectum (Fov) and the appearance of Fov chlamydospores in the tissues of the susceptible strain was observed, which was associated with significant degradation of the host tissue [
50]. Likewise, work with a Foc ST4 line labeled with Green Fluorescent Protein (GFP) verified the path of colonization after infection of susceptible cultivars and results showed the production of chlamydospores as soon as the host was invaded and prior to expression of external symptoms [
51,
52].
Overcoming of the genetic resistance is a constant threat to field management of Fusarium wilt and this problem can be closely related to the high variability of the fungus [
26,
27,
53,
54,
55]. According to Buddenhagen [
11], tropical race 4 (TR4) evolution was associated with the emergence of the ability to cause disease in banana cultivars of the Cavendish subgroup from spontaneous mutations that generated variants with potential to colonize and cause damage to this subgroup. Thus, the different strains of the pathogen can adapt to different conditions, especially during the extensive use of resistant cultivars since the presence of resistance genes in the plant lead to a directional seal in populations of pathogens, culminating in the increase in the frequency of individuals capable of overcoming plant resistance [
26,
56].
Based on our findings, it is possible to infer that the most aggressive isolates, identified in the cultivar ‘Grand Naine’, probably are new variants of race 1 isolates, a fact that requires extensive studies, in order to confirm this hypothesis. Similar results were found by Ribeiro et al. [
28] and Rocha et al. [
30] for the diversity among Foc isolates, where the 218A isolate was the most virulent, being selected as representative of most aggressiveness of the analyzed isolates compared to the standard ‘0801 isolate’. Rocha et al. [
43], evaluating a set of genotypes through histological and histochemical analyses, and gene expression, observed that the isolate 229A is associated with the suppression of defense responses in cultivars resistant to isolates 218A and 0801, indicating it is therefore more virulent.
Among the isolates evaluated, there was an extensive variation considering virulence and aggressiveness in the cultivars ‘Prata-Anã’ and ‘Grand Naine’ (
Figure 2). In a similar study, Cunha et al. [
57] collected Foc isolates in banana production areas in the State of Santa Catarina with cultivars of the Cavendish and Prata subgroups, aiming to estimate the diversity of the pathogen and stated that 78% of the isolates considered as more aggressive were collected in ‘Grand Naine’ (Cavendish subgroup). In the present study it is also shown that the majority of more aggressive/virulent isolates (68%) were also collected in Cavendish cultivars, especially in the State of São Paulo (
Table S1).
In addition to the emergence of new and more aggressive variants of
Fusarium oxysporum f. sp.
cubense (Foc), adverse environmental conditions also contribute to the incidence of Foc in banana-producing areas in Brazil. In the latest study on the population structure of Foc in Brazil, the presence of TR4 was investigated in a collection of isolates, and the genetic structure was determined based on an expanded set of markers [
27]. In this study, it was found that the Foc population appears to have low to moderate genetic variability, given the high clonal fraction.
The presented data reinforce the existence of variations in the virulence profile of some isolates and may be associated with the occurrence and widespread distribution of ST4 among
Fusarium oxysporum f. sp.
cubense (Foc) isolates in Brazil under favorable environmental conditions in different regions of the country [
27]. Furthermore, these data bring new perspectives on the characterization of SIX genes and their role in pathogenicity of Foc in bananas. Additional studies need to be carried out aiming at the more complete identification of these isolates, as well as a more accurate definition of the races, and the presence of pathogenicity and avirulence genes, for example via sequencing of the SIX regions or even sequencing the genome of these isolates.
Author Contributions
Conceptualization, K.V.C.V., A.d.J.R. and E.P.A.; methodology, A.d.J.R.; K.V.C.V., and M.d.S.F.; software, K.V.C.V., A.d.J.R., and M.d.S.F.; validation, K.V.C.V., A.d.J.R., M.d.S.F., K.N.P., V.B.d.O.A., S.A.S.d.O., C.F.F., F.H. and E.P.A.; formal analysis, A.d.J.R.; investigation, K.V.C.V., A.d.J.R., F.H., C.F.F., S.A.S.d.O., E.P.A.; resources, E.P.A.; data curation, K.V.C.V., A.d.J.R., M.d.S.F. and V.B.d.O.A.; writing—original draft preparation, K.V.C.V., A.d.J.R., E.P.A.; writing—review and editing, K.V.C.V., A.d.J.R., F.H., C.F.F., S.A.S.d.O., E.P.A.; visualization, K.V.C.V., A.d.J.R., M.d.S.F., V.B.d.O.A., K.N.P., V.B.d.O.A., C.F.F., F.H., S.A.S.d.O and E.P.A.; supervision, E.P.A.; project administration, E.P.A.; funding acquisition, E.P.A. All authors have read and agreed to the published version of the manuscript.