Homotypic clustering of OsMYB4 binding site motifs in promoters of the rice genome and cellular-level implications on sheath blight disease resistance

Rice WRKY13 TF protein binds to motifs in the promoter region to regulate downstream disease resistance-related genes

In plants, pathogen resistance is brought about by the binding of certain transcription factor (TF) proteins to the cis-elements of certain target genes. These cis-elements are present upstream in the motif of the promoters of each gene. This ensures the binding of a specific TF to a specific promoter, therefore regulating the expression of that gene. Therefore, the study of each promoter sequence of all the rice genes would help identify the target genes of a specific TF. Rice 1 kb upstream promoter sequences of 55,986 annotated genes were analyzed using the Perl program algorithm to detect WRKY13 binding motifs (bm). The resulting genes were grouped using Gene Ontology and gene set enrichment analysis. A gene with more than 4 TF bm in their promoter was selected. Ten genes reported to have a role in rice disease resistance were selected for further analysis. Cis-acting regulatory element analysis was carried out to find the cis-elements and confirm the presence of the corresponding motifs in the promoter sequences of these genes. The 3D structure of WRKY13 TF and the corresponding ten genes were built, and the interacting residues were determined. The binding capacity of WRKY13 to the promoter of these selected genes was analyzed using docking studies. WRKY13 was considered for docking analysis based on the prior reports of autoregulation. Molecular dynamic simulations provided more details regarding the interactions. Expression data revealed the expression of the genes that helped provide the mechanism of interaction. Further co-expression network helped to characterize the interaction of these selected disease resistance-related genes with the WRKY13 TF protein. This study suggests downstream target genes that are regulated by the WRKY13 TF. The molecular mechanism involving the gene network regulated by WRKY13 TF in disease resistance against rice fungal pathogens is explored.

Assessment of Rice Sheath Blight Resistance Including Associations with Plant Architecture, as Revealed by Genome-Wide Association Studies

BACKGROUND

Sheath blight (ShB) disease caused by Rhizoctonia solani Kühn, is one of the most economically damaging rice (Oryza sativa L.) diseases worldwide. There are no known major resistance genes, leaving only partial resistance from small-effect QTL to deploy for cultivar improvement. Many ShB-QTL are associated with plant architectural traits detrimental to yield, including tall plants, late maturity, or open canopy from few or procumbent tillers, which confound detection of physiological resistance.

RESULTS

To identify QTL for ShB resistance, 417 accessions from the Rice Diversity Panel 1 (RDP1), developed for association mapping studies, were evaluated for ShB resistance, plant height and days to heading in inoculated field plots in Arkansas, USA (AR) and Nanning, China (NC). Inoculated greenhouse-grown plants were used to evaluate ShB using a seedling-stage method to eliminate effects from height or maturity, and tiller (TN) and panicle number (PN) per plant. Potted plants were used to evaluate the RDP1 for TN and PN. Genome-wide association (GWA) mapping with over 3.4 million SNPs identified 21 targeted SNP markers associated with ShB which tagged 18 ShB-QTL not associated with undesirable plant architecture traits. Ten SNPs were associated with ShB among accessions of the Indica subspecies, ten among Japonica subspecies accessions, and one among all RDP1 accessions. Across the 18 ShB QTL, only qShB4-1 was not previously reported in biparental mapping studies and qShB9 was not reported in the GWA ShB studies. All 14 PN QTL overlapped with TN QTL, with 15 total TN QTL identified. Allele effects at the five TN QTL co-located with ShB QTL indicated that increased TN does not inevitably increase disease development; in fact, for four ShB QTL that overlapped TN QTL, the alleles increasing resistance were associated with increased TN and PN, suggesting a desirable coupling of alleles at linked genes.

CONCLUSIONS

Nineteen accessions identified as containing the most SNP alleles associated with ShB resistance for each subpopulation were resistant in both AR and NC field trials. Rice breeders can utilize these accessions and SNPs to develop cultivars with enhanced ShB resistance along with increased TN and PN for improved yield potential.

Enhancing rice grain production by manipulating the naturally evolved cis-regulatory element-containing inverted repeat sequence of OsREM20

Grain number per panicle (GNP) is an important agronomic trait that contributes to rice grain yield. Despite its importance in rice breeding, the molecular mechanism underlying GNP regulation remains largely unknown. In this study, we identified a previously unrecognized regulatory gene that controls GNP in rice, Oryza sativa REPRODUCTIVE MERISTEM 20 (OsREM20), which encodes a B3 domain transcription factor. Through genetic analysis and transgenic validation we found that genetic variation in the CArG box-containing inverted repeat (IR) sequence of the OsREM20 promoter alters its expression level and contributes to GNP variation among rice varieties. Furthermore, we revealed that the IR sequence regulates OsREM20 expression by affecting the direct binding of OsMADS34 to the CArG box within the IR sequence. Interestingly, the divergent pOsREM20^IR and pOsREM20^ΔIR alleles were found to originate from different Oryza rufipogon accessions, and were independently inherited into the japonica and indica subspecies, respectively, during domestication. Importantly, we demonstrated that IR sequence variations in the OsREM20 promoter can be utilized for germplasm improvement through either genome editing or traditional breeding. Taken together, our study characterizes novel genetic variations responsible for GNP diversity in rice, reveals the underlying molecular mechanism in the regulation of agronomically important gene expression, and provides a promising strategy for improving rice production by manipulating the cis-regulatory element-containing IR sequence.

Comparative Transcriptome Analysis of Rice Resistant and Susceptible Genotypes to Xanthomonas oryzae pv. oryzae Identifies Novel Genes to Control Bacterial Leaf Blight

The bacterial leaf blight in rice caused by Xanthomonas oryzae pv. oryzae (Xoo) affects crop losses worldwide. In spite of developing resistant varieties by introgressing different Xa genes, the occurrence of diseases is evident. Here we report identification of several genes that are associated with improved plant immunity against Xoo in a resistant genotype BPT-5204 in comparison with susceptible genotype TN-1. The RNA sequencing information was developed to identify the genes that could provide durable resistance in rice. Xoo-resistant rice genotype BPT-5204 with Xa 5, 13 and 21 genes is compared with sensitive Taichung Native 1 (TN-1) to identify the genetic pathways and gene networks involved in resistance mechanisms. The higher levels of salicylic acid resulted in upregulation of many pathogenesis-related (PR) and redox protein encoding transcripts which resulted in higher hypersensitive response in BPT-5204. Many Serine/threonine protein kinase, leucine-rich repeat (LRR) transmembrane protein kinase, protein kinase family genes, Wall-associated kinase (WAK) were upregulated that resulted in activation of bZIP, WRKY, MYB, DOF and HSFs transcription factors that are associated with improved plant immunity. The study provided roles of many genes and their associated plant immunity pathways that can be used for developing resistant rice cultivars.

Sheath blight of rice: a review and identification of priorities for future research

Rice sheath blight research should prioritise optimising biological control approaches, identification of resistance gene mechanisms and application in genetic improvement and smart farming for early disease detection. Rice sheath blight, caused by Rhizoctonia solani AG1-1A, is one of the most devasting diseases of the crop. To move forward with effective crop protection against sheath blight, it is important to review the published information related to pathogenicity and disease management and to determine areas of research that require deeper study. While progress has been made in the identification of pathogenesis-related genes both in rice and in the pathogen, the mechanisms remain unclear. Research related to disease management practices has addressed the use of agronomic practices, chemical control, biological control and genetic improvement: Optimising nitrogen fertiliser use in conjunction with plant spacing can reduce spread of infection while smart agriculture technologies such as crop monitoring with Unmanned Aerial Systems assist in early detection and management of sheath blight disease. Replacing older fungicides with natural fungicides and use of biological agents can provide effective sheath blight control, also minimising environmental impact. Genetic approaches that show promise for the control of sheath blight include treatment with exogenous dsRNA to silence pathogen gene expression, genome editing to develop rice lines with lower susceptibility to sheath blight and development of transgenic rice lines overexpressing or silencing pathogenesis related genes. The main challenges that were identified for effective crop protection against sheath blight are the adaptive flexibility of the pathogen, lack of resistant rice varieties, abscence of single resistance genes for use in breeding and low access of farmers to awareness programmes for optimal management practices.

Transcriptome analysis of a rice cultivar reveals the differentially expressed genes in response to wild and mutant strains of Xanthomonas oryzae pv. oryzae

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is a devastating disease in most of the rice growing regions worldwide. Among the 42 BB resistance (R) genes, Xa23 is an executor R gene, conferring broad-spectrum disease resistance to all naturally occurring biotypes of Xoo. In this study, CBB23, a rice line carrying Xa23 gene, was inoculated with wild PXO99^A and its mutant, P99M2, to retrieve the differentially expressed genes (DEGs). RNA-Seq analysis retrieved 1,235 DEGs (p-value ≤ 0.05) at 12, 24, 36, and 48 hours of post inoculation (hpi). Gene ontology (GO) analysis classified the DEGs functionally into biological process, cellular component and molecular function. KEGG pathway analysis categorized the DEGs into 11 different pathways, and the ribosome is a prominent pathway followed by biosynthesis of phenylpropanoids. Gene co-expression network analysis identified the clusters of transcription factors (TFs) which may be involved in PXO99^A resistance. Additionally, we retrieved 67 differentially expressed TFs and 26 peroxidase responsive genes which may be involved in disease resistance mechanism. DEGs involved in the host-pathogen interaction, e.g., signaling mechanism, cell wall and plant hormones were identified. This data would be a valuable resource for researchers to identify the candidate genes associated with Xoo resistance.

Gene network mediated by WRKY13 to regulate resistance against sheath infecting fungi in rice (Oryza sativa L.)

OsWRKY13 TF gene is known to play a regulatory role of signaling in physiological pathways related to either development or disease resistance in rice plants. Rice cultivars IR 50 and TRY 3, resistant and susceptible respectively to sheath blight, TRY 3 and CO 43 resistant and susceptible respectively to sheath rot were challenged with fungal pathogens and disease scoring was carried out. Percent Disease Index (PDI) was significantly higher in susceptible varieties than resistant varieties. RT-PCR and qPCR analyses of WRKY13 using RNA extracted from the plant tissues revealed higher WRKY13 expression in resistant varieties (both diseases) upon pathogen challenge compared to uninfected control and also the susceptible varieties. To compute and evaluate the possible molecular mechanism for observed resistance correlated to WRKY13 gene expression, rice gene expression profiles against bacterial leaf blight and leaf blast disease from ROAD database were used to prioritize the locus IDs that were used as input in RiceNet v2 tool. The expression of WRKY13-regulated TIFY9 gene was predicted and validated using RT-PCR and qRT-PCR along with WRKY12 and PR2. All three genes showed induced expression in R. solani challenged sheath blight resistant variety. WRKY12 and PR2 expression in S. oryzae challenged sheath rot resistant variety was higher. Agrobacterium mediated transformation was carried out in rice plants using overexpression construct of WRKY13 (agroinfection in seeds of varieties susceptible to sheath blight and sheath rot, followed by selection in antibiotic media, germinating and hardening of putative transgenic lines). Based on qPCR analysis, the expression level of WRKY13 and the co-expression levels of WRKY12, TIFY9 and PR2 were found higher in PCR-positive T₁ plants compared to wild-type. Infection bioassays in the transgenic plants of both varieties revealed enhanced resistance to the pathogens. A mechanism by which WRKY13 would influence the MAPK cascade with TIFY9 acting as a mediator, is proposed.

Comparative Transcriptome Profiling of Rice Near-Isogenic Line Carrying Xa23 under Infection of Xanthomonas oryzae pv. oryzae

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is an overwhelming disease in rice-growing regions worldwide. Our previous studies revealed that the executor R gene Xa23 confers broad-spectrum disease resistance to all naturally occurring biotypes of Xoo. In this study, comparative transcriptomic profiling of two near-isogenic lines (NILs), CBB23 (harboring Xa23) and JG30 (without Xa23), before and after infection of the Xoo strain, PXO99^A, was done by RNA sequencing, to identify genes associated with the resistance. After high throughput sequencing, 1645 differentially expressed genes (DEGs) were identified between CBB23 and JG30 at different time points. Gene Ontlogy (GO) analysis categorized the DEGs into biological process, molecular function, and cellular component. KEGG analysis categorized the DEGs into different pathways, and phenylpropanoid biosynthesis was the most prominent pathway, followed by biosynthesis of plant hormones, flavonoid biosynthesis, and glycolysis/gluconeogenesis. Further analysis led to the identification of differentially expressed transcription factors (TFs) and different kinase responsive genes in CBB23, than that in JG30. Besides TFs and kinase responsive genes, DEGs related to ethylene, jasmonic acid, and secondary metabolites were also identified in both genotypes after PXO99^A infection. The data of DEGs are a precious resource for further clarifying the network of Xa23-mediated resistance.