Analysis of RNA-Seq-Based Expression Profiles ofΔznf1 Mutants inMagnaporthe oryzae

YUE Xiao-feng,QUE Ya-wei,WANG Zheng-yi
DOI: https://doi.org/10.3864/j.issn.0578-1752.2016.17.009
2016-01-01
Abstract:Objective]Magnaportheoryzae is the causal agent of rice blast, which is one of the most important disease threatening the production of cultivated rice worldwide. Appressorium-mediated penetration is a key step in the disease cycle of the fungus. Previously, it was reported that a C2H2 zinc finger transcription factor encoded byZNF1 is essential for appressorium development, penetration and pathogenicity in the rice blast fungus. The objective of this study is to understand the regulatory mechanism of Znf1 and reveal the genes transcriptionally regulated by Znf1, thus providing new clues for further investigating molecular mechanism of pathogenicity in this fungus.[Method]The transcriptome profiles of vegetative mycelia of the wild-type strain Guy11 and aΔznf1mutant were assayed with the RNA-Seq technique. The gene expression levels were calculated using the FPKM method. The criteria of false discovery rate (FDR)≤0.001 and absolute value of log2 ratio≥1 were used to identify differentially expressed genes (DEGs). The sequences of the DEGs were subjected to BLAST queries against the gene ontology (GO) database and KEGG pathway database to predict their biological function and pathways. In order to define in more detail about the sub-set of genes regulated by Znf1, transcriptome profiles of a mutant lacking thePMK1 MAP kinase-encoding gene was also analyzed based on the RNA-Seq technique. To identify the genes regulated by both Znf1 and Pmk1, the DEGs betweenΔznf1andΔpmk1 were compared. In addition, the genes highly expressed during appressorium formation but down-regulated in eitherΔznf1orΔpmk1 were obtained by comparison with the previous transcriptional profile data.[Result] Totally, 709 DEGs in theΔznf1mutant, including 299 up-regulated and 410 down-regulated genes, were identified by comparison with the wild-type strain Guy11. Gene ontology enrichment analysis showed that 118, 299 and 308 DEGs were classified into cellular component, molecular function and biological process, respectively. KEGG pathway enrichment analysis revealed that the DEGs were mainly involved in metabolic pathways, biosynthesis of secondary metabolites and glycerophospholipid metabolism. Several known pathogenicity-related genes, includingLPP3,HOX7,PBS2andMPG1, were found down-regulated inΔznf1. The comparison of DEGs showed that about 56% DEGs inΔznf1 shared identical to those inΔpmk1. Three isotrichodermin C-15 hydroxylase encoding genes,MGG_03825, MGG_02329 andMGG_08498, were significantly down-regulated in bothΔznf1andΔpmk1. In addition, 48 genes up-regulated during appressorium formation were down-regulated in the two mutants, indicating that these putative appressorium-associated genes were regulated directly or indirectly by Znf1 and Pmk1. To confirm the reliability of the RNA-Seq data, 10 DEGs were randomly selected for qRT-PCR. The results showed that the expression patterns in qRT-PCR were consistent with those in RNA-Seq.[Conclusion] The expression profiling data and predicted molecular function of Znf1-dependent DEGs were obtained by RNA-Seq technique. Several pathogenicity-associated genes were regulated by Znf1. Additionally, several genes highly expressed during appressoria formation were also regulated by Znf1 as well as Pmk1. This study provided valuable information for further research on Znf1 downstream gene regulatory network.
What problem does this paper attempt to address?