Systematic Analysis of Zn2Cys6 Transcription Factors Required for Development and Pathogenicity by High-Throughput Gene Knockout in the Rice Blast Fungus

Because of great challenges and workload in deleting genes on a large scale, the functions of most genes in pathogenic fungi are still unclear. In this study, we developed a high-throughput gene knockout system using a novel yeast-Escherichia-Agrobacterium shuttle vector, pKO1B, in the rice blast fungus Magnaporthe oryzae. Using this method, we deleted 104 fungal-specific Zn2Cys6 transcription factor (TF) genes in M. oryzae. We then analyzed the phenotypes of these mutants with regard to growth, asexual and infection-related development, pathogenesis, and 9 abiotic stresses. The resulting data provide new insights into how this rice pathogen of global significance regulates important traits in the infection cycle through Zn2Cys6TF genes. A large variation in biological functions of Zn2Cys6TF genes was observed under the conditions tested. Sixty-one of 104 Zn2Cys6 TF genes were found to be required for fungal development. In-depth analysis of TF genes revealed that TF genes involved in pathogenicity frequently tend to function in multiple development stages, and disclosed many highly conserved but unidentified functional TF genes of importance in the fungal kingdom. We further found that the virulence-required TF genes GPF1 and CNF2 have similar regulation mechanisms in the gene expression involved in pathogenicity. These experimental validations clearly demonstrated the value of a high-throughput gene knockout system in understanding the biological functions of genes on a genome scale in fungi, and provided a solid foundation for elucidating the gene expression network that regulates the development and pathogenicity of M. oryzae.


Introduction
Magnaporthe oryzae is the best-studied phytopathogenic fungus, which was voted first in the top 10 list of fungal plant pathogens by an international community of molecular plant pathologists [1].The importance of this filamentous ascomycete fungus is not only owing to the fact that the rice blast disease caused by the fungus is the most destructive disease of rice throughout the world, which typically leads to 20-30% losses and even complete loss in grain production during regional epidemics [1,2], but also to its being a primary model in the study of hostfungal pathogen interactions [3].The rice blast fungus has a complicated life cycle including hyphal growth, conidiogenesis, conidial germination, appressorium formation and plant infection, which provides substantial biological information in eukaryotic development and pathogenesis.The rice blast fungus is highly amenable to molecular genetic manipulation, and the functions of numerous genes are identified by gene knockout or ectopic insertion [4][5][6].
Several loss-of-function techniques, such as homologous recombination [4], insertional mutagenesis [6] and RNA interference [26], have been used to investigate gene functions in fungi.Although insertional mutagenesis and RNA interference have been effectively used to inactivate gene functions, they all have limitations that prevent them from becoming high-throughput gene knockout techniques, including 1) inability to cover the genome and the heavy workload to clone genes in the insertional mutagenesis method [6], and 2) the silencing of potential unintended targets by the RNA interference method [26].Gene knockout through DNA homologous recombination is a primary technique used to inactivate gene functions.More than 650 putative TF genes were disrupted in Fusarium graminearum using homologous recombination [27].In that study, the double-joint PCR method [28] was used to generate gene knockout constructs, which were then transformed into fungal protoplasts.In Neurospora crassa, homologous recombination was developed into a high-throughput gene knockout technique [29], in which gene knockout constructs generated by making use of a yeast recombinational cloning method were transformed into conidia by electroporation.Yeast recombinational cloning is a more suitable method than double-joint PCR for high-throughput gene knockout study.However, protoplast transformation is an elaborate, time-consuming and inefficient method, and electroporation of germinating conidia has not been built in the rice blast fungus and many other filamentous fungi.The frequency of homologous recombination after transformation is often quite low in filamentous fungi, typically ,10% in M. oryzae [30].In addition, mutant screening is a tedious and inefficient step in highthroughput gene knockout studies.The implementation of KU70/ KU80 (mus-51/mus-52) mutations as gene knockout background strains greatly increases the frequency of homologous recombination in N. crassa, Aspergilli and other filamentous fungi, and reduces the workload in screening [29,[31][32][33][34].However, the null mutants obtained from KU70/KU80 deletion background strains need to be complemented with native KU70/KU80 genes before analyzing their mutant phenotypes [29].These shortcomings limit the application of the above methods in research in fungal functional genomics.Until now, there is still a surge in interest in functional genomics research through the systematic mutagenesis of identified genes sequenced in the genomes of a large number of fungi (http://www.ncbi.nlm.nih.gov/genome/).
To learn the biological functions of TFs at the genome level, we constructed a high-throughput gene knockout method that enables the rapid knockout of large numbers of genes in the rice blast fungus.In this method, gene knockout vectors are built by a yeast recombinational cloning method using a high-throughput way, DNA transformations are performed by Agrobacterium tumefaciens-mediated transformation (ATMT), and the null mutants are identified by the negative (GFP)/positive (a resistant gene) screening system through a novel yeast-Escherichia-Agrobacterium shuttle vector, pKO1B.With this system, we deleted 104 putative fungal-specific Zn 2 Cys 6 TF genes in M. oryzae.The null mutants were then examined for their phenotypes in development, pathogenicity and responses to stress conditions.In particular, many Zn 2 Cys 6 TF genes required for fungal growth, asexual development, conidial germination and appressorium formation, pathogenicity, and response to stress were identified.We further identified the genes regulated by GPF1 and CNF2 via RNAsequencing (RNA-seq) and found that GPF1 and CNF2 have similar mechanisms in the regulation of gene expression related to fungal pathogenicity.Our findings will provide new insights into the transcriptional regulation of fungal development and pathogenicity, and a way to study the biological functions of genes in fungi at the genome level.

Results
Strategy of the high-throughput gene knockout system via pKO1B vector in M. oryzae Three major challenges in high-throughput gene knockout of fungi are 1) how to build gene-deletion cassettes quickly, 2) how to transfer DNA into fungal cells easily, and 3) how to identify null mutants from the large number of transformants efficiently.To solve these problems, we developed a high-throughput gene knockout system using a newly designed yeast-Escherichia-Agrobacterium shuttle vector, pKO1B, (Figure 1A).Through this plasmid, three techniques suitable for high-throughput manipulation (yeast recombinational cloning to build gene-deletion cassettes, ATMT to transfer DNA into fungal cells and dual selection system to identify null mutants) were combined into the high-throughput gene knockout system.
pKO1B is a vector built on the framework of the binary vector pCAMBIA1300 (www.cambia.org).It contains the URA3-2mi-cro2_origin sequence from a yeast plasmid pYES2 (Invitrogen, USA) and eGFP gene under the control of a strong promoter of M. oryzae H3 histone gene [35] (Figure 1A).pKO1B has a new characteristic: the ability to be replicated in three organisms, Saccharomyces cerevisiae, Escherichia coli and A. tumefaciens.Through the gapped vector pKO1B, three DNA fragments (59 and 39 flanking fragments of the targeted gene and a resistant gene fragment) could be merged into a gene-deletion cassette in one step by yeast recombinational cloning (Figure 1B).The knockout vectors obtained were transformed into A. tumefaciens without transferring the gene-deletion cassettes from a yeast plasmid to another binary Agrobacterium plasmid.The gene-deletion cassettes in pKO1B were transformed into fungal cells using the ATMT method.When the gene-deletion cassettes were ectopically integrated into fungal genomic DNA, GFP was activated in ectopic insertional transformants and was easily observed under a

Author Summary
Magnaporthe oryzae is not only the fungus causing the rice blast disease, which leads to 20-30% losses in rice production, but also a primary model pathosystem for understanding host-pathogen interactions.However, there is no high-throughput gene knockout system constructed, and little is known about most of the genes in this fungus.We developed a high-throughput gene knockout system, and using this system, we obtained null mutants of 104 fungal-specific Zn 2 Cys 6 transcription factor (TF) genes by screening 8741 primary transformants in M. oryzae.We analyzed the functions of these TF genes in development, pathogenesis, and stress responses under 9 conditions.We found that 61 Zn 2 Cys 6 TF genes play indispensable and diversified roles in fungal development and pathogenicity.CNF1 is the first reported TF gene that strongly and negatively regulates asexual development in the rice blast fungus, and CCA1, CNF1, CNF2, CONx1, GPF1, GTA1, MoCOD1 and PCF1 are required for pathogenicity.We further found via RNA-seq that GPF1 and CNF2 have similar mechanisms in gene expression regulation related to pathogenicity.The resulting data provide new insights into how Zn 2 Cys 6 TF genes regulate important traits during the infection cycle of this rice blast pathogen.
fluorescence microscope as a negative selective marker, but not in null mutants (Figure 1C and 1D).
To screen null mutants efficiently, the genomic DNAs of the transformants without green fluorescence were isolated using an improved CTAB method in a high-throughput way (shown in Materials and Methods) and then were detected for the targeted gene and b-tubulin gene by double PCR.If the targeted gene was deleted in a mutant, only one band for b-tubulin appeared as a positive control on an electrophoretic gel; otherwise, there were two bands with one for the targeted gene and the other for btubulin in ectopic insertion transformants (Figure 1C and 1E).For those null mutants identified in the above negative screening PCR, we then continued to search for the unique recombinational DNA fragment that indicated a knockout event and only appeared in the null mutants by PCR.In this PCR, one primer was limited in the genomic DNA outside of the 59 or 39 flanking fragment of the targeted gene, and the other primer was limited in the resistant gene of the gene-deletion cassettes (primers p1+p2 or p3+p4 in Figure 1C).There was one band 1.2-2.0kb in length on the electrophoretic gel that appeared in null mutants.Otherwise, there was no band appearing in the ectopic transformants and wild-type strain (Figure 1F).Gene knockout and ectopic insertion maybe happened coincidentally (Figure 1C).To identify those mutants that were undetected in previous screening steps, copies of transformed gene deletion cassettes in null mutants were confirmed by qPCR after comparison with the wild-type strain using b-tubulin as a control.The null mutants had one copy of the gene deletion cassette, while the null+ectopic transformants had more than two copies and the wild-type strain none.Finally, the mutants containing a single copy of the gene deletion cassette were considered null mutants.
TFs containing fungal-specific Zn 2 Cys 6 zinc finger and fungal_TF_MHR Domain are the largest group of TFs in the rice blast fungus.The 163 putative fungal-specific Zn 2 Cys 6 genes, previously annotated at the Fungal Transcription Factor Database (http://ftfd.snu.ac.kr/index.php)[36] or obtained by BLAST searches at NCBI under organism item (rice blast fungus, taxid: 318829) with the Conserved Domain program [37], were selected to generate gene deletion mutants (Table S1 in Text S1).As a result, 104 TF genes were identified to be deleted in M. oryzae (Table S2 and Table S3 in Text S1).The knockout rate of 104 genes was 15.03612.02%(from 0.68 to 62.07%).The causes for tens of TF genes not deleted in this experiment lie in the following.1) Knockout vectors of twenty-one genes had not been constructed.2) For those genes with a knockout rate less than 1%, more transformants need to be screened.3) Because the SUR gene was selected as a resistance gene, the mutants must be screened on DCM medium.But the mutants of TF genes that are involved with amino acid metabolism cannot grow on DCM medium.4) Some genes may be essential genes.

Phenotypic analyses of fungal-specific Zn 2 Cys 6 transcription factor deletion mutants at developmental stages
We analyzed the phenotypes of the 104 Zn 2 Cys 6 TF null mutants at multiple developmental stages that M. oryzae likely encounters during its infection cycle in rice.These phenotypes analyzed included: developmental characteristics (mycelial growth, colony color and mycelial shape; conidiation; conidial germination; and appressorium formation) and virulence to plants (barley leaf explants, rice leaf explants and rice seedlings).A substantial fraction of the mutants (58.7%, 61/104) clearly displayed visible phenotypes (Figure 2A).The results (Figure 2B) revealed that 27 of the 104 TF genes studied were involved in mycelial growth (Figure S1), 25 TF genes in conidial production, 12 TF genes in conidial germination and 10 TF genes in appressorium formation, and 5 TF genes were involved in pathogenicity in barley (Figure S2), while 7 TF genes were involved in pathogenicity in rice (Figure S3).Table 1 and Table S4 in Text S1 provide a complete list of the phenotypic analyses.Among these mutants displaying clearly visible mutant phenotype changes, 42.6% (26/61) of mutants exhibited multiple mutant phenotypes and 57.4% (35/61) exhibited a single mutant phenotype (Figure 2A).The Zn 2 Cys 6 TFs regulating asexual reproduction also often controlled vegetative growth (conidial germination, colony growth, pigmentation and mycelial appearance) of the fungus (Figure 2C).Interestingly, 7 TF genes required for pathogenicity had multiple mutant phenotypes (Figure 2C, Table 1), and this phenomenon implied that the fungal pathogenic process is complex, where it is closely linked with many fungal development stages, and that the pathogenicity-related Zn 2 Cys 6 TF genes are also involved in other developmental processes.
We identified the homologs of the 104 M. oryzae Zn 2 Cys 6 TF genes in N. crassa and F. graminearum by Blastp (www.broadinstitute.org),and then compared the phenotypes of the knockout mutants with the homologs studied in N. crassa [29] or F. graminearum [27] (Table S5 in Text S1).Deletion of members of the Zn 2 Cys 6 TF family in F. graminearum and N. crassa resulted in mutant phenotypes at 16% (46/296) and 42% (30/72) of the mutants [27,29].However, more mutants (58.7%) in M. oryzae displayed visible mutant phenotypes.But only two mutants showed defects in growth and asexual and/or sexual development in three species simultaneously (Table S5 in Text S1).This comparison of mutant phenotypes between the three fungi suggested that the large majority of Zn 2 Cys 6 TF genes evolved to have a unique, not conserved function in regulating fungal development.

Zn 2 Cys 6 transcription factors involved in fungal growth
The null mutants of 27 Zn 2 Cys 6 TFs showed significant differences in fungal growth with the wild-type strain (Figure 2B, Table 1).Among these null mutants, three mutants (01C9-1/    S4 in Text S1).The TF genes deleted in these mutants were designated GCC1 (growth, conidiation and cell wall regulatory factor 1, MGG_07063), GPF1 (growth and pathogenicity regulatory factor 1, MGG_17841) and GTA1 (growth and tolerance to acidic stress regulatory factor 1, MGG_07149), respectively.To confirm the defects of the mutants in growth were caused by the knockout of the TF genes, we complemented the three mutants Dgcc1, Dgpf1 and Dgta1 with their respective native copies from wild-type strain 70-15 (Figure S4C).The phenotypic analyses showed the GCC1-rescued strain (gcc1-c) and GTA1-rescued strain (gta1-c) all recovered from their defects in growth, and GPF1-rescued strain (gpf1-c) also recovered most defects in growth when compared to the mutants and the wild-type strain (Table 2, Figure 3A).These results implied that Zn 2 Cys 6 TFs GCC1, GPF1 and GTA1 were required for the regulation of fungal growth in the rice blast fungus.
CNF1 is a strong conidial production negative regulatory factor in rice blast fungus, and its deletion may result in greatly increased conidial production (usually 20-to 40-fold of the wild-type strain).Moreover, the mutant Dcnf1 had significant changes in mycelial appearance, such as deep dark colony color from the gray-white of the wild-type strain, short spore-bearing aerial hyphae from long aerial hyphae of the wild-type strain, and spore-bearing aerial hyphae formed at the rim of mycelial colony earlier than in the wild-type strain (Figure 4A).Although the conidiophore development of Dcnf1 seemed similar to that of the wild-type strain and CNF1-rescued mutant cnf1-c (Figure 3B, Figure 4C), more aerial hyphae differentiated into conidiaophores in Dcnf1 (Figure 4B).The increase in conidial production of Dcnf1 was possibly partly due to the spore-bearing hyphae differentiating earlier and more from the aerial hyphae than in the wild-type strain.Dcca1 produced few and odd, long and vacuolated conidia (Table 2, Figure 5).Single or two tandem spores differentiated from the apex of conidiophores in the mutant Dcca1, and the CCA1-rescued mutant cca1-c produced normal conidiophores and conidia similar as the wild-type strain (Figure 4C, Figure 5).These results suggested that CCA1 is required for the differentiation of conidiophores and the formation of conidia.

Zn 2 Cys 6 transcription factors involved in conidial germination and appressorium formation
When the null mutants of 101 Zn 2 Cys 6 TFs were assayed for their conidial germination rate at 4 hpi (hour post inoculation) and their appressorium formation rate at 24 hpi, 12 mutants and 10 mutants displayed significant differences in conidial germination and appressorium formation, respectively (Figure 2B, Table 1).We assayed conidial germination and appressorium formation in the complemented strains of 6 null mutants (Dfzc16, Dgpf1, Dgta1, Dpcf1, Dcnf1 and Dcca1) (Figure S4C) and found that these complemented strains all recovered from their defects in conidial germination or in appressorium formation (Table 2).These experiments implied that Zn 2 Cys 6 TFs GPF1, PCF1, CNF1, and CCA1 are required for conidial germination and appressorium formation and that FZC16 and GTA1 are required for conidial germination in the rice blast fungus.
Among these 6 mutants, the mutant Dcca1 showed marked defects not only in conidial germination and appressorium formation, but also in conidiation (Table 2).The conidial germination rate of the mutant Dcca1 was 3.6%, and appressorium formation rate was 10.2% (Table 2).The structure and shape of conidium, germinated conidium and appressorium of the mutant Dcca1 displayed notable differences with the wild-type strain (Figure 5).Therefore, the mutant Dcca1 produced few multifarious and vacuolated conidia, which germinated and formed not fully melaninized appressoria at a low ratio.
After complemented with native genes (Figure S4C), we retested the pathogenicity of the 7 mutants with their complemented strains on rice seedlings by spraying with conidial suspensions (Dgpf1, Dcnf1, Dcnf2, Dgta1, Dcca1 and gpf1-c, cnf1-c, cnf2-c, gta1-c, and cca1-c) (Figure 6A) or on rice leaf explants by inoculating with mycelial plugs (DMocod1, Dconx1 and Mocod1-c,  conx1-c) (Figure 6B).The facts that the complemented strains recovered from their defects in virulence on rice implied that the loss or weakening of virulence in the mutants was caused by the deletion of the TF genes GPF1, CNF1, CNF2, GTA1, CCA1, CONx1 and MoCOD1, respectively.To assay the roles of these TF genes in plant penetration, onion cuticle and barley leaves were inoculated with a conidial suspension of the mutants and their complemented strains.Compared with the wild-type strain, the ability of Dgpf1, Dgta1, Dcnf2, Dcnf1 and Dcca1 to penetrate the onion cuticle and barley leaf cuticle was weakened to some extent (Figure 6C).Notably, the mutants Dgpf1 and Dcca1 were unable to penetrate barley cuticle, even after barley cuticle was abraded (for Dgpf1), while the wild-type strain penetrated into the epidermal cells of barley leaves and grew invasively after 24 hpi (Figure 6C, Figure 6D).These results suggested that GPF1 is required by the rice blast fungus for penetration into the plant cuticle and possibly for invasive growth.Although the conidia of Dcca1 were unable to penetrate barley cuticle and could not cause lesions in rice, the mycelia of Dcca1 could cause lesions in barley and rice.The difference in virulence between conidia and mycelia may lie in the conidial defects of Dcca1 in germination, appressorium formation and host penetration ability.

Phenotypic analyses of Zn 2 Cys 6 transcription factor deletion mutants under stress conditions
There are various stress conditions encountered by M. oryzae when it completes its infection cycle in rice.It is possible to predict the biological functions of these deleted Zn 2 Cys 6 TF genes in cell activity regulation by assaying the resistance ability to stress conditions of the mutants.We screened and assayed in a rapid way the mycelial growth of the mutants under 9 abiotic stress conditions (0.5 mM H 2 O 2 , 0.005% SDS, 0.3 M CaCl 2 , 0.5 M NaCl, 1 M sorbitol, olive oil as a single organic carbon source, pH 5.0 and pH 9.0, and minimal medium).Four mycelial inoculation blocks were placed in a 9-cm plate by needle point inoculation and the diameter of mycelial colonies were measured after 6 days culture.The results displayed that many mutants had significant visible phenotypes under those stress conditions when compared with the wild-type strain (Table S4 in Text S1): 14 mutants showed mutant phenotypes in MM medium (nitrogen metabolism), 16 mutants in CM-C medium containing olive oil (carbon metabolism), 22 mutants in CM medium with SDS (cell wall stress), 17 mutants or 19 mutants in CM medium with high concentration of Ca 2+ or Na + (ionic regulation), 18 mutants in CM medium with 1 M sorbitol (hypertonic pressure), 30 mutants in CM medium with 0.5 mM H 2 O 2 (oxidative stress), and 30 mutants or 14 mutants in CM medium at pH 5.0 or pH 9.0 (ambient pH regulation).
After comparison of the relative growth rate of mutants on CM medium, we selected 5 mutants (01A7-1, 01B11-1, 01C9-1, 01F1-1 and 02G4-1) that had significant changes in mycelial growth under stress conditions to be complemented with native genes to see if these phenotypes were caused by the deletion of TF genes (Figure S4).The TF genes in these 5 mutants were: MoNIT4 in mutant 01A7-1 (MGG_01518), GPF1 in mutant 01F1-1 (MGG_17841), GCC1 in mutant 01C9-1 (MGG_07063), TAS1 (tolerance to acidic stress regulatory factor 1) in mutant 01B11-1 (MGG_04108) and GTA1 in mutant 02G4-1 (MGG_07149).The mycelial growth of DMonit4 in minimal medium, Dgpf1 in CM-C medium containing olive oil, and Dtas1 at pH 5.0 was reduced (Figure 7A).However, the mycelial growth of Dgcc1 in CM medium containing SDS and Dgta1 at pH 5.0 was increased.When rescued with native genes, these mutants all recovered their growth under stress conditions (Figure 7A).Interestingly, the mutant DMonit4 grew normally as the wild-type strain on CM medium (Table S4 in Text S1, Figure S1), but grew slowly as a very sparse and thin mycelium layer on MM medium (Figure 7B).As previously reported, MoNIT4 is necessary for the regulation of nitrogen metabolism [38,41].These data also suggested that TF gene GPF1 is necessary for gene expression regulation relative to fat metabolism, that GCC1 is necessary for gene expression regulation relative to cell wall synthesis, and that TAS1 and GTA1 are necessary for gene expression regulation in response to acidic stress.

Genes regulated by two virulence-required genes GPF1 and CNF2
Eight Zn 2 Cys 6 TF genes (CCA1, CNF1, CNF2, CONx1, GPF1, GTA1, MoCOD1 and PCF1) are required for fungal pathogenicity to rice or barley.Dgpf1 and Dcnf2 are two mutants that have similar phenotypes on virulence with relatively limited differences in other fields, because DMocod1 and Dconx1 do not produce conidia, Dcca1 and Dgta1 have virulence when inoculated on rice leaves with mycelial plugs, Dpcf1 shows reduced virulence on barley inoculated with mycelial plugs, but normal when sprayed on rice, and Dcnf1 produces greatly increased conidia.To understand how these TF genes affect fungal pathogenicity, the two genes GPF1 and CNF2, whose mutants lose (Dgpf1) or have reduced (Dcnf2) virulence to both barley and rice, were selected to analyze their functions in gene expression regulation.Since Dgpf1 could not penetrate rice cuticle, we could not obtain its mRNA during infection.The rice blast fungus experiences temporal nitrogen starvation when it infects the plant, and the gene expression pattern of the mycelia under starvation is similar to that at plant infection [42][43][44][45], and the hyphae enduring starvation are similar to those inoculated as mycelial plugs in plant infection.We therefore analyzed the genome-wide gene expression of the mycelia after 4 h starvation of the wild-type strain, Dgpf1 and Dcnf2 by RNA-seq.The RNA-seq experiments were performed in biologic triplicate for each strain.
The transcripts of 10,864, 11,048 and 10,681 genes were identified in the starved mycelia of the wild-type strain, Dgpf1 and Dcnf2, respectively.Compared with the wild-type strain, 2641 genes were differentially expressed significantly in Dgpf1 (FDR, 0.05), with 1406 genes up-regulated and 1235 down-regulated (Table S6 in Text S1), and 3144 genes in Dcnf2 (FDR,0.05),with 1668 genes up-regulated and 1476 down-regulated (Table S7 in Text S1).The similarities and differences between genes regulated by GPF1 and CNF2 were analyzed in detail by comparing the DEGs in Dgpf1 or Dcnf2.The comparison assays showed that the expression of 611 genes or 1114 genes were independently regulated by GPF1 or CNF2, while 2030 genes were regulated by both GPF1 and CNF2 together (Figure 8).Among the 2030 genes commonly regulated by GPF1 and CNF2, 916 genes were downregulated and 1105 genes up-regulated simultaneously in both Dgpf1 and Dcnf2, while only 9 genes were regulated in the opposite direction in both mutants (Figure 8, Table S8 in Text S1).It is very surprising that the expression of so many DEGs were regulated in the same direction in Dgpf1 and Dcnf2 (Figure 8).The correlation coefficient between DEGs of Dgpf1 and Dcnf2 was determined by the linear trend model, and the result showed R 2 = 0.88 with p,0.0001.These similar patterns were not mainly caused by the starvation, as we also assayed the mutants of five TF genes containing non-Zn 2 Cys 6 domains by RNA-seq with the same treatment at the same time, but the mutants displayed specific DEG patterns different from Dgpf1 and Dcnf2.Therefore, these similar patterns in gene expression regulation in two mutants were possibly caused by the fact that TFs Gpf1 and Cnf2 had the same DNA-binding domain (Zn 2 Cys 6 domain), and Dgpf1 or Dcnf2 had similar phenotypes in pathogenicity.However, GPF1 and CNF2 still showed great differences in gene regulation since 1152 DEGs appeared in Dcnf2 when compared with Dgpf1 (Table S9 in Text S1).

Discussion
Since the genomes of a large number of fungi have been sequenced and are being sequenced, there is a surge of interest in functional genomics research through the systematic mutagenesis of identified genes.The construction of a genome-wide gene deletion mutant set of fungi based on the homologous recombinational gene knockout procedure is a valuable resource for the analysis of fungal development, pathogenicity, and many aspects of cell biology and biochemistry, such as those done in the budding yeast S. cerevisiae [65], the fission yeast Schizosaccharomyces pombe [66], and the saprobe filamentous fungus N. crassa [29].However, except for N. crassa, F. graminearum and Aspergilli which are highly efficient in gene deletion experiments, no efficient system has been available to perform high-throughput gene knockout in filamentous fungi.In this study, we present an approach to knockout large numbers of genes that utilizes several methods suitable to high-throughput manipulation.Several timeconsuming steps of the gene knockout procedure could be performed in a high-throughput way.In this procedure, the gene-deletion cassettes were built using a yeast homologous recombination method [29] with a yeast-Escherichia-Agrobacterium shuttle vector pKO1B, which is the first reported artificial plasmid that could be replicated in yeast, E. coli and Agrobacterium cells.The gene-deletion cassettes obtained in a binary vector pKO1B could be directly used to transform fungal cells through the ATMT method [67].This advantage avoids the tedious and inefficient work to transfer the gene-deletion cassettes from a yeast plasmid to another binary Agrobacterium plasmid.The use of GFP fluorescence as a negative marker to eliminate most ectopic insertion transformants reduced the number of transformants one has to identify by PCR or Southern blot.The employment of the negative (for the targeted gene)/positive (for the unique recombinational DNA) identification PCR and qPCR to identify null mutants instead of Southern blot, which is not high throughput, makes the identification of null mutants a high-throughput procedure.The copies of the gene-deletion cassette in mutant genomic DNA quantified by qPCR has been widely performed in animals and plants [68,69].The reliability and stability of PCR or qPCR were guaranteed by the extraction of high quality genomic DNA with CTAB [70] performed in a high-throughput way in this study.The result of the gene knockout events (Figure S4A) were also reconfirmed by Southern blot at the DNA level in 10 mutants (Figure S4B) and at the transcript level in 16 mutants (Figure S4C).As a result, two researchers were able to complete one cycle of the gene knockout experiment for 96 genes in a month, including preparing knockout vectors, performing ATMT transformation, and identifying null mutants, which is a major improvement over the individual gene knockout protocol.More importantly, this procedure can be adapted for knocking out genes in other fungi without any modification or after the substitution of organismspecific promoter to drive the GFP reporter.Besides the GFP gene, a second drug resistance gene (HPH/hygromycin B phosphotransferase gene, NEO/neomycin phosphotransferase II gene or BAR/glufosinate resistance gene, etc.) or a herpes simplex virus thymidine kinase (HSVtk) gene, which converts 5-fluoro-29deoxyuridine to a toxic compound [71], could be used as an alternative negative selection marker against ectopic transformants.As an alternative vector, pKO1B-HPH (shown in Materials and Methods) is another yeast-Escherichia-Agrobacterium shuttle which could be used in this high-throughput gene knockout system besides pKO1B.To increase knockout efficiency, M. oryzae KU80 null mutants, which are defective for non-homologous end joining (NHEJ) DNA repair [31], could be alternative strains used in targeted gene replacement instead of the wild-type strain.
In this study, we generated 104 fungal-specific Zn 2 Cys 6 TF gene-deleted mutants in M. oryzae by a high-throughput gene knockout procedure and analyzed the phenotypes of individual TF mutants.The deletion of Zn 2 Cys 6 TF genes resulted in phenotype changes in fungal development and pathogenicity in 58.7% mutants compared to the wild-type strain, while 26 mutants were defective in pleiotropic phenotype.Colony growth and asexual reproduction (conidiation) were the two phenotypic categories most observed in mutants, and most mutants defective in conidiation were often defective in vegetative growth (conidial germination, colony growth, pigmentation and mycelial appearance).Seven Zn 2 Cys 6 TF genes functional for pathogenicity in M. oryzae are also especially required for vegetative growth, conidiation or appressorium formation.Comparison of phenotypes between M. oryzae, F. graminearum and N. crassa Zn 2 Cys 6 TF gene orthologs seemed to agree with the previously reported viewpoint that Zn 2 Cys 6 TFs evolved divergently in how to regulate fungal growth and asexual development, rather than keeping the same function in different fungi [27].These divergences in functions may be due to the differences in the life cycles between the three fungi.N. crassa is an obligate saprophyte which lives on dead organic material and cannot attack a living host, while M. oryzae and F. graminearum are facultative saprophytic plant pathogens that additionally need specialized structures (such as an appressorium) and functions to infect and obtain nutrients from living plants.
Until now, 6 Zn 2 Cys 6 TF genes (MoCOD1, MoCOD2, MoNIT4, PIG1, TRA1 and XLR1) have been identified in the rice blast fungus [17,39,41,72,73], and 4 of them were also in our Zn 2 Cys 6 TF gene knockout mutant set.The expression of MoNIT4 and MoCOD1 were up-regulated during conidiation [41], and their mutants DMonit4 and DMocod1 had reduced conidiation [39,41], and DMocod1 also had lower pathogenicity [39].In our study, the mutants DMonit4 and DMocod1 were also defective in conidiation or pathogenicity; furthermore, our mutant DMocod1 did not produce any conidia on CM medium.Similar to a previous report [74], our DTra1 mutant also showed defects in conidial germination.PIG1 was identified to be involved in melanin biosynthesis, but not confirmed by knockout [73]; however, the Dpig1 mutant did not show notable phenotypes in fungal development in our study.These small discrepancies in mutant phenotypes with previously reported data may be due to the differences in wild-type strains and experimental conditions.TFs Mnh6, Moatf1, Mocrz1 and Mstu1 were reported to regulate hyphal growth in the rice blast fungus, and the colony growth of the mutants Dmnh6, DMoatf1, DMocrz1 and DMstu1 was about 70-90% of the wild-type strain [7][8][9]24].We found 27 Zn 2 Cys 6 TF genes involved in fungal growth, and 3 genes GCC1, GPF1 and GTA1 were required for normal colony growth.Interestingly, the colony growth of Dgta1 was only about 56% of that of wild-type strain, and GTA1 is a gene which affects colony growth the most among known TF genes in the rice blast fungus.GTA1 also functions in conidiation and pathogenicity to plants.
The change in asexual reproduction in a mutant is a phenotype that mostly happens when TF genes or other genes are deleted in the rice blast fungus.The TF gene-deleted mutants (Dcom1, Dcon7, Dcos1, Dmnh6, DMohox2) showed reduced conidiation [7,[11][12][13]15,16].In our study, 25 Zn 2 Cys 6 TF genes were identified as being involved in fungal conidiation.Of them, CCA1, GCC1, MoCOD1 and CONx1 are necessary for the differentiation of conidiophores.Unexpectedly, the deletion of four Zn 2 Cys 6 TF genes (CNF1, CNF2, CNF3 and CNF4) led to increased conidial production.In particular, Cnf1 is the strongest negative regulatory factor of conidial production identified in the rice blast fungus until now.The deletion of CNF1 led to earlier and more differentiation of spore-bearing hyphae than in the wildtype strain.Impressively, the mutants Dcnf1 and Dcnf2 produced more conidia, but their conidia had reduced virulence to barley and rice.Conidia are the main way to spread the rice blast disease.It is clear that a balance exists between the ability to produce conidia and the pathogenicity to plant in the phytopathogenic fungus.
We identified 2641 genes regulated by GPF1 and 3144 genes by CNF2 through RNA-seq.Interestingly, 2021 genes were regulated in the same direction in Dgpf1 and Dcnf2.This fact suggested GPF1 and CNF2 have similar mechanisms in the regulation of fungal pathogenicity.Nearly sixty DEGs in Dgpf1 and Dcnf2 were confirmed to be required for fungal pathogenicity (Table 3).In particular, the expression of two pathogenicity-required genes, an APSES TF gene MSTU1 [8] and a hydrophobin gene MPG1 [75], were greatly down-regulated in Dgpf1 (7.80-and 12.49-fold, respectively) and Dcnf2 (11.85-and 15.83-fold, respectively).About half of Zn 2 Cys 6 TF genes and many other DNA-binding domain TF genes were regulated by GPF1 and CNF2.These data primarily revealed the gene expression network in the regulation of pathogenicity controlled by GPF1 and CNF2.It is necessary to keep the well-balanced expression of TF genes and pathogenicityrequired genes for maintenance of pathogenicity in M. oryzae.
The high-throughput analysis of mutant phenotypes showed that some TF gene mutants shared highly overlapping phenotypes, such as mycelial growth in Dgcc1, Dgpf1 and Dgta1, and pathogenicity in Dgpf1, Dgta1, Dcnf1 and Dcnf2.On the other hand, several mutant phenotypes were shared by one TF gene mutant.Also, it is fascinating to know how to regulate similar phenotypes by different TF genes.In general, every TF gene regulates the expression of many downstream genes and the expression of every gene is regulated by many TF genes, and therefore each mutant phenotype is controlled by the changes in expression of several genes.When two TF genes show similar mutant phenotype, there are at least 3 types of phenotype regulation mechanisms between them.First, TF gene A regulates another TF gene B, and TF gene B continues to regulate a group of downstream genes.Second, TF genes A and B commonly regulate a group of downstream genes.Third, TF genes A and B regulate two different groups of downstream genes independently, but they lead to similar phenotypes.However, two or three of these gene regulation types often occur simultaneously.For example, GPF1 and other TF genes were down-regulated in Dcnf2, while CNF2 and other TF genes were also down-regulated in Dgpf1 (Table S11 in Text S1).A common group of genes (Figure 8) were regulated both by CNF1 and GPF1 which led to similar pathogenicity phenotypes in mutants, while two other different groups of genes (Figure 8) were regulated by CNF1 and GPF1 independently which led to different phenotypes in conidiation in mutants.
In conclusion, this study represents a major advance in a highthroughput gene knockout system suitable for filamentous fungi, provides the functional characterization of 104 Zn 2 Cys 6 TF genes in the rice blast fungus, and reveals gene expression patterns of two virulence-required TF genes GPF1 and CNF2.These studies will help us build more fungal gene-deletion mutant libraries in an economical way and uncover the transcriptional network in fungi and fungal pathogenic mechanisms.

Strains
M. oryzae strain 70-15 and its mutants, S. cerevisiae strain FY834, and E. coli strain DH5a and A. tumefaciens strain AGL1 were used in this study.

Building of a yeast-Escherichi-Agrobacterium shuttle vector pKO1B
A 1184-bp promoter fragment of M. oryzae H3 histone gene amplified from pKD5 [35] with primers H3sF and H3SR (Table S1 in Text S1) and 720-bp eGFP CDS fragment amplified from pEGFP (Clontech, USA) with primers KoGFPF and KoGFPR (Table S1 in Text S1) were inserted into the XhoI/EcoR1 sites and the EcoR1/SacI sites of the binary vector pCAMBIA1300 (Cambia, USA) to produce the vector pKO1.pKO1B was then produced by inserting a 2.9-kb URA3-2micro2_origin fragment amplified from pYES2 (Invitrogen, USA) with primers uraf and urar (Table S1 in Text S1) into the SacII site of pKO1 by a yeast recombinational cloning method using S. c. EasyComp Transformation Kit (Invitrogen, USA).

Primer design and construction of gene-deletion cassettes
A list of 2-kb regions on both sides of each ORF (12775 ORF totally) was retrieved from the rice blast fungus database MG8 (www.broadinstitute.org)and saved in Excel type files by a program written by Mr. Tan Cheng.23-nt primers (Table S1 in Text S1) of each gene-specific flank (1000 bp-1500 bp in length and 1200 bp in optimum length) were designed by the BatchPrimer3 program [76].For each gene, primers were designed and synthesized with the following common 33-nt 59 regions: 5f: GCTGTACAAGTAAGAGCTCGGTACCCGGGGATC… (Homologous to pKO1B) 5r: CCGGGAGATGTGGGGCACTGTGGCGTTGGCACA… (Homologous to SUR gene) 3f: TTGATTATTGCACGGGAATTGCATGCTCTCACA… (Homologous to SUR gene) 3r: TTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCA… (Homologous to pKO1B) The flank fragments were produced from genomic DNA of M. oryzae strain 70-15 using Primer Star or EX Taq (TaKaRa, China) in 96-well PCR plates.The SUR cassette fragment was generated by PCR with primers surf and surr (Table S1 in Text S1) from pBS-SUR [35].All PCR products were verified by agarose gel electrophoresis analysis and then used in subsequent steps without further purification.Among the 163 TF genes, both flanking fragments were successfully amplified for 142 genes.However, one or sometimes both flanking fragments failed to be amplified for 21 genes even after trying several reaction conditions.
The yeast transformation procedure was conducted following a small-scale yeast transformation protocol in the pYES2 user manual (Invitrogen, USA) in 96-well deep well plates as described by Colot et al [29].The cocktail mixture was made by adding 1.8 ml competent yeast cells, 100 ml linearized pKO1B by XbaI and HindIII (100 ng/ml), 210 ml PCR production of SUR cassette, 60 ml denatured salmon sperm DNA (Sangon, China) and 2.6 ml DMSO to 20.7 ml freshly prepared 16LiAc/40%PEG-3350/ 16TE solution (100 mM lithium acetate, pH 7.5; 50% PEG-3350; 10 mM Tris-HCl, pH 7.5; 1 mM EDTA).Next, 200 ml of the mixture, followed by 4 ml of 59 and 39 flank PCR fragments of the targeted gene, were pipetted into each well of a 96-deep-well plate with a multichannel pipette.The plates were sealed and processed further following the user manual.
After 3 days, the yeast cells cultured in SC-Ura liquid medium were collected and disrupted in Fastprep-24 homogenizer (MP, USA) at 4.0 m/s for 2 min.The plasmids were extracted with TIANprep yeast plasmid DNA kit (Tiangen Biotech, China) following the user manual.Next, the yeast plasmids were transformed into the competent cells of E. coli strain DH5a prepared following the Inoue method for ''ultra-competent'' cells [77].Four bacterial colonies on each LB plate were placed in the wells of 96-deep-well plates, with each well containing 1 ml of LB liquid medium with 50 mg/ml kanamycin.The plates were shaken at 220 rpm overnight at 37uC.
The correctness of homologous recombinational cloning for gene-deleted cassettes was confirmed by bacterial double PCR.Common primers were designed to amplify the 59 and 39 flanking fragment of each target gene-deleted cassette.Forward primer KO1Bf1 for the 59 flanking fragment and reverse primer KO1Br2 for the 39 flanking fragment were located inside pKO1B, and reverse primer SURr1 for the 59 flanking fragment and forward primer SURf2 for the 39 flanking fragment were located in SUR cassette (Table S1 in Text S1).The resulting PCR products in knockout plasmids constructed correctly would be 1.0-1.5 kb for 59 flanking fragment and 1.5-2.0kb for 39 flanking fragment.The gene deletion cassettes were further confirmed by sequencing with the primers SURf2 and SURr1.PCR and DNA sequencing results showed a success rate of 100% for the recombinant plasmids (Table S2 in Text S1).We found that it was enough to confirm the recombinant plasmids by PCR and not necessary to confirm by sequencing plasmids.The correctly built knockout plasmids were extracted from the bacterial cultures using the AxyPrep-96 plasmid purification kit (Axygen, China) following the protocol in the user manual.

Agrobacterium tumefaciens-mediated transformation with knockout vectors
The knockout plasmids were transformed into the competent cells of A. tumefaciens strain AGL1 using the freeze/thaw shock transformation method following the procedure described elsewhere [67].The plasmids in Agrobacterium cell were confirmed by culture PCR with primer set SurP1 and SurP2 (Table S1 in Text S1) for a 368-bp SUR gene fragment.
M. oryzae strain 70-15 was grown on CM medium for 12 days at 25uC under constant fluorescent light, and the conidia were harvested and transformed with the knockout plasmids mediated by A. tumefaciens according to a previously reported procedure [67], whereas performed in groups, usually 24 or 48 genes in a batch each time.The nitrocellulose membrane strips containing the conidia co-cultivated with A. tumefaciens were placed on AIM medium in the dark at 22uC for 2 days and then transferred onto the selection defined complex medium (DCM; 0.17% yeast nitrogen base without amino acids, 0.2% asparagine, 0.1% ammonium nitrate and 1% glucose, pH 6.0 with Na 2 HPO 4 ) plates containing 100 mg/ml sulfonylurea, 50 mg/ml kanamycin, 400 mg/ml cefotaxime and 60 mg/ml streptomycin.Sulfonylurearesistant transformants grown on selection medium were individually transferred onto new selection DCM plates using sterile toothpicks.A total of 8741 primary transformants corresponding to 133 Zn 2 Cys 6 TF genes were selected on the basis of sulfonylurea resistance (Table S2 in Text S1).No transformant was obtained for the other nine genes.

Identification of gene-deleted mutants by GFP fluorescence, double PCR and qPCR
After culture for 2 days, a little mycelium of each transformant was picked out and placed onto a glass slide.We usually placed the mycelia of six transformants on a slide in proper order.The green fluorescence emitted by transformants was then observed one by one under a fluorescence microscope.Ectopic transformants emitted green fluorescence and null mutants did not when excited under the fluorescence microscope (Figure 1D).The transformants without green fluorescence were picked out and inoculated on a new selective plate and grown for 3 days.In total, 3191 transformants without green fluorescence were screened from the primary transformants.
Extraction of genomic DNA was performed following the CTAB protocol of Rogers and Bendich [70] with modifications.A small piece (.9 mm 2 ) of mycelium along with some medium in the selective plate was transferred to a 2-ml round-bottom tube, and 400 ml ddH 2 O were added along with 0.1 g porcelain beads.The mycelial cells were disrupted in a Fastprep-24 homogenizer at 4.0 m/s for 2 min; 400 ml 46CATB buffer (4% CTAB, 100 mM Trisma base, 20 mM EDTA, 1.4 M NaCl) were then added to each tube and the tubes incubated at 65uC for 30 min.Next, genomic DNA was extracted following the normal CTAB protocol [70].In total, genomic DNA of the 3191 transformants was extracted with the CTAB method in a high-throughput way.
The gene-deleted mutants from the transformants without green fluorescence were identified by negative screening double PCR.Double PCR was performed using primers CKF and CKR internal to the targeted gene and primers Tbl-gF and Tbl-gR for the b-tubulin gene (Table S1 in Text S1).PCR was performed in 25-ml reaction mixtures in 96-well PCR plates: 0.15 ml primers Tbl-gF and Tbl-gR (20 mM), 0.5 ml primers CKF and CKR (20 mM), 2.5 ml 106PCR buffer, 0.4 ml dNTP mix (25 mM), 0.3 ml Taq (5 U/ml), 19.5 ml ddH 2 O and 1 ml genomic DNA.The PCR program was: 94uC 3 min followed by 35 cycles of 94uC for 30 s, 57uC for 30 s and 72uC for 30 s, and a final extension at 72uC for 10 min.The PCR products were detected by 1.0% agarose gel electrophoresis.If the targeted gene was deleted in a transformant, there was only one band for b-tubulin with 554 bp in length, appearing as a positive control; otherwise, in ectopic transformants and the wild-type strain there were two bands with one of 300-400 bp (the targeted gene) and another of 554 bp (b-tubulin) (Figure 1E).After checking 2280 genomic DNA samples, we identified 477 transformants with deletions in a total of 104 TF genes (Table S2 in Text S1, Figure S4A).No null mutants were found for the other 29 genes after screening all DNA samples.
For the transformants identified as null mutants in the negative screening double PCR, a second PCR was performed to verify the gene deletion event.One primer p1 or p4 (Table S1 in Text S1) was limited in the genomic DNA outside the 59 or 39 flanking fragment in gene-deletion cassettes, and another primer p2 or p3 (Table S1 in Text S1) was limited in the SUR gene in genedeletion cassettes.If the targeted gene was deleted, there was a band with 1.2-2.0kb in length appearing in the gel.Otherwise, there was no band for the ectopic transformant and wild-type strain 70-15 (Figure 1F).A total of 477 transformants of 104 genes identified by the preceding negative screening PCR were also identified as null mutants in this positive screening PCR (Table S2 in Text S1, Figure S4A).
The copies of transformed gene-deletion cassettes in null mutants were identified by qPCR.The concentration of DNA samples of null mutants was standardized to 25 ng/ml by DNA fluorometry.The primers for the genomic DNA of SUR gene are qSurF and qSurR (Table S1 in Text S1).The genomic DNA fragment coding for b-tubulin gene (one copy in the genome) was selected as a control.The primers of b-tubulin gene are qtblF and qtblR (Table S1 in Text S1).PCR mixture of 25 ml was prepared: 12.5 ml 26RT-PCR buffer (SYBR Green, premix Ex Taq, TaKaRa), 0.5 ml forward primer and 0.5 ml reverse primer (20 mM), 2 ml genomic DNA, and 9.5 ml ddH 2 O. Real-time PCR was performed in a Mastercycler (Eppendorf, USA) with the following program: 95uC for 2 min, 40 cycles (95uC for 10 s, 60uC for 20 s), and ending with a melting curve step.Each sample was repeated three times.The copies of transformed gene-deletion cassettes in null mutants were calculated by comparing with the data of the b-tubulin gene and wild-type strain.The exogenous gene was inserted into the mutant's genomic DNA in an integral multiple number.If the copy number of the selective marker gene was 1.060.2times the b-tubulin gene in genomic DNA, the mutant was regarded as containing a single insertion of the selective marker gene.As a result, 477 null mutants of 104 genes were identified as single insertion null mutants (Table S2 and Table S3 in Text S1).
The gene deletion events in mutants of ten randomly selected TF genes were reconfirmed by Southern blot (Figure S4B), which was performed according to a previously reported procedure [7].

Complementation of null mutants with native genes
The mutants were complemented with native gene copies of the wild-type strain 70-15.First, the pKO1B-HPH was built by the replacement of Ph3-GFP cassette with a HPH gene from pCB1003 [78] in pKO1B.The copies of the complementation genes were then cloned from the genomic DNA of the wild-type strain with the primers listed in Table S1 in Text S1 and were inserted into the XbaI/HindIII sites of pKO1B-HPH by the yeast recombinational cloning method.The constructed complementation plasmids were transformed into the mutants using the ATMT method, and the transformants were screened on selective medium containing 200 mg/ml hygromycin B. The gene-rescued transformants were identified by RT-PCR at the mRNA level (Figure S4C) with primers specific for the targeted genes (Table S1 in Text S1).

Phenotypic screening analyses of knockout mutants at development stages
The phenotypic screening analyses were performed by testing 20-30 knockout mutants and one control (the wild-type strain) once in a batch according to previously reported protocols [7,79,80].Mutant phenotypes were assayed in triplicate with five replicates for each strain.The values of colony growth, conidiation, conidial germination and appressorium formation of mutants in the different experimental groups were compared after normalization with the wild-type strain 70-15 in the same group.
Colony growth, morphology and conidiation.Mycelial blocks of 5 mm from 9-day-old wild-type strain 70-15 or mutants were inoculated in the center of CM solid medium in 6-cm plates followed by culture at 25uC under constant fluorescent light.For colony growth, the diameter of the mycelial colony was recorded and colony images were captured at 6 days post inoculation (dpi).For conidiation, the whole conidia of a 7-day-old CM-grown culture in 6-cm plates were harvested and counted using a hemacytometer.

Conidial
germination and appressorium formation.Conidia were harvested from a 10-day-old colony growing on a CM plate by adding 0.5 ml 25 ppm Tween-20 solution to the surface of the culture.The suspension was pipetted up and down once more to ensure enough conidia for the assay.A 40-ml aliquot of spore suspension was dropped onto a sterilized plastic coverslip, which was incubated in a moist chamber at 25uC for 4 and 24 h under dark conditions.A total of 300-500 conidia were counted for conidial germination at 4 hpi and for appressorium formation at 24 hpi.

Phenotypic screening analyses of knockout mutants under stress conditions
The phenotypes of the mutants were also screened under different stress conditions (Table S4 in Text S1) with 20-30 knockout mutants once in a batch and the wild-type strain in an experiment.Four pinpoint-like mycelia of 9-day-old strains were inoculated in the solid medium of different stress conditions in a 9cm plate with a space interval between each other and then incubated at 25uC.Colony images were captured and the diameters of the mycelial colonies were recorded at 6 dpi.The growth rates under stress conditions were compared between strains after normalization with the wild-type strain.
The phenotypes of the mutants of 5 TF genes and their complementation strains under stress conditions were reassessed.Mycelial blocks of 5 mm from 9-day-old strains were inoculated in the center of solid medium of different stress conditions in 6-cm plates followed by culture at 25uC under constant fluorescent light, along with the strains grown on CM medium as controls.Colony images were captured and the diameters of the mycelial colonies were recorded at 8 dpi.The experiments were performed in triplicate with five replicates for each strain.Growth inhibition rate of each mutant was calculated as the growth inhibition rate = (colony diameter on CM medium2colony diameter under stress condition)4colony diameter on CM medium6100.

Pathogenicity tests on barley and rice
The virulence tests on barley (Hordeum vulgare) and rice (Oryza sativa cv CO39) were performed following previously reported protocols [7].For assays with leaf explants of barley or rice, 5-mm mycelium blocks of mutants (along with the wild-type strain 70-15 and mock) were inoculated on the leaf, followed by incubation in a wet box at 25uC for 4 days.For spraying assay on rice seedlings, 4 ml conidial suspension (1610 5 conidia/ml) containing 0.2% (w/ v) gelatin were sprayed onto 15-20 rice seedlings between the third and fourth leaf stages using an artist's airbrush.Inoculated plants were placed in a wet box at 25 C for 2 days and then allowed to grow in controlled environment chambers with a photoperiod of 12 h using fluorescent lights for 5 days.
Plant penetration by the wild-type strain or mutant appressoria was assayed on leaf explants of barley or onion cuticles according to previously reported protocols [7].A droplet of conidia (5610 4 conidia/ml) was inoculated on barley leaf cuticles or onion cuticles and incubated at 25uC for 24 or 48 h.The barley leaves were then treated in methanol (overnight, room temperature) to remove chlorophyll, and fixed in alcoholic lactophenol (1 h, 95% alcohol/ lactophenol = 2:1).Appressorium penetration on barley leaf cuticle or onion cuticles was assessed under a microscope.

RNA-sequencing analysis
The wild-type strain, Dgpf1 and Dcnf2 were grown in liquid CM medium at 25uC with shaking at 180 rpm for 2 days.The cultures were then collected and incubated in H 2 O at 25uC for 4 h.The treatments were repeated in triplicate for each strain.Total RNA was extracted from ground mycelia in liquid nitrogen with the RNeasy Plant Mini Kit (QIAGEN), and mRNA was isolated using AMPure XP beads (Beckman).RNA-seq libraries were constructed using NEBNext RNA sample preparation kit (NEB) in accordance with the standard low-throughput protocol.Samples were sequenced in a 16100 nt way on an Illumina Hiseq2500 instrument using the TruSeq PE Cluster Kit v3 -cBot -HS (Illumina) and TruSeq SBS Kit v3-HS (Illumina).The clean reads were generated by removing adaptor sequences, tags with.10% ''N'', and low quality tags by FastQC (http://www.bioinformatics.babraham.ac.uk), and were mapped to the M. oryzae genome database (MG8) (www.broadinstitute.org)using Tophat software [81].The data from triple biological replicates were then analyzed using Cufflinks software and resulted in quantified genome-wide transcript levels of genes (expressed in fragments per kilobase of exon model per million mapped fragments -FPKM) [82].The significant differences in FPKM between different samples were assayed by the Cuffdiff component of the Cufflinks package [82].

List of genes deleted in this study
The functions of 104 Zn 2 Cys 6 TF genes were studied by knockout and phenotypic analysis in this study (Table 1).

Accession numbers
RNA-sequencing data were deposited in NCBI's Gene Expression Omnibus (GEO accession number GSE57146).The mycelial agar plugs of the mutants of 3 TF genes and the wildtype strain were placed on intact rice leaves for 4 days.(PDF)

Supporting Information
Figure S4 Knockout and complementation of Zn 2 Cys 6 transcription factor genes in M. oryzae.(A) Knockout event of null mutants confirmed by PCR.One mutant of each TF gene (which was assayed in mutant phenotype) was selected as a representative to show the PCR identification results.The size of DNA standards are indicated on the right of lanes (M).WT, M. oryzae strain 70-15; a, bands for b-tubulin; b, bands for the targeted genes; c, bands for unique recombinational DNA fragments relative to the target gene-deletion event.(B) Null mutants of ten randomly selected TF genes confirmed by Southern blot.Genomic DNAs were digested with restriction enzymes shown in Figure S4B and separated on 0.7% agarose gels.The DNAs were individually hybridized with the probes (indicated in Figure S4B).Only one band was detected in mutants and its size was different from that in the wild-type strain, indicating that homologous recombination occurred at a single site.(C) Complementation of 16 Zn 2 Cys 6 transcription factor gene-deleted mutants.The mutants were rescued with their native copy of gene in M. oryzae strain 70-15.RT-PCR results were shown after amplification with 35 cycles.RNA was isolated from the mycelia of the wild-type strain, mutants and complemented strains grown on CM medium.a, the targeted genes; b, b-tubulin gene; 1, mutants; 2, complemented strains; and 3, wild-type strain 70-15.(PDF) Text S1 Supplemental tables.Table S1.Primers used in this study.Table S2.Summary of the gene knockout of 163 Zn 2 Cys 6 transcription factor genes in the rice blast fungus.Table S3.List of fungal-specific Zn 2 Cys 6 zinc finger genes being knocked out and identification of copies of SUR gene in the mutants in this study.Table S4.Phenotypic analyses of the null mutants of 104 Zn 2 Cys 6 transcription factor genes.Table S5.Comparison of phenotypes of previously characterized transcription factors between three fungal mutants (M. oryzae, F. graminearum and N. crassa).Table S6.Differentially expressed genes in Dgpf1 when compared with the wild-type strain (FDR,0.05).Table S7.Differentially expressed genes in Dcnf2 when compared with the wild-type strain (FDR,0.05).Table S8.Differentially expressed genes shared by both Dgpf1 and Dcnf2.Table S9.Differentially expressed genes in Dcnf2 when compared with Dgpf1 (FDR, 0.05).Table S10.GPF1or CNF2-dependent genes that are dispensable for pathogenicity reported in previous studies.Table S11.Zn2Cys6 transcription factor genes that were listed as DEGs of Dgpf1 or Dcnf2.(XLSX)

Figure 1 .
Figure 1.Overview of the high-throughput gene knockout system in fungus.(A) Features of a new binary yeast-Escherichia-Agrobacterium shuttle vector, pKO1B.(B) Building of gene-deletion cassettes in pKO1B by yeast recombinational cloning.The 59 and 39 flanking fragments of the targeted genes were separately amplified from genomic DNA with primers 5f/5r and 3f/3r.Primers 5r and 3f have 59 tails homologous to the SUR cassette, whereas those for 5f and 3r are homologous to the vector.The two flanks were cotransformed into yeast along with the SUR cassette and gapped pKO1B.Homologous recombination created the circular knockout vector, and the final knockout vector was subsequently transformed into A. tumefaciens.(C) Deletion of the targeted gene.The gene-deletion cassette was transformed into the fungal cells via ATMT.Homologous recombination created three types of transformants: null mutants, ectopic insertion transformants, and null and ectopic insertion mutants.The GFP gene was discarded in the null mutants.Primers p1/p2 or p3/p4 were used to identify the unique recombinational DNA fragment, indicating a knockout event.(D) The transformants are screened for GFP fluorescence under a microscope.Putative null mutants do have not GFP fluorescence, but ectopic transformants do.(E) The transformants are screened by double PCR for the targeted gene using the b-tubulin gene as a positive control.The wild-type strain or ectopic transformants produced a characteristic band, indicating the targeted gene, while the null mutants did not.(F) The transformants are screened by PCR for a unique recombinational DNA fragment marked as a knockout event.The null mutants have a 1.2-2.0kb band on an electrophoretic gel, while the wild-type strain and the ectopic transformants do not.doi:10.1371/journal.ppat.1004432.g001

Figure 2 .
Figure 2. Analysis of Zn 2 Cys 6 transcription factor mutant phenotypes in fungal development stages.(A) Number of TFs showing multiple mutant phenotypes.(B) Number of mutants showing mutant phenotypes in each developmental stage.(C) Venn diagram showing the number of mutant phenotypes.The phenotypes included vegetative growth (conidial germination, colony growth, pigmentation and mycelial appearance), conidiation (asexual reproduction), appressorium formation and pathogenicity to rice and barley.doi:10.1371/journal.ppat.1004432.g002 phenotypes of the mutants in colony growth, conidiation, conidial germination, appressorium formation, virulence and response to stress were compared with the wild-type strain 70-15.The genes are named as FZC1,FZC87 (Fungal-specific Zn 2 Cys 6 transcription factor 1 to 87) if not designated in the succeeding text.The detailed phenotypic information of mutants was shown in TableS4in Text S1. doi:10.1371/journal.ppat.1004432.t001DMGG_07149) showed reduced colony growth on CM medium at 24.7, 26.6 and 43.5%, respectively (Table

Figure 4 .
Figure 4. Mycelial appearance and spore-bearing aerial hyphae of the M. oryzae strains.(A) The mycelial appearance of the wild-type strain, Dcnf1 and its complemented strain cnf1-c.Stars indicate nonsporulating hyphae.Bar = 1 cm.(B) Spore-bearing aerial hyphae of the wild-type strain, Dcnf1 and its complemented strain cnf1-c on CM medium.Bar = 20 mm.(C) Conidial development on conidiophores of the M. oryzae strains Dcnf1 and Dcca1.The conidiophore pictures of the wild-type, mutants (Dcnf1 and Dcca1), and their complemented strains (cnf1-c and cca1-c) were taken on the strains grown on CM medium on microscope slides for 1 day.Bar = 40 mm.doi:10.1371/journal.ppat.1004432.g004

Figure 8 .
Figure 8.Comparison of differentially expressed genes between Dgpf1 and Dcnf2.Green or red arrows are the directions of gene expression regulation in Dgpf1 or Dcnf2.doi:10.1371/journal.ppat.1004432.g008

Figure
Figure S1 The colonies of M. oryzae strains on CM medium.The mutants of 104 Zn 2 Cys 6 transcription factor genes and the wild-type strain 70-15 were cultured at 25uC for 6 days.(PDF) Figure S2 Pathogenicity screening assay of the mutants on barley leaf explants.The mycelial agar plugs of the mutants of 104 Zn 2 Cys 6 transcription factor genes and the wildtype strain 70-15 were placed on intact barley leaves for 4 days.(PDF) Figure S3 Pathogenicity screening assay of the mutants of 104 Zn 2 Cys 6 transcription factor genes on rice.(A) The rice seedlings were sprayed with conidial suspension (1610 5 spores/ml) of 101 M. oryzae mutants and cultured for 7 days.(B)The mycelial agar plugs of the mutants of 3 TF genes and the wildtype strain were placed on intact rice leaves for 4 days.(PDF)

Table 1 .
Phenotypic summary of 104 Zn

Table 2 .
Developmental characteristics of M. oryzae strains.
Notes: The strains (5-mm mycelial blocks) were grown on CM medium for 8 days, and the diameter of colonies was then measured (growth) or conidia were collected and counted (conidiation).40-ml (1610 5 conidia/ml) conidial suspensions were incubated on plastic slides for 4 hpi (conidial germination assay) and 24 hpi (appressorium formation assay).*Same capital letters indicate non-significant difference estimated by Duncan's test in each developmental item (P#0.

Table 3 .
Cont.The expression of genes increased or decreased in Dgpf1 and Dcnf2.The number means the fold change, and the symbol ''2'' means no significant change (FDR,0.05)afterwards compared with the wild-type strain.doi:10.1371/journal.ppat.1004432.t003 a)