Next Article in Journal
Biocontrol of Diseases Caused by Phytophthora capsici and P. parasitica in Pepper Plants
Next Article in Special Issue
Genomic and Transcriptomic Survey Provides Insights into Molecular Basis of Pathogenicity of the Sunflower Pathogen Phoma macdonaldii
Previous Article in Journal
Development of a Multiplex qPCR Assay for Fast Detection and Differentiation of Paracoccidioidomycosis Agents
Previous Article in Special Issue
Assembling Quality Genomes of Flax Fungal Pathogens from Oxford Nanopore Technologies Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conservation and Expansion of Transcriptional Factor Repertoire in the Fusarium oxysporum Species Complex

1
Department of Biochemistry and Molecular Biology, University of Massachusetts Amherst, Amherst, MA 01003, USA
2
Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, University of California Berkeley, Berkeley, CA 94720, USA
3
Department of Plant and Microbial Biology, University of California Berkeley, Berkeley, CA 94598, USA
*
Author to whom correspondence should be addressed.
Current address: Peking University Institute of Advanced Agricultural Sciences, Weifang 261000, China.
J. Fungi 2023, 9(3), 359; https://doi.org/10.3390/jof9030359
Submission received: 7 February 2023 / Revised: 11 March 2023 / Accepted: 13 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Genomics of Fungal Plant Pathogens, 2nd Edition)

Abstract

:
The Fusarium oxysporum species complex (FOSC) includes both plant and human pathogens that cause devastating plant vascular wilt diseases and threaten public health. Each F. oxysporum genome comprises core chromosomes (CCs) for housekeeping functions and accessory chromosomes (ACs) that contribute to host-specific adaptation. This study inspects global transcription factor profiles (TFomes) and their potential roles in coordinating CC and AC functions to accomplish host-specific interactions. Remarkably, we found a clear positive correlation between the sizes of TFomes and the proteomes of an organism. With the acquisition of ACs, the FOSC TFomes were larger than the other fungal genomes included in this study. Among a total of 48 classified TF families, 14 families involved in transcription/translation regulations and cell cycle controls were highly conserved. Among the 30 FOSC expanded families, Zn2-C6 and Znf_C2H2 were most significantly expanded to 671 and 167 genes per family including well-characterized homologs of Ftf1 (Zn2-C6) and PacC (Znf_C2H2) that are involved in host-specific interactions. Manual curation of characterized TFs increased the TFome repertoires by 3% including a disordered protein Ren1. RNA-Seq revealed a steady pattern of expression for conserved TF families and specific activation for AC TFs. Functional characterization of these TFs could enhance our understanding of transcriptional regulation involved in FOSC cross-kingdom interactions, disentangle species-specific adaptation, and identify targets to combat diverse diseases caused by this group of fungal pathogens.

1. Introduction

The fungal species complex of Fusarium oxysporum (FOSC) has been used as a model to study cross-kingdom fungal pathogenesis. Members within FOSC can cause devastating fusarium wilt diseases among economically important crops [1,2,3,4,5,6,7,8,9,10,11,12] and is listed among the top five most important plant pathogens [8]. With strong host specificity, plant pathogenic F. oxysporum strains are further grouped as formae speciales [13]. For instance, tomato pathogens are named F. oxysporum f. sp. lycopersici; cotton pathogens are F. oxysporum f. sp. vasinfectum [11], and banana pathogens are F. oxysporum f. sp. cubense [10]. Recently, members within FOSC have also been reported to be responsible for fusariosis, the top emerging opportunistic mycosis [1,7,12], and fusarium keratitis, one of the major causes of cornea infections in the developing world and the leading cause of blindness among fungal keratitis patients [14,15].
Comparative genomics studies on this cross-kingdom pathogen revealed that the FOSC genomes, both human and plant pathogens, are compartmentalized into two components: the core chromosomes (CCs) and accessory chromosomes (ACs). While CCs are conserved and vertically inherited to execute essential housekeeping functions, horizontally transmitted ACs are lineage- or strain-specific and are related to fungal adaptation and pathogenicity [1,7,12].
To coexist and function within the same genome, ACs and CCs coordinate their gene expression. One intriguing cross-regulation example, reported in the reference genome of F. oxysporum f. sp. lycopersici Fol4287, includes transcription factors Sge1 (SIX Gene Expression 1), Ftfs, and virulent factors SIX (Secreted in Xylem) proteins. Sge1 is a highly conserved CC-encoding TF. By name definition, Sge1 regulates the expression of SIX proteins [16,17]. The Fol4287 genome encodes an AC-encoding Ftf1 protein and one CC-encoding Ftf2 (Ftf1 CC homolog) [17]. Constitutive expression of either Ftf1 or Ftf2 induced the expression of effector genes [17]. Furthermore, it has been documented that DNA binding sites of Sge1 and Ftf1 are enriched among the cis-regulatory elements of in planta transcriptionally upregulated genes [17]. Another example of CC and AC cross-talking is the alkaline pH-responsive transcription factor PacC/Rim1p reported in F. oxysporum clinical strains [18]. In addition to the full-length PacC ortholog (PacC_O), located on a CC, the clinical isolate NRRL32931 genome encodes three truncated PacC homologs, named PacC_a, PacC_b, and PacC_c in ACs [18].
To thoroughly understand the coordination of the crosstalk between genome compartments and their contribution to the cross-kingdom fungal pathogenesis, this study compared the repertoire of TFs (i.e., TFome) among 15 F. oxysporum and 15 other ascomycete fungal genomes. Remarkably, we discovered a strong positive correlation (y = 0.07264x − 190.9, r2 = 0.9361) between the number of genes (x) and TFome size (y) of an organism. Primarily due to the acquisition of ACs, we observed increased TFome sizes among the FOSC genomes. All TFs were organized into 48 families based on the InterPro classification of proteins. Fourteen families involved in transcription/translation regulations and cell cycle controls were highly conserved. Thirty families, accounting for ¾ of all families, were expanded to various degrees among the FOSC genomes. Unique TF expansions driven by ACs include members of the Zn2-C6 fungal-type (Zn2-C6) and Zinc Finger C2H2 (Znf_C2H2) families. This comparative study highlighted conserved regulatory mechanisms. The signature of conservation established the foundation to study the various impacts of additional AC TFs on existing regulatory pathways. In combination with the existing expression data, this study provides insights into the fine-tuning of environmental adaptation performed by this group of diverse organisms in order to engage in cross-kingdom interactions with different hosts.

2. Materials and Methods

2.1. Generation of Fungal TFomes

The annotation pipeline is briefly summarized in Figure S1A,B. The fungal proteomes of 30 strains were downloaded from the JGI MycoCosm portal [19]. Protein annotation was performed using InterProScan/5.38–76.0 (https://www.ebi.ac.uk/interpro/search/sequence/, accessed on 7 February 2023) [20]. Annotations of proteins that putatively serve as TFs were filtered out using a table containing InterPro terms related to transcriptional regulatory functions summarized in the literature [21,22], with further addition by manual curation (Table S1). Orthologous analysis was conducted with OrthoFinder 2.5.4 (https://github.com/davidemms/OrthoFinder, accessed on 7 February 2023) [23] to probe the orthologs of functionally validated TFs (Tables S3 and S4) in Fusarium.

2.2. RNA-Seq Analysis

The RNA-Seq datasets were previously described [24,25] and deposited by those authors to the NCBI Short Read Archive with accession number GSE87352 and to the ArrayExpress database at EMBL-EBI (www.ebi.ac.uk/arrayexpress, accessed on 7 February 2023) under accession number E-MTAB-10597, respectively. For data reprocessing, reads were mapped to reference genomes of Arabidopsis [annotation version Araport11 [26]], Fo5176 [27], Fo47 [28], and Fol4287 [29] using HISAT2 version 2.0.5 [30]. Mapped reads were used to quantify the transcriptome by StringTie version 1.3.4 [31], at which step TPM (transcript per million) normalization was applied. Normalized read counts were first averaged per condition, transformed by log2 (normalized read count + 1) and Z-scaled. This was then visualized in pheatmap (version 1.0.12).

2.3. Genome Partition

The genome partition results for chromosome-level assemblies were retrieved from previous reports for Fol4287 [29], FoII5 [32], Fo5176 [27], and Fo47 [28]. Fo47 has a clear genome partition with 11 core chromosomes and one accessory chromosome, therefore serving as the reference for the genome partition of other F. oxysporum genomes. MUMmer/3.22 was applied to align scaffolds of genome assemblies against 11 core chromosomes of the reference genome Fo47 using default parameters. The scaffolds aligned to the core chromosomes of Fo47 with a coverage larger than 5% were annotated as core scaffolds. The rest of the scaffolds were partitioned as accessory scaffolds. Genes residing on the core and accessory scaffolds were annotated as the core and accessory genes, respectively.

2.4. Phylogenetics Analysis

Protein sequences were aligned via MAFFT/7.313 [33]. The iqtree/1.6.3 [34,35] was run on the sequence alignment to generate the phylogeny (by maximum likelihood method and bootstrapped using 1000 replicates) [36]. Visualization was conducted via the Interactive Tree of Life [37] to produce the phylogram. OrthoFinder 2.5.4 [23] was used for orthogroup determination. To build a species phylogram, 500 randomly selected conserved proteins (single-copy orthologs) were aligned first. The alignment was then concatenated, and the phylogeny was determined and visualized using the above methods.

2.5. Expansion Index Calculation

To understand genome regulation among FOSC, we developed two expansion index scores. The first uses two yeast lineages as the baseline (EIy):
E I y = A v e r a g e   n u m b e r   o f   T F s   i n   F O S C + 1 A v e r a g e   n u m b e r   o f   T F s   i n   y e a s t s + 1
In the second index score (EIf), we directly compared F. oxysporum with its Fusarium relatives to calculate the expansion index as follows:
E I f = A v e r a g e   n u m b e r   o f   T F s   i n   F O S C + 1 A v e r a g e   n u m b e r   o f   T F s   i n   F O S C   s i s t e r   s p e c i e s + 1

3. Results

3.1. FOSC TFome Expansion Resulted from the Acquisition of ACs

We compared 30 ascomycete fungal genomes (Figure 1 and Table 1) including 15 strains within the FOSC, nine sister species close to F. oxysporum, two yeast genomes (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and four other filamentous fungal species (Neurospora crassa, Aspergillus nidulans, Aspergillus acristatulus, and Magnaporthe oryzae). To maintain consistency, the protein sequences for all of these genomes were retrieved from the MycoCosm portal [19].
For a comprehensive TFome annotation, we used InterProScan (IPR) terms associated with fungal transcriptional regulation [21,22] and curated a mapping with updated IPR classification (interproscan version: 5.38–76.0) [52]. In addition, we searched the IPR classification of the protein families and obtained all other terms related to the transcriptional regulation activity. This resulted in 234 TF-related IPR terms (Table S1). Since most terms were initially defined in the mammalian systems, fungal genomes included in this study were associated with 71 out of the 234 TF-related IPR terms (Table S1, Materials and Methods, and Figure S1A,B for the annotation pipeline). After removing 13 terms for redundancy (two terms describing the identical domain) and 10 terms for minimal presentation (<4 among the 30 genomes), this comparative TFome study focused on 48 IPR terms in 27,967 TFs (Tables S1 and S2). Notably, a quarter of these terms were not reported to be affiliated with fungal transcriptional regulation by either Park et al. (2008) [21] or Shelest (2017) [22] (Table S1).
Comparing the total number of protein-coding genes (x) and the total number of TFs (y) within the same genome, we observed a strong positive correlation (y = 0.07264x − 190.9, r2 = 0.9361) (Figure 2A). FOSC TFomes were larger than other genomes included in this study, with an average of 1144 TFs per genome (Figure 2A, Table 1). After partitioning each FOSC genome into core and accessory regions (see Section 2.3 for details), we observed a positive correlation between the number of TFs encoded in the accessory chromosomal region of each strain (defined as accessory TFs hereafter) with the size of the accessory genomes (Mb) (y = 17.239x + 3.553) (Figure S2). This suggests that accessory chromosomes contribute directly to the expanded TFome.
To understand the genome regulation among the FOSC, we developed an expansion index score using two yeast lineages as the baseline (EIy):
E I y = A v e r a g e   n u m b e r   o f   T F s   i n   F O S C + 1 A v e r a g e   n u m b e r   o f   T F s   i n   y e a s t s + 1
Based on this index value, we classified TF families into three major groups (Table 2 and Table S1). Group 1 contained 14 TF families with an expansion score of 1, indicating high conservation. Group 2 included four families with an index score below 1, reflecting some level of gene family contraction. Group 3 contained 30 families with an expansion index greater than 1, indicating gene expansion.

3.2. Conserved TF Families That Are Primarily Associated with General/Global Transcription Factors

Fourteen TF families accounted for 30% of our annotated TF families. Most of these fourteen families had a single ortholog in all genomes included in this study (Figure 2B; Table 2 and Table S1), suggesting their functional conservation across the Ascomycota. These 30% conserved TF families accounted for less than 2% of the total TFomes. Annotation based on S. cerevisiae and other model organisms suggested their involvement in transcription/translation regulation and cell cycle control.

3.2.1. Transcription/Translation Regulation

Nine TF families were annotated to be related to transcription and translational regulation including TATA box-binding protein (TBP), TBP-associated factors (TAFs), RNA polymerase II elongation regulator Vps25, and CCAAT-binding factors (CBFs) related to ribosomal biogenesis.
One of the most conserved TF families, transcription initiation TBP binds directly to the TATA box to define the transcription start and initiate transcription facilitated by all three RNA polymerases. In fact, the function of TBP is so conserved that the yeast homolog can complement TBP mutations in humans [53,54]. Seven conserved TF families are classified as transcription positive/negative regulators, and transcription elongation factors. TAF12 and TAFII28 are parts of the transcription factor TFIID complex. Interacting with TBP, TAFs form the TFIID complex and positively participate in the assembly of the transcription preinitiation complex [55]. Similarly, TFIIH works synergistically with TFIID to promote the transcription [56]. In contrast, negative cofactor 2 (Ncb2) inhibits the preinitiation complex assembly [57]. Other factors include the CNOT1, a global regulator involved in transcription initiation and RNA degradation [58], and Vps72/YL1, which contributes to transcriptional regulation through chromatin remodeling, as reported in yeast [52,59]. Vps25 is a subunit of the ESCRT-II complex, which binds to the RNA polymerase II elongation factor to exert transcriptional control in mammalian systems [60]. One TF family is suggested to be involved in translational regulation. CCAAT box is a common cis-acting element found in the promoter and enhancer regions of genes in the eukaryotes [61,62]. CBFs are necessary for 60S ribosomal subunit biogenesis and are therefore involved in translational control [63,64,65]. This family including Noc3, Noc4, and Mak21 in S. cerevisiae had three members in each genome, and a clear single-copy orthologous relationship could be observed for each member (Figure S3A).

3.2.2. Cell Cycle Control

Five TF families were related to cell cycle control including cell cycle progression, DNA repair, and machinery/cell integrity maintenance.
One conserved TF family, APSES-type HTH, was reported to be involved in cell cycle control and crucial to development [66]. Every genome included in this study encoded four copies of the APSES-type HTH gene (Figure S3B) that formed four clades of single-copy orthologs in all genomes except yeasts. Genes in Clade 1 included StuA homologs. As a target of the cyclic AMP (cAMP)-dependent protein kinase A (PKA) signal transduction pathway, StuA was reported to be involved in dimorphic switch [67,68], fungal spore development, and the production of secondary metabolites [69]. Genes in Clade 2 and Clade 3 included S. cerevisiae Swi4 and Swi6, which were reported to form a protein complex regulating cell cycle progression from G1 to the S phase [70] as well as meiosis [71]. Genes in Clade 4 included homologs of S. pombe Bqt4, anchoring telomeres to the nuclear envelope [72].
The conserved TF family, DTT, represented by the S. cerevisiae homolog Itc1, is recognized as a subunit of the ATP-dependent Isw2p-Itc1p chromatin remodeling complex and is required for the repression of early meiotic gene expression during mitotic growth [73].
The other conserved TF family, RFX, was reported to be involved in DNA repair. Each strain encoded two orthologous copies, except for F. venenatum encoding two copies within the RFX1 clade (Figure S3C). Being a major transcriptional repressor of DNA-damage-regulated genes in S. cerevisiae, Rfx1 functions in DNA damage repair and replication checkpoint pathways [74]. In F. graminearum, Rfx1 was reported to be essential in maintaining the genome integrity [75]. The other copy, Rsc9 in S. cerevisiae, was reported to be a member of the chromatin structure-remodeling complex RSC, which is involved in transcription regulation and nucleosome positioning [76,77].
The conserved TF family NFYA was reported to bind to the CCAAT box. All strains maintained a single copy of this family. Its yeast homolog, Hap2, has been reported to induce the expression of mitochondrial electron transport genes [78] and its F. verticillioides homolog NFYA Hap2 was reported to be essential for fungal growth and the virulence on maize stalks [79].
The MADS MEF2-like TF family including S. cerevisiae Rlm1 was reported to be a component of the protein kinase C-mediated MAP kinase pathway involved in maintaining cell integrity [80]. Having a paralog from the whole genome duplication in S. cerevisiae, Rlm1 was detected as a single copy gene in all filamentous fungi included in this study. Its member in F. verticilioides, Mef2 has been reported to play a vital role in sexual development [81].

3.3. Gene Family Contractions in FOSC Partially Caused by Whole Genome Duplication in Yeast

We detected an expansion score of less than 1 for four TF families, MATalpha_HMGbox, NOT4, MADS_SRF-like, and HSF (heat shock factor), reflecting some level of gene family contraction among members of FOSC compared to the two yeast genomes (Figure S4).
TF family MATalpha_HMGbox including S. cerevisiae mating type protein alpha 1 has been reported to be a transcription activator that activates mating-type alpha-specific genes [82]. Reflecting the potential heterothallic mating strategy, all F. oxysporum Mat1-1 type strains contained this TF, but not the Mat1-2 strains, even though sexual reproduction has not been observed in FOSC [83].
TF family NOT4 was reported to be a component of the multifunctional CCR4-NOT complex, a global transcriptional repressor of the RNA polymerase II transcription [84]. Most genomes included in this study encoded one copy of this TF family, but some filamentous fungal genomes including A. nidulans, F. redolens, F. oxysporum strains NRRL26365, MRL8666, and PHW726 lost it. The functional implication of this loss remains to be discovered.
The contractions of the other two TF families, MADS SRF-like and HSF, were primarily caused by the whole genome duplication in yeast. In contrast to the contraction at the global scale, both TF families expanded among some FOSC strains when compared to other filamentous fungi (Figure S4).
As reported in M. oryzae, MADS SRF-like TF, essential for the transcriptional regulation of growth-factor-inducible genes [85], is important for microconidium production and virulence in host plants [86]. Due to the event of whole genome duplication, the S. cerevisiae genome contained two copies of this TF family, while all filamentous genomes encoded a single copy. However, we detected an average of 2.73 copies among the phytopathogenic FOSC strains. There were six copies in the genome of Fo5176, a pathogen of Brassicaceae plants including A. thaliana (Table S1).
TF family HSF has been reported to activate the production of heat shock proteins that prevent or mitigate protein misfolding under abiotic/biotic stresses [87]. The S. cerevisiae genome contained five copies of the HSF TF family, while all non-FOSC filamentous fungi had three copies. Members of FOSC exhibited some level of expansion to four or five copies (Fo47:4, Fol4287: 5, II5: 4, HDV274: 4, and Fo5176: 4), with 1–2 copies encoded in ACs. Phylogenetically, all HSF TFs were clustered into three major clades, named as Skn7, Sfl1, and Hsf1 (Figure 3A,B). All AC-encoding HSFs were phylogenetically close to Hsf1 (Figure 3A). Based on the study in M. oryzae, the family Sfl1 is essential for vegetative growth, conidiation, sexual reproduction, and pathogenesis [88]. Based on a study in F. graminearum, the family Skn7 is involved in regulating the oxidative stress response and is essential for pathogenicity [89]. Our expression data generated during the plant colonization [24] supported the involvement of all three core genes during plant colonization (Figure 3C). However, the Hsf1 accessory copies of these two strains were distinct, as the Fo47 AC-encoding Hsf1 was upregulated and the Fo5176 AC-encoding Hsf1 was downregulated, post inoculation (Figure 3C), suggesting that their distinct regulatory function is involved in these two distinct interactions.

3.4. Significant FOSC TFome Expansion Driven by a Few Exceedingly Expanded TF Families

3.4.1. Gain-of-Function among Filamentous Ascomycete Fungi

Three TF families, CP2 (EIy = 2.73), HTH_AraC (EIy = 2), and HTH_Psq (EIy = 3.53), were absent in both yeast genomes, suggesting a gain of function among filamentous ascomycete fungi (Table S1). TF family CP2 was studied in animal and fungal kingdoms with a function related to differentiation and development [90]. Both HTH_AraC and HTH_Psq are part of the helix-turn-helix (HTH) superfamily. HTH_AraC was first reported in bacteria as a positive regulator regulating the arabinose operon [91,92,93]. HTH_Psq, as part of the eukaryotic Pipsqueak protein family, has been reported in vertebrates, insects, nematodes, and fungi to regulate processes involved in cell death [94]. Most FOSC genomes encoded a single copy of HTH_AraC, while the number of HTH_Psq-containing proteins ranged from 0 to 9 in the FOSC and 0 to 3 in other Fusarium genomes. Since the HTH_Psq domain also exists in transposases [94], and ACs in FOSC are transposon-rich, it remains to be studied whether proteins containing the Psq domain are bona fide TFs.

3.4.2. Seven Exceedingly Expanded TF Families

Among the families containing minimally one yeast ortholog, seven TF families had expansion indices greater than 2 (Table 2 and Figure 2B) including Zn2-C6 (EIy = 15.09), bZIP (EIy = 5.80), Znf_C2H2 (EIy = 4.15), Homeobox (EIy = 2.28), PAI2 (EIy = 3.42), NDT80 (EIy = 3.47), and bHLH (EIy = 3.48). Based on the number of increments, the most significantly expanded TF families were Zn2-C6 (44 in yeasts versus 671 in FOSC) and Znf_C2H2 (40 in yeasts versus 167 in FOSC) (Figure 2C and Table S1). Because of their large expansion, these seven families accounted for more than 75% of the total TFome. All seven families exhibited a gradual expansion, following the pattern yeasts < non-Fusarium filamentous fungi < non-FOSC Fusarium < FOSC (Table S1). Annotating large TF families could be challenging. Here, we described some examples based on the literature.
Zn2-C6, a fungal TF family [95], was detected as the most significant expanded TF family, reaching over 600 members among the FOSC genomes and accounting for more than half of the total TFome. Able to form a homodimer, this group of TFs are also able to bind to the specific palindromic DNA sequence through direct contact with the major groove of the double-stranded DNA molecules [95]. The versatility of this group of TFs can be achieved by domain shuffling and by changing the nucleotide binding specificity. In addition to the well-documented Ftf1 [17,96,97,98,99], five additional TFs within this family were reported in F. oxysporum including Ctf1 [100], Ctf2 [100], Fow2 [101], XlnR [102] and Ebr1 [103]. Their functions were reported to be involved in the development, metabolism, stress response, and pathogenicity.
Znf_C2H2 was reported to be the most common DNA-binding motif found in the eukaryotic transcription factors [104]. Five reported F. oxysprum TFs are Czf1 [105], Con7-1 [106], PacC [18,107], ZafA [108], and St12 [109,110]. Classified in the Znf_C2H2 family, PacC has been linked to the fungal virulence in both plant and human hosts [18,107].
The other five families were bZIP, Homeobox, PAI2, Ndt80, and bHLH. The bZIP domain contains a region for sequence-specific DNA binding followed by a leucine zipper region required for dimerization [111]. Three F. oxysporum bZIP TFs were reported including Atf1 [112], Hapx [113], and MeaB [114], all of which are important for fungal pathogenicity. Homeobox is a DNA binding motif with a helix-turn-helix structure. In S. pombe, a homeobox-domain containing protein Phx1 was reported to be a transcriptional coactivator involved in yeast fission. In M. oryzae, the homeobox-domain containing protein Hox played roles in conidiation and appressorium development [115]. The TF family PAI2 is involved in the negative regulation of protease synthesis and sporulation of the Bacillus subtilis [116]. The TF family Ndt80 is essential for completing meiosis in S. cerevisiae [117,118] and Ustilago maydis [119] by promoting the expression of sporulation genes for the fulfillment of meiotic chromosome segregation [120]. The TF family bHLH forms a superfamily of transcriptional regulators found in almost all eukaryotes and is involved in diverse developmental processes [121]. In F. graminearum, a bHLH-domain containing protein Gra2 was reported to regulate the biosynthesis of phytotoxin gramillin [122], while a bHLH-domain containing protein SreA in Penicillium digitatum is required for anti-fungal resistance and full virulence in citrus fruits [123].

3.4.3. Other Families

The other 20 TF families (expanded but with EIy ≤ 2) accounted for 20% of the TFome, with an average 9.6 copies in each genome examined (Table S1).
Four TF families were functionally linked to chromatin remodeling including Bromodomain (EIy = 1.52), CBFA_NFYB (EIy = 1.35), Znf_RING-CH (EIy = 1.11), and ARID (EIy = 1.25). The TF family Bromodomain contained Spt7. As a crucial part of the SAGA complex in yeast, Spt7 has been reported to recognize the acetylated lysines of histones and eventually lead to chromatin unwinding [124]. The CBFA_NFYB domain was found in proteins (e.g., S. cerevisiae Dls1) that regulate RNA polymerase II transcription by controlling the chromatin accessibility (e.g., telomeric silencing) [125]. The TF family Znf_RING-CH also had a functional connection to chromatin modification (e.g., S. cerevisiae Rkr1) [126]. The domain ARID, a 100 amino acid motif, has been found in many eukaryotic TFs [127] such as Swi1 in S. cerevisiae, playing an important role in chromatin remodeling. This domain is also required to transcribe a diverse set of genes including some retrotransposons [128,129].
TFs belonging to the Ste12 family are only found in the fungal kingdom. Except for S. pombe, every genome encodes one copy. Binding to a DNA motif that mediates pheromone response, Ste12 TFs were reported to regulate fungal development and pathogenicity [130] and are involved in mating and pseudohyphae formation [131]. In F. oxysporum, Ste12, downstream of the Fmk1-mediated MAPK cascade, is involved in the control of invasive growth and fungal virulence [110].
Among the others, the Znf_NFX1 domain has been found in NK-X1, a repressor of the human disease-associated gene HLA-DRA [132]. The HMG_box (high mobility group box) in S. cerevisiae is seen in three proteins: Spp41, which is involved in the negative expression regulation of spliceosome components [133]; Nhp6a, which is required for the fidelity of some tRNA genes [134]; and Ixr1, a transcriptional repressor that regulates hypoxic genes [135]. Fep1, an example of Znf_GATA, was reported to be a transcriptional repressor involved in the regulation of some iron transporter genes under high iron concentrations [136]. S. cerevisiae Mbf1, belonging to Cro/C1-type HTH, is a transcriptional coactivator [137].

3.5. Orthologous Survey of TF Families That Were Manually Curated

To further understand expanded TFs and their impacts on transcriptional regulation, we curated a list of 102 TFs reported in the literature focusing on F. oxysporum, F. graminearum, and other phytopathogenic fungi (Table S3 and examples as described in the previous section). Compared to this list of curated TFs using Orthofinder, we defined 80 orthologous groups among the Fusarium genomes (Table S4). Sixty-two out of the 80 orthogroups were identified using the above IPR-annotated pipeline including 17 in Zn2C6, nine in Znf_C2H2, and one containing both the Znf_C2H2 and Zn2-C6 domains (Table S4). This helps add to the functional annotation of these large TF families while also adding additional annotation to 18 TF families (Table S4), accounting for 32 genes per genome (3% of average Fusarium TFome size). These newly annotated TFs include homologs of those without a domain annotation (e.g., disordered proteins F. oxysporum Ren1 [138], M. oryzae Som1 [139], and homologs of those with noncanonical TF domains such as Ankyrin_rpt and WD40_repeat).
We then directly compared F. oxysporum with its Fusarium relatives to calculate the expansion index as follows:
E I f = A v e r a g e   n u m b e r   o f   T F s   i n   F O S C + 1 A v e r a g e   n u m b e r   o f   T F s   i n   F O S C   s i s t e r   s p e c i e s + 1
The EIf ranged from the highest score of 3.54 (Fug, AreA_GATA) to the lowest score of 0.5 (Fox1, Fork_head) (Table S4). Among these 80 orthogroups, 36 groups were conserved (EIf = 1), with one gene per genome. Ten of these conserved groups were functionally validated in F. oxysporum (Table S4). Twenty four groups had scores less than 1, while 20 groups had a score greater than 1 (Table 3 and Table S4). Expanded groups included Fug1 (AreA_GATA, EIf = 3.54), Cos1 (Znf_C2H2, EIf = 2.8), Ftf1/Ftf2 (Zn2-C6, EIf = 2.7), Ebr1/Ebr2 (Zn2-C6, EIf = 2.5), and Ren1 (disordered, EIf = 2). We also identified PacC (EIf = 1.57) as the second most expanded group within the highly expanded Znf_C2H2 family. We will further discuss these six groups (Table 3).
Both Ftf1/Ftf2 and Ebr1/Ebr2, belonging to the Zn2-C6 family, contributes directly to fungal virulence [3,17,97]. The deletion of AC-encoding Ftf1 reduced the pathogenicity of F. oxysporum f. sp. phaseoli [97], highlighting the direct functional involvement of AC TF in virulence. In Fol, deleting either Ftf1 (AC encoding) or Ftf2 (CC encoding) reduced the fungal virulence [96,152]. Constitutive expression of either Ftf1 or Ftf2 induced the expression of effector genes [17]. The core copy Ftf2 was conserved among all Fusarium species, and the AC copy Ftf1 was only found in F. oxysporum and F. redolens (Figure 4). Ebr1 had multiple homologs in F. oxysporum, but a single copy in F. graminearum [103]. The F. oxysporum genome had three AC-encoding paralogs: Ebr2, Ebr3, and Ebr4. Interestingly, these AC-encoding paralogs are regulated by core copy Ebr1 [103]. It is worth noting that the Ebr2 coding sequence driven by the Ebr1 promoter was able to rescue the Ebr1 knockout mutation, indicating some functional redundancy of this family.
Both Cos1 and PacC are part of the Znf_C2H2 family. In M. oryzae, Cos1 was reported to be involved in conidiophore development [112] and functions as a negative regulator reducing fungal pathogenicity [153]. PacC has been reported as an important pH-responsive TF in F. oxysporum [18,107]. This TF family was expanded in clinical strains, showing an average accessory copy number of 3.7 of FOSC, whereas the non-clinical strains showed an average accessory copy number of 0.5. All of the Fusarium relatives’ genomes examined only contained a single copy of the core PacC. Our previous study using one F. oxysporum clinical isolate revealed that the expression of all PacC genes can be induced with a pH shift from 5.0 to 7.4 (the mammalian physiological pH), indicating a potential role in host adaptation [18]. Interestingly, the induction of AC-encoding PacC genes was CC-encoding PacC gene-dependent as the induction disappeared in the CC-encoding PacC knockout mutant, further supporting a cross-talk between the core and accessory TFs. Similar to EBR1, the expression of the AC-encoding PacC genes was much lower than that of the CC-encoding PacC gene, and knockouts of one AC PacC gene affected a small subset of genes compared with the CC PacC knockout, which had a broader effect on the cellular processes [49].
Fug1 has a role in pathogenicity (maize kernel colonization) and fumonisin biosynthesis in F. verticillioides [147]. The deletion of Fug1 increased the sensitivity to the antimicrobial compound 2-benzoxazolinone and to hydrogen peroxide, suggesting its role in mitigating stresses associated with the host defense [147]. Neither CC nor AC-encoding copies of these two genes were experimentally examined in FOSC. Ren1, a disordered protein without IPR functional domain, was expanded with a EIf score of 2 among the FOSC. However, the only reported study on its function is in F. oxysporum f. sp. melonis, regulating the development of conidiation [138].

3.6. Transcriptome Analysis to Probe the Essential TFs during Host Colonization

To understand the functional importance of FOSC TFs, we took advantage of two recently reported transcriptomics datasets [24,25] including pathogenic interactions (Fo5176 infecting Arabidopsis and Fol4287 infecting tomato) and endophytic interactions (Fo47 colonizing Arabidopsis) (Supplementary Dataset).
By examining patterns of expression (Table S5), we found that almost all genes within the conserved category (58 out of 60) (Group 1) were consistently expressed (TPM > 1 across all conditions), supporting their general roles in controlling life processes. Within the expanded category (Group 3), the proportion of genes being consistently expressed ranged from 41% to 59% for core TFs and only from 5% to 16% for AC-encoding TFs. With a less strict filter (TPM > 1 at minimum 1 condition), we found that all genes within the conserved category were expressed. Within the expanded category, 93% of core TFs and between 49% and 67% of AC-encoding TFs were expressed. The significant increase in the AC-encoding TFs with a lower stringency further supported their conditional involvement in niche adaptation.
We further reviewed the expression patterns of the reported TFs in Fol4287 (Table S6). Out of the 27 TFs encoded on the core genome, 18 were upregulated (defined by upregulation under at least three out of four in planta conditions compared to the axenic growth) during plant colonization, which is consistent with their reported roles in pathogenicity. The AC-encoding Ftf1 has been reported to play essential functions in fungal pathogenicity [96]. Of the ten accessory Ffs, eight were upregulated during plant colonization.
Using a higher stringency filter, selecting TFs upregulated under all in planta conditions, we searched for: (1) conserved core TFs that may be related to plant colonization and (2) expanded AC TFs that could be related to host-specific pathogenicity. In the Fol4287, Fo5176, and Fo47 genomes, 95, 62, and 44 core TFs were upregulated during plant colonization, respectively. Among them, ten copies were highly conserved (Table S7) as they were single-copy orthologs across all 15 F. oxysporum strains including Fow2 and Sfl1. Fow2, a Zn2C6 TF, is required for full virulence but not hyphal growth and conidiation in F. oxysporum f. sp. melonis [101]. Sfl1 was reported to be essential for vegetative growth, conidiation, sexual reproduction, and pathogenesis in M. oryzae [88].
Fol4287, Fo5176, and Fo47 contained 29, 34, and nine upregulated accessory TFs, respectively, including Ftf1 and Ren1 (Figure 4 and Table S8). Ftfs have been reported to play an essential role in the pathogenicity in Fol4287, although their involvement in other interactions has not. The Fol4287 genome encoded 10 accessory Ftfs and eight were upregulated during plant colonization. The Fo5176 genome encoded six accessory Ftfs, but only one copy was upregulated during plant colonization. Interestingly, eight upregulated Fol4287 and one upregulated Fo5176 Ftfs were clustered together (Figure 4). The unique expansion with regulatory adaptation (i.e., fine-tuned expression regulation) seemed to be restricted to Fol4287 and not the other pathogenic strain, Fo5176. In contrast, the Fo5176 genome encoded seven accessory Ren1 TF and two were upregulated during plant colonization (Figure 4), while the Fol4287 genome had only one accessory Ren1 not involved in host colonization. While functional validation is needed, strain-specific expansion followed by fine-tuned expression regulation when infecting host species exists and likely contributes to host-specific pathogenicity.

4. Discussion

For a soilborne pathogen with strong host specificity like FOSC, the adjustment of growth and cell cycle control in response to environmental cues is likely to be essential for survival. Expanded TF families likely contribute to the enhanced functions related to niche adaptation as these TF families play important roles in transmitting external and internal signals and regulating complex cellular signaling responses to the sensed stimuli. Therefore, it is not surprising that the genomes of FOSC had larger TFome sizes than the other fungi included in the study. The expansion of TFs among the FOSC resulted in a positive correlation between the total number of proteins and the size of the fungal TFome, which was also observed in other instances [22].
A total of 14 TF families that control the global transcriptional event such as TBP are highly conserved within the ascomycete fungal lineages. Conserved regulatory mechanisms revealed through this study suggest that the plant colonization process could be a common process among FOSC strains regardless of their host-specific pathogenesis. This notion is also supported by recent studies that highlighted the ability of FOSC strains as root colonizers regardless whether of they cause disease or function as endophytes [25,154].
In contrast to these stable TFs, 30 families were expanded to various degrees, and the most significant expansions occurred in the Zn2-C6 and Znf_C2H2 TF families among the FOSC genomes. The number of Zn2-C6 TFs increased significantly (with the highest expansion score) and made up most of the TFs (56.7%) found within the FOSC TFome. For example, Ftf1, a TF belonging to Zn2-C6 and is involved in tomato pathogenicity, was most significantly expanded to 10 copies of accessory Ftfs in the tomato pathogen Fol4287 genome. Eight out of 10 were induced during plant colonization.
The unique expansion of some TFs, driven by ACs, may provide clues as to host-specific interactions. Acquiring additional TFs will modify existing regulatory pathways, and this will require fine-tuning existing networks for this group of organisms to successfully adapt to different hosts under diverse environments. A previous survey of kinome (the complete set of protein kinases encoded in an organism’s genome) among the FOSC and other Ascomycetes also revealed a positive correlation between the size of the kinome and the size of the genome [48], identical to what we reported here for TFomes. As kinases and TFs are key regulators that modulate all important signaling pathways and are essential for the proper functions of almost all molecular and cellular processes, strong correlations between kinome and TFome suggest the ordered recruitment and establishment of ACs among FOSC genomes.
This realization further emphasizes the importance of additional functional studies. Reverse genetics is a powerful tool in defining the functional importance of a TF. For example, TF Ren1, a disordered protein, was identified by genetic and molecular characterization [138]. This TF is most significantly expanded (seven copies of accessory Rens) in the Arabidopsis pathogen Fo5176 genome and is involved in plant colonization. High throughput approaches such as chromatin immunoprecipitation sequencing (CHIP-Seq) and DNA affinity purification sequencing (DAP-Seq) [155] can be used to profile the cis-regulatory elements globally for a better understanding of transcriptional regulation in the fungal model F. oxysporum. Gene regulatory networks [156,157] can add more resolution to these complex regulatory processes. However, the ultimate understanding of the regulatory roles of each TF will come from careful molecular and biochemical characterization.
Our study offers a comprehensive look at the regulation from the evolutionary perspective while also providing an easily implemented computational pipeline to compare TFs and other functional groups in fungi. A better understanding of functions of TFs will not only inform Fusarium biology [158], but can also be extrapolated to other filamentous fungal systems.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof9030359/s1, Table S1. TF-type DNA-binding domains used to filter TFome. Table S2. TFome annotations across 30 genomes. Table S3. Curated TFs in phytopathogens. Table S4. Ortholog copy number of characterized TFs of phytopathogens across 30 genomes. Table S5. Number of TFs being expressed. Table S6. Orthologs of the reported TFs in Fol4287. Table S7. Highly conserved TFs that are constantly upregulated during plant colonization. Table S8. Strain-specific accessory TFs that are upregulated during plant colonization. Figure S1. TFome annotation pipelines. Figure S2. ACs contribute to the FOSC TFome expansion. Figure S3. Phylograms of three conserved families. Figure S4. Minimal gene family contractions in FOSC partially caused by whole genome duplication in yeast. Figure S5. The pipelines to probe the functional important TFs by RNA-Seq data. Supplementary Dataset. Normalized read counts of the RNA-Seq datasets

Author Contributions

Conceptualization, H.Y. (Houlin Yu) and L.-J.M.; Methodology, Analysis and Interpretation, All authors; Resources, S.H., R.D.H. and I.V.G.; Writing—H.Y. (Houlin Yu) and L.-J.M. with input from all authors; Supervision, L.-J.M.; Funding Acquisition, L.-J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the Natural Science Foundation (IOS-165241), the National Institutes of Health (R01EY030150), and the USDA National Institute of Food and Agriculture (MAS00496). The work (proposal:10.46936/10.25585/60001360) conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337, accessed on 7 February 2023), a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231. H.Y. is also supported by the Lotta M. Crabtree Fellowship and Constantine J. Gilgut Fellowship. S.M.-C. is also supported by the Vaadia-BARD Postdoctoral Fellowship.

Institutional Review Board Statement

No applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Unprocessed RNA-Seq data were retrieved from the public source as described in Section 2.2. All processed data supporting our results and conclusions are available in the Supplemental Materials.

Acknowledgments

We thank the Massachusetts Green High-Performance Computing Center for providing the high-performance computing capacity.

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

References

  1. Ma, L.-J.; Geiser, D.M.; Proctor, R.H.; Rooney, A.P.; O’Donnell, K.; Trail, F.; Gardiner, D.M.; Manners, J.M.; Kazan, K. Fusarium Pathogenomics. Annu. Rev. Microbiol. 2013, 67, 399–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ma, L.-J.; Shea, T.; Young, S.; Zeng, Q.; Kistler, H.C. Genome Sequence of Fusarium oxysporum f. sp. melonis Strain NRRL 26406, a Fungus Causing Wilt Disease on Melon. Genome Announc. 2014, 2, e00730-14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Michielse, C.B.; Rep, M. Pathogen profile update: Fusarium oxysporum. Mol. Plant Pathol. 2009, 10, 311–324. [Google Scholar] [CrossRef] [PubMed]
  4. Ploetz, R.C. Fusarium Wilt of Banana. Phytopathology 2015, 105, 1512–1521. [Google Scholar] [CrossRef] [Green Version]
  5. Edel-Hermann, V.; Lecomte, C. Current Status of Fusarium oxysporum Formae Speciales and Races. Phytopathology 2019, 109, 512–530. [Google Scholar] [CrossRef] [Green Version]
  6. Pegg, K.G.; Coates, L.M.; O’Neill, W.T.; Turner, D.W. The Epidemiology of Fusarium Wilt of Banana. Front. Plant Sci. 2019, 10, 1395. [Google Scholar] [CrossRef] [Green Version]
  7. Yang, H.; Yu, H.; Ma, L.-J. Accessory Chromosomes in Fusarium oxysporum. Phytopathology 2020, 110, 1488–1496. [Google Scholar] [CrossRef]
  8. Dean, R.; Van Kan, J.A.L.; Pretorius, Z.A.; Hammond-Kosack, K.E.; Di Pietro, A.; Spanu, P.D.; Rudd, J.J.; Dickman, M.; Kahmann, R.; Ellis, J.; et al. The Top 10 fungal pathogens in molecular plant pathology: Top 10 fungal pathogens. Mol. Plant Pathol. 2012, 13, 414–430. [Google Scholar] [CrossRef] [Green Version]
  9. Rahman, M.Z.; Ahmad, K.; Siddiqui, Y.; Saad, N.; Hun, T.G.; Mohd Hata, E.; Rashed, O.; Hossain, M.I.; Kutawa, A.B. First Report of Fusarium wilt disease on Watermelon Caused by Fusarium oxysporum f. sp. niveum (FON) in Malaysia. Plant Dis. 2021; Online ahead of print. [Google Scholar]
  10. Viljoen, A.; Mostert, D.; Chiconela, T.; Beukes, I.; Fraser, C.; Dwyer, J.; Murray, H.; Amisse, J.; Matabuana, E.L.; Tazan, G.; et al. Occurrence and spread of the banana fungus Fusarium oxysporum f. sp. cubense TR4 in Mozambique. S. Afr. J. Sci. 2020, 116, 1–11. [Google Scholar] [CrossRef]
  11. Halpern, H.C.; Bell, A.A.; Wagner, T.A.; Liu, J.; Nichols, R.L.; Olvey, J.; Woodward, J.E.; Sanogo, S.; Jones, C.A.; Chan, C.T.; et al. First Report of Fusarium Wilt of Cotton Caused by Fusarium oxysporum f. sp. vasinfectum Race 4 in Texas, USA. Plant Disease 2018, 102, 446. [Google Scholar] [CrossRef]
  12. Yu, H.; Ayhan, D.H.; Martínez-Soto, D.; Cochavi, S.M.; Ma, L.-J. Accessory Chromosomes of the Fusarium Oxysporum Species Complex and Their Contribution to Host Niche Adaptation. In Plant Relationships, The Mycota; Scott, B., Mesarich, C., Eds.; Springer International Publishing: Cham, Switzerland, 2023; pp. 371–388. [Google Scholar]
  13. Armstrong, G.M.; Armstrong, J.K. Formae speciales and races of Fusarium oxysporum causing wilt disease. In Fusarium: Disease, Biology, and Taxonomy; Nelson, P.E., Toussoun, T.A., Cook, R.J., Eds.; Pennsylvania State University Press: University Park, PA, USA, 1981; pp. 391–399. [Google Scholar]
  14. Kredics, L.; Narendran, V.; Shobana, C.S.; Vágvölgyi, C.; Manikandan, P.; Indo-Hungarian Fungal Keratitis Working Group. Filamentous fungal infections of the cornea: A global overview of epidemiology and drug sensitivity. Mycoses 2015, 58, 243–260. [Google Scholar] [CrossRef]
  15. Hassan, A.S.; Al-Hatmi AM, S.; Shobana, C.S.; van Diepeningen, A.D.; Kredics, L.; Vágvölgyi, C.; Homa, M.; Meis, J.F.; de Hoog, G.S.; Narendran, V.; et al. Antifungal Susceptibility and Phylogeny of Opportunistic Members of the Genus Fusarium Causing Human Keratomycosis in South India. Med. Myco. 2016, 54, 287–294. [Google Scholar] [CrossRef] [Green Version]
  16. Michielse, C.B.; van Wijk, R.; Reijnen, L.; Manders EM, M.; Boas, S.; Olivain, C.; Alabouvette, C.; Rep, M. The Nuclear Protein Sge1 of Fusarium oxysporum Is Required for Parasitic Growth. PLoS Pathog. 2009, 5, e1000637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. van der Does, H.C.; Fokkens, L.; Yang, A.; Schmidt, S.M.; Langereis, L.; Lukasiewicz, J.M.; Hughes, T.R.; Rep, M. Transcription Factors Encoded on Core and Accessory Chromosomes of Fusarium oxysporum Induce Expression of Effector Genes. PLoS Genet. 2016, 12, e1006401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zhang, Y.; Yang, H.; Turra, D.; Zhou, S.; Ayhan, D.H.; DeIulio, G.A.; Guo, L.; Broz, K.; Wiederhold, N.; Coleman, J.J.; et al. The genome of opportunistic fungal pathogen Fusarium oxysporum carries a unique set of lineage-specific chromosomes. Commun. Biol. 2020, 3, 50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Grigoriev, I.V.; Nikitin, R.; Haridas, S.; Kuo, A.; Ohm, R.; Otillar, R.; Riley, R.; Salamov, A.; Zhao, X.; Korzeniewski, F.; et al. MycoCosm portal: Gearing up for 1000 fungal genomes. Nucl. Acids Res. 2014, 42, D699–D704. [Google Scholar] [CrossRef] [PubMed]
  20. Jones, P.; Binns, D.; Chang, H.-Y.; Fraser, M.; Li, W.; McAnulla, C.; McWilliam, H.; Maslen, J.; Mitchell, A.; Nuka, G.; et al. InterProScan 5, genome-scale protein function classification. Bioinformatics 2014, 30, 1236–1240. [Google Scholar] [CrossRef] [Green Version]
  21. Park, J.; Park, J.; Jang, S.; Kim, S.; Kong, S.; Choi, J.; Ahn, K.; Kim, J.; Lee, S.; Kim, S.; et al. FTFD: An informatics pipeline supporting phylogenomic analysis of fungal transcription factors. Bioinformatics 2008, 24, 1024–1025. [Google Scholar] [CrossRef] [Green Version]
  22. Shelest, E. Transcription Factors in Fungi: TFome Dynamics, Three Major Families, and Dual-Specificity TFs. Front. Genet. 2017, 8, 53. [Google Scholar] [CrossRef] [Green Version]
  23. Emms, D.M.; Kelly, S. OrthoFinder: Phylogenetic orthology inference for comparative genomics. Genome Biol. 2019, 20, 238. [Google Scholar] [CrossRef] [Green Version]
  24. Guo, L.; Yu, H.; Wang, B.; Vescio, K.; DeIulio, G.A.; Yang, H.; Berg, A.; Zhang, L.; Edel-Hermann, V.; Steinberg, C.; et al. Metatranscriptomic Comparison of Endophytic and Pathogenic Fusarium–Arabidopsis Interactions Reveals Plant Transcriptional Plasticity. MPMI 2021, 34, 1071–1083. [Google Scholar] [CrossRef] [PubMed]
  25. Redkar, A.; Sabale, M.; Schudoma, C.; Zechmann, B.; Gupta, Y.K.; López-Berges, M.S.; Venturini, G.; Gimenez-Ibanez, S.; Turrà, D.; Solano, R.; et al. Conserved secreted effectors contribute to endophytic growth and multihost plant compatibility in a vascular wilt fungus. Plant Cell. 2022, 34, 3214–3232. [Google Scholar] [CrossRef]
  26. Cheng, C.-Y.; Krishnakumar, V.; Chan, A.P.; Thibaud-Nissen, F.; Schobel, S.; Town, C.D. Araport11, a complete reannotation of the Arabidopsis thaliana reference genome. Plant J. 2017, 89, 789–804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Fokkens, L.; Guo, L.; Dora, S.; Wang, B.; Ye, K.; Sánchez-Rodríguez, C.; Croll, D. A Chromosome-Scale Genome Assembly for the Fusarium oxysporum Strain Fo5176 To Establish a Model Arabidopsis-Fungal Pathosystem. G3 Genes|Genomes|Genet 2020, 7, 3549–3555. [Google Scholar] [CrossRef]
  28. Wang, B.; Yu, H.; Jia, Y.; Dong, Q.; Steinberg, C.; Alabouvette, C.; Edel-Hermann, V.; Kistler, H.C.; Ye, K.; Ma, L.-J.; et al. Chromosome-Scale Genome Assembly of Fusarium oxysporum Strain Fo47, a Fungal Endophyte and Biocontrol Agent. MPMI 2020, 33, 1108–1111. [Google Scholar] [CrossRef]
  29. Ma, L.-J.; van der Does, H.C.; Borkovich, K.A.; Coleman, J.J.; Daboussi, M.-J.; Di Pietro, A.; Dufresne, M.; Freitag, M.; Grabherr, M.; Henrissat, B.; et al. Comparative genomics reveals mobile pathogenicity chromosomes in Fusarium. Nature 2010, 464, 367–373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef]
  31. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, Y. Evolution of The Pathogenic Fusarium oxysporum through the Lens of Comparative Genomics. Ph.D. Thesis, University of Massachusetts Amherst, Amherst, MA, USA, 2019. [Google Scholar]
  33. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7, Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2, New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [Green Version]
  35. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  36. Hoang, D.T.; Chernomor, O.; von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2, Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef]
  37. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL) v5, an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  38. Goffeau, A.; Barrell, B.G.; Bussey, H.; Davis, R.W.; Dujon, B.; Feldmann, H.; Galibert, F.; Hoheisel, J.D.; Jacq, C.; Johnston, M.; et al. Life with 6000 Genes. Science 1996, 274, 546–567. [Google Scholar] [CrossRef] [Green Version]
  39. Wood, V.; Gwilliam, R.; Rajandream, M.-A.; Lyne, M.; Lyne, R.; Stewart, A.; Sgouros, J.; Peat, N.; Hayles, J.; Baker, S.; et al. The genome sequence of Schizosaccharomyces pombe. Nature 2002, 415, 871–880. [Google Scholar] [CrossRef] [Green Version]
  40. Galagan, J.E.; Calvo, S.E.; Cuomo, C.; Ma, L.-J.; Wortman, J.R.; Batzoglou, S.; Lee, S.-I.; Baştürkmen, M.; Spevak, C.C.; Clutterbuck, J.; et al. Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae. Nature 2005, 438, 1105–1115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Vesth, T.C.; Nybo, J.L.; Theobald, S.; Frisvad, J.C.; Larsen, T.O.; Nielsen, K.F.; Hoof, J.B.; Brandl, J.; Salamov, A.; Riley, R.; et al. Investigation of inter- and intraspecies variation through genome sequencing of Aspergillus section Nigri. Nat. Genet. 2018, 50, 1688–1695. [Google Scholar] [CrossRef] [Green Version]
  42. Galagan, J.E.; Calvo, S.E.; Borkovich, K.A.; Selker, E.U.; Read, N.D.; Jaffe, D.; FitzHugh, W.; Ma, L.-J.; Smirnov, S.; Purcell, S.; et al. The genome sequence of the filamentous fungus Neurospora crassa. Nature 2003, 422, 859–868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Dean, R.A.; Talbot, N.J.; Ebbole, D.J.; Farman, M.L.; Mitchell, T.K.; Orbach, M.J.; Thon, M.; Kulkarni, R.; Xu, J.-R.; Pan, H.; et al. The genome sequence of the rice blast fungus Magnaporthe grisea. Nature 2005, 434, 980–986. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Mesny, F.; Miyauchi, S.; Thiergart, T.; Pickel, B.; Atanasova, L.; Karlsson, M.; Hüttel, B.; Barry, K.W.; Haridas, S.; Chen, C.; et al. Genetic determinants of endophytism in the Arabidopsis root mycobiome. Nat. Commun. 2021, 12, 7227. [Google Scholar] [CrossRef]
  45. Gardiner, D.M.; McDonald, M.C.; Covarelli, L.; Solomon, P.S.; Rusu, A.G.; Marshall, M.; Kazan, K.; Chakraborty, S.; McDonald, B.A.; Manners, J.M. Comparative Pathogenomics Reveals Horizontally Acquired Novel Virulence Genes in Fungi Infecting Cereal Hosts. PLoS Pathog. 2012, 8, e1002952. [Google Scholar] [CrossRef] [Green Version]
  46. Cuomo, C.A.; Guldener, U.; Xu, J.-R.; Trail, F.; Turgeon, B.G.; Di Pietro, A.; Walton, J.D.; Ma, L.-J.; Baker, S.E.; Rep, M.; et al. The Fusarium graminearum Genome Reveals a Link Between Localized Polymorphism and Pathogen Specialization. Science 2007, 317, 1400–1402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Wiemann, P.; Sieber CM, K.; Bargen KW von Studt, L.; Niehaus, E.-M.; Espino, J.J.; Huß, K.; Michielse, C.B.; Albermann, S.; Wagner, D.; Bergner, S.V.; 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, e1003475. [Google Scholar] [CrossRef] [Green Version]
  48. DeIulio, G.A.; Guo, L.; Zhang, Y.; Goldberg, J.M.; Kistler, H.C.; Ma, L.-J. Kinome Expansion in the Fusarium oxysporum Species Complex Driven by Accessory Chromosomes; Mitchell, A.P., Ed.; mSphere: Washington, DC, USA, 2018; Volume 3, p. e00231-18. [Google Scholar]
  49. Yang, H. Accessory Genes Contribute to Rewiring the Transcriptional Network in Fusarium oxysporum. Ph.D. Thesis, University of Massachusetts Amherst, Amhert, MA, USA, 2020. [Google Scholar]
  50. Yu, H.; Ayhan, D.H.; Diener, A.C.; Ma, L.-J. Genome Sequence of Fusarium oxysporum f. sp. matthiolae, a Brassicaceae Pathogen. MPMI 2020, 33, 569–572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Williams, A.H.; Sharma, M.; Thatcher, L.F.; Azam, S.; Hane, J.K.; Sperschneider, J.; Kidd, B.N.; Anderson, J.P.; Ghosh, R.; Garg, G.; et al. Comparative genomics and prediction of conditionally dispensable sequences in legume–infecting Fusarium oxysporum formae speciales facilitates identification of candidate effectors. BMC Genom. 2016, 17, 191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Latrick, C.M.; Marek, M.; Ouararhni, K.; Papin, C.; Stoll, I.; Ignatyeva, M.; Obri, A.; Ennifar, E.; Dimitrov, S.; Romier, C.; et al. Molecular basis and specificity of H2A.Z–H2B recognition and deposition by the histone chaperone YL1. Nat. Struct. Mol. Biol. 2016, 23, 309–316. [Google Scholar] [CrossRef]
  53. Yamaguchi, Y.; Narita, T.; Inukai, N.; Wada, T.; Handa, H. SPT Genes: Key Players in the Regulation of Transcription, Chromatin Structure and Other Cellular Processes. J. Biochem. 2001, 129, 185–191. [Google Scholar] [CrossRef] [PubMed]
  54. Roberts, S.M.; Winston, F. SPT20/ADA5 encodes a novel protein functionally related to the TATA-binding protein and important for transcription in Saccharomyces cerevisiae. Mol. Cell Biol. 1996, 16, 3206–3213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Green, M.R. TBP-associated factors (TAF II s): Multiple, selective transcriptional mediators in common complexes. Trends Biochem. Sci. 2000, 25, 59–63. [Google Scholar] [CrossRef]
  56. Fribourg, S.; Kellenberger, E.; Rogniaux, H.; Poterszman, A.; Van Dorsselaer, A.; Thierry, J.-C.; Egly, J.-M.; Moras, D.; Kieffer, B. Structural Characterization of the Cysteine-rich Domain of TFIIH p44 Subunit. J. Biol. Chem. 2000, 275, 31963–31971. [Google Scholar] [CrossRef] [Green Version]
  57. Goppelt, A.; Stelzer, G.; Lottspeich, F.; Meisterernst, M. A mechanism for repression of class II gene transcription through specific binding of NC2 to TBP-promoter complexes via heterodimeric histone fold domains. EMBO J. 1996, 15, 3105–3116. [Google Scholar] [CrossRef]
  58. Chalabi Hagkarim, N.; Grand, R.J. The Regulatory Properties of the Ccr4–Not Complex. Cells 2020, 9, 2379. [Google Scholar] [CrossRef]
  59. Liang, X.; Shan, S.; Pan, L.; Zhao, J.; Ranjan, A.; Wang, F.; Zhang, Z.; Huang, Y.; Feng, H.; Wei, D.; et al. Structural basis of H2A.Z recognition by SRCAP chromatin-remodeling subunit YL1. Nat. Struct. Mol. Biol. 2016, 23, 317–323. [Google Scholar] [CrossRef]
  60. Kamura, T.; Burian, D.; Khalili, H.; Schmidt, S.L.; Sato, S.; Liu, W.J.; Conrad, M.N.; Conaway, R.C.; Conaway, J.W.; Shilatifard, A. Cloning and characterization of ELL-associated proteins EAP45 and EAP20. a role for yeast EAP-like proteins in regulation of gene expression by glucose. J. Biol. Chem. 2001, 276, 16528–16533. [Google Scholar] [CrossRef]
  61. Vuorio, T.; Maity, S.N.; de Crombrugghe, B. Purification and molecular cloning of the “A” chain of a rat heteromeric CCAAT-binding protein. Sequence identity with the yeast HAP3 transcription factor. J. Biol. Chem. 1990, 265, 22480–22486. [Google Scholar] [CrossRef] [PubMed]
  62. Becker, D.M.; Fikes, J.D.; Guarente, L. A cDNA encoding a human CCAAT-binding protein cloned by functional complementation in yeast. Proc. Natl. Acad. Sci. USA 1991, 88, 1968–1972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Milkereit, P.; Gadal, O.; Podtelejnikov, A.; Trumtel, S.; Gas, N.; Petfalski, E.; Tollervey, D.; Mann, M.; Hurt, E.; Tschochner, H. Maturation and Intranuclear Transport of Pre-Ribosomes Requires Noc Proteins. Cell 2001, 105, 499–509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Fromont-Racine, M.; Senger, B.; Saveanu, C.; Fasiolo, F. Ribosome assembly in eukaryotes. Gene 2003, 313, 17–42. [Google Scholar] [CrossRef]
  65. Edskes, H.K.; Ohtake, Y.; Wickner, R.B. Mak21p of Saccharomyces cerevisiae, a Homolog of Human CAATT-binding Protein, Is Essential for 60 S Ribosomal Subunit Biogenesis. J. Biol. Chem. 1998, 273, 28912–28920. [Google Scholar] [CrossRef] [Green Version]
  66. Xin, C.; Zhang, J.; Nian, S.; Wang, G.; Wang, Z.; Song, Z.; Ren, G. Analogous and Diverse Functions of APSES-Type Transcription Factors in the Morphogenesis of the Entomopathogenic Fungus Metarhizium rileyi. Appl. Environ. Microbiol. 2020, 86, e02928-19. [Google Scholar] [CrossRef] [PubMed]
  67. Pan, X.; Heitman, J. Sok2 Regulates Yeast Pseudohyphal Differentiation via a Transcription Factor Cascade That Regulates Cell-Cell Adhesion. Mol. Cell Biol. 2000, 20, 8364–8372. [Google Scholar] [CrossRef] [Green Version]
  68. Gimeno, C.J.; Fink, G.R. Induction of Pseudohyphal Growth by Overexpression of PHD1, a Saccharomyces cerevisiae Gene Related to Transcriptional Regulators of Fungal Developmentt. Mol. Cell. Biol. 1994, 14, 13. [Google Scholar]
  69. Lysøe, E.; Pasquali, M.; Breakspear, A.; Kistler, H.C. The Transcription Factor FgStuAp Influences Spore Development, Pathogenicity, and Secondary Metabolism in Fusarium graminearum. MPMI 2011, 24, 54–67. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Koch, C.; Moll, T.; Neuberg, M.; Ahorn, H.; Nasmyth, K. A Role for the Transcription Factors Mbp1 and Swi4 in Progression from G1 to S Phase. Science 1993, 261, 1551–1557. [Google Scholar] [CrossRef]
  71. Son, M.; Lee, Y.; Kim, K.-H. The Transcription Cofactor Swi6 of the Fusarium graminearum Is Involved in Fusarium Graminearum Virus 1 Infection-Induced Phenotypic Alterations. Plant Pathol. J. 2016, 32, 281–289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Chikashige, Y.; Yamane, M.; Okamasa, K.; Tsutsumi, C.; Kojidani, T.; Sato, M.; Haraguchi, T.; Hiraoka, Y. Membrane proteins Bqt3 and -4 anchor telomeres to the nuclear envelope to ensure chromosomal bouquet formation. J. Cell Biol. 2009, 187, 413–427. [Google Scholar] [CrossRef] [Green Version]
  73. Sugiyama, M.; Nikawa, J.-I. The Saccharomyces cerevisiae Isw2p-Itc1p Complex Represses INO1 Expression and Maintains Cell Morphology. J. Bacteriol. 2001, 183, 4985–4993. [Google Scholar] [CrossRef] [Green Version]
  74. Lubelsky, Y.; Reuven, N.; Shaul, Y. Autorepression of Rfx1 Gene Expression: Functional Conservation from Yeast to Humans in Response to DNA Replication Arrest. Mol. Cell. Biol. 2005, 25, 10665–10673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Min, K.; Son, H.; Lim, J.Y.; Choi, G.J.; Kim, J.-C.; Harris, S.D.; Lee, Y.-W. Transcription Factor RFX1 Is Crucial for Maintenance of Genome Integrity in Fusarium graminearum. Eukaryot. Cell. 2014, 13, 427–436. [Google Scholar] [CrossRef] [Green Version]
  76. Cairns, B.R.; Lorch, Y.; Li, Y.; Zhang, M.; Lacomis, L.; Erdjument-Bromage, H.; Tempst, P.; Du, J.; Laurent, B.; Kornberg, R.D. RSC, an Essential, Abundant Chromatin-Remodeling Complex. Cell 1996, 87, 1249–1260. [Google Scholar] [CrossRef] [Green Version]
  77. Hsu, J.; Huang, J.; Meluh, P.B.; Laurent, B.C. The Yeast RSC Chromatin-Remodeling Complex Is Required for Kinetochore Function in Chromosome Segregation. Mol. Cell Biol. 2003, 23, 3202–3215. [Google Scholar] [CrossRef] [Green Version]
  78. Olesen, J.T.; Fikes, J.D.; Guarente, L. The Schizosaccharomyces pombe homolog of Saccharomyces cerevisiae HAP2 reveals selective and stringent conservation of the small essential core protein domain. Mol. Cell Biol. 1991, 11, 9. [Google Scholar]
  79. Ridenour, J.B.; Bluhm, B.H. The HAP complex in Fusarium verticillioides is a key regulator of growth, morphogenesis, secondary metabolism, and pathogenesis. Fungal Genet. Biol. 2014, 69, 52–64. [Google Scholar] [CrossRef]
  80. Jung, U.S.; Sobering, A.K.; Romeo, M.J.; Levin, D.E. Regulation of the yeast Rlm1 transcription factor by the Mpk1 cell wall integrity MAP kinase: Reporters for cell wall integrity signalling. Mol. Microbiol. 2002, 46, 781–789. [Google Scholar] [CrossRef] [PubMed]
  81. Ortiz, C.S.; Shim, W.-B. The role of MADS-box transcription factors in secondary metabolism and sexual development in the maize pathogen Fusarium verticillioides. Microbiology 2013, 159, 2259–2268. [Google Scholar] [CrossRef] [Green Version]
  82. Martin, T.; Lu, S.-W.; Tilbeurgh, H.; van Ripoll, D.R.; Dixelius, C.; Turgeon, B.G.; Debuchy, R. Tracing the Origin of the Fungal α1 Domain Places Its Ancestor in the HMG-Box Superfamily: Implication for Fungal Mating-Type Evolution. PLoS ONE 2010, 5, e15199. [Google Scholar] [CrossRef] [Green Version]
  83. Arie, T.; Kaneko, I.; Yoshida, T.; Noguchi, M.; Nomura, Y.; Yamaguchi, I. Mating-type genes from asexual phytopathogenic ascomycetes Fusarium oxysporum and Alternaria alternata. Mol. Plant Microbe. Interact. 2000, 13, 1330–1339. [Google Scholar] [CrossRef] [Green Version]
  84. Albert, T.K.; Hanzawa, H.; Legtenberg, Y.I.; de Ruwe, M.J.; Heuvel, F.A.V.D.; Collart, M.A.; Boelens, R.; Timmers, H. Identification of a ubiquitin-protein ligase subunit within the CCR4-NOT transcription repressor complex. EMBO J. 2002, 21, 355–364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Messenguy, F.; Dubois, E. Role of MADS box proteins and their cofactors in combinatorial control of gene expression and cell development. Gene 2003, 316, 1–21. [Google Scholar] [CrossRef]
  86. Ding, Z.; Xu, T.; Zhu, W.; Li, L.; Fu, Q. A MADS-box transcription factor FoRlm1 regulates aerial hyphal growth, oxidative stress, cell wall biosynthesis and virulence in Fusarium oxysporum f. sp. cubense. Fungal. Biol. 2020, 124, 183–193. [Google Scholar] [CrossRef]
  87. Feder, M.E.; Hofmann, G.E. Heat-Shock Proteins, Molecular Chaperones, and the Stress Response: Evolutionary and Ecological Physiology. Annu. Rev. Physiol. 1999, 61, 243–282. [Google Scholar] [CrossRef] [Green Version]
  88. Li, G.; Zhou, X.; Kong, L.; Wang, Y.; Zhang, H.; Zhu, H.; Mitchell, T.K.; Dean, R.A.; Xu, J.-R. MoSfl1 Is Important for Virulence and Heat Tolerance in Magnaporthe oryzae. PLoS ONE 2011, 6, e19951. [Google Scholar] [CrossRef] [Green Version]
  89. Jiang, C.; Zhang, S.; Zhang, Q.; Tao, Y.; Wang, C.; Xu, J.-R. FgSKN 7 and FgATF 1 have overlapping functions in ascosporogenesis, pathogenesis and stress responses in Fusarium graminearum: Overlapping functions between FgSKN7 and FgATF1. Environ. Microbiol. 2015, 17, 1245–1260. [Google Scholar] [CrossRef]
  90. Paré, A.; Kim, M.; Juarez, M.T.; Brody, S.; McGinnis, W. The Functions of Grainy Head-Like Proteins in Animals and Fungi and the Evolution of Apical Extracellular Barriers. PLoS ONE 2012, 7, e36254. [Google Scholar] [CrossRef] [Green Version]
  91. Schleif, R. AraC protein, regulation of the l-arabinose operon in Escherichia coli, and the light switch mechanism of AraC action. FEMS Microbiol. Rev. 2010, 34, 779–796. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Gallegos, M.-T.; Michán, C.; Ramos, J.L. The XylS/AraC family of regulators. Nucl. Acids Res. 1993, 21, 807–810. [Google Scholar] [CrossRef] [Green Version]
  93. Bustos, S.A.; Schleif, R.F. Functional domains of the AraC protein. Proc. Natl. Acad. Sci. USA 1993, 90, 5638–5642. [Google Scholar] [CrossRef] [Green Version]
  94. Siegmund, T.; Lehmann, M. The Drosophila Pipsqueak protein defines a new family of helix-turn-helix DNA-binding proteins. Dev. Genes Evol. 2002, 212, 152–157. [Google Scholar] [CrossRef] [PubMed]
  95. MacPherson, S.; Larochelle, M.; Turcotte, B. A Fungal Family of Transcriptional Regulators: The Zinc Cluster Proteins. Microbiol. Mol. Biol. Rev. 2006, 70, 583–604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Niño-Sánchez, J.; Casado-Del Castillo, V.; Tello, V.; De Vega-Bartol, J.J.; Ramos, B.; Sukno, S.A.; Díaz Mínguez, J.M. The FTF gene family regulates virulence and expression of SIX effectors in Fusarium oxysporum. Mol. Plant Pathol. 2016, 17, 1124–1139. [Google Scholar] [CrossRef] [Green Version]
  97. Ramos, B.; Alves-Santos, F.M.; García-Sánchez, M.A.; Martín-Rodrigues, N.; Eslava, A.P.; Díaz-Mínguez, J.M. The gene coding for a new transcription factor (ftf1) of Fusarium oxysporum is only expressed during infection of common bean. Fungal Genet. Biol. 2007, 44, 864–876. [Google Scholar] [CrossRef]
  98. Zuriegat, Q.; Zheng, Y.; Liu, H.; Wang, Z.; Yun, Y. Current progress on pathogenicity-related transcription factors in Fusarium oxysporum. Mol. Plant Pathol. 2021, 22, 882–895. [Google Scholar] [CrossRef] [PubMed]
  99. Zhao, S.; An, B.; Guo, Y.; Hou, X.; Luo, H.; He, C.; Wang, Q. Label free proteomics and systematic analysis of secretome reveals effector candidates regulated by SGE1 and FTF1 in the plant pathogen Fusarium oxysporum f. sp. cubense tropical race 4. BMC Genom. 2020, 21, 275. [Google Scholar] [CrossRef] [PubMed]
  100. Rocha, A.L.M.; Di Pietro, A.; Ruiz-Roldán, C.; Roncero, M.I.G. Ctf1, a transcriptional activator of cutinase and lipase genes in Fusarium oxysporum is dispensable for virulence. Mol. Plant Pathol. 2008, 9, 293–304. [Google Scholar] [CrossRef]
  101. Imazaki, I.; Kurahashi, M.; Iida, Y.; Tsuge, T. Fow2, a Zn(II)2Cys6-type transcription regulator, controls plant infection of the vascular wilt fungus Fusarium oxysporum. Mol. Microbiol. 2007, 63, 737–753. [Google Scholar] [CrossRef] [PubMed]
  102. Calero-Nieto, F.; Di Pietro, A.; Roncero MI, G.; Hera, C. Role of the Transcriptional Activator XlnR of Fusarium oxysporum in Regulation of Xylanase Genes and Virulence. MPMI 2007, 20, 977–985. [Google Scholar] [CrossRef] [Green Version]
  103. Jonkers, W.; Xayamongkhon, H.; Haas, M.; Olivain, C.; van der Does, H.C.; Broz, K.; Rep, M.; Alabouvette, C.; Steinberg, C.; Kistler, H.C. EBR1 genomic expansion and its role in virulence of Fusarium species: EBR1 in virulence of Fusarium species. Environ. Microbiol. 2014, 16, 1982–2003. [Google Scholar] [CrossRef]
  104. Fedotova, A.A.; Bonchuk, A.N.; Mogila, V.A.; Georgiev, P.G. C2H2 Zinc Finger Proteins: The Largest but Poorly Explored Family of Higher Eukaryotic Transcription Factors. Acta Nat. 2017, 9, 47–58. [Google Scholar] [CrossRef] [Green Version]
  105. Yun, Y.; Zhou, X.; Yang, S.; Wen, Y.; You, H.; Zheng, Y.; Norvienyeku, J.; Shim, W.-B.; Wang, Z. Fusarium oxysporum f. sp. lycopersici C2H2 transcription factor FolCzf1 is required for conidiation, fusaric acid production, and early host infection. Curr. Genet. 2019, 65, 773–783. [Google Scholar] [CrossRef]
  106. Ruiz-Roldán, C.; Pareja-Jaime, Y.; González-Reyes, J.A.; GRoncero, M.I. The Transcription Factor Con7-1 Is a Master Regulator of Morphogenesis and Virulence in Fusarium oxysporum. MPMI 2015, 28, 55–68. [Google Scholar] [CrossRef] [Green Version]
  107. Caracuel, Z.; Roncero, M.I.G.; Espeso, E.A.; González-Verdejo, C.I.; García-Maceira, F.I.; Di Pietro, A. The pH signalling transcription factor PacC controls virulence in the plant pathogen Fusarium oxysporum: PacC controls virulence in Fusarium. Mol. Microbiol. 2003, 48, 765–779. [Google Scholar] [CrossRef]
  108. López-Berges, M.S. ZafA-mediated regulation of zinc homeostasis is required for virulence in the plant pathogen Fusarium oxysporum. Mol. Plant Pathol. 2020, 21, 244–249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Asunción García-Sánchez, M.; Martín-Rodrigues, N.; Ramos, B.; de Vega-Bartol, J.J.; Perlin, M.H.; Díaz-Mínguez, J.M. fost12, the Fusarium oxysporum homolog of the transcription factor Ste12, is upregulated during plant infection and required for virulence. Fungal Genet. Biol. 2010, 47, 216–225. [Google Scholar] [CrossRef] [PubMed]
  110. Rispail, N.; Di Pietro, A. Fusarium oxysporum Ste12 Controls Invasive Growth and Virulence Downstream of the Fmk1 MAPK Cascade. MPMI 2009, 22, 830–839. [Google Scholar] [CrossRef] [Green Version]
  111. Bader, A.G.; Vogt, P.K. Leucine Zipper Transcription Factors: bZIP Proteins. In Encyclopedic Reference of Genomics and Proteomics in Molecular Medicine; Springer: Berlin/Heidelberg, Germany, 2006; pp. 964–967. [Google Scholar]
  112. Li, X.; Han, X.; Liu, Z.; He, C. The function and properties of the transcriptional regulator COS1 in Magnaporthe oryzae. Fungal. Biol. 2013, 117, 239–249. [Google Scholar] [CrossRef]
  113. López-Berges, M.S.; Capilla, J.; Turrà, D.; Schafferer, L.; Matthijs, S.; Jöchl, C.; Cornelis, P.; Guarro, J.; Haas, H.; Di Pietro, A. HapX-Mediated Iron Homeostasis Is Essential for Rhizosphere Competence and Virulence of the Soilborne Pathogen Fusarium oxysporum. Plant Cell 2012, 24, 3805–3822. [Google Scholar] [CrossRef] [Green Version]
  114. López-Berges, M.S.; Rispail, N.; Prados-Rosales, R.C.; Di Pietro, A. A Nitrogen Response Pathway Regulates Virulence Functions in Fusarium oxysporum via the Protein Kinase TOR and the bZIP Protein MeaB. Plant Cell 2010, 22, 2459–2475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Kim, S.; Park, S.-Y.; Kim, K.S.; Rho, H.-S.; Chi, M.-H.; Choi, J.; Park, J.; Kong, S.; Park, J.; Goh, J.; et al. Homeobox Transcription Factors Are Required for Conidiation and Appressorium Development in the Rice Blast Fungus Magnaporthe oryzae. PLoS Genet. 2009, 5, e1000757. [Google Scholar] [CrossRef] [Green Version]
  116. Honjo, M.; Nakayama, A.; Fukazawa, K.; Kawamura, K.; Ando, K.; Hori, M.; Furutani, Y. A novel Bacillus subtilis gene involved in negative control of sporulation and degradative-enzyme production. J. Bacteriol. 1990, 172, 1783–1790. [Google Scholar] [CrossRef] [Green Version]
  117. Pierce, M.; Benjamin, K.R.; Montano, S.P.; Georgiadis, M.M.; Winter, E.; Vershon, A.K. Sum1 and Ndt80 Proteins Compete for Binding to Middle Sporulation Element Sequences That Control Meiotic Gene Expression. Mol. Cell Biol. 2003, 23, 4814–4825. [Google Scholar] [CrossRef] [Green Version]
  118. Tsuchiya, D.; Yang, Y.; Lacefield, S. Positive Feedback of NDT80 Expression Ensures Irreversible Meiotic Commitment in Budding Yeast. PLoS Genet. 2014, 10, e1004398. [Google Scholar] [CrossRef] [Green Version]
  119. Doyle, C.E.; Kitty Cheung, H.Y.; Spence, K.L.; Saville, B.J. Unh1, an Ustilago maydis Ndt80-like protein, controls completion of tumor maturation, teliospore development, and meiosis. Fungal. Genet. Biol. 2016, 94, 54–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  120. Hepworth, S.R.; Friesen, H.; Segall, J. NDT80 and the Meiotic Recombination Checkpoint Regulate Expression of Middle Sporulation-Specific Genes in Saccharomyces cerevisiae. Mol. Cell Biol. 1998, 18, 5750–5761. [Google Scholar] [CrossRef] [Green Version]
  121. Jones, S. An overview of the basic helix-loop-helix proteins. Genome Biol. 2004, 6, 226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Bahadoor, A.; Brauer, E.K.; Bosnich, W.; Schneiderman, D.; Johnston, A.; Aubin, Y.; Blackwell, B.; Melanson, J.E.; Harris, L.J. Gramillin A and B: Cyclic Lipopeptides Identified as the Nonribosomal Biosynthetic Products of Fusarium graminearum. J. Am. Chem. Soc. 2018, 140, 16783–16791. [Google Scholar] [CrossRef]
  123. Liu, J.; Yuan, Y.; Wu, Z.; Li, N.; Chen, Y.; Qin, T.; Geng, H.; Xiong, L.; Liu, D. A Novel Sterol Regulatory Element-Binding Protein Gene (sreA) Identified in Penicillium digitatum Is Required for Prochloraz Resistance, Full Virulence and erg11 (cyp51) Regulation. PLoS ONE 2015, 10, e0117115. [Google Scholar] [CrossRef] [Green Version]
  124. Donczew, R.; Warfield, L.; Pacheco, D.; Erijman, A.; Hahn, S. Two roles for the yeast transcription coactivator SAGA and a set of genes redundantly regulated by TFIID and SAGA. eLife 2020, 9, e50109. [Google Scholar] [CrossRef] [PubMed]
  125. Iida, T.; Araki, H. Noncompetitive Counteractions of DNA Polymerase ε and ISW2/yCHRAC for Epigenetic Inheritance of Telomere Position Effect in Saccharomyces cerevisiae. Mol. Cell Biol. 2004, 24, 217–227. [Google Scholar] [CrossRef] [Green Version]
  126. Braun, M.A.; Costa, P.J.; Crisucci, E.M.; Arndt, K.M. Identification of Rkr1, a Nuclear RING Domain Protein with Functional Connections to Chromatin Modification in Saccharomyces cerevisiae. Mol. Cell Biol. 2007, 27, 2800–2811. [Google Scholar] [CrossRef] [Green Version]
  127. Iwahara, J. The structure of the Dead ringer-DNA complex reveals how AT-rich interaction domains (ARIDs) recognize DNA. EMBO J. 2002, 21, 1197–1209. [Google Scholar] [CrossRef] [Green Version]
  128. Breeden, L.; Nasmyth, K. Cell cycle control of the yeast HO gene: Cis- and Trans-acting regulators. Cell 1987, 48, 389–397. [Google Scholar] [CrossRef] [PubMed]
  129. Hirschhorn, J.N.; Brown, S.A.; Clark, C.D.; Winston, F. Evidence that SNF2/SWI2 and SNF5 activate transcription in yeast by altering chromatin structure. Genes Dev. 1992, 6, 2288–2298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  130. Rispail, N.; Di Pietro, A. The homeodomain transcription factor Ste12, connecting fungal MAPK signaling to plant pathogenicity. Commun. Integr. Biol. 2010, 3, 327–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  131. Gancedo, J.M. Control of pseudohyphae formation in Saccharomyces cerevisiae. FEMS Microbiol. Rev. 2001, 25, 107–123. [Google Scholar] [CrossRef]
  132. Song, Z.; Krishna, S.; Thanos, D.; Strominger, J.L.; Ono, S.J. A novel cysteine-rich sequence-specific DNA-binding protein interacts with the conserved X-box motif of the human major histocompatibility complex class II genes via a repeated Cys-His domain and functions as a transcriptional repressor. J. Exp. Med. 1994, 180, 1763–1774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Maddock, J.R.; Weidenhammer, E.M.; Adams, C.C.; Lunz, R.L.; Woolford, J.L. Extragenic suppressors of Saccharomyces cerevisiae prp4 mutations identify a negative regulator of PRP genes. Genetics 1994, 136, 833–847. [Google Scholar] [CrossRef]
  134. Braglia, P.; Dugas, S.L.; Donze, D.; Dieci, G. Requirement of Nhp6 Proteins for Transcription of a Subset of tRNA Genes and Heterochromatin Barrier Function in Saccharomyces cerevisiae. Mol. Cell Biol. 2007, 27, 1545–1557. [Google Scholar] [CrossRef] [Green Version]
  135. Vizoso-Vázquez, Á.; Lamas-Maceiras, M.; Becerra, M.; González-Siso, M.I.; Rodríguez-Belmonte, E.; Cerdán, M.E. Ixr1p and the control of the Saccharomyces cerevisiae hypoxic response. Appl. Microbiol. Biotechnol. 2012, 94, 173–184. [Google Scholar] [CrossRef]
  136. Kim, H.-J.; Lee, K.-L.; Kim, K.-D.; Roe, J.-H. The iron uptake repressor Fep1 in the fission yeast binds Fe-S cluster through conserved cysteines. Biochem. Biophys. Res. Commun. 2016, 478, 187–192. [Google Scholar] [CrossRef]
  137. Takemaru, K.; Li, F.-Q.; Ueda, H.; Hirose, S. Multiprotein bridging factor 1 (MBF1) is an evolutionarily conserved transcriptional coactivator that connects a regulatory factor and TATA element-binding protein. Proc. Natl. Acad. Sci. USA 1997, 94, 7251–7256. [Google Scholar] [CrossRef] [Green Version]
  138. 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, 113–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Yan, X.; Li, Y.; Yue, X.; Wang, C.; Que, Y.; Kong, D.; Ma, Z.; Talbot, N.J.; Wang, Z. Two Novel Transcriptional Regulators Are Essential for Infection-related Morphogenesis and Pathogenicity of the Rice Blast Fungus Magnaporthe oryzae. PLoS Pathog. 2011, 7, e1002385. [Google Scholar] [CrossRef] [Green Version]
  140. Yue, X.; Que, Y.; Xu, L.; Deng, S.; Peng, Y.; Talbot, N.J.; Wang, Z. ZNF1 Encodes a Putative C2H2 Zinc-Finger Protein Essential for Appressorium Differentiation by the Rice Blast Fungus Magnaporthe oryzae. Mol. Plant Microbe. Interact. 2016, 29, 22–35. [Google Scholar] [CrossRef] [Green Version]
  141. Bravo-Ruiz, G.; Ruiz-Roldán, C.; Roncero, M.I.G. Lipolytic system of the tomato pathogen Fusarium oxysporum f. sp. lycopersici. Mol. Plant Microbe Interact. 2013, 26, 1054–1067. [Google Scholar] [CrossRef] [Green Version]
  142. Wight, W.D.; Kim, K.-H.; Lawrence, C.B.; Walton, J.D. Biosynthesis and role in virulence of the histone deacetylase inhibitor depudecin from Alternaria brassicicola. Mol. Plant Microbe. Interact. 2009, 22, 1258–1267. [Google Scholar] [CrossRef] [Green Version]
  143. Jones DA, B.; John, E.; Rybak, K.; Phan HT, T.; Singh, K.B.; Lin, S.-Y.; Solomon, P.S.; Oliver, R.P.; Tan, K.-C. A specific fungal transcription factor controls effector gene expression and orchestrates the establishment of the necrotrophic pathogen lifestyle on wheat. Sci. Rep. 2019, 9, 15884. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  144. Oh, M.; Son, H.; Choi, G.J.; Lee, C.; Kim, J.-C.; Kim, H.; Lee, Y.-W. Transcription factor ART1 mediates starch hydrolysis and mycotoxin production in Fusarium graminearum and F. verticillioides. Mol. Plant Pathol. 2016, 17, 755–768. [Google Scholar] [CrossRef] [Green Version]
  145. Dufresne, M.; Perfect, S.; Pellier, A.-L.; Bailey, J.A.; Langin, T. A GAL4-like Protein Is Involved in the Switch between Biotrophic and Necrotrophic Phases of the Infection Process of Colletotrichum lindemuthianum on Common Bean. Plant Cell 2000, 12, 1579–1590. [Google Scholar] [CrossRef] [Green Version]
  146. Son, H.; Fu, M.; Lee, Y.; Lim, J.Y.; Min, K.; Kim, J.-C.; Choi, G.J.; Lee, Y.-W. A novel transcription factor gene FHS1 is involved in the DNA damage response in Fusarium graminearum. Sci. Rep. 2016, 6, 21572. [Google Scholar] [CrossRef] [Green Version]
  147. Ridenour, J.B.; Bluhm, B.H. The novel fungal-specific gene FUG1 has a role in pathogenicity and fumonisin biosynthesis in Fusarium verticillioides. Mol. Plant Pathol. 2017, 18, 513–528. [Google Scholar] [CrossRef] [PubMed]
  148. Jiang, C.; Zhang, C.; Wu, C.; Sun, P.; Hou, R.; Liu, H.; Wang, C.; Xu, J.-R. TRI6 and TRI10 play different roles in the regulation of deoxynivalenol (DON) production by cAMP signalling in Fusarium graminearum. Environ. Microbiol. 2016, 18, 3689–3701. [Google Scholar] [CrossRef]
  149. Schumacher, J.; Simon, A.; Cohrs, K.C.; Viaud, M.; Tudzynski, P. The Transcription Factor BcLTF1 Regulates Virulence and Light Responses in the Necrotrophic Plant Pathogen Botrytis cinerea. PLoS Genet. 2014, 10, e1004040. [Google Scholar] [CrossRef] [Green Version]
  150. Wang, Q.; Pokhrel, A.; Coleman, J.J. The Extracellular Superoxide Dismutase Sod5 from Fusarium oxysporum Is Localized in Response to External Stimuli and Contributes to Fungal Pathogenicity. Front. Plant Sci. 2021, 12, 608861. [Google Scholar] [CrossRef] [PubMed]
  151. Mendoza-Mendoza, A.; Eskova, A.; Weise, C.; Czajkowski, R.; Kahmann, R. Hap2 regulates the pheromone response transcription factor prf1 in Ustilago maydis. Mol. Microbiol. 2009, 72, 683–698. [Google Scholar] [CrossRef] [PubMed]
  152. de Vega-Bartol, J.J.; Martín-Dominguez, R.; Ramos, B.; García-Sánchez, M.-A.; Díaz-Mínguez, J.M. New Virulence Groups in Fusarium oxysporum f. sp. phaseoli: The Expression of the Gene Coding for the Transcription Factor ftf1 Correlates with Virulence. Phytopathology 2011, 101, 470–479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  153. Zhou, Z.; Li, G.; Lin, C.; He, C. Conidiophore Stalk-less1 Encodes a Putative Zinc-Finger Protein Involved in the Early Stage of Conidiation and Mycelial Infection in Magnaporthe oryzae. MPMI 2009, 22, 402–410. [Google Scholar] [CrossRef] [Green Version]
  154. Martínez-Soto, D.; Yu, H.; Allen, K.S.; Ma, L.-J. Differential colonization of the plant vasculature between endophytic versus pathogenic Fusarium oxysporum strains. MPMI 2022, 36, 4–13. [Google Scholar] [CrossRef]
  155. Baumgart, L.A.; Lee, J.E.; Salamov, A.; Dilworth, D.J.; Na, H.; Mingay, M.; Blow, M.J.; Zhang, Y.; Yoshinaga, Y.; Daum, C.G.; et al. Persistence and plasticity in bacterial gene regulation. Nat. Methods 2021, 18, 1499–1505. [Google Scholar] [CrossRef]
  156. Guo, L.; Ji, M.; Ye, K. Dynamic network inference and association computation discover gene modules regulating virulence, mycotoxin and sexual reproduction in Fusarium graminearum. BMC Genom. 2020, 21, 179. [Google Scholar] [CrossRef] [Green Version]
  157. Guo, L.; Zhao, G.; Xu, J.; Kistler, H.C.; Gao, L.; Ma, L. Compartmentalized gene regulatory network of the pathogenic fungus Fusarium graminearum. New Phytol. 2016, 211, 527–541. [Google Scholar] [CrossRef]
  158. Gordon, T.R. Fusarium oxysporum and the Fusarium Wilt Syndrome. Annu. Rev. Phytopathol. 2017, 55, 23–39. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogeny of fungal genomes included in this study. Both left and right phylograms were constructed by the concatenated alignment of randomly selected 500 single-copy orthologous proteins, followed by the maximum likelihood method with 1000 bootstraps. Left shows a phylogram of FOSC (represented by the reference genome Fol4287) together with the other 15 ascomycetes. The right shows a phylogram of members within FOSC, rooted by F. verticillioides (not shown). Only bootstrap values not equal to 100 are shown.
Figure 1. Phylogeny of fungal genomes included in this study. Both left and right phylograms were constructed by the concatenated alignment of randomly selected 500 single-copy orthologous proteins, followed by the maximum likelihood method with 1000 bootstraps. Left shows a phylogram of FOSC (represented by the reference genome Fol4287) together with the other 15 ascomycetes. The right shows a phylogram of members within FOSC, rooted by F. verticillioides (not shown). Only bootstrap values not equal to 100 are shown.
Jof 09 00359 g001
Figure 2. TFome conservation and variation among ascomycete fungi: baseline description. (A) There is a positive correlation between the number of protein-coding genes and TFome size of an organism. JGI fungal genome identifiers were used as labels. (B) Histogram illustrates the distribution of expansion indices among different families. (C) Average number of TFs of the two most drastically expanded families (Znf_C2H2 and Zn2-C6) within each genome set. Genome Set 1 (G1) includes two yeast genomes (S. cerevisiae and S. pombe). Genome Set 2 (G2) includes four filamentous fungal species (N. crassa, A. nidulans, A. acristatulus, and M. oryzae). Genome Set 3 (G3) includes nine sister species close to F. oxysporum. Genome Set 4 (G4) includes 15 FOSC genomes.
Figure 2. TFome conservation and variation among ascomycete fungi: baseline description. (A) There is a positive correlation between the number of protein-coding genes and TFome size of an organism. JGI fungal genome identifiers were used as labels. (B) Histogram illustrates the distribution of expansion indices among different families. (C) Average number of TFs of the two most drastically expanded families (Znf_C2H2 and Zn2-C6) within each genome set. Genome Set 1 (G1) includes two yeast genomes (S. cerevisiae and S. pombe). Genome Set 2 (G2) includes four filamentous fungal species (N. crassa, A. nidulans, A. acristatulus, and M. oryzae). Genome Set 3 (G3) includes nine sister species close to F. oxysporum. Genome Set 4 (G4) includes 15 FOSC genomes.
Jof 09 00359 g002
Figure 3. Evolutional trajectory of heat shock factors (HSFs) suggests genome expansion and adaptation. (A) Phylograms of HSFs were constructed by the maximum likelihood method with 1000 bootstraps. Branches of Fusarium HSFs are colored in yellow. Accessory HSFs of FOSC are shaded in red. (B) Number of accessory HSFs in some FOSC genomes. (C) Expression of HSF genes during plant colonization (hpi indicates hours post inoculation) compared to axenic growth. Transcriptome data were previously described in Guo et al. 2021 [24]. See Section 2 for details of the data reprocessing and visualization.
Figure 3. Evolutional trajectory of heat shock factors (HSFs) suggests genome expansion and adaptation. (A) Phylograms of HSFs were constructed by the maximum likelihood method with 1000 bootstraps. Branches of Fusarium HSFs are colored in yellow. Accessory HSFs of FOSC are shaded in red. (B) Number of accessory HSFs in some FOSC genomes. (C) Expression of HSF genes during plant colonization (hpi indicates hours post inoculation) compared to axenic growth. Transcriptome data were previously described in Guo et al. 2021 [24]. See Section 2 for details of the data reprocessing and visualization.
Jof 09 00359 g003
Figure 4. Unique expansion of some TFs, driven by ACs, may provide clues to host-specific adaptation. RNA-Seq data were previously described [24,25]. (A) Ftf1, the TF involved in tomato pathogenicity was most significantly expanded (10 copies of accessory FTFs) in the tomato pathogen Fol4287 genome and the expression of eight out of 10 were induced during plant colonization. (B) Ren1 was the most significantly expanded (seven copies of accessory RENs) in the Arabidopsis pathogen Fo5176 genome, and two of them were induced during plant colonization.
Figure 4. Unique expansion of some TFs, driven by ACs, may provide clues to host-specific adaptation. RNA-Seq data were previously described [24,25]. (A) Ftf1, the TF involved in tomato pathogenicity was most significantly expanded (10 copies of accessory FTFs) in the tomato pathogen Fol4287 genome and the expression of eight out of 10 were induced during plant colonization. (B) Ren1 was the most significantly expanded (seven copies of accessory RENs) in the Arabidopsis pathogen Fo5176 genome, and two of them were induced during plant colonization.
Jof 09 00359 g004aJof 09 00359 g004b
Table 1. Fungal genomes used in this study.
Table 1. Fungal genomes used in this study.
Fungal Species or StrainsMycoCosm IdentifierGenome Size (MB)No. of GenesTFome SizeHostReference
Saccharomyces cerevisiaeSacce112.076575284 [38]
Schizosaccharomyces pombeSchpo112.615134228 [39]
Aspergillus nidulansAspnid130.4810,680635 [40]
Aspergillus acristatulusAspacri132.5911,221666 [41]
Neurospora crassaNeucr241.049730447 [42]
Magnaporthe oryzaeMagor140.4912,673520Rice[43]
Fusarium solaniFusso152.9317,6561137broad hosts[44]
F. pseudograminearumFusps136.3312,395627Wheat[45]
F. graminearumFusgr136.4513,321608Wheat[46]
F. venenatumFusven137.4512,845802 [44]
F. tricinctumFustr143.6914,106925Broad hosts[44]
F. verticillioidesFusve241.7815,869917Corn[29]
F. fujikuroiFusfu143.8314,813901Broad hosts[47]
F. redolensFusre152.5617,0511098Broad hosts[44]
F. communeFusco148.3715,7311012Broad hosts[44]
F. oxysporum f. sp. cubense (II5)FoxII549.4316,0481047Banana[32]
F. oxysporum f. sp. radicis-lycopersici (CL57)Fusoxrad149.3618,2381151Tomato[48]
F. oxysporum Fo47 (Fo47)FusoxFo47_250.3616,2071082 [28]
F. oxysporum f. sp. lycopersici (MN25)Fusoxlyc148.6417,9311119Tomato[48]
F. oxysporum NRRL26365 (NRRL26365)Fox26365_148.4616,0471036Human[49]
F. oxysporum f. sp. melonis (FoMelon)Fusoxmel154.0319,6611219Melon[2]
F. oxysporum f. sp. lycopersici (Fol4287) Fusox261.3620,9251292Tomato[29]
F. oxysporum NRRL32931 (NRRL32931)Fusox3293147.9117,2801072Human[18]
F. oxysporum MRL8996 (MRL8996)FoxMRL899650.0716,6311057Human[18]
F. oxysporum f. sp. matthiolae (PHW726)FoxPHW726_157.2217,9961157Brassica[50]
F. oxysporum f. sp. vasinfectum (FoCotton)Fusoxvas152.9119,1431189Cotton[48]
F. oxysporum f. sp. pisi (HDV247)Fusoxpis155.1919,6231229Pea[51]
F. oxysporum f. sp. raphani (PHW815)Fusoxrap153.519,3061132Brassica[48]
F. oxysporum f. sp. conglutinans (PHW808)Fusoxcon153.5819,8541142Brassica[48]
F. oxysporum Fo5176 (Fo5176)FoxFo517667.9819,1301236Arabidopsis[27]
Table 2. Expansion index (EIy) of 48 TF families.
Table 2. Expansion index (EIy) of 48 TF families.
IPRTermEIy
Group 1
IPR000814TBP1
IPR003228TFIID_TAF121
IPR004595TFIIH_C1-like1
IPR006809TAFII281
IPR042225Ncb21
IPR008570Vps251
IPR008895Vps72/YL11
IPR007196CNOT11
IPR005612CBF1
IPR001289NFYA1
IPR018004APSES-type HTH1
IPR003150RFX1
IPR033896MADS_MEF2-like1
IPR018501DDT1
Group 2
IPR006856MATalpha_HMGbox0.8
IPR039515NOT40.9
IPR033897MADS_SRF-like0.95
IPR000232HSF0.98
Group 3
IPR003163Tscrpt_reg_HTH_APSES-type1.04
IPR001766Fork_head1.05
IPR011016Znf_RING-CH1.11
IPR001965Znf_PHD1.11
IPR009071HMG_box1.12
IPR004181Znf_MIZ1.24
IPR001606ARID1.25
IPR000679Znf_GATA1.3
IPR001005SANT/Myb1.32
IPR000818TEA/ATTS1.33
IPR003120Ste121.33
IPR003958CBFA_NFYB1.35
IPR001083Cu_fist1.37
IPR000967Znf_NFX11.4
IPR006565Bromodomain1.52
IPR001387Cro/C1-type_HTH1.6
IPR001841Znf_RING1.64
IPR000571Znf_CCCH1.74
IPR001878Znf_CCHC1.83
IPR010666Znf_GRF2
IPR018060HTH_AraC *2
IPR001356Homeobox2.28
IPR007604CP2 *2.73
IPR007396PAI23.42
IPR024061NDT803.47
IPR011598bHLH3.48
IPR007889HTH_Psq *3.53
IPR013087Znf_C2H24.15
IPR004827bZIP5.8
IPR001138Zn2-C615.09
* Indicates the families without a presence in yeasts.
Table 3. Ortholog copy number and expansion index (EIf) of the characterized and expanded TFs in F. oxysporum.
Table 3. Ortholog copy number and expansion index (EIf) of the characterized and expanded TFs in F. oxysporum.
TFReported SpeciesReferencesFamilyOverlap *Average_FoAverage_
Non-Fo
EIf
Ftf1/Ftf2F. oxysporum[96]Zn2-C6Yes4.801.112.75
Ebr1/Ebr2F. oxysporum[103]Zn2-C6Yes5.271.562.45
Znf1M. oryzae[140]Zn2-C6Yes6.472.781.98
Ctf2F. oxysporum[141]Zn2-C6Yes2.931.331.69
Fow2F. oxysporum[101]Zn2-C6Yes2.071.001.53
Dep6A. brassicicola[142]Zn2-C6Yes0.930.671.16
Pf2A. brassicicola[143]Zn2-C6Yes1.201.001.10
Art1F. verticilioides[144]Zn2-C6Yes1.000.891.06
Clta1C. lindemuthianum[145]Zn2-C6Yes1.071.001.03
Fhs1F. graminearum[146]Zn2-C6Yes1.071.001.03
Cos1M. oryzae[112]Znf_C2H2Yes1.800.002.80
PacCF. oxysporum[107]Znf_C2H2Yes2.131.001.57
Fug1F. verticillioides[147]AreA_GATANo7.271.333.54
Ren1F. oxysporum[138]disorderedNo3.001.002.00
Tri10F. graminearum[148]Fun_TFNo1.130.331.60
Ltf1B. cinerea[149]Znf_GATAYes4.002.441.45
Ndt80U. maydis[119]NDT80Yes1.731.111.29
Hap3pF. verticillioides[147]CBFA_NFYBYes1.331.001.17
Sod1F. oxysporum[150]SOD_Cu_ZnNo1.471.221.11
Prf1F. oxysporum[151]HMG_boxYes1.071.001.03
* Column ‘Overlap’ indicates whether orthologous mapping probed families were already included in our domain-based TF annotation.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yu, H.; Yang, H.; Haridas, S.; Hayes, R.D.; Lynch, H.; Andersen, S.; Newman, M.; Li, G.; Martínez-Soto, D.; Milo-Cochavi, S.; et al. Conservation and Expansion of Transcriptional Factor Repertoire in the Fusarium oxysporum Species Complex. J. Fungi 2023, 9, 359. https://doi.org/10.3390/jof9030359

AMA Style

Yu H, Yang H, Haridas S, Hayes RD, Lynch H, Andersen S, Newman M, Li G, Martínez-Soto D, Milo-Cochavi S, et al. Conservation and Expansion of Transcriptional Factor Repertoire in the Fusarium oxysporum Species Complex. Journal of Fungi. 2023; 9(3):359. https://doi.org/10.3390/jof9030359

Chicago/Turabian Style

Yu, Houlin, He Yang, Sajeet Haridas, Richard D. Hayes, Hunter Lynch, Sawyer Andersen, Madison Newman, Gengtan Li, Domingo Martínez-Soto, Shira Milo-Cochavi, and et al. 2023. "Conservation and Expansion of Transcriptional Factor Repertoire in the Fusarium oxysporum Species Complex" Journal of Fungi 9, no. 3: 359. https://doi.org/10.3390/jof9030359

APA Style

Yu, H., Yang, H., Haridas, S., Hayes, R. D., Lynch, H., Andersen, S., Newman, M., Li, G., Martínez-Soto, D., Milo-Cochavi, S., Hazal Ayhan, D., Zhang, Y., Grigoriev, I. V., & Ma, L. -J. (2023). Conservation and Expansion of Transcriptional Factor Repertoire in the Fusarium oxysporum Species Complex. Journal of Fungi, 9(3), 359. https://doi.org/10.3390/jof9030359

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop