Abstract:Breeding high-yield crop cultivars has improved agronomic performance for key grain crops (Hirano et al., 2017). However traditional breeding takes years. Molecular genetic studies have identified beneficial trait-associated alleles in elite cultivars providing a platform for gene-editing techniques such as the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) system (Rothan et al., 2019). Most current 'super rice' varieties with high yields have several beneficial agronomic traits. For example the variety Liang-You-Pei-Jiu has strong culms for lodging resistance and large panicles for high yield. Our previous study using recombinant inbred lines (RILs) from a cross between 93-11 and PA64 identified 43 quantitative trait loci (QTLs) associated with many agronomic traits including heading date spikelet number per panicle and grain shape (Gao et al., 2013). However no QTL for culm strength or panicle size was detected. Here we studied culm strength by measuring stem cross-section area (SCSA) of the fourth internode (Figure 1a). By using 127 chromosome substitution segment lines (CSSLs) from a cross between 93-11 and Nipponbare (NPB) we detected three QTLs associated with SCSA (Figure 1b,c,d) (Xu et al., 2010). The position of qSCSA1-1 (Chr1: 42290679-43345963) overlapped with STRONG CULM 1 (SCM1) (Mulsanti et al., 2018) and SEMI DWARF-1 (SD-1) and qSCSA6-1 (Chr6: 28410389-28439277) overlapped with SCM2 (Figure 1b,d). For both QTLs the alleles from NPB positively contributed to SCSA. The region of qSCSA3-1 on the long arm of Chr3 (28406301-30831094) included SCM3/RICE TEOSINTE BRANCH1 (OsTB1)/FINE CULM1 (FC1) (Figure 1b,c,d) (Yano et al., 2015) and the 93-11 allele had a positive effect on culm strength. To test whether OsTB1 is the causal gene for qSCSA3-1 fine mapped qSCSA3-1 using near-isogenic lines (NILs). qSCSA3-1 mapped to a 51-kb region between markers CT11 and CT17 a region that includes OsTB1 (Figure 1e). Quantitative-PCR (qPCR) showed that the expression of OsTB1 in 93-11 inflorescences was higher than that of NPB and NIL-93-11 carrying the NPB allele NIL-qSCSA3-1NPB (Figure 1f). In addition to affecting SCSA introgression of the NPB allele resulted in a decrease in bending stress stem area and number of spikelets per panicle and an increase in tiller number but had no effect on plant height (Figure 1f). These observations indicated that OsTB1 has pleiotropic effects on different traits. Therefore we examined the effect of OsTB1 using mutant and complemented plants (Figure 1g). The null mutant of OsTB1 fc1-2 showed lower bending stress smaller panicles and increased tiller number compared with rescued plants expressing OsTB1 from 93-11 or Kasalath (Figure 1h). These observations confirmed that qSCSA3-1 is allelic to SCM3/OsTB1 and has pleiotropic effects. We compared the sequence of OsTB1 among rice cultivars NPB (japonica) Koshihikari (japonica) and four varieties with large SCSA: 93-11 (indica) Kasalath (indica) Zhongchao123 (japonica) and Chugoku117 (indica) and found polymorphisms in the 5'-flanking coding sequence (CDS) (Figure 1i,j). Interestingly the four varieties with large SCSA contained a TGTG insertion at + 219 in the 5'-noncoding region suggesting that this insertion might be important for OsTB1 expression and consequently affect SCSA. Based on this prediction we introduced mutations in the promoter and 5'-noncoding regions of SCM3/OsTB1 in NPB using CRISPR-Cas9. We designed six single-guide RNAs (sgRNAs) four targeting the promoter sequence and two sgRNA5 and 6 targeting the region proximal to the TGTG insertion in the 5'-flanking CDS. These six sgRNAs were integrated into one plasmid and transformed into NPB using Agrobacterium. We obtained 23 plants carrying 9 different mutations (Figure 1j). Relative to wild-type (WT) NPB Type 1 had more tillers and smaller culms and panicles; Type 2 had similar phenotypes to WT; Type 3 had fewer tillers and larger culms and panicles (Figure 1k). The Type 1 plants 67-1-1 63-1-6 71-3-2 and 66-2-3 were similar to the null mutant of OsTB1 fc1-2 (Minakuchi et al., 2010). Indeed 71-3-2 and 66-2-3 contained deletions in the 5'-noncoding region of OsTB1. OsSPL14/IPA1 positively regulates OsTB1 expression through binding to a GTAC motif (Lu et al., 2013) and this motif was deleted in these plants likely resulting in the null phenotypes. The Type 3 plants 83-4-1 and 76-1-2 mimicked the phenotypes of NILSCM3 (Yano et al., 2015) and 93-11 suggesting that mutations in these plants enhance OsTB1 expression. Indeed these plants showed higher OsTB1 expression than the control whereas the Type 1 and Type 2 plants showed reduced and unchanged expression levels respectively (Figure 1l). 76-1-2 had a TGTG insertion which also occurs in varieties with large SCSA (Figure 1i j) confirming that the TGTG insertion enhances OsTB1 expression. Although 83-4-1 had mutations at six sites in the 5'-flanking CDS region it showed a gain-of-function phenotype (Figure 1k). Thus different sequence(s) in the 5'-flanking CDS region in WT may repress or enhance transcription. In 84-5-1 a 36-bp deletion around the TGTG insertion did not increase OsTB1 expression suggesting that the TGTG insertion did not disrupt a repressive site. Instead the TGTG insertion may create a binding site for transcription activator(s) hypotheses that will require further validation. For Type 2 plants which showed similar phenotypes to the control plants (Figure 1k) OsTB1 expression did not differ from the control although they had various deletions and single-nucleotide polymorphisms upstream of the gene (Figure 1j). These results demonstrated that targeted editing of OsTB1 cis-regulatory elements could produce alleles having different expression levels and phenotypes (Figure 1l m). The application of CRISPR-Cas9 for studying cis-elements has been discussed but a detailed strategy has not been established in rice (Rodriguez-Leal et al., 2017). In this study we introduced mutations in cis-regulatory sequences which led to alterations in gene expression and phenotypes and reproduced beneficial alleles in new rice varieties. QTL studies have identified beneficial alleles in rice elite germ plasm including alleles affecting gene expression. For example the molecular mechanism for major QTLs controlling plant structure such as GN1A IPA1/WFP NUMBER OF GRAINS 1 (NOG1) FZP1/qSRN7 and also SCM2/APO1 as mentioned above depends on the difference in the expression level of these genes (Huo et al., 2017). These spontaneous alleles have been introgressed into elite varieties by breeding which is time-consuming and labour-intensive. CRISPR-Cas9 however easily produced alleles with similar expression levels as shown for OsTB1. Furthermore we may be able to isolate desirable alleles with more appropriate levels of expression for optimizing breeding targets. In this regard saturation editing of known genes associated with agronomic traits of interest may accelerate molecular design of desirable traits. This research was financially supported by the Natural Science Foundation of China (Grant No: 31461143014). Xingming Hu, Qian Qian and Makoto Matsuoka conceived the project and designed the research strategies. Yongtao Cui, Xingming Hu, Lan Shen, Fanmiao Wang, Shuang Ruan and Anhui Feng cloned the CRISPR-Cas9 plasmid, identified mutants and performed qRT-PCR. Guohuang Liang, Bin Zhang, Yongjun Lin, Dongdong Chen, Li Zhu and Longbiao Guo performed the SASA QTL analysis and mapping. Qian Qian, Xingming Hu, Guojun Dong and Jiang Hu performed phenotype observations and genetic tests. Xingming Hu, Makoto Matsuoka and Qian Qian wrote the article. The authors declare no competing interests. Figure S1 Mapping of SCSA (Stem Cross-Section Area) QTLs CSSLs. Figure S2 Identification of hygromycin in different editedplants. Table S1 Primer sequences used in this study. Data S1 Materials and Methods. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

Abstract:Rice blast, caused by Magnaporthe oryzae, is a serious disease in rice (Oryzae sativa L.) worldwide. Because the pathogen mutates rapidly, rice blast has become one of the major threats for global rice production (Wang and Valent, 2017). Using host resistance is the most effective method for blast control. However most rice resistance genes are race-specific (Li et al., 2017). Discovering race non-specific resistance genes will enhance blast resistance in rice cultivars. Fusarium graminearum is a hemibiotrophic fungal pathogen for wheat Fusarium head blight (FHB) and shares a similar lifestyle with Magnaporthe oryzae for rice blast. Previously, we have cloned Fhb1 encoding a histidine-rich calcium-binding protein (TaHRC) from wheat for race non-specific resistance to Fusarium species, which has been widely deployed in many wheat cultivars worldwide (Bai et al., 2018; Su et al., 2019). TaHRC serves as a susceptibility factor to promote FHB spread and loss-of-function mutation in TaHRC confers FHB resistance (Chen et al., 2022). HRC is highly conservative with the closest TaHRC homologue from rice (OsHRC), however, the functions of OsHRC in rice blast resistance remains to be investigated. We cloned OsHRC (LOC_Os01g03060) from ‘Zhonghua11' (ZH11) and the OsHRC DNA sequence is 2390 bp with three exons and two introns. OsHRC encodes a putative 33.5 kDa protein of 285-amino acids and shares 63.5%–64.5% amino acid identity with TaHRC from three wheat subgenomes (A, B and D). The TaHRC and OsHRC proteins were predicted to have similar biological features and both contain a histidine and arginine enriched region without any known functional domain except for a motif of a nuclear localization signal (NLS), indicating the conserved HRC function between rice and wheat. The OsHRC protein was localized to the nucleus in rice leaf protoplasts (Figure 1a). Transactivation activity assay showed that yeast transiently expressing BD-OsHRC and BD-TaHRC did not grow in a deficient medium (SD/-Trp-His) (Figure 1b), suggesting that both nuclear proteins do not have transactivational activity and share the same biological functions. To verify the function of OsHRC on rice blast resistance, we designed a guide RNA to target a 20-nt sequence in the 59 bp downstream of the translation start codon of OsHRC and conducted CRISPR/Cas9-mediated gene editing to knock out (KO) OsHRC in ZH11. Among the five OsHRC KO mutants generated, three (KO-1, KO-2 and KO-3) showed premature translation termination due to frameshifts, and were predicted to generate truncated proteins of 39, 236, and 39 amino acids, respectively (Figure 1c). The three mutants and ZH11 were inoculated with RB22, a highly virulent M. oryzae strain from China, 60 days after planting. All the three mutants showed significantly smaller lesions and less fungal biomass than their non-edited controls at 14 days post inoculation (DPI) (Figure 1d–f). To examine whether the resistance is race non-specific, the three mutants were inoculated with another virulant strain S005 and showed the same level of resistance as to RB22 (Figure S1a,b), demonstrating that the loss-of-function of OsHRC confers race non-specific resistance to M. oryzae. To explore the molecular basis of OsHRC resistance to rice blast, RNA-seq was conducted to identify differentially expressed genes (DEGs) between an OsHRC KO mutant and ZH11. At 14 DPI with RB22, a total of 1333 DEGs were detected at the fold-change of 2 and false discovery rate (FDR) < 0.05 using the R package DESeq2. KEGG analysis indicated that the DEGs were mainly enriched in the pathways of photosynthesis, biosynthesis of secondary metabolites, and phenylpropanoid biosynthesis (Figure 1g). Phenylpropanoids play important roles in plant responses to biotic and abiotic stresses (Li et al., 2020). Quantitative RT-PCR results demonstrated that the genes in the phenylpropanoid biosynthesis pathway including genes encoding two phenylalanine ammonia-lyase (PAL6, Os04g0518400; PAL4, Os02g0627100), a 4-coumarate-CoA ligase (4CL, Os01g0901600), a cinnamyl-CoA reductase (CCR, Os02g0180700) and a cinnamyl-alcohol dehydrogenase (CAD, Os04g0229100) were significantly upregulated in the three OsHRC mutants compared with ZH11 (Figures 1hand S2a), consistent with Vanholme et al. (2019). In addition, the transcription factors, MYB30 and MYB55 involved in the regulation of the phenylpropanoid pathway (Li et al., 2020) were upregulated in the three OsHRC mutants (Figures 1h and S2a). Lignin is one of the most important phenylpropanoid metabolites, and its accumulation in leaf inhibits the penetration of M. oryzae (Li et al., 2020; Zhou et al., 2018). Since all the seven genes are pivotal for the lignin biosynthesis (Vanholme et al., 2019), we compared the lignin content between all mutants and the non-edited control, and found that the leaves of three mutants contained significantly higher lignin content than the control at 14 DPI with the isolate RB22 (Figures 1i and S2b), suggesting that OsHRC-mediated blast resistance is most likely due to the activation of the lignin synthesis pathway after M. oryzae infection. HRC is a conserved gene in cereal crops and the loss-of-function mutation in HRC confers resistance to wheat FHB (Su et al., 2019) and maize ear rot without significant penalty on yield (Liu et al., 2022). To evaluate effects of OsHRC on major rice agronomic traits, we compared heading date, plant height, and grain size between the mutants and WT control in greenhouse experiments and found that the mutants were slightly shorter, but with larger spikes, more spikelets per spike, and earlier heading date than the control (Figure 1j,k). These data indicate that OsHRC mutants showed significant improvement in blast resistance without obvious yield reduction in greenhouses, therefore, can be used as a new source of blast resistance in rice breeding. This work was supported by the National Key Research and Development Program of China (2022YFD1201502), Talent Funds of China Agricultural University (2021RC009) and US Wheat and Barley Scab Initiative. Authors thank Dr. H. Fang from Chinese Academy of Agricultural Sciences (CAAS) for technical help with rice blast phenotyping and Prof. Y. Ning from CAAS and Prof. Y. Peng from China Agricultural University for providing M. oryzae isolates. USDA is an equal opportunity provider and employer. ZS and GB designed the project and wrote the manuscript. ZS, YD, FZ, ZZ, XH, FS, JL performed experiments. All authors revised and approved the manuscript. The authors declare no conflicts of interest. Data S1 Material and methods. Figure S1 Comparison of OsHRC of rice for responses to the inoculation of M. oryzae isolate S005. (a) Lesions caused by M. oryzae on the non-edited control (WT) and three knockout (KO) mutant lines (KO-1, KO-2 and KO-3). (b) Relative lesion areas (means ± standard errors) in the three mutants were significantly smaller than in non-edited control at 14 days post inoculation (DPI) with the M. oryzae isolate S005. (c) Relative fungal biomass (means ± standard errors) in the three mutants was significantly lower than in non-edited ZH11 control at 14 DPI. Figure S2 Comparisons of transcript levels of lignin-related genes and lignin content in the leaves between ZH11 and the mutants after different time points of inoculation with the M. oryzae isolates RB22. (a) Comparisons of seven differentially expressed lignin-related genes between three mutants and ZH11 (WT control) at 0, 5, 10 and 14 days post inoculation (dpi). Data were normalized using the OsActin gene control. Error bars represent the standard error (SE) of three mutants. (b) Comparisons of lignin content between the inoculated and uninoculated leaves of the three mutants and ZH11 (WT control) at 14 DPI with isolate RB22. CK represent uninoculated control. Treatment represents 14 days post inoculation with isolate RB22. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

An Expanded Cysteine-Rich Receptor-Like Kinase Gene Cluster Functionally Differentiates in Drought, Cold, Heat, and Pathogen Stress Responses in Rice

Rice Gene,osckx2-2, Regulates Inflorescence and Grain Size by Increasing Endogenous Cytokinin Content

Loss-function mutants of OsCKX gene family based on CRISPR-Cas systems revealed their diversified roles in rice

A Rice Calcium-Dependent Protein Kinase OsCPK9 Positively Regulates Drought Stress Tolerance and Spikelet Fertility

Cytokinin oxidase/dehydrogenase family genes exhibit functional divergence and overlap in ricegrowth and development, especially in control of tillering

Identification, Characterization, and Mutational Analysis of a Probable KEAP1 Ortholog in Rice (oryza Sativa L.).

[Identification and expression profile analysis of rice CRF gene family]

Production of Novel Beneficial Alleles of a Rice Yield‐related QTL by CRISPR/Cas9

Genome-wide identification of CRF gene family members in four rice subspecies and expression analysis of OsCRF members in response to cold stress at seedling stage

OsWRKY70 Plays Opposite Roles in Blast Resistance and Cold Stress Tolerance in Rice

Loss of OsHRC function confers blast resistance without yield penalty in rice.

Overexpression of OsRRK1 Changes Leaf Morphology and Defense to Insect in Rice

Combined Gene Family Characterization and RNA-Seq to Study the Response of Β-Ketoacyl-coa Synthase to Abiotic Stress in Rice (oryza Sativa L.)

Overexpression of the receptor-like kinase gene OsNRRB enhances drought-stress tolerance in rice

The OsWRKY63-OsWRKY76-OsDREB1B Module Regulates Chilling Tolerance in Rice.

Cytokinin oxidase2-deficient mutants improve panicle and grain architecture through cytokinin accumulation and enhance drought tolerance in indica rice

Integrative transcriptomic analysis deciphering the role of rice bHLH transcription factor Os04g0301500 in mediating responses to biotic and abiotic stresses

Isolation of a Novel Leucine-Rich Repeat Receptor-Like Kinase (oslrr2) Gene from Rice and Analysis of Its Relation to Abiotic Stress Responses

The Suppressed Expression of a Stress Responsive Gene 'osdsr2' Enhances Rice Tolerance in Drought and Salt Stress.

OsCKq1 Regulates Heading Date and Grain Weight in Rice in Response to Day Length

Genome-Wide Identification and Expression Analysis of RCC1 Gene Family under Abiotic Stresses in Rice (Oryza sativa L.)