Transcriptome changes in Fusarium verticillioides caused by mutation in the transporter-like gene FST1
© Niu et al.; licensee BioMed Central. 2015
Received: 8 December 2014
Accepted: 19 February 2015
Published: 25 April 2015
Fusarium verticillioides causes an important seed disease on maize and produces the fumonisin group of mycotoxins, which are toxic to humans and livestock. A previous study discovered that a gene (FST1) in the pathogen affects fumonisin production and virulence. Although the predicted amino acid sequence of FST1 is similar to hexose transporters, previous experimental evidence failed to prove function.
Three new phenotypes were identified that are associated with the FST1 mutant of F. verticillioides (Δfst1), namely reduction in macroconidia production, increased sensitivity to hydrogen peroxide, and reduced mycelial hydrophobicity. A transcriptome comparison of the wild type and strain Δfst1 grown on autoclaved maize kernels for six days identified 2677 genes that were differentially expressed. Through gene ontology analysis, 961 genes were assigned to one of 12 molecular function categories. Sets of down-regulated genes in strain Δfst1 were identified that could account for each of the mutant phenotypes.
The study provides evidence that disruption of FST1 causes several metabolic and developmental defects in F. verticillioides. FST1 appears to connect the expression of several gene networks, including those involved in secondary metabolism, cell wall structure, conidiogenesis, virulence, and resistance to reactive oxygen species. The results support our hypothesis that FST1 functions within the framework of environmental sensing.
Fusarium verticillioides (telemorph, Gibberella moniliformis), which is present in most maize fields, can be an asymptomatic endophyte or the causal agent of seedling, stalk, ear, and kernel diseases . The pathogen produces fumonisins, a group of structurally related polyketide mycotoxins, during colonization of maize kernels. Ingestion of fumonisin B1 (FB1), the most predominant fumonisin analog, can result in leukoencephalomalacia in horses and pulmonary edema in swine. The mycotoxin also has been implicated in human diseases, including cancer and birth defects . Guidelines for maximum fumonisin levels in human food and animal feeds have been established worldwide . Furthermore, economic losses associated with fumonisin contamination in maize exports by the three major maize-exporting nations (US, China and Argentina) were estimated at $100 million annually  with the US losses alone at nearly $40 million annually .
Recent publications describe the complexity of genes that influence regulation of fumonisin biosynthesis. Pathway-specific activator FUM21 (FVEG_14633), which controls transcription of the cluster of FUM genes , was shown to increase when F. verticillioides was treated with the histone deacetylase inhibitor chostatin A . These results support evidence that histone modification plays an important role in the epigenetic regulation of fumonisin production [7,8]. There are several intriguing reports indicating that environmental conditions (nutrients and pH) also affect the transcription of FUM genes and FB1 production. Expression of the nitrogen utilization gene AREA (FVEG_02033) was found to be responsible for repression of FB1 production by ammonium . Under repression conditions, AREA is hypothesized to bind to GATA sequences in the promoters of the FUM genes. Generally, acidic conditions favor FB1 production . Experimental evidence indicates that PACC (FVEG_05393), which has homology to the alkaline-activator gene PACC in A. nidulans , inhibits FB1 production and FUM1 (FVEG_00316) transcription at pH 8 . Finally, carbon source and availability, especially amylopectin, greatly affect FB1 biosynthesis . Studies on carbon utilization have led to the identification of two genes, HXK1 (FVEG_00957) and FST1 (FVEG_08441). HXK1, a putative hexose kinase was shown to be required for fructose metabolism . Strains without a functional HXK1 also produced less FB1 and were less virulent on maize kernel than the wild type (WT). The function of FST1 is the focus of the current study.
FST1 was identified through a comparative analysis of genes expressed in colonized maize germ and endosperm tissues . Of 50 putative sugar transporter genes represented on a microarray, FST1 was one of six genes identified as highly expressed during fungal growth in endosperm tissue compared to germ tissue . Expression of FST1 was also reduced in a F. verticillioides strain with a disrupted ZFR1 gene, a putative Zn2Cys6 transcription factor . FST1 encodes a 574-amino-acid protein with 12 putative transmembrane domains. Heterologous expression of FST1 in yeast system failed to show hexose transporter activity . Disruption of FST1 in F. verticillioides resulted in reduced virulence and FB1 production [15,16]. The reduced virulence phenotype in inoculated kernels was manifested as slower growth and rot symptoms when compared to the WT . When inoculated onto autoclaved kernels or synthetic media, mutant growth was the same as WT . In contrast, the mutant failed to produce FB1 on either living or dead kernels.
In the current study, we describe three new phenotypes attributed to a non-functional FST1. Furthermore, we describe the effects of FST1 on whole genome expression by comparing the transcriptomes of the WT and ∆fst1 strains of F. verticillioides grown on autoclaved maize kernels. The results support our hypothesis that FST1 has a regulatory function that globally impacts gene expression.
Macroconidia production and sensitivity to H2O2
Effect of Δfst1 on conidiation
90 ± 7
1,252 ± 118
9 ± 3*
1,052 ± 123
75 ± 7
1,603 ± 163
Analysis of transcriptome
Summary of RNAseq data from Illumina sequencing a
Expression of tubulin and elongation factor (EF) genes during colonization of autoclaved maize kernels by strains Δfst1 and WT a
Tubulin alpha chain
Tubulin alpha chain
Tubulin beta chain
Tubulin gamma chain
Molecular function ontology of differentially expressed genes in WT and Δfst1 during colonization of autoclaved maize kernels a
Molecular Function b
Up in ∆fst1
Down in ∆fst1
Not expressed in WT
Not expressed in ∆fst1
Peptidases and Proteases
Integral Membrane Proteins
FUM gene cluster
Comparison of expression of FUM genes in wild type (WT) and Δfst1 a
FVEG ID b
P value c
Log 2 fold change d
Expression of selected genes in strain Δfst1 relative to expression in wild type (WT) of F. verticillioides
FVEG ID a
Relative Expression b
−33.3 (−25.8, − 48.5)
−9.4 (−7.6, − 11.7)
−3.9 (−3.7, − 4.1)
−2.0 (−2.0, − 2.0)
−86.7 (−73.3, − 102.6)
−102.1 (−71.0, −147.0)
−119.1 (−115.8, − 122.5)
−4.1 (−4.1, − 4.2)
−12.0 (−10.1, − 14.2)
−32.4 (−21.8, − 47.1)
3.7 (3.3, 4.0)
2.2 (1.9, 2.4)
1.5 (1.5, 1.6)
4.9 (4.8, 4.9)
−76.0 (−73.3, − 78.7)
33.7 (31.8, 35.7)
−34.4 (−23.5, − 50.4)
−2.9 (−2.8, − 3.0)
−2.6 (−2.4, − 2.8)
−22.1 (−20.5, − 24.0)
Comparison of hydrophobin ( HYD ) genes during colonization of autoclaved maize kernels by strains Δfst1 and WT a
FVEG ID b
P value c
Log 2 fold change d
Ma et al.  predicted 683 putative transcription factor (TF) genes in F. verticillioides and Wiemann et al.  predicted 640. Of these predicted TF, our analysis identified 115 differentially expressed (Table 4). Transcription factors in fungi have been classified into 61 families , and we found that 108 of the differentially expressed TF genes were in 12 of the 61 families. Most (80%) of the TFs were C2H2 zinc finger (16 genes) and Zn2Cys6 (76 genes). FUM21 is classified in the Zn2Cys6 family and its expression in Δfst1 was 4.6-fold less compared to that of the WT (Table 5). However, its P-value (0.012) was just outside the threshold we selected for statistical testing.
Classification of putative transporter genes differentially expressed during colonization of autoclaved maize kernels by wild type (WT) and strain Δfst1 a
Transporter Type b
Up in Δfst1
Down in Δfst1
Amino acid related
Major facilitator superfamily
A total of 189 of the differentially expressed genes were categorized with putative oxidase functions (Table 4). Compared to the WT, two-thirds of these genes exhibited reduced expression in strain Δfst1 and the other third were expressed at higher levels. Additionally, the expression of 21 oxidase genes was only measured in the WT and seven only in Δfst1. A word-search of the F. verticillioides genome database identified 30 putative peroxidase and seven catalase genes, and ten of these genes were differentially expressed. The peroxidase genes (POD1 FVEG_10866; POD3 FVEG_12884; POD4 FVEG_12465; FVEG_04790) were all down-regulated as much as 100-fold in strain Δfst1 compared to WT. Four catalase genes (CAT1 FVEG_05529; FVEG_05976; FVEG_03348; FVEG_05591) also were down-regulated in Δfst1. Expression of the putative catalases CAT2 (FVEG_12611) and CAT3 (FVEG_11955) was up-regulated 4-fold and 2-fold, respectively, in strain Δfst1. We used qPCR analysis to measure the expression of peroxidases and catalases in both autoclaved kernels and infected living kernels. In autoclaved kernels, expression of three peroxidase genes (POD1, POD3 and POD4) and three catalase genes (CAT1, CAT2 and CAT3) were found to be similar to expression indicated by the RNAseq results (Table 6). qPCR analysis of the inoculated living kernels indicated similar effects on expression of the peroxidases and catalases (Table 6).
Differences in expression of putative, secreted, cell wall-degradation genes during colonization of autoclaved maize kernels by wild type (WT) and strain Δfst1 a
FVEG ID b
P value c
Log 2 fold change d
Previous studies indicated that deletion of FST1 in F. verticillioides results in reduced fumonisin production and virulence [15,16]. Here we have linked the mutation to increased sensitivity to H2O2, reduced macroconidia production and reduced hydrophobicity. Considering these diverse phenotypes, the goal of this research was to characterize the effects of FST1 on genome-wide expression during colonization of maize kernels. Autoclaved kernels were chosen to eliminate the effects associated with reduced biomass and fungal development caused by the slower growth of the FST1 mutant when inoculated to living kernels. Even without a living host environment, significant changes in transcription were found in the mutant, many of which may contribute to the observed phenotypes.
For our comparison of the transcriptomes of WT and Δfst1, we relied on the F. verticillioides reference genome at the Broad Institute. Recent updates in the annotation of the genome created changes in gene reference numbers and gene identifications. Two changes were important to our study. First, the FUM8 gene (originally: FVEG_00318, GenBank Accession No AAG27130) was separated into two genes: FVEG_14634 and 14635. In the original annotations, FUM8 contained a 2532-bp open reading frame encoding a 839 amino acid protein described as the aminotransferase responsible for the condensation of alanine to the polyketide backbone of B-series fumonisins . The disruption of FUM8 in F. verticillioides, which blocks fumonisin production and mycotoxin production, was recovered in the mutant by complementation with the WT FUM8 gene . In the latest annotation of the genome, the sequence encoding the first 279 amino acids of FVEG_00318 plus 11 additional amino acids was designated as FVEG_14635, and the sequence encoding the last 554 amino acids of FVEG_00318, which contains aminotransferase domain, was designated as FVEG_14634. Regardless of this particular annotation error, expression of FUM8 is significantly reduced in strain Δfst1 along with most of the other FUM genes, confirming the role of FST1 in fumonisin production.
The second peculiar annotation change in the reference genome was that for FST1 (FVEG_08441). Originally listed as a “hypothetical protein”, with similarity to hexose transporters, the gene is now listed as a “myo-inositol transporter”. Inositol is a polyol that functions as an essential constituent of cell membranes as derivatives of phosphatidylinositol and as important cell signaling molecules of inositol phosphates . Two myo-inositol transporter genes have been described in S. cerevisiae by complementation of a strain defective in myo-inositol uptake . A BLAST analysis of the F. verticillioides genome with the yeast ITR1p sequence identified eight genes with high sequence similarity (FVEG_01519, FVEG_01638, FVEG_02081, FVEG_03992, FVEG_06504, FVEG_07757, FVEG_11293, and FVEG_12687). The sequence of FST1 was not identified by the search. Among the eight identified genes, expression was significantly down-regulated in Δfst1 for FVEG_06504 (named ITR1) (19-fold) and FVEG_03992 (5-fold), while the expression of FVEG_12687 was significantly up-regulated (12-fold). We measured the expression of ITR1 by qPCR and verified that its expression was significantly reduced (Table 6). In light of these observations, the assignment of the functional role of myo-inositol transporter to FST1 is premature.
Kim and Woloshuk  described the phenotype of Δfst1 as having slower growth and symptom development, and thus reduced virulence, compared to WT on wound-inoculated maize kernels. This growth inhibition was not observed on autoclaved kernels . We hypothesized that the reduced virulence of Δfst1 resulted from an increased sensitivity to the effects of reactive oxygen species (ROS), which includes H2O2 produced by the living kernel [26,27]. The greater inhibition of the growth of strain Δfst1 by H2O2 compared to the WT and fst1-comp strains supports this hypothesis.
During pathogenesis, F. verticillioides could encounter ROS produced in maize kernels through several independent pathways. Kim and Woloshuk  inoculated the crown of maize kernels at the R4 (dough) stage of development, a period when the endosperm tissue is undergoing program cell death (PCD) . ROS molecules, including H2O2, are produced during PCD in plants  and likely during endosperm development . ROS production is also a characterized response of plants to pathogen invasion and plays a major role in host defense . Most pathogens respond to ROS by the production of peroxidases and catalases . Our RNAseq analysis of F. verticillioides grown on autoclaved kernels identified several putative catalases and peroxidases whose expression was changed in Δfst1 mutants. Four putative peroxidase genes were down-regulated in Δfst1, as were four of the six putative catalases. We also found that these oxidases were similarly affected in living kernels infected with the F. verticillioides strains.
To gain greater insight into a possible function of the catalases and how they may affect virulence, we examined their function in other plant pathogens. Catalases have been separated by phylogenetic analysis into four clades: peroxisomal, cytoplasmic, spore-specific, and secreted . We found sequence similarity in the five differentially expressed catalases from our study when compared to the catalases assigned to the four clades in Giles . FVEG_11955 was most similar to XP324526 in Neurospora crassa and FG02881 in Gibberella zeae, both of which belong to the peroxisomal catalase (clade P). FVEG_05976 was similar to FG05695 in G. zeae, which belongs to the cytoplasmic catalase (clade C). FVEG_05591 was similar to AAK15808 in N. crassa and FG06554 in G. zeae, which belong to the spore-specific catalase (clade A). Sequence analysis of the N-termini of the five predicted catalase proteins indicated that none are secreted.
As mentioned, catalases also have an important role in fungal development, including conidiogenesis. The sequences of the five differentially expressed, putative catalases in F. verticillioides are highly similar to CATB in Magnaporthe grisea, CATA and CATB in A. nidulans, CAT1 and CAT3 in N. crassa, and CATB in Blumeria graminis. In M. grisea, CATB is up regulated during infection of rice . A strain disrupted in CATB was reduced in virulence with increased sensitivity to hydrogen peroxide, and was severely affected in conidia production. In addition, CATA mutants in A. nidulans exhibited reduction in conidiation and increased sensitivity to hydrogen peroxide . The vast majority of conidial produced by F. verticillioides are microconidia. Although the number of macroconidia produced by the WT used in our study comprised only about 7% of the total conidia population, the reduction of macroconidia was consistently observed in strain Δfst1. From our study, it is not possible to determine if the altered expression of the five catalases in strain Δfst1 is responsible for the reduced production of macroconidia.
Aside from the role of catalases in conidial development, transcription factors are known to impact conidiation in fungi, and the expression of several putative TF genes were down-regulated in strain Δfst1. These genes include FVEG_16516 similar to REN1 of Fusarium oxysporum, FVEG_09661 and FVEG_00646 similar to BRLA and ABAA of A. nidulans, respectively, FVEG_12826 similar to FL (fluffy) in N. crassa, and FVEG_06118 similar to FGSG_06160 in F. graminearum. Mutants of REN1 and ABAA fail to produce normal conidia because of developmental malfunctions associated with phialides, the conidiogenous cells [35,36]. Mutants of BRLA fail to produce conidiophores  and FL mutants fail to produce conidia in chains . Furthermore, expression of the conidiation-specific gene CON-10 is not induced in FL mutants of N. crassa. In strain Δfst1, a gene (FVEG_00227) with high sequence identity to CON-10 was down-regulated 14-fold compared to the WT. In F. graminearum, Son et al.  reported that deletion of FGSG_06160 results in a reduction in conidia production but no effect on virulence. We measured the expression of FL-like gene (FLF1 FVEG_12826) by qPCR (Table 6). The expression was 2.9-fold of WT, which is near the 2.2-fold reduction obtained from the RNAseq analysis. These results indicate that reduced expression of one or more of these TFs may impact production of macroconidia but not microconidia.
Hydrophobins are another family of proteins that are associated with conidiogenesis as well as aerial hypha formation and have been shown to be involved in virulence [17,40-42]. Hydrophobins are separated into two classes based on spacing of cysteine residues and physical characteristics. Class I hydrophobins are highly insoluble proteins that form rodlets, and class II are more soluble and do not form rodlets. Fuchs et al.  predicted that hydrophobin genes in F. verticillioides encode three class I proteins (HYD1 FVEG_03689, HYD2 FVEG_03685 and HYD3 FVEG_06538) and two class II proteins (HYD4 FVEG_01575 and HYD5 FVEG_07695). Examination of the protein sequences derived from HYD6 (FVEG_01573) and HYD8 (FVEG_10008) suggests they are class II and class I hydrophobins, respectively. We could not discern the class of HYD7 (FVEG_09843) based on sequence alignments. Mutants of F. verticillioides with deletions of HYD1 or HYD2 are not defective in radial growth, conidial numbers, or corn seedling infection. However, these mutants fail to form microconidial chains . Expression of these two genes was unaffected in Δfst1 and the strain produced normal microconidial chains. We observed the spreading of droplets of detergent solution placed on the surface of strain Δfst1, suggesting a deficiency in the more soluble class II hydrophobins . The down-regulated expression of HYD4, HYD5 and HYD6 in strain Δfst1 is likely associated with this phenotype.
Previous studies have shown that fumonisin production and FST1 expression are higher in the endosperm than in germ tissues [15,16]. These observations suggest that components within the endosperm provide an environment conducive for the pathogen. Strain Δfst1 grows as well as the WT on autoclaved maize, implying that it produces the secreted enzymes needed to breakdown macromolecules in the kernel and transporters to move nutrients into growing hyphae. However, our transcriptome results indicate that the mutation in FST1 greatly impacts the expression of several genes that encode secreted enzymes. We found that the expression of genes encoding enzymes that degrade complex carbohydrate polymers, which make up host cell walls, was altered in strain Δfst1, but not uniformly. The lack of a growth phenotype in the mutant when grown on autoclaved maize and culture media may reflect functional redundancy in these large gene families . For example, the expression of the alpha-amylase gene FVEG_12957  was not affected in strain Δfst1. Expression of this gene would likely mask the potential effects caused by the down regulation of the two starch degradation genes (FVEG_12681 and FVEG_14136).
In this study, we described three new phenotypes associated with a mutation in FST1 that may contribute to the reduced virulence phenotype, namely the increased sensitivity to hydrogen peroxide, reduction of macroconidia production, and changes in mycelial hydrophobicity associated with Δfst1 mutants. We propose that reduced resistance to H2O2 in Δfst1 may impede the strain’s ability to respond to ROS encountered during pathogenesis. Our analysis of the transcriptomes of WT and Δfst1 indicated that the mutation of FST1 affects the expression of 17% of the genes in F. verticillioides. Among the genes affected were many that impact mycotoxin biosynthesis, virulence, resistance to H2O2, and conidiogenesis. Our study supports the hypothesis that FST1 has a role other than sugar transport. Other researchers have described putative sugar transporters that appear to have broader functions. Mutants of RCO-3 in N. crassa displayed altered responses to increasing glucose concentrations in culture media . The authors suggested that RCO-3 functions as a sugar sensor and a regulator of conidia production. In Magnaporthe oryzae, mutations affecting MOST1 result in reduced conidiation and production of the secondary metabolite melanin . The authors were not able to complement the defects by expression of other sugar transporter genes. Further studies are needed to determine how these genes (including FST1) regulate the function of multiple cell processes.
Fungal strains and culture conditions
Fusarium verticillioides strain 7600 (wild type, WT) is deposited in the Fungal Genetics Stock Center, University of Kansas Medical School, Kansas City, KS, USA. The mutant strain ∆fst1 and corresponding complemented stain fst1-comp were previously described by Bluhm et al. . Cultures were stored long-term in 50% glycerol at −80°C and maintained as working stock on PDA medium (B&D, Sparks, MD).
To assess conidiation, strains were inoculated onto Petri dishes containing 1.5% water agar with six to eight gamma-irradiated carnation leaves (average size 18 mm2) on the agar surface . For each fungal strain, nine carnation leaves were sampled after 7 days of incubation. Individual carnation leaves were placed into 1.5 ml microcentrifuge tubes containing 0.3 ml of water and vortexed briefly. Conidial number was determined with a hemacytometer . Macroconidia and microconidia were recorded as the number of conidia per carnation leaf.
Resistance to hydrogen peroxide was measured as described by Lessing  and Ridenour  with some modifications. Conidia (1 ml of 1 × 106 conidia) were mixed with 20 ml of molten PDA and poured into a Petri plate. After incubation for 24 hours at room temperature, a well was cut into the center of the plate with a cork borer (1 cm). To each well, 200 μl of 15% H2O2 (v/v) was added. Plates were incubated for another 24 hours at room temperature in the dark. Inhibition of growth appeared as a clear zone around the well. The area of the inhibition zone was determined. Test on each fungal strain was replicated three times.
Mycelial hydrophobicity was tested by placing droplets (10 μl) of water or a detergent solution (0.2% SDS, 50 mM EDTA) on the colony surface of strains grown on PDA medium for six days in the dark at room temperature. After 30 minutes, we determined whether or not the droplets maintained their spherical shape on the surface of the mycelium [18,49].
Next-generation sequencing methods were used to obtain transcriptome data from the WT and strain ∆fst1 grown on autoclaved maize kernels. Kernels of maize inbred B73 were submerged in deionized water and autoclaved for 15 min. Afterwards, the kernels were crushed slightly to disrupt the pericarp, and approximately 7 g of kernels (10–12 kernels) were placed in glass vials (20 ml) and autoclaved for 30 min. Four replicate vials of the WT and ∆fst1 were inoculated with 100 μl of 106 conidia/ml. Vials were incubated at 28°C for 6 days, then flash frozen in liquid nitrogen and stored at −80°C.
Total RNA was isolated from the content of each vial as described by Bluhm et al.  and purified with the RNeasy Mini Kit (Qiagen, Valencia, CA, USA). Further purification was achieved by treatment with the DNA-Free RNA kit (Zymo Research, Irvine, CA, USA). The Purdue Genetic Core Facility conducted quality assessment, processing, and sequencing of the RNA. The RNA samples had a RIN (RNA Integrity Number) over 7.0 as determined with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Paired-end sequences were obtained with an Illumina HiSeq 2500 sequencer (Illumina, San Diego, CA, USA). Sequence data were trimmed of adapters and filtered to remove low quality sequence and reads less than 30 nt.
RNA sequence data from each sample were mapped to the reference genome of Gibberella moniliformis, which was downloaded (June 2014) from the Broad Institute Fusarium Comparative Database (http://www.broadinstitute.org). Sequence data were mapped to the reference genome was done with CLC Genomics Workbench (version 7.0.4, CLC Bio, Boston), and gene expression was quantified as reads per kilobase per million mapped reads (RPKM) . Statistical analysis (pairwise t-testing) was also conducted with the CLC Genomics software. Differentially expressed genes between WT and ∆fst1 were sorted to identify those with absolute fold change values of > 2.0 and P value < 0.01. Genes expressed uniquely in each fungal strain were identified also. The selected genes from the differentially expressed and those in the uniquely expressed groups were analyzed for gene ontology (GO). For each gene, the translated sequence was analyzed with Blast2GO (version 2.7.2, Blast2Go.com). Results were sorted with respect to molecular function of the top BLAST descriptors.
Quantitative real time-PCR
Quantitative PCR (qPCR) analysis was conducted on RNA isolated from both autoclaved and living maize kernels. For autoclaved kernels, equal amounts of purified RNA were pooled from the four biological replicates of WT and strain ∆fst1 used in the RNAseq anlaysis. To obtain living kernels, maize B73 was greenhouse-grown and ears were inoculated with the F. verticillioides strains as described by Kim and Woloshuk . Six days after inoculation, infected kernels were collected from three ears (biological replicates) and total RNA was isolated. As with the autoclaved kernels, purified RNA were pooled from the three biological replicates of WT and strain ∆fst1.
PCR primers used in this study
FVEG ID a
Primer Sequence (5’ -> 3’)
CTT CTG ATG CTC TTC TCT TCC TCG C
TCT GGT ATA TCT CAC CAA TGA ACG CGA T
ACA CCA AAG CCT CTA CAA GTG A
AGG TAT CGG GCA CCG CT
TTG CGA GGA TCT GTT CTT CTA TC
TAT TAC CGA GCT TGC GCT ATA C
TCA TTG ACC GTG CTC AAC TCC TCA
TGT CGA GTT GAC GAA GAA GT
TCC TGG AAC AAC TGG AAT GG
CAA TCA AGA CAG ACA GGA GAG G
GGC TAG CTA CAT CCA AGA AGA C
GTA CCA TCA GCC ATG ATC TCA A
GAT CTT CTG GAC CAA CCT CAA T
CCT GAA CTT GGG CTC CTT ATA C
AGA AGA AGG CTG GTG CTA ATG
GGC TCC ATG ACC TGA ACA TAC
GAG CGA CAC GCA AAC CAT TGA AGT
ACC ACC AAC AGT CGA GAT TCG TGT
TGC TCA TTT CCA AGA TCC GCG
GCG CAT GCA GAT ATC GTA GAG G
TTG CTC CAC CAA CTC TTA CTG
GCG TTG ATG TTG ATG AGA GCA
AGC TCT CCG CCA TCT TCT A
GCT CAA TGT CTC TCT CCT CAA C
GTC TCT CCC GTT CAT GAT TCT C
GGG TTG ACT TGG GTG GTA TT
AGC GAT GCT TCT TGT CCT TAC
AAC CAA GCT CAC GAC CTA TTT
GGG ACC TGT TGC CAT TAA GA
TCA TCC TCC GGC ATT TCA TAG
CGT ATG GCA GAG TGG GTA AAT
CAT CAG GAT TCG GAC GGT ATA TG
Availability of supporting data
Supporting sequence data are available in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE66044 (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE66044).
We thank Gregory R. OBrian and Dr. Larry Dunkle for their critical review of the manuscript. We thank Aparna Natarajan, who identified the increased sensitivity to H2O2 phenotype in strain ∆fst1. Funding for this research was provided by USDA/NIFA/AFRI, award number 10-65108-20567 to CPW and GAP.
- Oren L, Ezrati S, Cohen D, Sharon A. Early events in the Fusarium verticillioides–maize interaction characterized by using a green fluorescent protein–expressing transgenic isolate. Appl Environ Microbiol. 2003;69(3):1695–701.View ArticlePubMed CentralPubMedGoogle Scholar
- Hendricks K. Fumonisins and neural tube defects in South Texas. Epidemiology. 1999;10(2):198–200.View ArticlePubMedGoogle Scholar
- Wu F. Mycotoxin reduction in Bt corn: potential economic, health, and regulatory impacts. Transgenic Res. 2006;15(3):277–89.View ArticlePubMedGoogle Scholar
- Wu F. Mycotoxin risk assessment for the purpose of setting international regulatory standards. Environ Sci Technol. 2004;38(15):4049–55.View ArticlePubMedGoogle Scholar
- Brown DW, Butchko RA, Busman M, Proctor RH. The Fusarium verticillioides FUM gene cluster encodes a Zn(II)2Cys6 protein that affects FUM gene expression and fumonisin production. Eukaryot Cell. 2007;6(7):1210–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Visentin I, Montis V, Doll K, Alabouvette C, Tamietti G, Karlovsky P, et al. Transcription of genes in the biosynthetic pathway for fumonisin mycotoxins is epigenetically and differentially regulated in the fungal maize pathogen Fusarium verticillioides. Eukaryot Cell. 2012;11(3):252–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Liu SY, Lin JQ, Wu HL, Wang CC, Huang SJ, Luo YF, et al. Bisulfite sequencing reveals that Aspergillus flavus holds a hollow in DNA methylation. PloS One 2012, 7(1).Google Scholar
- Woloshuk CP, Shim WB. Aflatoxins, fumonisins, and trichothecenes: a convergence of knowledge. FEMS Microbiol Rev. 2013;37(1):94–109.View ArticlePubMedGoogle Scholar
- Kim H, Woloshuk CP. Role of AREA, a regulator of nitrogen metabolism, during colonization of maize kernels and fumonisin biosynthesis in Fusarium verticillioides. Fungal Genet Biol. 2008;45(6):947–53.View ArticlePubMedGoogle Scholar
- Shim WB, Woloshuk CP. Regulation of fumonisin B–1 biosynthesis and conidiation in Fusarium verticillioides by a cyclin–like (C–type) gene, FCC1. Appl and Environ Microbiol. 2001;67(4):1607–12.View ArticleGoogle Scholar
- Tilburn J, Arst HN, Penalva MA. Regulation of gene expression by ambient pH. Cell Mol Biol Filamentous Fungi 2010:480--487.Google Scholar
- Flaherty JE, Pirttila AM, Bluhm BH, Woloshuk CP. PAC1, a pH–regulatory gene from Fusarium verticillioides. Appl and Environ Microbiol. 2003;69(9):5222–7.View ArticleGoogle Scholar
- Bluhm BH, Woloshuk CP. Amylopectin induces fumonisin B–1 production by Fusarium verticillioides during colonization of maize kernels. Mol Plant Microbe In. 2005;18(12):1333–9.View ArticleGoogle Scholar
- Kim H, Smith JE, Ridenour JB, Woloshuk CP, Bluhm BH. HXK1 regulates carbon catabolism, sporulation, fumonisin B(1) production and pathogenesis in Fusarium verticillioides. Microbiol--Sgm. 2011;157:2658–69.View ArticleGoogle Scholar
- Bluhm BH, Kim H, Butchko RAE, Woloshuk CP. Involvement of ZFR1 of Fusarium verticillioides in kernel colonization and the regulation of FST1, a putative sugar transporter gene required for fumonisin biosynthesis on maize kernels. Mol Plant Pathol. 2008;9(2):203–11.View ArticlePubMedGoogle Scholar
- Kim H, Woloshuk CP. Functional characterization of fst1 in Fusarium verticillioides during colonization of maize kernels. Mol Plant Microbe In. 2011;24(1):18–24.View ArticleGoogle Scholar
- Fuchs U, Czymmek KJ, Sweigard JA. Five hydrophobin genes in Fusarium verticillioides include two required for microconidial chain formation. Fungal Genet Biol. 2004;41(9):852–64.View ArticlePubMedGoogle Scholar
- Ridenour JB, Bluhm BH. The HAP complex in Fusarium verticillioides is a key regulator of growth, morphogenesis, secondary metabolism, and pathogenesis. Fungal Genet Biol. 2014;69:52–64.View ArticlePubMedGoogle Scholar
- Ma LJ, van der Does HC, Borkovich KA, Coleman JJ, Daboussi MJ, Di Pietro A, et al. Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010;464(7287):367–73.View ArticlePubMed CentralPubMedGoogle Scholar
- Wiemann P, Sieber CMK, Von Bargen KW, Studt L, Niehaus EM, Espino JJ, et al.: Deciphering the cryptic genome: genome--wide analyses of the rice pathogen Fusarium fujikuroi reveal complex regulation of secondary metabolism and novel metabolites. Plos Pathog 2013, 9(6).Google Scholar
- Park J, Park J, Jang S, Kim S, Kong S, Choi J, et al. FTFD: an informatics pipeline supporting phylogenomic analysis of fungal transcription factors. Bioinformatics. 2008;24(7):1024–5.View ArticlePubMedGoogle Scholar
- Ravalason H, Grisel S, Chevret D, Favel A, Berrin JG, Sigoillot JC, et al. Fusarium verticillioides secretome as a source of auxiliary enzymes to enhance saccharification of wheat straw. Bioresource Technol. 2012;114:589–96.View ArticleGoogle Scholar
- Seo JA, Proctor RH, Plattner RD. Characterization of four clustered and coregulated genes associated with fumonisin biosynthesis in Fusarium verticillioides. Fungal Genet Biol. 2001;34(3):155–65.View ArticlePubMedGoogle Scholar
- Barker CJ, Illies C, Gaboardi GC, Berggren PO. Inositol pyrophosphates: structure, enzymology and function. Cell Mol Life Sci. 2009;66(24):3851–71.View ArticlePubMedGoogle Scholar
- Nikawa J, Tsukagoshi Y, Yamashita S. Isolation and characterization of two distinct myo–inositol transporter genes of Saccharomyces cerevisiae. J Biol Chem. 1991;266(17):11184–91.PubMedGoogle Scholar
- Torres MA, Jones JD, Dangl JL. Reactive oxygen species signaling in response to pathogens. Plant Physiol. 2006;141(2):373–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Heller J, Tudzynski P. Reactive oxygen species in phytopathogenic fungi: signaling, development, and disease. Annu Rev Phytopathol. 2011;49:369–90.View ArticlePubMedGoogle Scholar
- Young TE, Gallie DR. Programmed cell death during endosperm development. Plant Mol Biol. 2000;44(3):283–301.View ArticlePubMedGoogle Scholar
- Van Breusegem F, Dat JF. Reactive oxygen species in plant cell death. Plant Physiol. 2006;141(2):384–90.View ArticlePubMed CentralPubMedGoogle Scholar
- Sabelli PA. Replicate and die for your own good: endoreduplication and cell death in the cereal endosperm. J Cereal Sci. 2012;56(1):9–20.View ArticleGoogle Scholar
- Torres MA. ROS in biotic interactions. Physiol Plant. 2010;138(4):414–29.View ArticlePubMedGoogle Scholar
- Giles SS, Stajich JE, Nichols C, Gerrald QD, Alspaugh JA, Dietrich F, et al. The Cryptococcus neoformans catalase gene family and its role in antioxidant defense. Eukaryot Cell. 2006;5(9):1447–59.View ArticlePubMed CentralPubMedGoogle Scholar
- Skamnioti P, Henderson C, Zhang ZG, Robinson Z, Gurr SJ. A novel role for catalase B in the maintenance of fungal cell–wall integrity during host invasion in the rice blast fungus Magnaporthe grisea. Mol Plant Microbe In. 2007;20(5):568–80.View ArticleGoogle Scholar
- Navarro RE, Stringer MA, Hansberg W, Timberlake WE, Aguirre J. catA, a new Aspergillus nidulans gene encoding a developmentally regulated catalase. Cur Genet. 1996;29(4):352–9.Google Scholar
- Ohara T, Inoue I, Namiki F, Kunoh H, Tsuge T. REN1 is required for development of microconidia and macroconidia, but not of chlamydospores, in the plant pathogenic fungus Fusarium oxysporum. Genetics. 2004;166(1):113–24.View ArticlePubMed CentralPubMedGoogle Scholar
- Sewall TC, Mims CW, Timberlake WE. Abaa controls phialide differentiation in Aspergillus–nidulans. Plant Cell. 1990;2(8):731–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Adams TH, Boylan MT, Timberlake WE. brlA is necessary and sufficient to direct conidiophore development in Aspergillus nidulans. Cell. 1988;54(3):353–62.View ArticlePubMedGoogle Scholar
- Bailey LA, Ebbole DJ. The fluffy gene of Neurospora crassa encodes a Gal4p–type C6 zinc cluster protein required for conidial development. Genetics. 1998;148(4):1813–20.PubMed CentralPubMedGoogle Scholar
- Son H, Seo YS, Min K, Park AR, Lee J, Jin JM, et al. A phenome--based functional analysis of transcription factors in the cereal head blight fungus, Fusarium graminearum. Plos Pathog 2011;7(10):e1002310Google Scholar
- Talbot NJ, Kershaw MJ, Wakley GE, de Vries OMH, Wessels JGH, Hamer JE. MPG1 encodes a fungal hydrophobin involved in surface interactions during infection–related development of Magnaporthe grisea. Plant Cell. 1996;8(6):985–99.View ArticlePubMed CentralPubMedGoogle Scholar
- Wosten HAB. Hydrophobins: multipurpose proteins. Annu Rev Microbiol. 2001;55:625–46.View ArticlePubMedGoogle Scholar
- Wosten HAB, Devries OMH, Wessels JGH. Interfacial self–assembly of a fungal hydrophobin into a hydrophobic rodlet layer. Plant Cell. 1993;5(11):1567–74.View ArticlePubMed CentralPubMedGoogle Scholar
- Wagner A. Distributed robustness versus redundancy as causes of mutational robustness. Bioessays. 2005;27(2):176–88.View ArticlePubMedGoogle Scholar
- Madi L, McBride SA, Bailey LA, Ebbole DJ. rco–3, a gene involved in glucose transport and conidiation in Neurospora crassa. Genetics. 1997;146(2):499–508.PubMed CentralPubMedGoogle Scholar
- Saitoh H, Hirabuchi A, Fujisawa S, Mitsuoka C, Terauchi R, Takano Y. MoST1 encoding a hexose transporter–like protein is involved in both conidiation and mycelial melanization of Magnaporthe oryzae. FEMS Microbiol Let. 2014;352(1):104–13.View ArticleGoogle Scholar
- Fisher NL, Burgess LW, Toussoun TA, Nelson PE. Carnation leaves as a substrate and for preserving cultures of Fusarium species. Phytopathology. 1982;72(1):151–3.View ArticleGoogle Scholar
- Aberkane A, Cuenca–Estrella M, Gomez–Lopez A, Petrikkou E, Mellado E, Monzon A, et al. Comparative evaluation of two different methods of inoculum preparation for antifungal susceptibility testing of filamentous fungi. J Antimicrob Chemother. 2002;50(5):719–22.View ArticlePubMedGoogle Scholar
- Lessing F, Kniemeyer O, Wozniok I, Loeffler J, Kurzai O, Haertl A, et al. The Aspergillus fumigatus transcriptional regulator AfYap1 represents the major regulator for defense against reactive oxygen intermediates but is dispensable for pathogenicity in an intranasal mouse infection model. Eukaryot Cell. 2007;6(12):2290–302.View ArticlePubMed CentralPubMedGoogle Scholar
- Yan X, Li Y, Yue XF, Wang CC, Que YW, Kong DD, et al. Two novel transcriptional regulators are essential for infection--related morphogenesis and pathogenicity of the rice blast fungus Magnaporthe oryzae. Plos Pathog 2011;7(12):e1002385Google Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA–Seq. Nat Methods. 2008;5(7):621–8.View ArticlePubMedGoogle Scholar
- Reese BN, Payne GA, Nielsen DM, Woloshuk CP. Gene expression profile and response to maize kernels by Aspergillus flavus. Phytopathology. 2011;101(7):797–804.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real–time quantitative PCR and the 2(Delta Delta C(T)) Method. Methods. 2001;25(4):402–8.View ArticlePubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.