Abstract
Introduction:
Fusarium verticillioides (Fv) is a major phytopathogen responsible for maize root rot, affecting crop productivity globally. A probable infection mechanism has been suggested in Fusarium, involving the disruption, and partial degradation of the plant cell wall by the colonizing fungal hyphae.
Methods:
In this study, the highly virulent Fv DA42 strain was subjected to whole-genome sequencing, in silico secretome analysis and SEM structural analysis to elucidate its pathogenic mechanism.
Results:
The assembled genome comprised 175 contigs (=200 pb) totaling 42.27 Mb, with an N50 of 1.24 Mb and GC content of 48.51%. A total of 14,198 protein-coding genes were predicted, of which 997 (7.03%) correspond to classical secreted proteins. The predicted secretome includes 262 carbohydrate-active enzymes (CAZymes), 62 proteases, 400 effectors, 481 virulence factors and 288 uncharacterized proteins. Functional annotation revealed enrichment in enzymatic activities such as pectinesterases, feruloyl esterases, and glucosidases, highlighting their role in host cell wall degradation. Chromosomal distribution showed secretome genes concentrated on chromosomes 4 and 8, with the highest density (49.2 genes/Mb) on chromosome 10. Scanning electron microscopy confirmed degradation of maize root hairs and epidermis seven days post-infection, this degradation may have occurred days prior to the observation. STRING analysis identified key proteins like FVEG_10795 (pectinesterase) and FVEG_09361 (feruloyl esterase) as central to coordinated enzymatic attacks.
Discussion:
This integrative analysis offers crucial insights into Fv pathogenicity and provides a molecular basis for targeted antifungal strategies and resistance breeding in maize.
1 Introduction
Maize root rot is an important fungal disease, which is a constant menace to maize agriculture worldwide. The main fungal phytopathogen associated with root rot are Fusarium spp. However, the Fusarium incidence statistics in maize rots is limited and it varies by geographic location, climate (humidity and temperature) during the agriculture cycle, agricultural practices, and the maize genetic background (). Additionally, visual identification of root infection is not easy. Reports show prevalence of this phytopathogen in corn field from diverse agriculture regions around the world, for example, Serbia, Croatia (), Portugal (), Poland (), China (), Germany (), and Mexico ().
Fusarium species also cause other types of rot and blight, affecting seedling health, stalk resistance, and grain quality (; ; ). Fusarium maize rots are caused by several species (; ), with the most predominant being F. graminearum Schwabe and F. verticillioides Sacc. Nirenberg (; ; ; ; ; ).
Fusarium verticillioides (Fv) like many fungal pathogens, possesses a remarkable ability to infect maize plants at every developmental stage, from seed and seedling through the productive stage (). The infectious potential of Fusarium is evident in its ability to infiltrate not only in damaged roots but also in stalks, and even developing grains through the stigmas and silk channels (; ; ), suggesting that this phytopathogen could utilize diverse infection strategies (). In soil, the root cell wall is the first barrier that soilborne pathogens, such as Fv, must overcome to infect plants. The mechanism of infection in tomato and pea roots by F. oxysporum occurs mainly through cell-to-cell spaces, lateral root emergence, or cell damage caused by rocks or insects. The hyphae enter by disrupting, and causing partial degradation of the cell wall (). These hyphae form appressorium-like infection structures or penetration pegs. These structure have also been visualized in wheat coleoptiles () and floral organs () colonized by F. graminearum and F. culmorum (), as well as by other fungal necrotrophic phytopathogens, such as Botrytis cinerea () and Rhizoctonia solani (). Using SEM and red or green fluorescent protein tags, F. verticillioides has been shown colonizing maize by penetration of the lateral root breakage zones, and root epidermis through appressorium-like structures (; ; ). During the first days of Fusarium root rot infection, symptoms start as small pale-yellow root discolorations that progress to dark brown color and eventually form necrotic spots (; ). Simultaneously, root hairs collapse, providing a visual cue that suggests events associated with cell wall degradation in the epidermal cells caused possibly by Fv secreted proteins (). This intricate process of invasion and colonization may shed light on the dynamic interplay between Fv and maize, highlighting the importance of understanding the physiological and molecular aspects of this bipartite interaction.
The plant cell wall is a complex structure primarily composed of polysaccharides such as cellulose, hemicellulose, pectin, and lignin as well as proteins and phenolic compounds, such as ferulic acid. These components provide structural support by crosslinking different polysaccharides (). Fungal phytopathogens, such as Fv, secrete an arsenal of lytic enzymes and effector proteins, collectively known as the protein secretome. These proteins act synergistically to degrade cell wall components, which is crucial for nutrient acquisition, fungal colonization, and dispersion within plant tissues (; ). The secretome consists of proteases, lipases, protein effectors, and carbohydrate-active enzymes (CAZymes), such as cellulases, xylanases, and pectinases (). Understanding the mechanisms that Fv and other phytopathogens use to degrade the cell wall and overcome host defenses will provide novel targets for disease control.
Whole-genome sequencing (WGS) has revolutionized fungal pathogen research, enabling advancements in disease diagnostics and the determination of evolutionary relationships between species (). Additionally, WGS offers a comprehensive view of potential genetic factors. Next-generation sequencing (NGS) technologies facilitate the identification of numerous genes, including those involved in pathogenesis and secretome gene products. In silico analysis of the secretome, which predicts secreted gene products from gene sequences, allows for the identification of key lytic enzymes and protein effectors involved in plant cell wall degradation and pathogenic mechanisms. Six Fv genome assemblies have been deposited at NCBI, ranging from contigs to chromosome level, with genome sizes ranging from 41.84 to 44.65 Mb and containing between 15,053 to 20,574 protein-coding genes (; ). Recently, a report was published on the gapless genome assembly of Fv strain 7600 using PacBio HiFi technology. The resulting final genome assembly is 41.994 Mb and contains 15,230 protein-coding genes (). This genome assembly represents an improvement over the first version, which had a size of 41.79 Mb and contained 14,335 protein-coding genes, along with other genome statistics ().
The Fv secretome could plays an important role in plant pathogenesis () and human Fusariosis (), as well as in biotechnological potential uses. Despite the importance of this fungal species, to our knowledge, there are only a few reports on the Fv secretome. One of these reports shows the potential for commercial use of Fv secretome as a complement to commercial cellulases cocktail in the saccharification (lignocellulose hydrolysis) of wheat straw, 166 Fv secreted proteins were identified by proteomic analysis, of which 57 are involved in the degradation of lignocellulosic material, these include CAZymes that degrade cellulose, xylans, pectins, lignins; and other carbohydrate-modifying enzymes (). Another report involved the in silico comparative analysis of the Fv BIONCL4 strain, which identified 2,058 proteins related to the secretome. Of these, 676 corresponded to the classical secretion pathway, while 569 were identified as pathogenesis-related proteins ().
Roots act as the first line of defense against soil pathogens such as F. verticillioides. The pathogen has evolved mechanisms to degrade the root cell wall, allowing its entry and subsequent colonization of the plant. This occurs mainly through the secretion of hydrolytic enzymes and effectors that turn off plant immune system () and break down structural components () such as proteins, cellulose, hemicellulose, lignin and pectin, facilitating its dispersion inside the plant tissues and the establishment of infection. In this study, we used SEM, WGS and in silico secretome analysis to investigate the molecular arsenal of Fv and its significant role in maize root rot. Our aim was to identify potentially key secreted proteins involved in pathogenesis and cell wall degradation. This study could provide the basis for the development of more effective control strategies for the Fv-maize interaction and contribute to the protection of this important crop, as well as other crops affected by Fusarium species.
2 Materials and methods
2.1 Seed infection
To evaluate the response to F. verticillioides (Fv) infection, a maize inbred line highly susceptible (IL09) to root rot was used and previously characterized by and . The virulent Fv DA42 strain, described by was used for infection assays. Seed preparation followed protocols previously described by and , with minor modifications. Seeds were superficially disinfected by sonication for 5 minutes in sterile distilled water with Tween 20 (5 drops/100 mL) using an ultrasonic bath (2.8 L, Fisher Scientific). This was followed by immersion in 1.5% sodium hypochlorite at 52°C for 20 minutes in a thermostatic water bath (FE-377, Felisa), and then rinsed thrice with sterile distilled water under aseptic conditions in a biological safety cabinet (Herasafe KS, Thermo Scientific).
The Fv strain was cultured on Spezieller Nährstoffarmer agar (SNA) medium () supplemented with a 1 cm² filter paper and incubated at 25 ± 2°C for 14 days. Conidia were harvested by adding 5 mL of sterile 0.8% NaCl solution to the culture, followed by gentle agitation. A working conidial suspension was prepared at 1 × 106 conidia/mL, quantified using a Neubauer counting chamber (cat. No. 3110, Hausser Scientific, USA) and a light microscope (B-383-M11, Optika, Italy). Disinfected seeds were immersed for 5 minutes in the conidial suspension, while control seeds were treated with sterile water. Ten seeds were placed every 2 cm along sterile, moistened Kraft paper sheets (40 × 20 cm), which were rolled and placed in plastic bags. Seed germination was carried out at 25°C under a 16:8 h light:dark photoperiod for 7 days. Moisture was maintained by daily irrigation with 15 mL of sterile water. Daily visual observations and photographs were taken using a stereo microscope (M205FA, Leica, Germany).
2.2 Microscopy analyses
Primary root tissues with and without necrosis symptoms were harvested for scanning electronic microscopy (SEM) studies. These samples were treated as previously described by and observed using an FEI-Quanta250 FEG microscope (Czech Republic) operated at an acceleration voltage of 5 kV.
2.3 DNA isolation and sequencing
The Fv DA42 strain culture () was grown in Petri dishes containing potato dextrose broth (Cat. 7041, MCD Lab) and incubated at 25 ± 2°C. After seven days of incubation, the mycelium was harvested by scraping it with a spatula. Approximately 100 mg of mycelium was used for genomic DNA (gDNA) extraction using the CTAB protocol (). The quality of gDNA was evaluated on a 0.8% agarose gel stained with GelRed (cat. 41003, Biotum) and visualized using a Gel Doc XR+ (BioRad). Genomic DNA was quantified with a Nanodrop spectrophotometer (Nanodrop 8000, Thermo Fisher) and stored at -20°C. The DNA samples were prepared according to the NGS library preparation workflow and sequenced using the Illumina platform by Macrogen (whole genome de novo sequencing, paired-end reads with a read length of 151).
2.4 Assembly of raw data and gene and secretome prediction
Raw reads quality was evaluated using FastQC () and trimmed to remove low quality bases using Trimmomatic () with default parameters for paired-end data. These trimmed raw reads were filtered by performing a BLASTn search against the mitochondrial genome of the teleomorph of Fv (Gibberella moniliformis) (NCBI RefSeq: NC_016687.1). High quality reads were then assembled de novo using the Shovill tool () with assemblers, where the Velvet assembler gave the best metrics. The quality of the assembled genome was verified by QUAST (), which provides key metrics such as genome size, N50, L50, the number of contigs, and the largest contig. As a second filter to remove potential mitochondrial contaminants, the contigs were BLASTed against a NR mitochondrial database (https://ftp.ncbi.nlm.nih.gov/genomes/refseq/mitochondrion), and those showing ≥90% identity was manually removed. Gene (nt) and protein (aa) predictions were performed using AUGUSTUS through gene models (), which were built from the Fv 7600 strain (GCA_000149555.1) using Train Augustus. These programs were implemented on the Galaxy USA platform (). The predicted proteins were functionally annotated by eggNOG-mapper databases (). RefSeq protein annotation was supported by local blast (rblast Version 0-99.4) using the GCF_000149555.1_ASM14955v1_protein.faa database with the argument “max_target_seqs 1” in Rstudio version 4.3.3. The UniProt database (), considering Fv 7600 (Taxon ID 334819) was used to support the FVEG gene name of the selected protein dataset.
Subcellular localization was predicted for this dataset using Wolf PSORT (). Signal peptides were predicted using DeepSig (). These programs were implemented on the Galaxy Europe platform (). Transmembrane proteins were predicted using Phobius by forcing the predictor to choose between two submodels (). Proteins lacking glycosylphosphatidylinositol (GPI) were predicted using the PredGPI web tool (). Carbohydrate-active enzyme (CAZymes) were predicted in the DBCAN3 web server (). Effector proteins were determined using EffectorP-fungi 3.0 (), while virulence factors were predicted using the pathogen–host interactions database PHI-base () with NCBI BLAST+ blastp () on the Galaxy platform. Phylogenetic analysis of the filtered proteins was conducted in MEGA11 () using the UPGMA method.
2.5 Genome representation
Marker positions of the filtered protein dataset were obtained from the GFF file of the Fv 7600 strain (GCA_000149555.1) and then imported into R version 4.3.3. The qtl and LinkageMapView packages () were loaded for linkage map analysis and visualization. The maximum position (maxpos) across chromosomes 1 to 11 was calculated, and an axis-tick vector (at.axis) from 0 to maxpos in 1 cM increments was generated, along with a label vector (axlab) at 20 cM intervals. Marker density was estimated using lmvdencolor with an RColorBrewer color palette. Finally, all chromosomes were plotted with marker labels (dupnbr=TRUE) and a density gradient (denmap=TRUE) using lmv.linkage.plot with “Chromosome %” as the unit.
2.6 GO and KEGG analysis
The online DAVID software () was used to obtain the gene ontology (Biological processes) of the filtered proteins. GO terms and biological processes were summarized and consolidated to build a specific functional database. This database was used with the GGPLOT2 package to build a GO bar plot indicating the most abundant GO terms and its functions. GGPLOT2 in Rstudio version 4.3.3. was used to build a bar plot showing the number of proteins identified in each pathway.
2.7 String network analysis
Based on their reported function, a subset of the selected amino acid sequences with assigned gene names was mapped in the STRINGdb (https://string-db.org/) using the rba_string_map_ids function from rbioapi package, specifying the specie ID: 334819 (Fv) in Rstudio version 4.3.3. A selection of the STRING IDs corresponding to “high identity sequences to the PHI-base” was used to retrieve the interaction network for each protein using the rba_string_interaction_partners function. The top ten protein interactions were selected based on their scores for each chosen protein, and a subset of proteins was visualized using the rba_string_network_image function.
2.8 Phylogenetic analysis of Fusarium genomes
Fusarium verticillioides genomes (GCA013759275.1, GCA020882315.1, GCA003316975.1, GCA017309895.1, GCA_033110985.1, and GCA000149555.1) were retrieved from GenBank and analyzed alongside the Fv DA42 genome using whole-genome phylogeny with the ParSNP tool (), with the GCA000149555.1 genome as reference. The resulting Newick file was then visualized using the iTOLL web server ().
2.9 Model of bipartite interaction
A basic version of the figure was initially generated using an AI-assisted tool (ChatGPT, OpenAI), and it was adapted and refined to accurately represent the data and findings.
3 Results
3.1 Fusarium verticillioides degrades epidermis and root hairs
To determine whether the mechanism of root infection involves secretome-mediated degradation of the epidermis and root hair cell walls, seeds were infected with Fv DA42 strain. The primary roots, including apparently healthy tissue (without visible necrosis) and necrotic tissue from 7-day-old seedlings, were observed by SEM (Figure 1). In the root zone without visible necrosis, F. verticillioides hyphae grew over the surface of the epidermis and among the root hairs. Root hairs were abundant, and many of them in the basal zone was collapsed and sunken, suggesting a loss of cell wall firmness (Figure 1A). In contrast, in the necrotic root tissues, there were no roots hairs, and the epidermal surface exhibited extensive degradation, even showing holes (Figure 1B). Fv hyphae grew profusely over the damaged tissue. The extensive degradation of the root epidermis and root hairs, along with the loss of cell wall integrity, is consistent with the action of cell wall-degrading enzymes present in the fungal secretome, supporting the hypothesis that secretome-mediated degradation is a key mechanism in the Fv root infection process.
Figure 1
3.2 Fusarium verticillioides genome assembly
Previous evidence indicates that Fv secretes enzymes that play a crucial role during the root infection process. Therefore, it is essential to characterize the complete repertoire of Fv secreted proteins (the secretome). The Fv genome was assembled de novo from paired-end Illumina reads. The final assembly comprises 175 contigs, with the largest contig measuring 2.911 Mb, yielding a total assembly size of 42.27 Mb, an N50 of 1.244 Mb, and an N90 of 0.431 Mb. Predicted gene lengths ranged from 201 bp to 22,734 bp, with a mean gene length of 1,492 bp, and an overall GC content of 48.51%. Mapping to the reference genome was performed. Raw paired-end reads were aligned to the F. verticillioides reference gap-less genome (GCA_027571605.1) using Bowtie2 (v2.5.3). Mapping quality was evaluated with Qualimap BamQC. Of 37,630,594 total reads, 36,178,219 (96.14%) mapped to the reference, leaving 452,375 reads (3.86%) unmapped. Properly paired reads accounted for 36,178,219 (96.14%), while singleton mappings represented 307,027 reads (0.82%). No duplicates were flagged in the BAM file, but an estimated 13,628,671 reads were marked as duplicates, yielding a duplication rate of 63.22%. The mean read length was 117 bp, and 1.8% of read pairs overlapped. The mean mapping quality score was 40.21. Coverage analysis showed an average depth of 100.93× (SD = 55.68), with per‐chromosome mean coverages ranging from 95.68× (chr11: CM000588) to 102.92× (chr3: CM000580), and a GC content of 48.83%. Insert sizes were broadly distributed (median = 486 bp; 1st–3rd quartiles = 419–567 bp), and the overall error rate (mismatches plus indels) was 0.53%, comprising 16,648,075 mismatches, 770,398 insertions (2.02% of reads), and 352,021 deletions (0.8% of reads). This assembly quality provides a robust foundation for downstream annotation and functional analyses.
3.3 Secretome prediction
To identify secreted proteins potentially involved in root cell wall degradation and pathogenesis an in silico prediction of the secretome was performed (Figure 2). A total of 14,198 protein-coding genes were predicted from the F. verticillioides genome. The secretome genes and proteins functionally annotated in the eggNOG-mapper database identified 3,181 extracellular proteins; 1,279 proteins with signal peptide; 1,108 without transmembrane domains, and 111 proteins with glycosylphosphatidylinositol (GPI) anchors.
Figure 2
Out of the 14,198 proteins, 997 were predicted as classical secreted proteins, representing 7.03% of the protein-coding genes (Supplementary Table 1). The secretome of Fv revealed 16 lipases (1.6%), 62 proteases (6.7%), and 288 unknow or uncharacterized proteins (28.88%). The analyzed subset includes 262 CAZymes, 400 effectors, and 481 virulence factors. The use of diverse databases for functional prediction resulted in redundant function or annotation. In general, 26% of the secreted proteins corresponded to unclassified metabolism, with three major groups related to carbohydrate and glycan metabolism (31%), 11% to amino acid metabolism, and 4% related to lipid metabolism, as identified by BlastKoala (). Smaller proportions were associated with genetic information processing (4%), protein families related to metabolism (6%), and other minor categories (<1%). This distribution highlights the abundance of carbohydrate and amino acid metabolism related-proteins in the secretome, suggesting their relevance in host-pathogen interactions and cell wall degradation mechanism mediated by the secretome.
3.4 Genomic distribution and functional annotation of the Fusarium verticillioides secretome
The chromosomal distribution of secretome genes varied widely across chromosomes (Figure 3), with the highest number found on chromosomes 4 (137 genes) and 8 (123 genes). In contrast, the lowest number of protein-coding genes were observed on chromosomes 7 (62 genes) and 9 (63 genes). However, the secretome gene density was higher on chromosome 10, with 49.2 genes/Mb, while chromosome 1 had the lowest density with 14 genes/Mb. Eleven genes, including FVEG_14091 (ricin B lectin), FVEG_13989 (gluconolactonase), and FVEG_14136 (alpha-glucosidase), were not located on any chromosome, but were instead found in the scaffolds NW_017387866.1, NW_017387870.1, and NW_017387871.1 of the reference genome. Gene ontology (GO) terms for the secretome proteins were identified by functional category (Figure 3B). The molecular function annotation term (green bars) for the secretome revealed the presence of enzymatic activities associated with plant cell wall degradation. The analysis of the most significantly enriched terms identified proteins involved in rhamnogalacturonan endolyase activity (GO:0102210), mannan endo-1,4-β-mannosidase activity (GO:0016985), and glucan exo-1,3-β-glucosidase activity (GO:0004338). Additional activities included galactose oxidase (GO:0045480), endo-1,3(4)-β-glucanase (GO:0052861), and carbon-oxygen lyase acting on polysaccharides (GO:0016837). Furthermore, proteolytic activities such as tripeptidyl-peptidase (GO:0008240) and metallocarboxypeptidase (GO:0004181) were also detected, indicating potential roles in protein degradation during host colonization. The presence of endo-1,4-β-xylanase (GO:0031176) and α-L-arabinofuranosidase (GO:0046556), along with other relevant and abundant activities for root cell wall degradation were identified, such as pectin lyase activity (GO:0047490) and feruloyl esterase activity (GO:0030600).
Figure 3
The cellular component prominent GO terms (blue bars) include apoplastic space (GO:0042597), extracellular region (GO:0005615 and GO:0005576), and cell‐wall, which are consistent with the secretory nature of these proteins. The biological processes GO terms (red bars) include catabolic processes to degrade the principal component of cell wall, including polysaccharides (simple and complex) (GO:0000272, GO:0046355, GO:0045493, and GO:0046373), proteins (GO:0009063 and GO:0006508), and lipids (GO:0016042). In general, the extensive diversity and genomic distribution of these genes, along with their functional profiles, suggest that many secreted proteins participate in host interactions and extracellular enzymatic activities characteristic of virulence factors and effectors in pathogens.
3.5 Phylogenetic clustering of CAZymes and proteases in the Fusarium verticillioides secretome
A phylogenetic analysis was performed based on protein evolutionary similarity within the two more abundant functional groups of the Fv DA42 secretome: CAZymes (carbohydrate-active enzymes) and proteases (Figure 4). The phylogenetic relationships among various CAZymes involved in the degradation, modification, or biosynthesis of carbohydrates are evident (Figure 4A). Branches group proteins with similar sequences, indicating families that are functionally or evolutionarily related. Several well-defined clades (e.g., highlighted in red, orange, and purple) suggest that certain CAZymes share conserved domains or common activities, such as glucosidases, cellulases, chitinases, and glucosyl hydrolases. Additionally, a phylogenetic analysis of secreted proteases, enzymes potentially involved in degrading host proteins and contributing to virulence, reveals clustering by sequence similarity (Figure 4B). Colored clusters reflect specific protease families, such as serine proteases or metalloproteases, and the presence of multiple distinct clades underscores the functional diversification among the secreted proteases. Overall, the phylogenetic analysis highlights the structural and functional differences and similarities between the main secreted protein groups in Fv DA42 strain. Highly conserved proteins within each group suggest key roles in host interaction, either by degrading plant cell wall components or by interfering with host defenses through proteolytic activity.
Figure 4
3.6 Protein-protein interaction of key enzymes and effectors identified through STRING
Based on a global STRING analysis of 14 selected genes related to cell wall degradation and infection process and their high identity in the PHI-base database. (Supplementary Table 2), several genes (FVEG_03246, FVEG_10795, and FVEG_05849) with high functional connectivity (scores > 0.7) were identified, exhibiting very strong interactions related to cell cycle, plant cell wall degradation, and putative adhesion functions, respectively. We also detected polysaccharide degradation networks involving multiple genes, such as FVEG_10795 (putative pectinesterase), FVEG_09361 (putative ferulic acid esterase), FVEG_05642 (putative chitin binding protein), and FVEG_09702 (putative pectate lyase), suggesting their coordinated role in modifying the plant cell wall, potentially during fungal infection. In parallel, FVEG_09149 (Mg peroxidase) and FVEG_04647 (putative necrosis inducing protein) appear to be involved in cell wall degradation through the depolymerization of lignin and necrosis processes. Based on their high STRING score analysis and relevant literature, a subset of genes was selected for visualization (Figure 5), each linked to a specific enzyme or effector: FVEG_09361 (feruloyl esterase, EC 3.1.1.73), FVEG_10795 (pectinesterase, EC 3.1.1.11), FVEG_13183 (cell wall glycosyl hydrolase), FVEG_05642 (chitin-binding type-4 domain protein), FVEG_09149 (peroxidase, EC 1.11.1), and FVEG_04647 (necrosis-inducing protein). This analysis confirms the roles for these proteins in plant cell wall degradation (pectinesterases, polygalacturonases, and esterases), as virulence effectors (necrosis inducers or chitin-binding camouflage proteins), and in oxidative responses (peroxidases). Specifically, FVEG_10795 (pectinesterase) catalyzes the demethylation of pectin, a prerequisite for its breakdown. It interacts most strongly with FVEG_08451 (endo-polygalacturonase; score 0.933). Together, these enzymes likely act sequentially; FVEG_10795 demethylates pectin, and FVEG_08451 cleaves it into monomers. FVEG_09361 (feruloyl esterase) contributes to the breakdown of ferulic acid cross-links in the plant cell wall. Its main partner is FVEG_16566 (score 0.372), an AB hydrolase-1 domain protein that may participate in related hydrolytic pathways. FVEG_13183 (cell wall glycosyl hydrolase) probably targets glucans or other complex polysaccharides. Its principal interactor is FVEG_12180 (score 0.366), an SGNH-type esterase involved in polysaccharide ester modification, suggesting a shared depolymerization pathway. FVEG_05642 (chitin-binding type-4 domain protein) acts as a fungal effector, masking chitin from plant defenses. Its main interactor is FVEG_12434 (pectate lyase; score 0.224); although the score is low, this suggests potential collaborative mechanisms for attacking chitin and pectate. FVEG_09149 (peroxidase) is involved in oxidative processes, likely related to defense or stress adaptations, and interacts primarily with FVEG_04068 (ribonuclease H2 subunit B; score 0.723), indicating a possible co-expression or joint participation in damage-response or oxidative-stress pathways. FVEG_04647 (necrosis-inducing protein) is a putative fungal effector that triggers host cell death. Its main interaction is with FVEG_04937 (cytochrome c oxidase subunit; score 0.490), implying a link to mitochondrial processes that may be tied to necrosis or stress signaling in the host.
Figure 5
3.7 Phylogenetic relation of draft genome analysis of Fusarium verticillioides DA42 strain
The phylogenetic analysis of the genomes included strains isolated from the United States, Australia, India, and Mexico. The phylogenetic tree (Figure 6) showed that the Fv DA42 strain genome from Mexico clustered with the Fv 7600 reference genome (GCA000149555.1) isolate in USA and an Indian strain (GCA033110985.1), however, Fv DA42 strain genome was more closely related to the reference strain than to the Indian strain. In a different tree branch another Fv isolate from USA (GCA013759275.1) groups with the isolate from Italy (GCA020882315.1) and both are more distant to the isolate from Australia (GCA003316975.1). The phylogenetic analysis supports a close evolutionary relationship between the DA42 strain and the reference genome 7600 strain, while revealing partial genetic distances among isolates from different geographic regions of origin.
Figure 6
4 Discussion
The extended agricultural life cycle of maize (6–8 months) exacerbates the challenges associated with Fusarium infection. From the vulnerable seedling stage to the critical phases of vegetative and reproductive growth, the potential for Fusarium to disrupt and delay normal development of maize remains a persistent concern for both farmers and researchers (). The fungus’s ability to penetrate various plant tissues, including roots, stalks, and seeds, underscores the need for a comprehensive understanding of the mechanisms that facilitate its entry and subsequent colonization within the host (; ; ; ). and describe the early manifestations of Fusarium infection, which include yellow-brown discoloration and the subsequent formation of necrotic spots on the roots. These symptoms not only serve as visible indicators of infection but also suggests a complex interaction between the pathogen and the host plant at the cellular level. The collapse of root hairs further emphasizes the severity of this interaction, suggesting the active degradation of cell walls in the epidermal cells during a necrotrophic phase (Figure 1). Although, the SEM observations were made at 7 days post infection and thus correspond to an advanced stage of root necrosis, it is highly likely that the degradative enzymatic machinery and signaling are strongly activated several days earlier (). In the current study, the de novo high-quality genome assembly of Fv DA42 strain allowed for downstream annotation and functional analyses. The DA42 strain showed a closest phylogenetic similarity to the genome reference Fv 7600 (GCA_000149555.1) reported by and the improved gapless genome in , than to the BIONCL4 (GCA_033110985.1) strain (). This indicates a close evolutionary relationship (Figure 6) and potentially similar phenotypic behavior. Similar to the findings of , the genome assembly from Italy (GCA020882315.1) showed greater distance from the DA42 strain. The separation of strains from Mexico (FvDA42) and U.S. (Fv USA), and as well as the convergence of Mexican (FvDA42), U.S. (Fv7600 USA) and Indian (Fv India) strains into different clades or shared clades respectively, may reflect historical gene flow associated with the global movement and exchange of commercial maize germplasm, local adaptation and long-term cultivation of genetically distinct maize varieties in each region could have driven the divergence of associated Fv populations. Although our coverage analysis was over 100.9X and the mapping to the first version reference genome was 96.14%, these differences could be related either to sequence technology or intra-species variability. Therefore, the genome of Fv DA42 strain, which clusters closest to the reference Fv 7600 strain, may offer valuable insights into infection mechanisms within the Fusarium genus, and represents a promising candidate for comparative future pathogenicity assays.
Despite the importance of the Fv secretome, a few studies have been conducted to identify secreted proteins. The results of this study reinforce the key role of the secretome for Fv fungal infection and colonization processes in maize roots. According to our characterization, approximately 7% of the Fv DA42 proteome corresponds to secreted proteins via classical pathways, primarily consisting of CAZymes, proteases, effectors, and virulence factors, which play important roles in cell wall degradation and overcoming maize defenses. In accordance with the relevance of the integrity of the host cell wall, recent finding on the maize glycosyltransferase (ZmXYXT2) showed that reinforcement of the cell walls restricted F. verticillioides infection, whereas thinner cell walls in the zmxyxt2 mutant facilitated the colonization and fumonisin accumulation (). Studies on plant pathogens emphasize the significance of secretomes in plant–fungus interactions, the proportion of these secreted proteins ranges from 4% to 14% for necrotrophs compared to biotrophs (). reported the first proteome of Fv 7600, identifying a total of 166 secreted proteins, which represented only 1.17% of the predicted proteome, a relatively low proportion of the total protein-coding genes. This limited detection may be attributed to the proteomic technique used at the time, which likely had lower sensitivity and coverage compared to other advanced or complementary techniques available today. A few years later, performed whole-genome sequencing of F. verticillioides BIONCL4 strain isolated from maize grains using Illumina technology. They predicted a total of 15,053 protein-coding-genes, of which 2,058 were secreted proteins and 676 were associated with the classical secretory system. These represent 13.6% and 4.49% of the total of secreted proteins and classical secretory pathways, respectively. The predicted secretome of F. graminearum, is approximately 4.11% (). In the in silico secretome of the palm dieback-causing agent F. oxysporum f. sp. albedinis, 1,464 out of 16,887 genes (8.6%) were predicted as secreted proteins (). A recent study in F. graminearum provides the first proteome secreted into the apoplast of wheat coleoptiles, identifying the metalloprotease effector Fg28, which induces cell death and reactive oxygen species at the first day of infection (). Similarly, in F. oxysporum, a novel secreted cysteine-rich protein, FolSCP1, was identified as a virulence factor that promotes infection by binding and inhibiting PR-5, a positive regulator of tomato immunity (). In F. verticillioides, a RNA-seq analysis revealed that FvLcp1, a secreted LysM protein required for fumonisin production, which accumulates in appressoria and contributes to host recognition, suppression of host cell death, and promotion of fumonisin biosynthesis during maize kernel colonization (). These studies highlight the role of secreted proteins in Fusarium host interactions. In particular, these emphasize the biological relevance of the Fv secretome and the need for further studies to better understand its role in the molecular mechanisms of maize–Fusarium interactions.
The GO analysis of the DA42 strain secretome (Figure 3B) revealed enrichment in enzymatic activities such as endogalacturonases, exoglucosidases, and feruloyl esterases, which are essential for breaking down pectin, hemicellulose, and lignin. These biological processes were classified under the carbohydrate catabolism. Also, processes related to proteolysis, cell adhesion and cell wall organization were identified. These results are consistent with secretome analyses of necrotrophic fungi, such as Neofusicoccum parvum () and Botrytis cinerea (), where CAZymes, including pectinesterases, endopolygalacturonases, and proteases, are involved in host protein degradation.
Similarly, several effectors were identified in Fv DA42, including peroxidases and necrosis-inducing proteins, which are comparable to those known to induce host necrosis or suppress host immunity in other fungal pathogens (; ; ). A meta-analysis conducted by on 150 phytopathogenic fungi revealed that secretomes vary with fungal lifestyle, with necrotrophs tending to secrete a wide range of CAZymes targeting pectin and xylan. One study showed that the predicted secretome of F. graminearum, which includes 574 proteins, contains phytotoxic enzymes and effectors, many of which are transcriptionally upregulated during infection (). The interaction between CAZymes and proteases (Figure 4) and the connected nodes in the protein–protein interaction network (Figure 5) for FVEG_10795, FVEG_09361, and FVEG_09149 reinforce the hypothesis of a coordinated mechanism for fungal attack. This could involve sequential pectin demethylation, enzymatic breakdown, and oxidative stress, forming a tripartite virulence strategy. The collapse of root hairs and degradation of the cell wall epidermis observed by SEM microscopy (Figure 1) confirms that enzymatic attack is effective during the early stages of infection, as previously hypothesized by .
This research integrates scan electronic microscopy, genome sequencing, in silico secretome analysis, and protein–protein interaction analysis to decipher the infection strategy of F. verticillioides in maize roots. Previous studies have addressed the general pathogenicity of this species, primarily focusing on ear and stalk infections. Our study provides the first integrated correlation between in silico secretome composition and SEM observations of early-stage root tissue degradation, including root hairs collapse and breakage of epidermal cell walls. To our knowledge, this is the first report that combines multi-layered genomic and proteomic predictions with structural (SEM) evidence of infection in maize roots by Fv. This integrative approach opens new aspects to understand the molecular arsenal underlying fungal root pathogenesis and to identify specific secreted proteins as promising targets for resistance breeding and antifungal strategies, thereby contributing with valuable tools for sustainable crop protection. In summary, the secretome of F. verticillioides plays a central role in colonizing susceptible maize roots, making evident the enrichment of enzymes related to cell wall degradation (carbohydrate polymers, glucans and proteins), effectors, and virulence. Although Fv has been described as a hemibiotroph, our OMIC data and SEM observations (in this and our previous report) suggest a predominantly necrotrophic lifestyle when interacting with susceptible genotypes during the early infection stages, likely driven by the early and sustained activation of cell wall–degrading enzymes and necrosis-inducing factors through still unknown mechanisms (Figure 7).
Figure 7
5 Conclusion
This study provides a multidimensional analysis of Fv DA42, combining SEM, genome assembly, in silico secretome prediction, interaction network analyses, and phylogenetic. We identified that approximately 7% of the strain’s proteome consists of secreted proteins, which are rich in CAZymes, proteases, peroxidases, effectors and necrosis-inducing proteins. Through SEM, we detected early-stage damage, including the collapse of root hairs and degradation of root epidermal cell walls, directly linking molecular predictions with structural evidence. Significant interaction, such as the sequential action of FVEG_10795 (pectinesterase) and FVEG_08451 (polygalacturonase), support a model of coordinated enzymatic machinery that deconstructs pectin in a stepwise manner during the initial colonization. To our knowledge, this is the first integrated report that combines whole-genome secretome annotation, protein–protein interaction modeling, and root tissue-level SEM observations to elucidate the molecular mechanisms of F. verticillioides root infection. The specific secreted proteins identified emerge as promising targets for future genetic resistance breeding or/and antifungal control strategies. Our work thus lays a solid foundation for sustainable prevention methods against maize root rot and enhances the molecular understanding of Fusarium–maize interactions.
Statements
Data availability statement
The datasets presented in this study can be found online the BioProject PRJNA1294327.
Author contributions
JN-Z: Conceptualization, Data curation, Formal analysis, Methodology, Software, Supervision, Writing – original draft, Writing – review & editing. YA-C: Data curation, Formal analysis, Methodology, Writing – original draft, Writing – review & editing, Investigation. JC-M: Methodology, Writing – review & editing. ME-M: Methodology, Writing – review & editing, Data curation. IM-M: Methodology, Writing – review & editing, Investigation. RS: Investigation, Methodology, Writing – review & editing, Data curation. FQ-F: Data curation, Investigation, Methodology, Writing – review & editing, Conceptualization, Formal analysis, Funding acquisition, Project administration, Software, Supervision, Validation, Writing – original draft.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This research was supported by the Mexican National Council on Science and Technology (CONACyT) or Secretaría de Ciencia, Humanidades, Tecnología e Innovación (SECIHTI), with financial support from Infraestructura no. 250738, Fronteras de la Ciencia no. 1070, and CBF-2025-I-1850 and the Instituto Politécnico Nacional-SIP (grant numbers 20230074, 20240601, 20250878, and 20250006).
Acknowledgments
We thank Sabina Quiroz-García for adapting and improving the figure created with ChatGPT (OpenAI).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that Generative AI was used in the creation of this manuscript. A basic version of the figure was initially generated using an AI-assisted tool (ChatGPT, OpenAI), and it was adapted and refined to accurately represent the data and findings.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1672761/full#supplementary-material
References
1
AlbertsB.JohnsonA.LewisJ. (2002). Molecular biology of the cell (New York: Garland Science).
2
AndrewsS. (2023). FastQC A Quality Control tool for High Throughput Sequence Data. Available online at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (Accessed June 13, 2024]).
3
Atanasova-PenichonV.BernillonS.MarchegayG.LornacA.Pinson-GadaisL.PontsN.et al. (2014). Bioguided isolation, characterization, and biotransformation by fusarium verticillioides of maize kernel compounds that inhibit fumonisin production. Mol. Plant-Microbe Interact.27, 1148–1158. doi: 10.1094/Mpmi-04-14-0100-R
4
BaldwinT. T.ZitomerN. C.MitchellT. R.ZimeriA. M.BaconC. W.RileyR. T.et al. (2014). Maize Seedling Blight Induced by Fusarium verticillioides: Accumulation of Fumonisin B-1 in Leaves without Colonization of the Leaves. J. Agric. Food Chem.62, 2118–2125. doi: 10.1021/jf5001106
5
BeccariG.HaoG. X.LiuH. Q. (2022). Pathogenesis: Infection mechanisms and disease progression in host plants. Front. Plant Sci.13. doi: 10.3389/fpls.2022.1020404
6
BishopC. D.CooperR. M. (1983). An ultrastructural study of root invasion in three vascular wilt diseases. Physiol. Plant Pathol22, 15–IN13. doi: 10.1016/S0048-4059(83)81034-0
7
BoenischM. J.SchaferW. (2011). Fusarium graminearum forms mycotoxin producing infection structures on wheat. BMC Plant Biol.11, 110. doi: 10.1186/1471-2229-11-110
8
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, 2114–2120. doi: 10.1093/bioinformatics/btu170
9
BouqellahN. A.ElkadyN. A.FaragP. F. (2023). Secretome analysis for a new strain of the blackleg fungus plenodomus lingam reveals candidate proteins for effectors and virulence factors. J. Fungi (Basel)9, 740. doi: 10.3390/jof9070740
10
BradshawM. J.BartholomewH. P.FonsecaJ. M.GaskinsV. L.PruskyD.JurickW. M. (2021). Delivering the goods: Fungal secretion modulates virulence during host-pathogen interactions. Fungal Biol. Rev.36, 76–86. doi: 10.1016/j.fbr.2021.03.007
11
BrownN. A.AntoniwJ.Hammond-KosackK. E. (2012). The predicted secretome of the plant pathogenic fungus fusarium graminearum: A refined comparative analysis. PloS One7. doi: 10.1371/journal.pone.0033731
12
BrylaM.PierzgalskiA.ZapasnikA.UwinezaP. A.Ksieniewicz-WozniakE.ModrzewskaM.et al. (2022). Recent research on fusarium mycotoxins in maize-A review. Foods11, 3465. doi: 10.3390/foods11213465
13
ChethanaK. W. T.JayawardenaR. S.ChenY. J.KontaS.TibprommaS.AbeywickramaP. D.et al. (2021). Diversity and function of appressoria. Pathogens10, 746. doi: 10.3390/pathogens10060746
14
CheungN.TianL.LiuX.LiX. (2020). The destructive fungal pathogen botrytis cinerea-insights from genes studied with mutant analysis. Pathogens9, 923. doi: 10.3390/pathogens9110923
15
CockP. J. A.ChiltonJ. M.GrüningB.JohnsonJ. E.SoranzoN. (2015). NCBI BLAST+ integrated into Galaxy. GigaScience4, s13742–015–0080–7. doi: 10.1186/s13742-015-0080-7
16
CzarneckaD.CzubackaA.Agacka-MoldochM.Trojak-GoluchA.KsiezakJ. (2022). The occurrence of fungal diseases in maize in organic farming versus an integrated management system. Agronomy-Basel12, 558. doi: 10.3390/agronomy12030558
17
CzemborE.StepienL.WaskiewiczA. (2015). Effect of environmental factors on species and associated mycotoxins in maize grain grown in Poland. PloS One10, e0133644. doi: 10.1371/journal.pone.0133644
18
de SousaR. R.OsorioP. R. A.NoseN. P. E.de ArrudaG. L.FerreiraT. P. D.HaesbaertF. M.et al. (2022). Detection and transmission of Fusarium verticillioides in corn seeds according to the plant stage. Acta Scientiarum-Agronomy44, e53213. doi: 10.4025/actasciagron.v44i1.53213
19
DinhS. Q.JoyceD. C.IrvingD. E.WearingA. H. (2011). Histology of waxflower (Chamelaucium spp.) flower infection by Botrytis cinerea. Plant Pathol.60, 278–287. doi: 10.1111/j.1365-3059.2010.02347.x
20
DuncanK. E.HowardR. J. (2010). Biology of maize kernel infection by fusarium verticillioides. Mol. Plant-Microbe Interact.23, 6–16. doi: 10.1094/Mpmi-23-1-0006
21
GaiX.DongH.WangS.LiuB.ZhangZ.LiX.et al. (2018). Infection cycle of maize stalk rot and ear rot caused by Fusarium verticillioides. PloS One13, e0201588. doi: 10.1371/journal.pone.0201588
22
GalaxyC. (2024). The Galaxy platform for accessible, reproducible, and collaborative data analyses: 2024 update. Nucleic Acids Res.52, W83–W94. doi: 10.1093/nar/gkae410
23
GeorgiadouS. P.VelegrakiA.ArabatzisM.NeonakisI.ChatzipanagiotouS.DalekosG. N.et al. (2014). Cluster of Fusarium verticillioides bloodstream infections among immunocompetent patients in an internal medicine department after reconstruction works in Larissa, Central Greece. J. Hosp. Infection86, 267–271. doi: 10.1016/j.jhin.2014.01.011
24
GörtzA.OerkeE. C.SteinerU.WaalwijkC.de VriesI.DehneH. W. (2008). Biodiversity of Fusarium species causing ear rot of maize in Germany. Cereal Res. Commu.36, 617–622. doi: 10.1556/CRC.36.2008.Suppl.B.51
25
GuptaN. C.YadavS.AroraS.MishraD. C.BudhlakotiN.GaikwadK.et al. (2022). Draft genome sequencing and secretome profiling of Sclerotinia sclerotiorum revealed effector repertoire diversity and allied broad-host range necrotrophy. Sci. Rep.12, 21855. doi: 10.1038/s41598-022-22028-z
26
GurevichA.SavelievV.VyahhiN.TeslerG. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics29, 1072–1075. doi: 10.1093/bioinformatics/btt086
27
HajnalE. J.KosJ.RadicB.AnicM.RadovicR.KudumijaN.et al. (2023). Impact of climate changes on the natural prevalence of mycotoxins in maize harvested in Serbia and Croatia. Foods12, 1002. doi: 10.3390/foods12051002
28
HortonP.ParkK. J.ObayashiT.FujitaN.HaradaH.Adams-CollierC. J.et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res.35, W585–W587. doi: 10.1093/nar/gkm259
29
Huerta-CepasJ.SzklarczykD.HellerD.Hernandez-PlazaA.ForslundS. K.CookH.et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res.47, D309–D314. doi: 10.1093/nar/gky1085
30
HusseinH. M.ChristensenM. J.BaxterM. (2003). Occurrence and distribution of Fusarium species in maize fields in New Zealand. Mycopathologia156, 25–30. doi: 10.1023/a:1021307023039
31
JiaM.GongX.FanM.LiuH.ZhouH.GuS.et al. (2023). Identification and analysis of the secretome of plant pathogenic fungi reveals lifestyle adaptation. Front. Microbiol.14. doi: 10.3389/fmicb.2023.1171618
32
JonesJ. D. G.DanglJ. L. (2006). The plant immune system. Nat.444, 323–329. doi: 10.1038/nature05286
33
KallL.KroghA.SonnhammerE. L. (2007). Advantages of combined transmembrane topology and signal peptide prediction–the Phobius web server. Nucleic Acids Res.35, W429–W432. doi: 10.1093/nar/gkm256
34
KanehisaM.SatoY.MorishimaK. (2016). BlastKOALA and ghostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol.428, 726–731. doi: 10.1016/j.jmb.2015.11.006
35
KangZ.BuchenauerH. (2002). Studies on the infection process of Fusarium culmorum in wheat spikes: Degradation of host cell wall components and localization of trichothecene toxins in infected tissue. Eur. J. Plant Pathol.108, 653–660. doi: 10.1023/A:1020627013154
36
KilleB.NuteM. G.HuangV.KimE.PhillippyA. M.TreangenT. J. (2024). Parsnp 2.0: scalable core-genome alignment for massive microbial datasets. Bioinformatics40, btae311. doi: 10.1093/bioinformatics/btae311
37
LeslieJ. F.SummerellB. A. (2007). The fusarium laboratory manual (USA, UK, Australia: Blackwell Publishing Ltd).
38
LetunicI.BorkP. (2024). Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res.52, W78–W82. doi: 10.1093/nar/gkae268
39
Leyva-MadrigalK. Y.Larralde-CoronaC. P.Apodaca-SanchezM. A.Quiroz-FigueroaF. R.Mexia-BolanosP. A.Portillo-ValenzuelaS.et al. (2015). Fusarium species from the fusarium fujikuroi species complex involved in mixed infections of maize in northern sinaloa, Mexico. J. Phytopathol.163, 486–497. doi: 10.1111/jph.12346
40
LiP.ZhaoR.FangY.FanY.HuQ.HuangW.et al. (2025). Proteomic analysis of the fusarium graminearum secretory proteins in wheat apoplast reveals a cell-death-inducing M43 peptidase. J. Fungi (Basel)11, 240. doi: 10.3390/jof11040240
41
LoweR. G.HowlettB. J. (2012). Indifferent, affectionate, or deceitful: lifestyles and secretomes of fungi. PloS Pathog.8, e1002515. doi: 10.1371/journal.ppat.1002515
42
MaL. J.van der DoesH. C.BorkovichK. A.ColemanJ. J.DaboussiM. J.Di PietroA.et al. (2010). Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature464, 367–373. doi: 10.1038/nature08850
43
McCotterS. W.HorianopoulosL. C.KronstadJ. W. (2016). Regulation of the fungal secretome. Curr. Genet.62, 533–545. doi: 10.1007/s00294-016-0578-2
44
MurilloI.CavallarinL.San SegundoB. (1999). Cytology of infection of maize seedlings and immunolocalization of the pathogenesis-related PRms protein. Phytopathology89, 737–747. doi: 10.1094/Phyto.1999.89.9.737
45
NavaleV. D.SawantA. M.GowdaV. U.VamkudothK. R. (2022). Assembly, annotation, and comparative whole genome sequence of fusarium verticillioides isolated from stored maize grains. Pathogens11, 810. doi: 10.3390/pathogens11070810
46
Nazar PourF.PedrosaB.OliveiraM.FidalgoC.DevreeseB.DriesscheG. V.et al. (2022). Unveiling the secretome of the fungal plant pathogen neofusicoccum parvum induced by in vitro host mimicry. J. Fungi (Basel)8, 971. doi: 10.3390/jof8090971
47
OkelloP. N.PetrovićK.KontzB.MathewF. M. (2019). Eight species of fusarium cause root rot of corn (Zea mays) in south dakota. Plant Health Prog.20, 38–43. doi: 10.1094/php-11-18-0075-rs
48
Olivares-GarciaC. A.Mata-RosasM.Pena-MontesC.Quiroz-FigueroaF.Segura-CabreraA.ShannonL. M.et al. (2020). Phenylpropanoids are connected to cell wall fortification and stress tolerance in avocado somatic embryogenesis. Int. J. Mol. Sci.21, 5679. doi: 10.3390/Ijms21165679
49
OrenL.EzratiS.CohenD.SharonA. (2003). Early events in the Fusarium verticillioides-maize interaction characterized by using a green fluorescent protein-expressing transgenic isolate. Appl. Environ. Microbiol.69, 1695–1701. doi: 10.1128/AEM.69.3.1695-1701.2003
50
OuelletteL. A.ReidR. W.BlanchardS. G.BrouwerC. R. (2018). LinkageMapView-rendering high-resolution linkage and QTL maps. Bioinformatics34, 306–307. doi: 10.1093/bioinformatics/btx576
51
PfordtA.RomeroL. R.SchiwekS.KarlovskyP.von TiedemannA. (2020). Impact of environmental conditions and agronomic practices on the prevalence of species associated with ear- and stalk rot in maize. Pathogens9, 236. doi: 10.3390/pathogens9030236
52
PierleoniA.MartelliP. L.CasadioR. (2008). PredGPI: a GPI-anchor predictor. BMC Bioinf.9, 392. doi: 10.1186/1471-2105-9-392
53
PintoL.Soler-LópezL.SerranoA.Sánchez-RodríguezC. (2025). Between host and invaders: the subcellular cell wall dynamics at the plant–pathogen interface. Annu. Rev. Plant Biol.76, 255–284. doi: 10.1146/annurev-arplant-061824-115733
54
QianH.XiaoZ.ChengL.GengR.MaY.BiY.et al. (2025). A novel secreted protein of fusarium oxysporum promotes infection by inhibiting PR-5 protein in plant. Plant Cell Environ.48, 1021–1036. doi: 10.1111/pce.15200
55
QiuJ. B.XuJ. H.DongF.YinX. C.ShiJ. R. (2015). Isolation and characterization of Fusarium verticillioides from maize in eastern China. Eur. J. Plant Pathol.142, 791–800. doi: 10.1007/s10658-015-0652-5
56
QiuH.ZhaoX.FangW.WuH.AbubakarY. S.LuG.-d.et al. (2019). Spatiotemporal nature of Fusarium graminearum-wheat coleoptile interactions. Phytopathol. Res.1, 26. doi: 10.1186/s42483-019-0033-7
57
Quiroz-FigueroaF. R.Cruz-MendivilA.Ibarra-LacletteE.Garcia-PerezL. M.Gomez-PerazaR. L.Hanako-RosasG.et al. (2023). Cell wall-related genes and lignin accumulation contribute to the root resistance in different maize (Zea mays L.) genotypes to Fusarium verticillioides (Sacc.) Nirenberg infection. Front. Plant Sci.14, 1195794. doi: 10.3389/fpls.2023.1195794
58
RafiqiM.JelonekL.DioufA. M.MbayeA.RepM.DiarraA. (2022). Profile of the in silico secretome of the palm dieback pathogen, Fusarium oxysporum f. sp. albedinis, a fungus that puts natural oases at risk. PloS One17, e0260830. doi: 10.1371/journal.pone.0260830
59
RavalasonH.GriselS.ChevretD.FavelA.BerrinJ. G.SigoillotJ. C.et al. (2012). Fusarium verticillioides secretome as a source of auxiliary enzymes to enhance saccharification of wheat straw. Bioresource Technol.114, 589–596. doi: 10.1016/j.biortech.2012.03.009
60
RománS. G. (2017). “Caracterización de genotipos de maíz (Zea mays L.) a la infección de Fusarium verticillioides en diferentes fases del ciclo de vida de la planta y su correlación con marcadores moleculares de tipo SNPs,” in Matria en Recursos Naturales y Medio Ambiente (Instituto Politécnico Nacional, Guasave, Mexico).
61
RomanS. G.Quiroz-ChavezJ.VillalobosM.Urias-GutierrezV.Nava-PerezE.Ruiz-MayE.et al. (2020). A global screening assay to select for maize phenotypes with a high tolerance or resistance to fusarium verticillioides (Sacc.) nirenberg rots. Agronomy-Basel10, 1990. doi: 10.3390/agronomy10121990
62
SavojardoC.MartelliP. L.FariselliP.CasadioR. (2018). DeepSig: deep Learn. improves Signal Pept. detection proteins. Bioinf.34, 1690–1696. doi: 10.1093/bioinformatics/btx818
63
SeemannT. (2017). Shovill: Faster SPAdes assembly of Illumina reads. Available online at: https://github.com/tseemann/shovill (Accessed Jun 18, 2024).
64
ShermanB. T.HaoM.QiuJ.JiaoX.BaselerM. W.LaneH. C.et al. (2022). DAVID: a web server for functional enrichment analysis and functional annotation of gene list update). Nucleic Acids Res.50, W216–W221. doi: 10.1093/nar/gkac194
65
SimoesD.CarbasB.SoaresA.FreitasA.SilvaA. S.BritesC.et al. (2023). Assessment of agricultural practices for controlling and mycotoxins contamination on maize grains: exploratory study in maize farms. Toxins15, 136. doi: 10.3390/toxins15020136
66
SperschneiderJ.DoddsP. N. (2022). EffectorP 3.0: prediction of apoplastic and cytoplasmic effectors in fungi and oomycetes. Mol. Plant Microbe Interact.35, 146–156. doi: 10.1094/MPMI-08-21-0201-R
67
StankeM.WaackS. (2003). Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics19 Suppl 2, ii215–ii225. doi: 10.1093/bioinformatics/btg1080
68
TamuraK.StecherG.KumarS. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Mol. Biol. Evol.38, 3022–3027. doi: 10.1093/molbev/msab120
69
ThompsonM. E. H.RaizadaM. N. (2018). Fungal pathogens of maize gaining free passage along the silk road. Pathogens7, 81. doi: 10.3390/pathogens7040081
70
UniProtC. (2025). UniProt: the universal protein knowledgebase in 2025. Nucleic Acids Res.53, D609–D617. doi: 10.1093/nar/gkae1010
71
UrbanM.CuzickA.SeagerJ.WoodV.RutherfordK.VenkateshS. Y.et al. (2022). PHI-base in 2022: a multi-species phenotype database for Pathogen-Host Interactions. Nucleic Acids Res.50, D837–D847. doi: 10.1093/nar/gkab1037
72
Velarde FélixS.Valdez RubioN.Zamora GalvánF.López MolinaR.Melgoza VillagómezC. M.Garzón TiznadoJ. A. (2018). Molecular identification of Fusarium spp. isolated maize in Sinaloa. Rev. mexicana Cienc. agrícolas9, 15. doi: 10.29312/remexca.v9i8.756
73
WarhamE. J.ButlerL. D.SuttonB. C. (1996). Seed testing of maize and wheat: A laboratory guide (México: Cimmyt).
74
WuL.WangX.-M.XuR.-Q.LiH.-J. (2011). Root infection and systematic colonization of dsRed-labeled fusarium verticillioides in maize. Acta Agronomica Sin.37, 793–802. doi: 10.1016/S1875-2780(11)60023-0
75
XiaL. K.CaoY. Y.WangJ.ZhangJ.HanS. B.LiH. Y.et al. (2022). First report of fusarium culmorum causing maize stalk rot in China. Plant Dis.106, 1521–1521. doi: 10.1094/Pdis-07-21-1442-Pdn
76
XuY.FanY.LiuL.CaoJ.ZhouJ.LiuE.et al. (2025). Enhancing maize resistance to Fusarium verticillioides through modulation of cell wall structure and components by ZmXYXT2. J. Advanced Res. doi: 10.1016/j.jare.2025.02.023
77
YaoG.ChenW.SunJ.WangX.WangH.MengT.et al. (2023). Gapless genome assembly of Fusarium verticillioides, a filamentous fungus threatening plant and human health. Sci. Data10, 229. doi: 10.1038/s41597-023-02145-8
78
ZhangH.KimM. S.HuangJ.YanH.YangT.SongL.et al. (2022). Transcriptome analysis of maize pathogen Fusarium verticillioides revealed FvLcp1, a secreted protein with type-D fungal LysM and chitin-binding domains, that plays important roles in pathogenesis and mycotoxin production. Microbiol. Res.265, 127195. doi: 10.1016/j.micres.2022.127195
79
ZhangY.P.UyemotoJ. K.KirkpatrickB. C. (1998). A small-scale procedure for extracting nucleic acids from woody plants infected with various phytopathogens for PCR assay. J. Virological Methods71, 45–50. doi: 10.1016/S0166-0934(97)00190-0
80
ZhengJ.GeQ.YanY.ZhangX.HuangL.YinY. (2023). dbCAN3: automated carbohydrate-active enzyme and substrate annotation. Nucleic Acids Res.51, W115–W121. doi: 10.1093/nar/gkad328
Summary
Keywords
Fusariosis, genomic, lytic enzymes, omic, NGS, Zea mays, SEM
Citation
Narváez-Zapata JA, Alaníz-Cuevas Y, Chávez-Medina JA, Espinoza-Mancillas MG, Maldonado-Mendoza IE, Singh RK and Quiroz-Figueroa FR (2025) Unveiling the molecular arsenal: genome sequencing and in silico secretome analysis of Fusarium verticillioides provide insights into maize root rot pathogenesis. Front. Plant Sci. 16:1672761. doi: 10.3389/fpls.2025.1672761
Received
24 July 2025
Accepted
15 September 2025
Published
30 September 2025
Volume
16 - 2025
Edited by
Malkhan Singh Gurjar, Indian Agricultural Research Institute (ICAR), India
Reviewed by
Sanjay Kumar Goswami, Indian Institute of Sugarcane Research (ICAR), India
Deeksha Joshi, ICAR-Indian Agricultural Research Institute, India
Updates
Copyright
© 2025 Narváez-Zapata, Alaníz-Cuevas, Chávez-Medina, Espinoza-Mancillas, Maldonado-Mendoza, Singh and Quiroz-Figueroa.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ignacio E. Maldonado-Mendoza, ignacioemaldonado@yahoo.com.mx; Francisco Roberto Quiroz-Figueroa, fquiroz@ipn.mx
†These authors have contributed equally to this work and share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.