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.