Weighted Gene Co-Expression Network Coupled with a Critical-Time-Point Analysis during Pathogenesis for Predicting the Molecular Mechanism Underlying Blast Resistance in Rice

Background Rice blast, caused by the ascomycete fungus M. oryzae, is one of the most important diseases of rice. Although many blast resistance (R) genes have been identified and deployed in rice varieties, the molecular mechanisms responsible for the R gene-mediated defense responses are yet not fully understood. Results In this study, we used comparative transcriptomic analysis to explore the molecular mechanism involved in Piz-t-mediated resistance in a transgenic line containing Piz-t (NPB-Piz-t) compared to Nipponbare (NPB). Clustering and principal component analysis (PCA) revealed that the time-point at 24-h post inoculation (hpi) was the most important factor distinguishing the four time-points, which consisted of four genes of mitogen-activated protein kinases (MAPKs) signaling pathway, one gene related to WRKY DNA-binding domain containing protein, five pathogenesis-related protein (OsPR1s) genes, and three genes of R proteins involving in the most significant protein-protein interaction (PPI) pathway. Using weighted gene co-expression network analysis (WGCNA) to investigate RNA-seq data across 0, 24, 48, and 72 hpi, nine modules with similar patterns expression pattern (SEP) and three modules with differential expression pattern (DEP) between NPB-Piz-t and NPB across 0, 24, 48, and 72 hpi with KJ201 (referred to as Piz-t-KJ201 and NPB-KJ201) were identified. Among these the most representative SEP green-yellow module is associated with photosynthesis, and DEP pink module comprised of two specific expressed nucleotide-binding domain and leucine-rich repeat (NLR) genes of LOC_Os06g17900 and LOC_Os06g17920 of Pi2/9 homologous, three NLR genes of LOC_Os11g11810, LOC_Os11g11770, and LOC_Os11g11920 which are putatively associated with important agronomic traits, and a B3 DNA binding domain containing protein related genes (LOC_Os10g39190). Knockout of LOC_Os10g39190 via CRISPR-Cas9 resulted in plant death at the seedling stage. Conclusions The research suggested that Piz-t and multiple NLR network might play important roles in the regulation of the resistance response in the Piz-t-KJ201 interaction system. The identified genes provide an NLR repository to study the rice-M. oryzae interaction system and facilitate the breeding of blast-resistant cultivars in the future.