- Research article
- Open Access
Environmental pH modulates transcriptomic responses in the fungus Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus
BMC Genomics volume 19, Article number: 721 (2018)
The Ambrosia Fusarium Clade phytopathogenic Fusarium fungi species have a symbiotic relationship with ambrosia beetles in the genus Euwallacea (Coleoptera: Curculionidae). Related beetle species referred to as Euwallacea sp. near fornicatus have been spread in California, USA and are recognized as the causal agents of Fusarium dieback, a disease that causes mortality of many plant species. Despite the importance of this fungi, no transcriptomic resources have been generated. The datasets described here represent the first ever transcripts available for these species. We focused our study on the isolated species of Fusarium that is associated with one of the cryptic species referred to as Kuroshio Shot Hole Borer (KSHB) Euwallacea sp. near fornicatus.
Hydrogen concentration is a critical signal in fungi for growth and host colonization, the aim of this study was to evaluate the effect of different pH conditions on growth and gene expression of the fungus Fusarium sp. associated with KSHB. An RNA-seq approach was used to compare the gene expression of the fungus grown for 2 weeks in liquid medium at three different pH levels (5.0, 6.0, and 7.0). An unbuffered treatment was included to evaluate the capability of the fungus to change the pH of its environment and the impact in gene expression. The results showed that the fungus can grow and modulate its genetic expression at different pH conditions; however, growth was stunted in acidic pH in comparison with neutral pH. The results showed a differential expression pattern in each pH condition even when acidic conditions prevailed at the end of the experiment. After comparing transcriptomics data from the three treatments, we found a total of 4,943 unique transcripts that were differentially expressed.
We identified transcripts related to pH signaling such as the conserved PAL/RIM pathway, some transcripts related to secondary metabolism and other transcripts that were differentially expressed. Our analysis suggests possible mechanisms involved in pathogenicity in this novel Fusarium species. This is the first report that shows transcriptomic data of this pathogen as well as the first report of genes and proteins involved in their metabolism identifying potential virulence factors.
Through the infection process, phytopathogenic fungi can sense and adapt to its host environment to ensure growth, reproduction, and virulence. Hydrogen ion concentration (pH) is a critical host environment variable which can limit or enhance the growth of the phytopathogenic fungi and its ability to compromise the host defense response through the secretion of pathogenicity factors such as hydrolytic enzymes [1,2,3] and the production of toxins [4, 5]. Because of this critical feature, some fungi have evolved to modify the host environment by actively secreting acids or alkali. For example, Fusarium oxysporum can alkalinize the environment by secreting ammonia , while the secretion of fusaric acid gives rise to extracellular acidification [4, 7, 8]. Interestingly, it has been reported that fusaric acid, through nitric oxide signaling, induces programmed cell death in their plant hosts . These studies show that some pathogens have the ability to modulate environmental pH .
Studies have also shown that filamentous fungi have the potential to sense the concentration of hydrogen ions and respond with regulatory systems such as RIM or PAL in yeasts or filamentous fungi, respectively . PAL proteins sense alkaline pH and transmit this external stimulus in two consecutive steps that produce proteolytic cleavage of a zinc finger transcription factor named PacC in filamentous fungi and RIM101 in yeasts . The active PacC is translocated to the nucleus where it can recognize the consensus DNA site 5’-GCCAAG-3′ and consequently up- or down-regulate the transcription of acidic- or alkaline-pH responsive genes, respectively [12, 13]. Some pathogenic factors regulated by the Pal/Rim system in fungi include: (i) pectate lyase (pelB) in Colletotrichum gloeosporioides , (ii) endoglucanases in Alternaria alternata , (iii) synthesis of extracellular polysaccharides and oxalic acid as well as production of polygalacturonase and laccase enzymes in Botrytis cinerea  and others .
Fusarium is a cosmopolitan genus that includes species of necrotrophic fungi that cause serious phytosanitary problems, and there is evidence that some relevant virulence factors are regulated by external pH, both in a PacC-dependent and -independent manner. For example, the homolog of PacC in Fusarium graminearum, FgPac1, regulates Tri gene cluster expression and the trichothecene production . In Fusarium verticillioides, the fumonisin production is enhanced in acidic pH and is negatively regulated by the PacC homolog . In addition, two relevant endo-polygalacturonase (pg1 and pg5) are downregulated by a mechanism independent of PacC in F. oxysporum at alkaline pH .
Recently, two Fusarium species belonging to the Ambrosia Fusarium Clade (AFC) have been discovered. These species have a symbiotic relationship with cryptic species of the Asian beetle referred to as Euwallacea sp. near fornicatus and have been spread in California, USA. Moreover, they have been recognized as the causal agents of Fusarium dieback in forests, landscape trees, and in the agricultural sector, causing significant losses in the avocado industry [18, 19]. Fusarium euwallaceae is known to be associated with the Polyphagous Shot Hole Borer (PSHB)  and Fusarium sp. is associated with the Kuroshio Shot Hole Borer (KSHB) that recently has been named F. kuroshium sp. . In 2016, KSHB was also discovered in Baja California, Mexico , and despite the considerable interest in Fusarium species associated with ambrosia beetles, no transcriptomic data are available to help understand the nature of these pathogens, as well as the virulence factors required for disease establishment and development. In this sense, an effective strategy to evaluate the effect of pH on the behavior of a fungus is to challenge it at different pH levels and to analyze their effect on conidiation, radial growth, or fresh weight , to identify specific responses at the gene level , as well as to evaluate large-scale responses .
In this study, we evaluated the ability of this fungus to growth at different pH levels and analyzed the associated transcriptomic responses. We used the Fusarium sp. associated with KSHB (strain HFEW-16-IV-019) and cultured it at distinct pH conditions (5.0, 6.0, and 7.0) for 2 weeks. Following that, we used an RNA-seq approach to identify pH-responsive genes including the signaling route Pal/Rim. In addition, we identified some known elements (genes/proteins) in the fungi virulence arsenal and evaluated the expression profile at different pH levels.
Besides, those nr-unigenes that showed the most notable change and homologous to those genes involved in the biosynthesis of polyketide-derived mycotoxins were identified. Additionally, to the canonical genes, we were able to identify some nr-unigenes related with pathogenicity that for the first time has been recognized in this fungus. The innovation of these data implies the first transcriptomic knowledge of this important phytopathogen that is a symbiont of the beetle Euwallacea sp. near fornicatus, giving interesting clues about the molecular mechanisms that govern its growth and adaptation to abiotic stress. Finally, this is the first report that documents the ability of this fungus to grow in different pH conditions and some transcriptomic responses associated with this adaptation.
Effect of pH on conidiation and growth of Fusarium sp. associated with KSHB
To characterize the influence of pH on fungal conidiation and growth, Fusarium sp. associated with KSHB was grown in liquid media (PDB) and microscopic observations were recorded, and fresh weight was measured. The same strain was grown on solid (PDA) culture media, and the radial growth and the colony phenotype differences were observed. After 9 days of incubation at 28 ± 2 °C in complete darkness, Fusarium sp. associated with KSHB showed markedly reduced radial growth at pH 5.0 in comparison with pH 6.0 and 7.0, whereas pH 7.0 is the condition with the highest observed colony diameter (Fig. 1). Colony color on PDA plates (obverse view) was white and pale yellow at pH 5.0 and 6.0, while at pH 7.0, the color was purplish-gray (Fig. 1a). The reverse view of the plates showed a colony-pigmentation that ranged from pale yellow (pH 5.0 and 6.0) to brownish orange (pH 7.0) (Fig. 1b).
The mycelial growth in liquid media at different pH levels was also evaluated by recording differences in fresh weight. The highest growth was at pH 7.0 (13.09 ± 0.49 g), whereas it was intermediate at pH 6.0 (10.69 ± 0.53 g) and was the lowest at pH 5.0 (8.29 ± 0.16 g). At the end of the experiment (2 weeks post-inoculation), the pH of the liquid medium was measured. Interestingly, regardless of the initial pH condition, in all the cases the pH was more acidic than the initial condition. The final values of the initial pH 5.0, 6.0, and 7.0 conditions were 3.8 ± 0.04, 4.35 ± 0.14, and 4.84 ± 0.20, respectively. The conidia analysis indicated that at pH 7.0 the production of the particular conidia of this species was induced (a dolphin-like appearance; ). These kinds of conidia are swollen in their upper half commonly with 4 or 5 septa. The conidia at pH 5.0 and 6.0 showed no differences and were small with 1 or 2 septa, and no dolphin-like conidia were detected (Fig. 1c).
Construction of the unigenes set for Fusarium sp. associated with KSHB species
A total of 151,009,935 paired-end reads were generated from nine distinct libraries that represent each of the evaluated conditions (pH 5.0, 6.0, and 7.0) which were sequenced in triplicate. On average, each of the three sequenced replicates contributed approximately 17 million paired reads. Before de novo assembly, raw reads were filtered to only include those with acceptable quality (Additional file 1: Table S1; see Methods for more details). Longer reads were generated by merging paired-end reads with overlapping regions. Both data sets (long reads and no merged paired-end reads) were used to assemble the transcriptome of Fusarium sp. associated with KSHB (Additional file 1: Table S1). A total of 37,550 unigenes were generated ranging from 200 to 24,657 bp, with an average length of 1,900.81 bp (Additional file 2).
Out-of-frame insertions/deletions in coding regions were corrected using the AlingWise pipeline  which correct the erroneous frame shifts through a homology-based method (see Methods section for more details). After correcting the frame-shifts, redundant sequences were eliminated. We considered a sequence redundant if presented an identity of at least 95% over 90% or more of the length of the sequences against which it was compared. The corrected and non-redundant (nr) unigenes set comprises a curated transcripts collection of 25,070 sequences representative of the Fusarium sp. associated with KSHB transcriptome with identified and translated open reading frames (ORFs) ranging between 15 to 5,692 amino acids (Additional file 1: Table S2). Only these sequences were considered for future analysis.
Homology search and functional annotation of Fusarium sp. associated with KSHB unigenes
The annotation process for nr-unigenes included a functional classification and a search for similar sequences in sets of predicted proteins of the filamentous fungi whose genomes have been completely sequenced and that we considered closely related species or representatives of specific classes. The species included in the comparisons which were conducted using the Basic Local Alignment Search Tool (BLAST; ) were: Botrytis cinerea [28, 29] [Class: Leotiomycetes]; Symbiotaphrina kochii [Class: Xylonomycetes]; Acremonium alcalophilum, Fusarium fujikuroi , Fusarium graminearum , Fusarium oxysporum , Nectria haematococca (; commonly referred to by its asexual name Fusarium solani), Fusarium verticillioides [31, 32], Neurospora crassa  and Verticillium dahliae . With only a few exceptions (S. kochii, A. alcalophilum, and N. crassa), most of the selected fungi are considered highly pathogenic in a wide range of crops. Species such as S. kochii were included because Symbiotaphrina is a yeast-like genus of endosymbionts of beetles which have been suggested to be implicated in B-vitamin biosynthesis, fatty acid metabolism, and detoxification of noxious plant compounds by insects . On the other hand, A. alcalophilum is a cellulolytic fungus whose optimal growth depends on alkaline conditions . N. crassa was included because it is the model system for the understanding of light responses in fungi . Also, Saccharomyces cerevisiae [Class: Saccharomycetes]; was included because it has the best-annotated genome from the Ascomycota phylum [39, 40].
As expected, once unigenes were corrected and filtered, most of the proteins produced by translated the ORFs, showed high similarity (e-value < 10− 5) with proteins from at least one of the ten species against which they were compared (Additional file 1: Table S2). In more than 90% of the Fusarium sp. associated with KSBH nr-unigenes set, homologous proteins were detected from other species of Fusarium genus. In minor proportion and decreasing order, homologous proteins were detected from V. dahliae, N. crassa, A. alcalophilum belonging to Sordariomycetes, B. cinerea belonging to Leotiomycetes, S. kochii belonging to Xylonomycetes, and S. cerevisiae belonging to Saccharomycetes (Fig. 2). Fusarium sp. associated with KSBH proteins produced by translated ORFs were classified as ‘complete’ or ‘partial’. Proteins were considered complete if they covered at least 70% the length of four or more identified homologous proteins. Based on this criterion, 9,392 nr-unigenes (37.5% of the total), may represent full-length cDNAs from Fusarium sp. associated with KSBH (see Additional file 1: Table S2).
Gene ontology (GO) terms , enzyme commission (EC) number , the KOG identification number (from Eukaryotic Orthologous Groups of proteins ) and the most representative InterPro (IPR) domains [44,45,46], were assigned to each nr-unigenes based on the information available from the homologous detected (Additional file 3: Table S3). Of the three major GO annotation categories, molecular functions comprised 66.52% of the total assigned annotations, whereas biological processes and cellular components comprised 22.50% and 10.99%, respectively. The GO terms with the largest number of assigned sequences in the molecular function category were zinc ion binding (2,157 nr-unigenes), catalytic activity (2,087), ATP binding (1,826), and transcription factor activity (1,310). Out of the biological processes, the terms with the most sequences were metabolic processes (2,228), regulation of transcription/DNA-dependent (1,417), transport (1,247), and electron transport (766). In the cellular components’ category, the terms with the most sequences were intracellular (1,210), integral to membrane (1,184), membrane (990), and nucleus (643). Figure 3 shows the top ten most represented GO terms from each of the three major categories: molecular functions, biological processes, and cellular components.
Metabolic pathway assignments
Regarding the representation or coverage of metabolic pathways represented in Fusarium sp. associated with KSBH transcriptome, EC numbers assigned to nr-unigenes were used to generate a metabolic network for this species. This Fusarium sp. associated with KSBH metabolic network was compared against a background metabolic network which was generated combining the total of EC numbers assigned to each of the proteins presents in the genomes of distinct species of the genus Fusarium (F. fujikuroi, F. graminearum, F. oxysporum, F. solani and F. verticillioides). Metabolic networks were reconstructed using “Search & Color Pathway” tool (http://www.genome.jp/kegg/tool/map_pathway2.html) available in the tools suite of KEGG database [47, 48]. The background network of the Fusarium genus comprises a total of 115 metabolic pathways in which 4,174 EC numbers were assigned. Although the metabolic network of Fusarium sp. associated with KSBH shows the same number of metabolic pathways to those present in the Fusarium background metabolic network, these pathways are represented by slightly fewer numbers of enzymes (3,372 EC numbers, 80% of the total; Additional file 4: Figure S1 and Additional file 3: Table S4). These data suggest that the presented transcriptome comprises a deep coverage of the protein coding-genes present in the genome of Fusarium sp. associated with KSBH.
Identification of orthologous groups across Fusarium species
OrthoMCL is an analysis pipeline that uses the reciprocal BLAST and the Markov CLuster (MCL) algorithm to infer and group orthologs (and paralogs) across multiple taxa (see Methods for more details). In addition to homolog searches presented and discussed before, an orthoMCL clustering analysis was performed to identify and group (putative) orthologs across Fusarium species. S. cerevisiae was also included in the OrthoMCL analysis because it remains the best available annotated genome in the Ascomycota phylum. A total of 17,908 ortholog groups were identified among compared species (Additional file 3: Table S5). These OrthoMCL-defined protein families group a total of 91,680 proteins (~ 82% of the total included in the analysis). 17,791 of Fusarium sp. associated with KSBH nr-unigenes (70.96% of the total), were assigned to at least one these groups. In total, 6,424 proteins (grouped in a total of 6,205 families) were common among all Fusarium species (Additional file 4: Figure S2 and Additional file 3: Table S5). The ortholog identification made possible the non-erroneous assignation of the functions and expression profiles to each of the nr-unigenes that were analyzed and discussed in sections described below.
Phylogenetic relationship with additional strains of the Ambrosia Fusarium Clade
Fusarium species that maintain mutualistic association with ambrosia beetles are grouped in the Ambrosia Fusarium Clade (AFC), a strongly supported monophyletic group composed of numbered species without a name [19, 49]. This monophyletic group (AFC) was resolved within the F. solani species complex (FSSC) using a four-locus phylogenetic analysis. The locus included in the analysis were: (i) a portion of the nuclear ribosomal RNA gene repeat (rDNA) comprising the internal transcribed spacer region (ITS) together with domains D1 and D2 of the nuclear ribosomal large subunit (LSU, 1,004 bp total alignment); (ii) translation elongation factor 1 alpha (EF-1α, 687 bp alignment); (iii) DNA-directed RNA polymerase II subunit 1 (RPB1, 1,588 bp alignment) and (iv) DNA-directed RNA polymerase subunit 2 (RPB2, 1,635 bp alignment) . To infer the phylogenetic relationship of Fusarium sp. associated with KSBH, strain HFEW-16-IV-019, which was used in the presented work, we first identify the corresponding nr-unigenes to those homologous used in the previously reported four-locus analysis [21, 49]. With only one exception (EF-1α, UN25284), the identified nr-unigenes represent complete coding sequences from each of the locus mentioned above (UN01385 for RB1, UN04406 for RB2, and UN02180 for LSU). The four loci were concatenated into a single matrix (4,703 bp aligned) and maximum likelihood phylogenetic tree (Additional file 4: Figure S3) was generated for a total of 49 sequences which were alignment (Additional file 5). The strain HFEW-16-IV-019 was brought together with the two sister clades which group the isolates of Fusarium sp. previously designated to particular genealogical lineage named AF-12. The later lineages represent Fusarium species that come from Euwallacea sp. and were isolated either from the beetle or tree galleries of the host species identified in California, USA. The strain HFEW-16-IV-019 recently included in the four-locus analysis and isolated from beetles or elm (Ulmus sp.) tree galleries from Tijuana, Baja California, Mexico was not only grouped with the AF-12 lineage but also with AF-13 lineage, which is represented by Fusarium species isolated from the Euwallacea sp. from Taiwan (Additional file 4: Figure S3).
Differential massive expression analysis of pH-responsive genes in Fusarium sp. associated with KSBH
To compare the expression levels of the nr-unigenes among samples (pH 5.0, 6.0 and 7.0) and identify those which are differentially expressed, sequences count-based differential expression analysis was performed. RSEM  and DEseq  were the programs used for these purposes. First, an expression profile matrix was created containing the number of expected reads for each of the 25,070 nr-unigenes (rows) in each of conditions analyzed (columns). Also, the normalized values TPM (transcripts per million) and FPKM (fragments per Kilobases of contigs/unigenes for per million mapped reads) were calculated (Additional file 3: Table S6). Considering that recently it has been reported that FPKM is an inconsistent measure to compare among samples [52, 53], in this work the TPM values were chosen as the representative value of the expression profiles. Differentially expressed genes were identified comparing the expected read counts normalized by the DEseq package between each pair of samples. In total, 4,785 genes (Additional file 3: Table S7) were identified as differentially expressed in at least one of the three comparisons that were performed (pH 5.0/pH 6.0, pH 5.0/pH 7.0, and pH 7.0/pH 6.0).
The Venn diagrams generated with the differentially expressed genes show that most of these genes are common between pH 5.0/pH 6.0 and pH 5.0/pH 7.0 comparisons (Fig. 4a). On the other hand, only a few differentially expressed genes are identified between the pH 7.0 and pH 6.0 treatments. Hierarchical clustering and heatmap analyses show similar expression patterns with subtle changes between the pH 6.0 or pH 7.0 treatments. Interestingly, both treatments show significant differences when they are compared with pH 5.0 (Fig. 4b and c).
As expected, after surveyed the list of differentially expressed genes, we found many orthologous genes related with the Pal/Rim system and involved in pH signaling. In this sense, Fig. 5 illustrates all these genes in the canonical Pal/Rim pathway [54, 55] and show the expression profiles in Fusarium sp. associated with KSBH at each pH condition. PacC, which is a transcriptional activator of alkaline-induced genes and a transcriptional repressor of acid-induced genes , was not selected as a differentially expressed unigene; however, considering that is a key player of the PAL/Rim system, we identified the corresponding ortholog in order to analyze its expression profile. The corresponding ortholog (UN05861) shows a similar expression profile in all pH conditions analyzed (Fig. 5 and Additional file 3: Table S8).
PHO85 (UN07212), an additional protein also related to the Pal/Rim pathway was also found in the list of differentially expressed genes showing higher transcript levels at pH 6.0 and 7.0 (Fig. 5 and Additional file 3: Table S8). PHO85 is a cyclin-dependent kinase which plays important roles in a number of cellular responses including carbon source utilization, glycogen metabolism, morphogenesis, cell cycle progression, phosphate metabolism and pH tolerance [56,57,58].
Additionally, we were interested in identifying genes that showed the greatest fold change value between different pHs, even if they were not the most abundant. In this sense, the nr-unigenes identified as UN01196, UN01849, UN03862, and UN02411 were the genes that most change their expression during the transition from pH 6.0 to 5.0 (upregulated at pH 5.0) with fold-change values (on the log2-scale) of 12.36, 12.22, 12.10 and 11.93, respectively. On the other hand, the genes that most change during the transition of pH 5.0 to 6.0 (downregulated at pH 5.0) were: UN01657, UN16350, UN13020 and UN04495 with values of log2 fold-change of 12.14, 11.07, 10.85 and 10.54, respectively.
For the nr-unigenes that were most changed during the transition of pH 7.0 to 6.0, the transcriptomic analysis revealed that UN12294, UN07751, and UN06565 are downregulated at pH 7.0 (log2 fold-change values of − 10.84, − 9.81 and − 9.08, respectively), while UN11374, UN12042, UN12867, UN02745, and UN10679 are highly upregulated (log2 fold-change values of 10.08, 9.83, 9.32, 9.00, and 8.36, respectively; see Additional file 3: Table S8 for more details).
To confirm the expression patterns determined by RNA-sequencing (RNA-seq) analysis, we used qRT-PCR to analyze the expression of eight randomly selected genes. Although the expression values of the eight genes showed slight variations as compared to the corresponding values from the qRT-PCR analyses, the expression profile obtained from the RNA-Seq analysis were highly similar to those obtained from qRT-PCR assays (Additional file 4: Figure S4), indicating that the deep sequencing and qRT-PCR data were consistent.
Pathogenicity related to environmental alkalization/acidification
It has been suggested that external pH regulates pathogenicity both in a dependent or independent PAL/RIM manner. Because of this, pH-responsive nr-unigenes were compared by BLASTp algorithm against the pathogen-host interaction database (PHI-base) v.3.5 . The PHI-base groups experimentally verified pathogenicity and virulence genes of different pathogens  and contain manually curated molecular and biological information, along with the descriptions about mutant phenotypes associated with gene information. Moreover, F. oxysporum and F. graminearum are comprised in this database . Interestingly, nearly 40% of Fusarium sp. associated with KSBH differentially expressed nr-unigenes (1,846 from 4,875), found significant similarity (e-value 10− 5) to some of the genes reported in PHI-base (Additional file 3: Table S9), suggesting an important role in pathogenesis (Additional file 4: Figure S5). For instance, transcript levels of LaeA protein (UN02482 and UN02483) increased in all treatments (Fig. 6a). LaeA is a methyltransferase considered a key master regulatory protein of secondary metabolism associated with the epigenetic control.
Both LaeA homologs from Fusarium sp. associated with KSBH (UN02482 and UN02483) were 100% identical to each other, and they shared 30% of identity with the protein FVEG_00539. This protein is represented in an F. verticillioides microarray which made possible the identification of a putative LaeA homolog of F. fujikuroi . Both candidates show a conserved protein methyltransferase s-adenosyl methionine (SAM) binding site, which has been previously identified in LaeA proteins from distinct ascomycetes species  (Fig. 6b).
Consistently, homolog genes to polyketide synthases PKS-1 (UN01158, UN03788, and UN04670), and to FUM1 (UN01196, UN06387, and UN13689) were identified as differentially expressed. Both PKS-1 and FUM1, are polyketide synthases (PKS) involved in melanin  and fumonisin  biosynthesis respectively and have been considered as relevant virulence factors [66, 67]. Regarding nonribosomal peptides, we also found some differentially expressed Fusarium sp. associated with KSBH nr-unigenes UN01784 and UN10282 which are homologs to different nonribosomal peptides synthases (NRPS). Expression profiles from PKS and NRPS at the different pH analyzed are shown in Fig. 6c.
The secondary metabolites produced by PKSs or NRPSs, such as the toxins, need to be exported to overcome the intercellular accumulation. Because of this, the fungi have developed efflux pumps powered by ATP-binding cassette transporters (ABC). We found 13 genes differentially expressed that are homologous to ABC-transporters (Additional file 3: Table S9) specifically, the nr-unigenes UN00536, UN00589 and UN03631 were induced at acidic pH.
The host pH modulation has also been considered as a pathogen’s strategy to promote: i) better hydrolytic activities, ii) permeability increment of the host membranes and/or iii) induce the PAL/RIM system described above. In this sense, one of the mechanisms that fungi use to alkalize the medium is through the secretion of ammonia by the glutamate dehydrogenase (GDH2) enzyme. The GDH2 (homolog to UN02168, UN02220, UN02850, and UN02918) is induced at acidic pH. Also, the MEPB transporter has a key role in ammonia secretion, and we found two nr-unigenes (UN11029 and UN12518) that were induced at acidic (pH 5.0 and 6.0, respectively).
Proteases are enzymes that have an important role during infection because these are capable of mediating many cellular functions such as the adhesion of host cells, the initial penetration of the plant cell wall, and colonization and the disease establishment. In this sense, using the PHI database analysis, we obtained 26 nr-unigenes of proteases that have different and particular gene expression among pH conditions (Additional file 3: Table S9).
Finally, an interesting finding was the induction of the putative homolog to the effector GIP2 of Phytophthora sojae (UN19266) at neutral pH. The GIP2 is an inhibitor protein that share structural similarity with trypsin class of serine proteases but are nonfunctional. Nevertheless, GIP2 inhibits the soybean endo-β-1,3-glucanase, an enzyme that function as a defense response of the plant. If this finding is biologically functional, it will be the second case that describes a counter defensive weapon used by plant pathogens to suppress a plant defense response .
Identification of Fusarium sp. associated with KSBH nr-unigenes with potential relevance to the biosynthetic pathway of fusaric acid
The FA biosynthetic pathway in Fusarium species was proposed by Brown et al.  and comprise a condensation of three acetate units to form a fully reduced 6-carbon polyketide chain which after assimilation of nitrogen from glutamine or oxaloacetate leads to FA formation (Fig. 7). From oxaloacetate to FA, several homologs/orthologs to genes involved in these pathways were identified in the Fusarium sp. associated with KSBH transcriptome.
Based on the EC numbers assigned during the annotation process, the seven key enzymes involved in the oxaloacetate biosynthetic pathway were identified in the nr-unigenes: (1) citrate synthase [EC: 18.104.22.168], (2) aconitate hydratase [EC: 22.214.171.124], (3) isocitrate dehydrogenase [EC: 126.96.36.199], (4) succinyl-CoA synthetase [EC: 188.8.131.52], (5) succinate dehydrogenase [EC: 184.108.40.206], (6) fumarate hydratase [EC: 220.127.116.11] and (7) NAD-dependent malate dehydrogenase [EC: 18.104.22.168] (Additional file 4: Figure S6 and Additional file 3: Table S10).
The OrthoMCL-defined protein family in which well characterized FUB1 proteins were grouped comprise a total of 76 orthologous proteins: 12 from Fusarium sp. associated with KSBH, 15 from F. fujikuroi, 15 from F. verticillioides, 12 from F. oxysporum, 10 from F. graminearum and 12 from F. solani, respectively (Additional file 3: Table S11). These data suggest that Fusarium species contain a family of PKS which varies slightly in the number of members that compose it. Considering that for Fusarium sp. associated with KSBH some of these PKS are represented by partial OFRs, we looked for their complete coding sequences in the genome of Fusarium sp. associated with KSBH which has been recently announced  and was generated from the same strain used in this study.
Figure 8 shows the phylogenetic tree constructed by maximum likelihood (ML) for all Fusarium PKS proteins which aligned (Additional file 6). Three major clades were recognized on the phylogenetic tree. Clades I and II brought together at least one PKS from most of Fusarium species (only F. solani was absent on these clades). The third clade (clade III) which included most of the PKS, were spliced in two well-recognized sub-clades (A and B, respectively). Orthologues to FUB1 are contained in clade III-B. These FUB1-like proteins (one for each Fusarium species) were independently aligned and represented in a new ML phylogenetic tree (Fig. 9a). The multiple sequence alignments show that FUB1 proteins across Fusarium species have identities which ranged from 30 to 95% and they have conserved common motifs which were identified and classified using the Pfam  database (Fig. 9b and c).
Regarding the expression profile of the FUB-like Fusarium sp. associated with KSBH genes which are orthologous to those reported in other Fusarium species, we noticed that none of them were identified as differentially expressed. Interestingly, fub1(UN19324), fub11(UN12405), and fub12(UN05786), show higher transcripts levels at neutral or alkaline pH (pH 6.0 and 7.0) while fub3, fub7, and fub9 present greater amounts of transcripts in the acidic treatment (pH 5.0) showing a negative correlation with pH level (Additional file 4: Figure S7). With only one exception (fub1), all FUB-like nr-unigenes from Fusarium sp. associated with KSBH show low transcript values (TPM values less than 20). Furthermore, the TPM values for fub1 are close to zero across all evaluated pH conditions. These results were expected, considering that FA contributes to acidification of the growth medium and the pH was similar in all the treatments when the tissue was harvested to isolate the RNA.
Identity of the strain used in this study
Based on the phylogenetic analysis, the strain HFEW-16-IV-019 used in this study is grouped with isolates corresponding to Fusarium sp. associated with KSBH that were collected in the locations of El Cajon, Fallbrook, Bonsall and Escondido, CA. It is important to mention that in a previous work  we named this strain as F. euwallaceae because until that time there were not enough genomic data available that allow us to infer the differences between the Fusarium strains associated with KSBH or PSHB .
Effect of pH on conidiation and growth of Fusarium sp. associated with KSHB
Acidic/alkaline requirement for fungi growth is quite broad, ranging from pH 3.0 to more than pH 8.0, with the optimum around pH 5.0 if nutritional requirements are satisfied . Some Fusarium species (e.g. F. oxysporum), have shown the ability to grow in acid conditions (pH 2.0) with an optimal value of pH 6.0 for growth and sporulation . In contrast, our results showed that a highly acid medium (pH 3.0 or less) completely inhibited the growth of Fusarium sp. associated with KSHB.
Differential massive expression analysis of pH-responsive genes in Fusarium sp. associated with KSBH
Hierarchical clustering and heatmap analyses show similar expression patterns with subtle changes between the pH 6.0 or pH 7.0 treatments. Interestingly, both treatments show significant differences when they are compared with pH 5.0 (Fig. 4 and Additional file 3: Table S7). These data suggest that molecular responses of Fusarium sp. associated with KSBH are similar when the fungus is growing at pH 7.0 or 6.0, while significant differences can be observed when the fungus is growing at pH 5.0.
The fungi have developed a signal transduction pathway in order to sense external pH designated as RIM or PAL for yeasts or filamentous fungi, respectively. We found many orthologous genes related with the Pal/Rim system and involved in pH signaling during the differential massive analysis in Fusarium sp. associated with KSBH. The activation of the PAL/Rim system requires firstly the perception of the external pH by a molecular complex located in the plasma membrane that includes PalH/Rim21 (UN01086) a 7-transmembrane domain receptor which senses alkaline pH  and facilitate ubiquitination of the arrestin-like protein PalF/Rim8 (UN03271) . These proteins, as well as the chaperone PalI/Rim9 protein (UN01473), are a component of the membrane sensing complex . PalF/Rim8 binding to the VPS23 protein (ESCRT-I) induces the assembly of the ESCRT-II and -III complexes as a scaffolding platform for Pal/Rim pathway activation [77, 78]. On this platform, the Pal/Rim proteolysis complex is assembled by PalA/Rim20 (UN07027), and the PalB/Rim13 protease (UN01037) . This proteolysis complex catalyzes cleavage and activation of its target, the PacC/Rim101 transcription factor (UN05739) which acts directly as transcriptional activator of “alkaline” genes . According to the expression profile of the orthologous genes identified in Fusarium sp. associated with KSBH (Fig. 5 and Additional file 3: Table S8), PALH/RIM21 receptor is upregulated in neutral/alkaline pH just like other proteins related to the endosomal sorting complexes required for transport (ESCRT I, II, and III, respectively). However, proteins from the proteolysis complex show different expression profile. On the one hand, PalA/Rim20 is also upregulated at neutral/alkaline pH, while PalB/Rim13 should be considered as ubiquitous because even when its expression profile is similar to others upregulated genes, the ortholog to PalB/Rim13 of Fusarium sp. associated with KSBH did not show significant values after DEseq normalization (Fig. 5 and Additional file 3: Table S8). This was expected since the pH of all the samples had decreased to an acidic level at the time of sample collection. PacC also shows lowest transcript levels in comparison with some other genes of the PAL/Rim system (e.g. PalH, Vsp4, Vsp24 and Snf7). This data suggest that the basal transcript levels detected for PacC could be related with some other functions performed in addition of perception and signaling to the alkaline conditions. For example, it has been reported that PacC null mutant is comprised in growth and conidiation in Aspergillus . The fact that Fusarium sp. associated with KSBH, grew at initially acidic conditions (media pH 5.0) and even after 2 weeks the pH decreased more, suggests two alternatives: first an active Pal/Rim route, and second, an alternative regulation and/or adaptation mechanism of the fungus to grow at high concentration of hydrogen ions but with some restrictions since the fungus was not able to grow at pH 3.0. Together, these data show that the PAL/RIM pathway is present in Fusarium sp. associated with KSBH.
Also, we identified genes that showed the greatest fold change value between different pHs.
The first group were the genes that most change their expression during the transition from pH 6.0 to 5.0 (UN01196, UN01849, UN03862, and UN02411). UN01196 gene encodes for a polyketide synthase protein which was homologous to FUB1, the importance of this protein and its role in pathogenicity is discussed below. UN01849 is a homolog of the vacuolar segregation 7 protein (Vac7), which in yeast functions as a Fab1 activator and together forms part of the Phosphatidylinositol 3,5-bisphosphate [PI(3,5)P2] regulation system. The PI(3,5)P2 system is critical for acidification of the vacuole since experiments with fab1Δ and vac7Δ mutants showed lower acidification of the vacuole compared to the wild-type . UN03862 is a hypothetical protein and has a multidrug and toxic compound extrusion (MATE) domain. The MATE transporters are implicated in the multidrug resistance through the exclusion of xenobiotics, and toxic metabolites from cells only have been studied in Saccharomyces cerevisiae and confers resistance to the methionine analog ethionine . UN02411 is predicted to be similar to a serine carboxypeptidase F from Fusarium solani F33, and the enzyme can release arginine and lysine from the carboxy terminus of peptides .
UN01657, UN16350, UN13020 and UN04495 were the genes that most change during the transition of pH 5.0 to 6.0 (downregulated at pH 5.0). UN01657 encodes a carnitine O-acetyltransferase, UN016350 a 2,3- dihydroxybenzoate decarboxylase, UN13020 a succinate-fumarate transporter and UN04495 an aldo-keto reductases. Both up- and down-regulated genes at pH 5.0 are included in the Additional file 3: Table S8.
The nr-unigenes UN12294, UN07751, and UN06565 were downregulated at pH 7.0 (during the transition of pH 7.0 to 6.0), while UN11374, UN12042, UN12867, UN02745, and UN10679 were highly upregulated.UN12294 encodes a 2-dehydropantoate 2-reductase, this enzyme participates in pantothenate (Vitamin B5) biosynthesis pathway. UN07751 is recognized as an amine oxidase such as a polyamine oxidase and UN06565 encodes for an unknown protein that present two domains assigned as fungal Zn binuclear cluster domain containing protein and viral late glycoprotein. With respect with the nr-unigens that are highly upregulated both UN11374 and UN10679 encode a type of permease, for example, the first one encodes an α-glucoside permease and the second one a related to neutral amino acid permease. UN12867 is an uncharacterized protein, UN02745 codes for a low-affinity potassium transporter protein, and UN12042 a coatomer subunit delta. Interestingly, all these genes encode proteins that are related to transport through the membranes and suggest their role during an efficient metabolites mobilization as well during the adjustment of pH condition of the medium. This observation is based in part on the existence of a readjustment of the initial conditions of the medium, which changed to acidic conditions after 2 weeks of growth of Fusarium sp. associated with KSBH. Finally, it is worth mentioning that the global transcriptomic analysis at each condition showed clear differences independent of the high concentration of hydrogen ions prevailing at the end of all the experiments. Interestingly, this suggests that we could use the transcriptional expression profile to know the immediate past ecologic niche of the fungus.
Pathogenicity related to environmental alkalization/acidification
In order to identify genes related with pathogenicity in Fusarium sp. associated with KSBH, pH-responsive nr-unigenes were compared by BLASTp algorithm against the pathogen-host interaction database (PHI-base) v.3.5 . In this sense various genes were recognized such as LaeA gene (Fig. 6). LaeA is a methyltransferase considered a key master regulatory protein of secondary metabolism and involved in an initial process that converts heterochromatin to euchromatin, perhaps by interfering with methylases or deacetylases associated with heterochromatin. Commonly, secondary metabolites are codified by gene clusters which are maintained as silenced heterochromatin, but, under certain stimulus, the chromatin is reprogrammed to result in increased transcription of these secondary metabolite genes . It has been proved that deletion of LaeA gene results in widespread reduction of most secondary metabolites, mainly polyketides and nonribosomal peptides including metabolites with toxic properties . In addition, in F. graminearum, it was revealed that the FgLaeA is a positive regulator of mycotoxin production such as trichothecene and zearalenone .
Likewise, homolog genes to PKS-1 and to FUM1 were differentially expressed. Both PKS-1 and FUM1, are polyketide synthases (PKS) involved in melanin  and fumonisin  biosynthesis respectively and have been considered as relevant virulence factors [66, 67]. It has been reported that disruption of pks1 gene reduced virulence in Cochliobolus heterostrophus - Zea mays interaction assays  while some fumonisins, mainly, fumonisin B1, is recognized as a virulence factor since inhibits three different plant targets: i) the sphinganine N-acyl transferases, disrupting sphingolipid metabolism (for review ), ii) the maize plasma membrane H + -ATPase  and iii) certain maize β-1,3-glucanase isoforms .
Regarding nonribosomal peptides, we identified some differentially expressed such as PES1 (UN10282), that is an NRP synthetase which confers protection against oxidative stress  and is essential for Fumigaclavine C production in Aspergillus fumigatus . Also, a homolog to NPS6 (UN01784) was also found as a remarkable NRP related with the pathogenesis since is a conserved virulence factor of plant pathogenic ascomycetes involved in siderophore-mediated iron metabolism . Also, the NRPS are responsible for beauvericin synthesis, a Fusarium mycotoxin that induces cell death in tomato protoplasts .
Moreover, we identified 13 nr-unigenes homologous to ABC-transporters. Many of these transporters are responsible also for exporting the host-derived antimicrobial compounds, so the transporters, undoubtedly, have a crucial role in the pathogenesis . Interestingly, it seems that UN00536 is homolog to the ABC transporter ABC4 from Magnaporthe grisea. The abc4 mutant did not form functional appressoria and is nonpathogenic . Also, we found two Fusarium sp. associated with KSBH nr-unigenes (UN00589, UN03631) that are homologs to ABC1 transporter of F. graminearum. It has been reported that FgABC1 and FgABC3 deletion mutants were impeded in virulence on wheat, barley, and maize .
We also recognized, 26 nr-unigenes of proteases that have differential expression among the treatments. Is already known de importance of this kind of enzymes during pathogenesis and we found that the nr-unigenes UN05400 and UN06203 have homology with BcNma, a protease of B. cinerea involved in apoptosis and viability, interestingly, the mutant has slightly reduced virulence . Also, we identified a protease that was induced at neutral pH, and it is the homolog of a vacuole-localized protease named SPM1 (UN12426) that is a protease required for the infection of the rice blast fungus Magnaporthe oryzae .
Identification of Fusarium sp. associated with KSBH nr-unigenes with potential relevance to the biosynthetic pathway of fusaric acid
Fusaric acid (FA) is a polyketide-derived secondary metabolite which is produced by multiple Fusarium species (revised by [69, 98, 99]). FA (or 5-butylpicolinic acid), acidifies plant surfaces , elicit plant defense responses , induces programmed cell death  and it has been strongly associated with plant wilt symptoms . Moreover, it has been reported that in several species such as F. verticillioides, F. fujikuroi and also different strains of F. oxysporum, genes related to FA production and fumonisins are clustered, and both clusters (FUB and FUM) are uniform in gene organization [98, 101].
Despite small variants in the genomic regions from distinct species, a total of 12 genes are present in FUB cluster. Even when FUB genes are coexpressed during specific growth conditions (e.g. GYAM medium at pH 3.0, ), specific roles from some FUB genes during fusaric acid biosynthesis  and some of the biochemical mechanism by which fusaric acid is biosynthesized from polyketide precursors remains unknown or are speculative [69, 98].
Regarding the FA biosynthetic pathway, the Fusarium sp. associated with KSBH orthologous proteins were identified by searching nr-unigenes grouped with FUB proteins that have been previously characterized in others Fusarium species (Additional file 3: Table S11). Orthologues to FUB1, FUB3, FUB7, FUB9, FUB11, and FUB12 were identified, while the remaining FUB orthologous proteins (FUB2, FUB4–6, FUB8, and FUB10) were absent in the Fusarium sp. associated with KSBH transcriptome. Orthologues to the same FUB proteins were also absent in F. solani and F. graminearum species (Additional file 3: Table S11). These data suggest that maybe in some Fusarium species, not all FUB genes are essential to produce this toxin. This is consistent with the observations about the FUB genes present in F. solani, a closely related species to Fusarium sp. associated with KSBH in which the ability to produce fusaric acid has been proved experimentally . Also, even when the FUB cluster consists of at least 12 genes, analysis of deletion mutants has been shown that only nine FUB genes are required for fusaric acid production and only the deletion of some of them can partially or completely arrest the synthesis of this metabolite . For instance, in F. fujikuroi, fub3 and fub5 deletion reduced fusaric acid production by 20% and 40%, respectively . Meanwhile, deletion of fub1 (which encodes a PKS) completely arrested the fusaric acid production in both F. verticillioides and F. fujikuroi [69, 103]. Recently, it has been reported that fub1 in F. oxysporum is positively controlled by LaeA and PacC through the modulation of chromatin accessibility at the fub1 locus .
The ability to acidify the external environment is shared by some pathogenic fungi and might contribute to generating an appropriate environment to produce some virulence factors which contribute to damaging the host tissues and promote the development of the disease. For instance, plant-necrotizing fungi such as Sclerotiorum clerotinia and Botrytis sp., secrete significant amounts of oxalic acid , while Penicillium sp. and Aspergillus sp. secrete mainly gluconic and citric acids [105, 106]. The produced acids not only acidify the tissues but can also lower the activity of reactive oxygen species produced by the host . The fusaric acid produced by some Fusarium species acidifies plant surfaces and activates the membrane H+ -ATPase, a pH-regulated process that leads to the expression of proteases and subsequent tissue invasion . In contrast, other fungi can alkalize the environment by the active secretion of ammonia. This capability has been recognized in Colletotrichum species, Fusarium oxsyporum, and Alternaria alternata . This work shows that Fusarium sp. associated with KSBH has the machinery to adapt to acidic environment. It will be necessary in the future to determine if this is useful for its virulence and to establish disease. Besides, the behavior of the fungus must be evaluated at alkaline conditions. Finally, Fusarium sp. associated with KSBH has the ability to induce transcription of genes related to virulence such as proteases, PKs, ABC transporters at both alkaline pH and acidic pH. The plasticity that is observed in this fungus to adapt to such broad pH conditions suggests that it can be a difficult pathogen to treat.
The ‘Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación (SAGARPA)’ through ‘Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria (SENASICA)’ and the ‘Dirección General de Sanidad Vegetal (DGSV)’ provided the biological material and the corresponding permission to work with the fungus in the ‘Centro Nacional de Referencia Fitosanitaria (CNRF)’ facilities under strict biosecurity conditions. The fungus Fusarium sp. associated with KSBH was isolated from the mycangia of beetles collected from internal galleries found in elm (Ulmus sp.) trees identified by Mexican phytosanitary authorities in urban areas of Tijuana municipality showing typical symptoms caused by the polyphagous shot hole borer and by their associated fungi [18, 20, 49, 108]. The beetle was taxonomically identified, and the species of the fungus was confirmed by searching for and aligning four common markers previously used in phylogenetic analysis to distinguish between different strains of Fusarium sp. associated with KSBH which were isolated from various sources and locations. For additional information see  and the corresponding section described above (Phylogenetic relationship with additional Fusarium sp. associated with KSBH strains). The four loci were: LSU (nuclear ribosomal large subunit), EF-1α (translation elongation factor 1a), RPB1 and RPB2 (subunits 1 and 2 of the DNA-directed RNA polymerase II).
Potato Dextrose Agar (PDA) medium was used for maintenance of stock cultures. For inoculum preparation, the fungus was initially grown at 28 °C on a PDA plate for 7 days. Then, a 0.7 cm2 plug from the outer zone of the colony was punched with a sterile well cutter and transferred to 100 mL Potato Dextrose Broth (PDB) medium contained in a 250 mL Erlenmeyer flask. Before autoclaving, the pH of the medium was adjusted to specific levels using either 10 N HCl or 10 N NaOH. The fungus was grown at 28 ± 2 °C on a rotary shaker at 200 rpm, and three replicate sets were grown in each case. After 2 weeks (14 days), the culture broth was filtered through a filter paper (Whatman No.1), and the collected fungus was pulverized cryogenically using a mortar and pestle.
RNA isolation, library preparation, and sequencing
Total RNA was isolated from 200 mg homogenized tissue according to the manufacturer’s protocol and using a Norgen RNA Purification Kit (Norgen Biotek Corporation). The RNA integrity was evaluated by chip-based capillary gel electrophoresis using a Bioanalyzer system (Agilent Technologies). The RNA concentration was determined by absorbance at 260 nm using a NanoDrop 2000 UV-Vis spectrophotometer (Thermo Fisher Scientific). A total amount of 500 ng of RNA was used as input material for each RNA sample preparation. Nine libraries (three independent biological replicates per analyzed conditions: pH 5.0, 6.0 and 7.0, respectively) were generated using the TruSeq RNA Sample Preparation Kit (Illumina, Cat. N°RS–122–2002) following manufacturer’s recommendations and index codes were added to identify each sample independently. The libraries were sequenced on the NextSeq 500 platform (Illumina) using the 2 × 150 bp paired-end sequencing protocol. The sequence images were transformed to bcl files with Illumina NextSeq control software (v1.4) and Real-Time Analysis software (v2), then bcl files generated were demultiplexed to fastq files with bcl2fastq (v2.0) conversion software. Files containing sequence reads and quality scores were deposited in the Short Read Archive (SRA) of the National Center for Biotechnology Information (NCBI). Accession number SRP156234.
Quality control of sequenced reads and de novo transcriptome assembly
Before assembly, raw data in fastq format were processed through a Python-based script available at GitHub repository (https://github.com/Czh3/NGSTools/blob/master/qualityControl.py). The stringency parameters employed to select High Quality (HQ) paired-end reads were: -q 20 (Minimum quality score to keep), −p 90 (Minimum percent of bases that must have [−q] quality) and -a 30 (the minimum average quality in the paired-end reads selected). SeqPrep v1.1 (https://github.com/jstjohn/SeqPrep) was used to remove adapters and merge HQ paired-end reads with overlap into a single longer reads. Adapter sequences were trimmed by searching for a ≥ 10 nt overlap between the end of each read and the adapter sequence allowing up to 2% mismatched bases of superior quality (≥ 20) and requiring a minimum of 87% matching bases in the overlapping region. Reads were discarded if they were shorter than 75 bp and orphan reads were also removed to keep pairs only. Forward and reverse reads were merged into single reads if they overlapped for at least 25 nt allowing up to 2% mismatching bases of superior quality and requiring a minimum of 90% matching bases in the overlapping region. Transcriptome de novo assembly was performed including adapters-free HQ paired-end reads and those longer sequences resulting from merging reads (R1 and R2) through their overlapping regions. Trinity assembler  with default parameters was used for this purpose in which all data set generated (libraries from pH 5.0, 6.0, and 7.0, and their corresponding replicates) were combined. Finally, the contigs resulting from the assembly process (also commonly named as unique transcripts or unigenes) were masked using the SeqClean software to eliminate sequence regions that could cause problems during the annotation process. Targets for masking include poly A/T tails and end rich with Ns (undetermined bases) and/or low complexity sequences. A unigenes set from Fusarium sp. associated with KSBH was generated considering only resulting contigs with a minimum length of 200 bp.
Identification of protein coding regions, and annotation
To generate an accurate and non-redundant data set representative from the Fusarium sp. associated with KSBH transcriptome, erroneous frame-shifts identified within the coding region were corrected based on the alignment of each assembled unigene against its homologs (at least three) identified by BLAST searches. AlignWise , a pipeline which drives several programs such as BLAST (; for similarity searches), Muscle ( to lead the alignments) and GeneWise ( to identify ORF and correct frame-shifts), was used for such purposes. Once the ORF nucleotide sequences were identified and corrected for potential frame shifts, the BLASTClust program  was used to create a non-redundant set of sequences. Assembled sequences 95% identical over an area covering at least 90% the length of another unigene, were removed to avoid redundancies.
BLAST searches were performed to identify homologs to the proteins translated from coding regions identified in the nr-unigenes for each set of reference proteins, those. BLASTp algorithm  was used for such purposes. An E-value threshold of 10− 5 was considered in all comparisons. The latest version for all reference proteins was downloaded from JGI/MycoCosm portal (http://genome.jgi.doe.gov/programs/fungi/index.jsf). The selected fungi species were mentioned in the Results section above. Additionally, in order to classify the nr-unigenes based on a functional annotation, InterPro (IPR) domains, the Gene Ontology (GO) terms, and the EuKaryotic Orthologous Groups (KOG) derived from Kyoto Encyclopedia of Genes and Genomes (KEGG), were assigned based on the information available for the homologous proteins which were identified in BLAST searches using the SBH (single-directional best hit) method.
Ortholog groups identification
The complete proteomes encoded by the genomes of S. cerevisiae, F. fujikuroi, F. oxysporum, F. verticillioides, F. graminearum, and F. solani as well as the proteins produced after translated the ORFs identified in the Fusarium sp. associated with KSBH nr-unigenes, were globally compared. The complete data set represents 111,688 protein-coding sequences (5,982 from S. cerevisiae, 14,813 from F. fujikuroi, 13,322 from F. graminearum, 20925 from F. oxysporum, 15,707 from F. solani, 15,869 from F. verticillioides, and 25,070 from Fusarium sp. associated with KSBH). Except for Fusarium sp. associated with KSBH, only proteins translated from representative gene models (without alternative isoforms) were included. An all-against-all comparison was performed using BLASTp (threshold E-value of ≤10− 10) followed by clustering with the v2.0.9 of OrthoMCL program . Default parameters, including an inflation value of 1.5, were used to define orthologous cluster structure.
Phylogenetic analyses of PKS proteins were performed in a maximum-likelihood (ML) framework using SeaView v2.4 software , which drives the Muscle  (for alignment) and PhyML  (for phylogenetic analysis) programs. The PhyML option was used under LG (Le and Gascuel) model . Equilibrium frequencies, topologies, and branch lengths were optimized, the starting tree was determined using BioNJ, and both nearest-neighbor interchange (NNI) and subtree pruning and regrafting (SPR) algorithms for tree searching were used. Branch robustness was analyzed by approximate likelihood-ratio test (aLRT) .
Identification of differentially expressed genes
The pH-responsive genes were identified by comparing global expression profiles obtained from each transcriptomic data set of Fusarium sp. associated with KSBH generated for the three analyzed conditions (pH 5.0, 6.0, and 7.0). First, the high-quality reads from each sample were mapped to the reference transcripts (unigenes) using RSEM software . RSEM is a software package that use short reads mappers such as Bowtie2  to independently assign the reads generated from each condition (distinct pH levels) to each of the unigenes resulting from the assembly process. The number of reads assigned to a specific unigene represents the abundance of mRNA produced by a particular gene in a given sample. The Expectation-maximization algorithm implemented in RSEM returns an expression profile matrix where the transcript abundance is represented by the expected reads counts (a non-normalized value) and their normalized values, TPM (transcripts per million) and FPKM (fragments per Kilobases of contigs/unigenes for per million mapped reads). Expected reads counts are then imported by tximport package  to be then used as input for DESeq, an R/Bioconductor package  which normalize the samples comparing them by pairs and using the relative log expression method to adjust the ratio values of the gene expression levels to baseline, giving the hypothesis that most genes are not differentially expressed. After comparisons of paired samples (pH 5.0/pH 6.0; pH 5.0/pH 7.0 and pH 7.0/pH 6.0), an adjusted P-value of < 0.05 was considered as threshold to identify differentially expressed unigenes.
Gene expression validation using real-time PCR
According to the manufacturer’s instructions, cDNA templates for PCR amplification were prepared from all samples by using reverse specific primers and SuperScript III reverse transcriptase (Invitrogen). Each quantitative PCR reaction contained cDNA template from ~ 1 μg total RNA, 1× SYBR Green® PCR Master Mix (Life Technologies) and 500 nM forward and reverse primers (Additional file 3: Table S12). Real-time PCR was performed on an STRATAGEN MX3005P (Agilent Technologies) detection system under the following thermal cycling conditions: 10 min at 95 °C followed by a total of 40 cycles of 30 s at 95 °C, 30 s at 55 °C and 1 min at 72 °C. For qRT-PCR, relative transcript abundance was calculated and normalized with respect to Pol II (DNA-directed RNA polymerase II core subunit RPB2; UN04406) a commonly used housekeeping gene which permits minimize variation in cDNA template levels. All calculations and analyses were performed using MxPro PCR software (Agilent Technologies) and the 2−ΔΔCt method . Amplification efficiency (85.49 to 101.56%) for the primer sets was determined by amplification of cDNA dilution series (1:5). Specificity of the RT-PCR products was followed by a melting curve analysis with continuous fluorescence data acquisition during the 55–95 °C melt.
ATP-binding cassette transporters
Ambrosia Fusarium Clade
Approximate Likelihood-Ratio Test
Basic Local Alignment Search Tool
Enzyme Commission (Number)
Fragments Per Kilobases of Contigs/Inigenes for Per Million Mapped Reads
F. solani Species Complex
Internal Transcribed Spacer (Region)
Kuroshio Shot Hole Borer
Markov CLuster (Algorithm)
Nearest Neighbour Interchange
Non-Ribosomal Peptides Synthases
Open Reading Frames
Potato Dextrose Agar
Potato Dextrose Broth
S-adenosyl methionine (binding site)
Single-Directional Best Hit (method)
Subtree pruning and regrafting (algorthims)
Short Read Archive
Transcripts Per Million
Eshel D, Miyara I, Ailing T, Dinoor A, Prusky D. pH regulates endoglucanase expression and virulence of Alternaria alternata in persimmon fruit. Mol Plant Microbe Interact : MPMI. 2002;15(8):774–9.
Drori N, Kramer-Haimovich H, Rollins J, Dinoor A, Okon Y, Pines O, Prusky D. External pH and nitrogen source affect secretion of pectate lyase by Colletotrichum gloeosporioides. Appl Environ Microbiol. 2003;69(6):3258–62.
Yakoby N, Kobiler I, Dinoor A, Prusky D. pH regulation of pectate lyase secretion modulates the attack of Colletotrichum gloeosporioides on avocado fruits. Appl Environ Microbiol. 2000;66(3):1026–30.
Caracuel Z, Roncero MI, Espeso EA, Gonzalez-Verdejo CI, Garcia-Maceira FI, Di Pietro A. The pH signalling transcription factor PacC controls virulence in the plant pathogen Fusarium oxysporum. Mol Microbiol. 2003;48(3):765–79.
Merhej J, Urban M, Dufresne M, Hammond-Kosack KE, Richard-Forget F, Barreau C. The velvet gene, FgVe1, affects fungal development and positively regulates trichothecene biosynthesis and pathogenicity in Fusarium graminearum. Mol Plant Pathol. 2012;13(4):363–74.
Miyara I, Shnaiderman C, Meng X, Vargas WA, Diaz-Minguez JM, Sherman A, Thon M, Prusky D. Role of nitrogen-metabolism genes expressed during pathogenicity of the alkalinizing Colletotrichum gloeosporioides and their differential expression in acidifying pathogens. Mol Plant Microbe Interact : MPMI. 2012;25(9):1251–63.
Brandao RL, Castro IM, Passos JB, Nicoli JR, Thevelein JM. Glucose-induced activation of the plasma membrane H(+)-ATPase in Fusarium oxysporum. J Gen Microbiol. 1992;138(Pt 8):1579–86.
Davis D. Fusaric acid in selective pathogenicity of Fusarium oxysporum. Phytopathology. 1969;59(10):1391–5.
Jiao J, Zhou B, Zhu X, Gao Z, Liang Y. Fusaric acid induction of programmed cell death modulated through nitric oxide signalling in tobacco suspension cells. Planta. 2013;238(4):727–37.
Alkan N, Espeso EA, Prusky D. Virulence regulation of phytopathogenic fungi by pH. Antioxid Redox Signal. 2013;19(9):1012–25.
Penalva MA, Arst HN Jr. Regulation of gene expression by ambient pH in filamentous fungi and yeasts. Microbiology and molecular biology reviews : MMBR. 2002;66(3):426–46 table of contents.
Tilburn J, Sarkar S, Widdick DA, Espeso EA, Orejas M, Mungroo J, Penalva MA, Arst HN Jr. The Aspergillus PacC zinc finger transcription factor mediates regulation of both acid- and alkaline-expressed genes by ambient pH. EMBO J. 1995;14(4):779–90.
Espeso EA, Penalva MA. Three binding sites for the Aspergillus nidulans PacC zinc-finger transcription factor are necessary and sufficient for regulation by ambient pH of the isopenicillin N synthase gene promoter. J Biol Chem. 1996;271(46):28825–30.
Manteau S, Abouna S, Lambert B, Legendre L. Differential regulation by ambient pH of putative virulence factor secretion by the phytopathogenic fungus Botrytis cinerea. FEMS Microbiol Ecol. 2003;43(3):359–66.
Velho AC, Mondino P, Stadnik MJ. Extracellular enzymes of Colletotrichum fructicola isolates associated to apple bitter rot and Glomerella leaf spot. Mycology. 2018;9(2):145–54.
Merhej J, Richard-Forget F, Barreau C. Regulation of trichothecene biosynthesis in Fusarium: recent advances and new insights. Appl Microbiol Biotechnol. 2011;91(3):519–28.
Flaherty JE, Pirttilä AM, Bluhm BH, Woloshuk CP. PAC1, a pH-regulatory gene from Fusarium verticillioides. Appl Environ Microbiol. 2003;69(9):5222–7.
Eskalen A, Stouthamer R, Lynch SC, Rugman-Jones PF, Twizeyimana M, Gonzalez A, Thibault T. Host range of Fusarium dieback and its ambrosia beetle (Coleoptera: Scolytinae) vector in southern California. Plant Dis. 2013;97(7):938–51.
O'Donnell K, Sink S, Libeskind-Hadas R, Hulcr J, Kasson MT, Ploetz RC, Konkol JL, Ploetz JN, Carrillo D, Campbell A, et al. Discordant phylogenies suggest repeated host shifts in the Fusarium-Euwallacea ambrosia beetle mutualism. Fungal Genet Biol : FG & B. 2015;82:277–90.
Freeman S, Sharon M, Maymon M, Mendel Z, Protasov A, Aoki T, Eskalen A, O'Donnell K. Fusarium euwallaceae sp. nov.--a symbiotic fungus of Euwallacea sp., an invasive ambrosia beetle in Israel and California. Mycologia. 2013;105(6):1595–606.
Na F, Carrillo JD, Mayorquin JS, Ndinga-Muniania C, Stajich JE, Stouthamer R, Huang Y-T, Lin Y-T, Chen CY, Eskalen A. Two novel fungal symbionts Fusarium kuroshium sp. nov. and Graphium kuroshium sp. nov. of Kuroshio shot hole borer (Euwallacea sp. nr. fornicatus) cause Fusarium dieback on woody host species in California. Plant Disease. 2017;102(6):1154–64.
García-Avila CDJ, Trujillo-Arriaga FJ, López-Buenfil JA, González-Gómez R, Carrillo D, Cruz LF, Ruiz-Galván I, Quezada-Salinas A, Acevedo-Reyes N. First report of Euwallacea nr. fornicatus (Coleoptera: Curculionidae) in Mexico. Fla Entomol. 2016;99(3):555–6.
Hossain MA, Ali MA, Moni ZR, Islam MS, Islam MR. Effect of temperature and pH on the growth and sporulation of Fusarium moniliforme: causing bakanae disease of rice. Sci Agric. 2015;11(3):151–4.
Rolland S, Bruel C, Rascle C, Girard V, Billon-Grand G, Poussereau N. pH controls both transcription and post-translational processing of the protease BcACP1 in the phytopathogenic fungus Botrytis cinerea. Microbiology (Reading, England). 2009;155(Pt 6):2097–105.
Li T, Gong L, Wang Y, Chen F, Gupta VK, Jian Q, Duan X, Jiang Y. Proteomics analysis of fusarium proliferatum under various initial pH during fumonisin production. J Proteome. 2017;164:59–72.
Evans T, Loose M. AlignWise: a tool for identifying protein-coding sequence and correcting frame-shifts. BMC bioinformatics. 2015;16:376.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Amselem J, Cuomo CA, van Kan JA, Viaud M, Benito EP, Couloux A, Coutinho PM, de Vries RP, Dyer PS, Fillinger S, et al. Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PLoS Genet. 2011;7(8):e1002230.
Staats M, van Kan JA. Genome update of Botrytis cinerea strains B05.10 and T4. Eukaryot Cell. 2012;11(11):1413–4.
Wiemann P, Sieber CM, von Bargen KW, Studt L, Niehaus EM, Espino JJ, Huss K, Michielse CB, Albermann S, Wagner D, 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):e1003475.
Cuomo CA, Guldener U, Xu JR, Trail F, Turgeon BG, Di Pietro A, Walton JD, Ma LJ, Baker SE, Rep M, et al. The Fusarium graminearum genome reveals a link between localized polymorphism and pathogen specialization. Science (New York, NY). 2007;317(5843):1400–2.
Ma LJ, van der Does HC, Borkovich KA, Coleman JJ, Daboussi MJ, Di Pietro A, Dufresne M, Freitag M, Grabherr M, Henrissat B, et al. Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature. 2010;464(7287):367–73.
Coleman JJ, Rounsley SD, Rodriguez-Carres M, Kuo A, Wasmann CC, Grimwood J, Schmutz J, Taga M, White GJ, Zhou S, et al. The genome of Nectria haematococca: contribution of supernumerary chromosomes to gene expansion. PLoS Genet. 2009;5(8):e1000618.
Galagan JE, Calvo SE, Borkovich KA, Selker EU, Read ND, Jaffe D, FitzHugh W, Ma LJ, Smirnov S, Purcell S, et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature. 2003;422(6934):859–68.
Klosterman SJ, Subbarao KV, Kang S, Veronese P, Gold SE, Thomma BP, Chen Z, Henrissat B, Lee YH, Park J, et al. Comparative genomics yields insights into niche adaptation of plant vascular wilt pathogens. PLoS Pathog. 2011;7(7):e1002137.
Shen SK, Dowd PF. Detoxification spectrum of the cigarette beetle symbiont Symbiotaphrina kochii in culture. Entomologia Experimentalis et Applicata. 1991;60(1):51–9.
Nagai K, Sakai T, Rantiatmodjo RM, Suzuki K, Gams W, Okada G. Studies on the distribution of alkalophilic and alkali-tolerant soil fungi I. Mycoscience. 1995;36(3):247–56.
Liu Y, He Q, Cheng P. Photoreception in Neurospora: a tale of two white collar proteins. Cell Mol Life Sci : CMLS. 2003;60(10):2131–8.
Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, et al. Life with 6000 genes. Science (New York, NY). 1996;274(5287):563–47 546.
Karathia H, Vilaprinyo E, Sorribas A, Alves R. Saccharomyces cerevisiae as a model organism: a comparative study. PLoS One. 2011;6(2):e16015.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.
Bairoch A. The ENZYME database in 2000. Nucleic Acids Res. 2000;28(1):304–5.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN, et al. The COG database: an updated version includes eukaryotes. BMC bioinformatics. 2003;4:41.
Finn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, Chang H-Y, Dosztányi Z, El-Gebali S, Fraser M, et al. InterPro in 2017—beyond protein family and domain annotations. Nucleic Acids Res. 2017;45(Database issue):D190–9.
Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, et al. InterPro: the integrative protein signature database. Nucleic Acids Res. 2009;37(Database issue):D211–5.
Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, Bernard T, Binns D, Bork P, Burge S, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2012;40:D306–12.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457–62.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.
Kasson MT, O’Donnell K, Rooney AP, Sink S, Ploetz RC, Ploetz JN, Konkol JL, Carrillo D, Freeman S, Mendel Z, et al. An inordinate fondness for Fusarium: phylogenetic diversity of fusaria cultivated by ambrosia beetles in the genus Euwallacea on avocado and other plant hosts. Fungal Genet Biol. 2013;56:147–57.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12:323.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26(4):493–500.
Wagner GP, Kin K, Lynch VJ. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci = Theorie in den Biowissenschaften. 2012;131(4):281–5.
Maeda T. The signaling mechanism of ambient pH sensing and adaptation in yeast and fungi. FEBS J. 2012;279(8):1407–13.
Bignell E. The molecular basis of pH sensing, signaling, and homeostasis in Fungi. Adv Appl Microbiol. 2012;79:1–18.
Moffat J, Huang D, Andrews B. Functions of Pho85 cyclin-dependent kinases in budding yeast. Prog Cell Cycle Res. 2000;4:97–106.
Carroll AS, O'Shea EK. Pho85 and signaling environmental conditions. Trends Biochem Sci. 2002;27(2):87–93.
Nishizawa M, Tanigawa M, Hayashi M, Maeda T, Yazaki Y, Saeki Y, Toh-e A. Pho85 kinase, a cyclin-dependent kinase, regulates nuclear accumulation of the Rim101 transcription factor in the stress response of Saccharomyces cerevisiae. Eukaryot Cell. 2010;9(6):943–51.
Winnenburg R, Baldwin TK, Urban M, Rawlings C, Kohler J, Hammond-Kosack KE. PHI-base: a new database for pathogen host interactions. Nucleic Acids Res. 2006;34(Database issue):D459–64.
Urban M, Pant R, Raghunath A, Irvine AG, Pedro H, Hammond-Kosack KE. The pathogen-host interactions database (PHI-base): additions and future developments. Nucleic Acids Res. 2015;43(Database issue):D645–55.
Urban M, Irvine A, Cuzick A, Hammond-Kosack K. Using the pathogen-host interactions database (PHI-base) to investigate plant pathogen genomes and genes implicated in virulence. Front Plant Sci. 2015;6(605):1–4.
Wiemann P, Brown DW, Kleigrewe K, Bok JW, Keller NP, Humpf H-U, Tudzynski B. FfVel1 and FfLae1, components of a velvet-like complex in Fusarium fujikuroi, affect differentiation, secondary metabolism and virulence. Mol Microbiol. 2010;77(4):972–94.
Bok JW, Keller NP. LaeA, a regulator of secondary metabolism in Aspergillus spp. Eukaryot Cell. 2004;3(2):527–35.
Takano Y, Kubo Y, Shimizu K, Mise K, Okuno T, Furusawa I. Structural analysis of PKS1, a polyketide synthase gene involved in melanin biosynthesis in Colletotrichum lagenarium. Mol Gen Genet. 1995;249(2):162–7.
Proctor RH, Plattner RD, Brown DW, Seo JA, Lee YW. Discontinuous distribution of fumonisin biosynthetic genes in the Gibberella fujikuroi species complex. Mycol Res. 2004;108(Pt 7):815–22.
Abbas HK, Smeda RJ, Gerwick BC, Shier WT. Fumonisin B1 from the fungus Fusarium moniliforme causes contact toxicity in plants: evidence from studies with biosynthetically labeled toxin. J Nat Toxins. 2000;9(1):85–100.
Taborda CP, da Silva MB, Nosanchuk JD, Travassos LR. Melanin as a virulence factor of Paracoccidioides brasiliensis and other dimorphic pathogenic fungi: a minireview. Mycopathologia. 2008;165(4–5):331–9.
Finkelshtein A, Shlezinger N, Bunis O, Sharon A. Botrytis cinerea BcNma is involved in apoptotic cell death but not in stress adaptation. Fungal Genet Biol : FG & B. 2011;48(6):621–30.
Brown DW, Butchko RA, Busman M, Proctor RH. Identification of gene clusters associated with fusaric acid, fusarin, and perithecial pigment production in Fusarium verticillioides. Fungal Genet Biol : FG & B. 2012;49(7):521–32.
Ibarra-Laclette E, Sánchez-Rangel D, Hernández-Domínguez E, Pérez-Torres CA, Ortiz-Castro R, Villafán E, Alonso-Sánchez A, Rodríguez-Hass B, López-Buenfil A, García-Avila C, et al. Draft genome sequence of the phytopathogenic fungus Fusarium euwallaceae, the causal agent of Fusarium dieback. Genome announcements. 2017;5(999):e00881–17.
Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44(D1):D279–85.
Pardo E, Marin S, Ramos AJ, Sanchis V. Ecophysiology of ochratoxigenic Aspergillus ochraceus and Penicillium verrucosum isolates. Predictive models for fungal spoilage prevention - a review. Food Addit Contam. 2006;23(4):398–410.
Tyagi S, Paudel R. Effect of different pH on the growth and sporulation of Fusarium oxysporum: the causal organism of wilt disease of tomato. International Journal of Basic and Applied Biology (IJBAB). 2014;2(1):103–6.
Obara K, Yamamoto H, Kihara A. Membrane protein Rim21 plays a central role in sensing ambient pH in Saccharomyces cerevisiae. J Biol Chem. 2012;287(46):38473–81.
Herranz S, Rodriguez JM, Bussink HJ, Sanchez-Ferrero JC, Arst HN Jr, Penalva MA, Vincent O. Arrestin-related proteins mediate pH signaling in fungi. Proc Natl Acad Sci U S A. 2005;102(34):12141–6.
Calcagno-Pizarelli AM, Negrete-Urtasun S, Denison SH, Rudnicka JD, Bussink HJ, Munera-Huertas T, Stanton L, Hervas-Aguilar A, Espeso EA, Tilburn J, et al. Establishment of the ambient pH signaling complex in Aspergillus nidulans: PalI assists plasma membrane localization of PalH. Eukaryot Cell. 2007;6(12):2365–75.
Obara K, Kihara A. Signaling events of the Rim101 pathway occur at the plasma membrane in a ubiquitination-dependent manner. Mol Cell Biol. 2014;34(18):3525–34.
Herrador A, Herranz S, Lara D, Vincent O. Recruitment of the ESCRT machinery to a putative seven-transmembrane-domain receptor is mediated by an arrestin-related protein. Mol Cell Biol. 2010;30(4):897–907.
Penalva MA, Lucena-Agell D, Arst HN Jr. Liaison alcaline: pals entice non-endosomal ESCRTs to the plasma membrane for pH signaling. Curr Opin Microbiol. 2014;22:49–59.
Diez E, Alvaro J, Espeso EA, Rainbow L, Suarez T, Tilburn J, Arst HN Jr, Penalva MA. Activation of the Aspergillus PacC zinc finger transcription factor requires two proteolytic steps. EMBO J. 2002;21(6):1350–9.
McCartney AJ, Zhang Y, Weisman LS. Phosphatidylinositol 3,5-bisphosphate: low abundance, high significance. BioEssays. 2014;36(1):52–64.
Omote H, Hiasa M, Matsumoto T, Otsuka M, Moriyama Y. The MATE proteins as fundamental transporters of metabolic and xenobiotic organic cations. Trends Pharmacol Sci. 2006;27(11):587–93.
Matsumura E, Sato T, Toyoda N. Isolation and characterization of a microbial Arg/Lys carboxypeptidase, carboxypeptidase F. Lett Appl Microbiol. 1995;20(3):157–9.
Keller NP, Turner G, Bennett JW. Fungal secondary metabolism - from biochemistry to genomics. Nat Rev Microbiol. 2005;3(12):937–47.
Bok JW, Balajee SA, Marr KA, Andes D, Nielsen KF, Frisvad JC, Keller NP. LaeA, a regulator of morphogenetic fungal virulence factors. Eukaryot Cell. 2005;4(9):1574–82.
Kim HK, Lee S, Jo SM, McCormick SP, Butchko RA, Proctor RH, Yun SH. Functional roles of FgLaeA in controlling secondary metabolism, sexual development, and virulence in Fusarium graminearum. PLoS One. 2013;8(7):e68441.
Yang G, Rose MS, Turgeon BG, Yoder OC. A polyketide synthase is required for fungal virulence and production of the polyketide T-toxin. Plant Cell. 1996;8(11):2139–50.
Sánchez-Rangel D, Plasencia J. The role of sphinganine analog mycotoxins on the virulence of plant pathogenic fungi. Toxin Rev. 2010;29(3–4):73–86.
Gutiérrez-Nájera N, Muñoz-Clares RA, Palacios-Bahena S, Ramírez J, Sánchez-Nieto S, Plasencia J, Gavilanes-Ruíz M. Fumonisin B1, a sphingoid toxin, is a potent inhibitor of the plasma membrane H+-ATPase. Planta. 2005;221(4):589–96.
Sanchez-Rangel D, Sanchez-Nieto S, Plasencia J. Fumonisin B1, a toxin produced by Fusarium verticillioides, modulates maize beta-1,3-glucanase activities involved in defense response. Planta. 2012;235(5):965–78.
Reeves EP, Reiber K, Neville C, Scheibner O, Kavanagh K, Doyle S. A nonribosomal peptide synthetase (Pes1) confers protection against oxidative stress in Aspergillus fumigatus. FEBS J. 2006;273(13):3038–53.
O’Hanlon KA, Gallagher L, Schrettl M, Jochl C, Kavanagh K, Larsen TO, Doyle S. Nonribosomal peptide synthetase genes pesL and pes1 are essential for Fumigaclavine C production in Aspergillus fumigatus. Appl Environ Microbiol. 2012;78(9):3166–76.
Oide S, Moeder W, Krasnoff S, Gibson D, Haas H, Yoshioka K, Turgeon BG. NPS6, encoding a nonribosomal peptide synthetase involved in siderophore-mediated iron metabolism, is a conserved virulence determinant of plant pathogenic ascomycetes. Plant Cell. 2006;18(10):2836–53.
Paciolla C, Dipierro N, Mulè G, Logrieco A, Dipierro S. The mycotoxins beauvericin and T-2 induce cell death and alteration to the ascorbate metabolism in tomato protoplasts. Physiol Mol Plant Pathol. 2004;65(1):49–56.
Gupta A, Chattoo BB. Functional analysis of a novel ABC transporter ABC4 from Magnaporthe grisea. FEMS Microbiol Lett. 2008;278(1):22–8.
Abou Ammar G, Tryono R, Doll K, Karlovsky P, Deising HB, Wirsel SG. Identification of ABC transporter genes of Fusarium graminearum with roles in azole tolerance and/or virulence. PLoS One. 2013;8(11):e79042.
Saitoh H, Fujisawa S, Ito A, Mitsuoka C, Berberich T, Tosa Y, Asakura M, Takano Y, Terauchi R. SPM1 encoding a vacuole-localized protease is required for infection-related autophagy of the rice blast fungus Magnaporthe oryzae. FEMS Microbiol Lett. 2009;300(1):115–21.
Brown DW, Lee SH, Kim LH, Ryu JG, Lee S, Seo Y, Kim YH, Busman M, Yun SH, Proctor RH, et al. Identification of a 12-gene fusaric acid biosynthetic gene cluster in Fusarium species through comparative and functional genomics. Molecular plant-microbe interactions : MPMI. 2015;28(3):319–32.
Desjardins AE. Fusarium mycotoxins. St. Paul, MN: APS Press, American Phytopathological Society; 2006.
Bouizgarne B, El-Maarouf-Bouteau H, Madiona K, Biligui B, Monestiez M, Pennarun AM, Amiar Z, Rona JP, Ouhdouch Y, El Hadrami I, et al. A putative role for fusaric acid in biocontrol of the parasitic angiosperm Orobanche ramosa. Molecular plant-microbe interactions : MPMI. 2006;19(5):550–6.
Proctor RH, Van Hove F, Susca A, Stea G, Busman M, van der Lee T, Waalwijk C, Moretti A, Ward TJ. Birth, death and horizontal transfer of the fumonisin biosynthetic gene cluster during the evolutionary diversification of Fusarium. Mol Microbiol. 2013;90(2):290–306.
Bacon CW, Porter JK, Norred WP, Leslie JF. Production of fusaric acid by Fusarium species. Appl Environ Microbiol. 1996;62(11):4039–43.
Niehaus EM, von Bargen KW, Espino JJ, Pfannmuller A, Humpf HU, Tudzynski B. Characterization of the fusaric acid gene cluster in Fusarium fujikuroi. Appl Microbiol Biotechnol. 2014;98(4):1749–62.
Lopez-Diaz C, Rahjoo V, Sulyok M, Ghionna V, Martin-Vicente A, Capilla J, Di Pietro A, Lopez-Berges MS. Fusaric acid contributes to virulence of Fusarium oxysporum on plant and mammalian hosts. Mol Plant Pathol. 2017.
Ruijter GJ, Visser J. Characterization of Aspergillus niger phosphoglucose isomerase. Use for quantitative determination of erythrose 4-phosphate. Biochimie. 1999;81(3):267–72.
Prusky D, Yakoby N. Pathogenic fungi: leading or led by ambient pH? Mol Plant Pathol. 2003;4(6):509–16.
Chen C, Dickman MB. Proline suppresses apoptosis in the fungal pathogen Colletotrichum trifolii. Proc Natl Acad Sci U S A. 2005;102(9):3459–64.
O’Donnell K, Libeskind-Hadas R, Hulcr J, Bateman C, Kasson MT, Ploetz RC, Konkol JL, Ploetz JN, Carrillo D, Campbell A, et al. Invasive Asian Fusarium – Euwallacea ambrosia beetle mutualists pose a serious threat to forests, urban landscapes and the avocado industry. Phytoparasitica. 2016;44(4):435–42.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Birney E, Clamp M, Durbin R. GeneWise and Genomewise. Genome Res. 2004;14(5):988–95.
Alva V, Nam SZ, Soding J, Lupas AN. The MPI bioinformatics toolkit as an integrative platform for advanced protein sequence and structure analysis. Nucleic Acids Res. 2016;44(W1):W410–5.
Li L, Stoeckert CJ Jr, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13(9):2178–89.
Gouy M, Guindon S, Gascuel O. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27(2):221–4.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52(5):696–704.
Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25(7):1307–20.
Anisimova M, Gascuel O. Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative. Syst Biol. 2006;55(4):539–52.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
Soneson C, Love M, Robinson M. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences [version 1; referees: 2 approved]. F1000Research. 2015;4(1521):1–18.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods (San Diego, Calif). 2001;25(4):402–8.
Ost KS, O’Meara TR, Huda N, Esher SK, Alspaugh JA. The Cryptococcus neoformans alkaline response pathway: identification of a novel rim pathway activator. PLoS Genet. 2015;11(4):e1005159.
Bardou P, Mariette J, Escudie F, Djemiel C, Klopp C: jvenn: an interactive Venn diagram viewer. BMC bioinformatics 2014, 15:293.
Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6(7):e21800.
We thank Martín Aluja Schuneman Hofer for the managing of funding with Mexican phytosanitary authorities. Also, we thank Antonio Cárcamo-Rodríguez and Karen Preuss-Ángeles from the Mycology Laboratory at CNRF. We are also thankful to our other colleagues from CNRF: Clemente de Jesús García Avila, Andrés Quezada Salinas and José Abrahán Ramírez Pool for their valuable suggestions which improved the quality of this manuscript. We are grateful for the group of colleagues from INECOL (most of them included in this work) who allowed us early access to the genome of Fusarium sp. associated with KSBH. We also thank the staff of mycology laboratory at CNRF for their technical assistance. We thank to Yifan Lii for editing the paper. Special thanks to Josué Pérez-Lira for his excellent technical assistance.
This work was supported by SAGARPA-SENASICA through an agreement between INECOL and Mexican phytosanitary authorities (SENASICA-INECOL, 2016). This work also was supported by FORDECYT-CONACYT project 292399.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Summary of sequencing data generated from Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus. Table S2. Annotation of Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus no-redundant unigenes. (XLSX 25088 kb)
The nucleotide sequence of the non-redundant unigenes from Fusarium sp. (ZIP 20501 kb)
Table S3. Functional categorization of the Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus non-redundant unigenes Table S4. Metabolomic pathways network. Table S5. OrthoMCL-defined protein families shared between some Ascomycota species. Table S6. Expression profile matrix. Table S7. pH-responsive genes identified in Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus. Table S8. pH-responsive genes from Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus involved in specific matabolic pathways. Table S9. pH-responsive genes from Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus with significant similarity (e-value 10-05) to those genes/proteins previosly reported in PHI-base (pathogen-host interaction database). Tabla S10. Non-redundnat unigenes from Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus involved in the tricarboxylic acid cycle (oxaloacetate biosynthetic pathway). Table S11. Genes involved in fusaric acid (FA) biosynthesis. Table S12. Primers used in RT-qPCR assays. (XLSX 14848 kb)
Figure S1. Metabolic network (and their corresponding pathways) represented in the nr-unigene set from Fusarium sp. associated with KSBH Euwallacea sp. near fornicatus. Figure S2. Comparison of Fusarium species proteomes. Figure S3. Maximum likelihood phylogenetic tree based on the concatenated sequences of four-locus analyzed (RB1, RB2, EF-1a, LSU). Figure S4. Results of gene expression validated by quantitative real-time PCR analysis. Figure S5. GO treemap for the differentially expressed genes and related to pathogen-host interaction. Figure S6. Oxaloacetate biosynthesis pathway reconstructed based on the de novo assembly and annotation of the Fusarium sp. associated with KSBH Euwallacea sp. near fornicatus transcriptome. Figure S7. Proposed biosynthetic pathway for fusaric acid biosynthesis in species of Fusarium sp. associated with KSBH Euwallacea sp. near fornicatus and the expression profile of genes involved in the pathway [122, 123]. (PPTX 5498 kb)
Concatenated multi-locus sequence alignment used in Maximum Likelihood phylogenetic tree. (FASTA 184 kb)
Alignment of the full-length sequences of the polyketide synthases (PKSs) from distinct Fusarium species. (FASTA 1114 kb)
About this article
Cite this article
Sánchez-Rangel, D., Hernández-Domínguez, EE., Pérez-Torres, CA. et al. Environmental pH modulates transcriptomic responses in the fungus Fusarium sp. associated with KSHB Euwallacea sp. near fornicatus. BMC Genomics 19, 721 (2018). https://doi.org/10.1186/s12864-018-5083-1