Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice.

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many QTLs and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five SNP variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity towards LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice, and suggest the potential of WLGhap.B in improving rice yield.

DEGAP: Dynamic elongation of a genome assembly path.

Genome assembly remains to be a major task in genomic research. Despite the development over the past decades of different assembly software programs and algorithms, it is still a great challenge to assemble a complete genome without any gaps. With the latest DNA circular consensus sequencing (CCS) technology, several assembly programs can now build a genome from raw sequencing data to contigs; however, some complex sequence regions remain as unresolved gaps. Here, we present a novel gap-filling software, DEGAP (Dynamic Elongation of a Genome Assembly Path), that resolves gap regions by utilizing the dual advantages of accuracy and length of high-fidelity (HiFi) reads. DEGAP identifies differences between reads and provides 'GapFiller' or 'CtgLinker' modes to eliminate or shorten gaps in genomes. DEGAP adopts an iterative elongation strategy that automatically and dynamically adjusts parameters according to three complexity factors affecting the genome to determine the optimal extension path. DEGAP has already been successfully applied to decipher complex genomic regions in several projects and may be widely employed to generate more gap-free genomes.

Fine Mapping of Five Grain Size QTLs Which Affect Grain Yield and Quality in Rice.

Grain size is a quantitative trait with a complex genetic mechanism, characterized by the combination of grain length (GL), grain width (GW), length to width ration (LWR), and grain thickness (GT). In this study, we conducted quantitative trait loci (QTL) analysis to investigate the genetic basis of grain size using BC1F2 and BC1F2:3 populations derived from two indica lines, Guangzhan 63-4S (GZ63-4S) and TGMS29 (core germplasm number W240). A total of twenty-four QTLs for grain size were identified, among which, three QTLs (qGW1, qGW7, and qGW12) controlling GL and two QTLs (qGW5 and qGL9) controlling GW were validated and subsequently fine mapped to regions ranging from 128 kb to 624 kb. Scanning electron microscopic (SEM) analysis and expression analysis revealed that qGW7 influences cell expansion, while qGL9 affects cell division. Conversely, qGW1, qGW5, and qGW12 promoted both cell division and expansion. Furthermore, negative correlations were observed between grain yield and quality for both qGW7 and qGW12. Nevertheless, qGW5 exhibited the potential to enhance quality without compromising yield. Importantly, we identified two promising QTLs, qGW1 and qGL9, which simultaneously improved both grain yield and quality. In summary, our results laid the foundation for cloning these five QTLs and provided valuable resources for breeding rice varieties with high yield and superior quality.

DNA methylation remodeling and the functional implication during male gametogenesis in rice.

BACKGROUND: Epigenetic marks are reprogrammed during sexual reproduction. In flowering plants, DNA methylation is only partially remodeled in the gametes and the zygote. However, the timing and functional significance of the remodeling during plant gametogenesis remain obscure. RESULTS: Here we show that DNA methylation remodeling starts after male meiosis in rice, with non-CG methylation, particularly at CHG sites, being first enhanced in the microspore and subsequently decreased in sperm. Functional analysis of rice CHG methyltransferase genes CMT3a and CMT3b indicates that CMT3a functions as the major CHG methyltransferase in rice meiocyte, while CMT3b is responsible for the increase of CHG methylation in microspore. The function of the two histone demethylases JMJ706 and JMJ707 that remove H3K9me2 may contribute to the decreased CHG methylation in sperm. During male gametogenesis CMT3a mainly silences TE and TE-related genes while CMT3b is required for repression of genes encoding factors involved in transcriptional and translational activities. In addition, CMT3b functions to repress zygotic gene expression in egg and participates in establishing the zygotic epigenome upon fertilization. CONCLUSION: Collectively, the results indicate that DNA methylation is dynamically remodeled during male gametogenesis, distinguish the function of CMT3a and CMT3b in sex cells, and underpin the functional significance of DNA methylation remodeling during rice reproduction.

A double-stranded RNA binding protein enhances drought resistance via protein phase separation in rice.

Drought stress significantly impacts global rice production, highlighting the critical need to understand the genetic basis of drought resistance in rice. Here, through a genome-wide association study, we reveal that natural variations in DROUGHT RESISTANCE GENE 9 (DRG9), encoding a double-stranded RNA (dsRNA) binding protein, contribute to drought resistance. Under drought stress, DRG9 condenses into stress granules (SGs) through liquid-liquid phase separation via a crucial α-helix. DRG9 recruits the mRNAs of OsNCED4, a key gene for the biosynthesis of abscisic acid, into SGs and protects them from degradation. In drought-resistant DRG9 allele, natural variations in the coding region, causing an amino acid substitution (G267F) within the zinc finger domain, increase DRG9's binding ability to OsNCED4 mRNA and enhance drought resistance. Introgression of the drought-resistant DRG9 allele into the elite rice Huanghuazhan significantly improves its drought resistance. Thus, our study underscores the role of a dsRNA-binding protein in drought resistance and its promising value in breeding drought-resistant rice.

Common and specific genetic basis of metabolite-mediated drought responses in rice.

Plants orchestrate drought responses at metabolic level but the genetic basis remains elusive in rice. In this study, 233 drought-responsive metabolites (DRMs) were quantified in a large rice population comprised of 510 diverse accessions at the reproductive stage. Large metabolic variations in drought responses were detected, and little correlation of metabolic levels between drought and normal conditions were observed. Interestingly, most of these DRMs could predict drought resistance in high accuracy. Genome-wide association study revealed 2522 significant association signals for 233 DRMs, and 98% (2471/2522) of the signals were co-localized with the association loci for drought-related phenotypic traits in the same population or the linkage-mapped QTLs for drought resistance in other populations. Totally, 10 candidate genes were efficiently identified for nine DRMs, seven of which harbored cis-eQTLs under drought condition. Based on comparative GWAS of common DRMs in rice and maize, representing irrigated and upland crops, we have identified three pairs of homologous genes associated with three DRMs between the two crops. Among the homologous genes, a transferase gene responsible for metabolic variation of N-feruloylputrescine was confirmed to confer enhanced drought resistance in rice. Our study provides not only genetic architecture of metabolic responses to drought stress in rice but also metabolic data resources to reveal the common and specific metabolite-mediated drought responses in different crops.

Paternal DNA methylation is remodeled to maternal levels in rice zygote.

Epigenetic reprogramming occurs during reproduction to reset the genome for early development. In flowering plants, mechanistic details of parental methylation remodeling in zygote remain elusive. Here we analyze allele-specific DNA methylation in rice hybrid zygotes and during early embryo development and show that paternal DNA methylation is predominantly remodeled to match maternal allelic levels upon fertilization, which persists after the first zygotic division. The DNA methylation remodeling pattern supports the predominantly maternal-biased gene expression during zygotic genome activation (ZGA) in rice. However, parental allelic-specific methylations are reestablished at the globular embryo stage and associate with allelic-specific histone modification patterns in hybrids. These results reveal that paternal DNA methylation is remodeled to match the maternal pattern during zygotic genome reprogramming and suggest existence of a chromatin memory allowing parental allelic-specific methylation to be maintained in the hybrid.

Control of rice ratooning ability by a nucleoredoxin that inhibits histidine kinase dimerization to attenuate cytokinin signaling in axillary buds.

Rice ratooning, the fast outgrowth of dormant buds on stubble, is an important cropping practice in rice production. However, the low ratooning ability (RA) of most rice varieties restricts the application of this cost-efficient system, and the genetic basis of RA remains unknown. In this study, we dissected the genetic architecture of RA by a genome-wide association study in a natural rice population. Rice ratooning ability 3 (RRA3), encoding a hitherto not characterized nucleoredoxin involved in reduction of disulfide bonds, was identified as the causal gene of a major locus controlling RA. Overexpression of RRA3 in rice significantly accelerated leaf senescence and reduced RA, whereas knockout of RRA3 significantly delayed leaf senescence and increased RA and ratoon yield. We demonstrated that RRA3 interacts with Oryza sativa histidine kinase 4 (OHK4), a cytokinin receptor, and inhibits the dimerization of OHK4 through disulfide bond reduction. This inhibition ultimately led to decreased cytokinin signaling and reduced RA. In addition, variations in the RRA3 promoter were identified to be associated with RA. Introgression of a superior haplotype with weak expression of RRA3 into the elite rice variety Guichao 2 significantly increased RA and ratoon yield by 23.8%. Collectively, this study not only uncovers an undocumented regulatory mechanism of cytokinin signaling through de-dimerization of a histidine kinase receptor-but also provides an eximious gene with promising value for ratoon rice breeding.

DELLA-mediated gene repression is maintained by chromatin modification in rice.

DELLA proteins are master regulators of gibberellic acid (GA) signaling through their effects on gene expression. Enhanced DELLA accumulation in rice and wheat varieties has greatly contributed to grain yield increases during the green revolution. However, the molecular basis of DELLA-mediated gene repression remains elusive. In this work, we show that the rice DELLA protein SLENDER RICE1 (SLR1) forms a tripartite complex with Polycomb-repressive complex 2 (PRC2) and the histone deacetylase HDA702 to repress downstream genes by establishing a silent chromatin state. The slr1 mutation and GA signaling resulted in dissociation of PRC2 and HDA702 from GA-inducible genes. Loss-of-function or downregulation of the chromatin regulators impaired SLR1-dependent histone modification and gene repression. Time-resolved analysis of GA signaling revealed that GA-induced transcriptional activation was associated with a rapid increase of H3K9ac followed by H3K27me3 removal. Collectively, these results establish a general epigenetic mechanism for DELLA-mediated gene repression and reveal details of the chromatin dynamics during transcriptional activation stimulated by GA signaling.

Serine protease NAL1 exerts pleiotropic functions through degradation of TOPLESS-related corepressor in rice.

NARROW LEAF 1 (NAL1) is a breeding-valuable pleiotropic gene that affects multiple agronomic traits in rice, although the molecular mechanism is largely unclear. Here, we report that NAL1 is a serine protease and displays a novel hexameric structure consisting of two ATP-mediated doughnut-shaped trimeric complexes. Moreover, we identified TOPLESS-related corepressor OsTPR2 involved in multiple growth and development processes as the substrate of NAL1. We found that NAL1 degraded OsTPR2, thus modulating the expression of downstream genes related to hormone signalling pathways, eventually achieving its pleiotropic physiological function. An elite allele, NAL1A, which may have originated from wild rice, could increase grain yield. Furthermore, the NAL1 homologues in different crops have a similar pleiotropic function to NAL1. Our study uncovers a NAL1-OsTPR2 regulatory module and provides gene resources for the design of high-yield crops.

Genome-wide profiling of rice Double-stranded RNA-Binding Protein 1-associated RNAs by targeted RNA editing.

RNA-binding proteins (RBPs) play essential roles in regulating gene expression. However, the RNA ligands of RBPs are poorly understood in plants, not least due to the lack of efficient tools for genome-wide identification of RBP-bound RNAs. An RBP-fused adenosine deaminase acting on RNA (ADAR) can edit RBP-bound RNAs, which allows efficient identification of RNA ligands of RBPs in vivo. Here, we report the RNA editing activities of the ADAR deaminase domain (ADARdd) in plants. Protoplast experiments indicated that RBP-ADARdd fusions efficiently edited adenosines within 41 nucleotides (nt) of their binding sites. We then engineered ADARdd to profile the RNA ligands of rice (Oryza sativa) Double-stranded RNA-Binding Protein 1 (OsDRB1). Overexpressing the OsDRB1-ADARdd fusion protein in rice introduced thousands of A-to-G and T-to-C RNA‒DNA variants (RDVs). We developed a stringent bioinformatic approach to identify A-to-I RNA edits from RDVs, which removed 99.7% to 100% of background single-nucleotide variants in RNA-seq data. This pipeline identified a total of 1,798 high-confidence RNA editing (HiCE) sites, which marked 799 transcripts as OsDRB1-binding RNAs, from the leaf and root samples of OsDRB1-ADARdd-overexpressing plants. These HiCE sites were predominantly located in repetitive elements, 3'-UTRs, and introns. Small RNA sequencing also identified 191 A-to-I RNA edits in miRNAs and other sRNAs, confirming that OsDRB1 is involved in sRNA biogenesis or function. Our study presents a valuable tool for genome-wide profiling of RNA ligands of RBPs in plants and provides a global view of OsDRB1-binding RNAs.

The CCT transcriptional activator Ghd2 constantly delays the heading date by upregulating CO3 in rice.

CONSTANS, CO-like, and TOC1 (CCT) family genes play important roles in regulating heading date, which exerts a large impact on the regional and seasonal adaptation of rice. Previous studies have shown that Grain number, plant height, and heading date2 (Ghd2) exhibit a negative response to drought stress by directly upregulating Rubisco activase and exerting a negative effect on heading date. However, the target gene of Ghd2 regulating heading date is still unknown. In this study, CO3 is identified by analyzing ChIP-seq data. Ghd2 activates CO3 expression by binding to the CO3 promoter through its CCT domain. EMSA experiments show that the motif CCACTA in the CO3 promoter was recognized by Ghd2. A comparison of the heading dates among plants with CO3 knocked out or overexpressed and double mutants overexpressing Ghd2 with CO3 knocked out shows that CO3 negatively and constantly regulates flowering by repressing the transcription of Ehd1, Hd3a, and RFT1. In addition, the target genes of CO3 are explored via a comprehensive analysis of DAP-seq data and RNA-seq data. Taken together, these results suggest that Ghd2 directly binds to the downstream gene CO3, and the Ghd2-CO3 module constantly delays heading date via the Ehd1-mediated pathway.

GLW7.1, a Strong Functional Allele of Ghd7, Enhances Grain Size in Rice.

Grain size is a key determinant of both grain weight and grain quality. Here, we report the map-based cloning of a novel quantitative trait locus (QTL), GLW7.1 (Grain Length, Width and Weight 7.1), which encodes the CCT motif family protein, GHD7. The QTL is located in a 53 kb deletion fragment in the cultivar Jin23B, compared with the cultivar CR071. Scanning electron microscopy analysis and expression analysis revealed that GLW7.1 promotes the transcription of several cell division and expansion genes, further resulting in a larger cell size and increased cell number, and finally enhancing the grain size as well as grain weight. GLW7.1 could also increase endogenous GA content by up-regulating the expression of GA biosynthesis genes. Yeast two-hybrid assays and split firefly luciferase complementation assays revealed the interactions of GHD7 with seven grain-size-related proteins and the rice DELLA protein SLR1. Haplotype analysis and transcription activation assay revealed the effect of six amino acid substitutions on GHD7 activation activity. Additionally, the NIL with GLW7.1 showed reduced chalkiness and improved cooking and eating quality. These findings provide a new insight into the role of Ghd7 and confirm the great potential of the GLW7.1 allele in simultaneously improving grain yield and quality.

Field rice panicle detection and counting based on deep learning.

Panicle number is directly related to rice yield, so panicle detection and counting has always been one of the most important scientific research topics. Panicle counting is a challenging task due to many factors such as high density, high occlusion, and large variation in size, shape, posture et.al. Deep learning provides state-of-the-art performance in object detection and counting. Generally, the large images need to be resized to fit for the video memory. However, small panicles would be missed if the image size of the original field rice image is extremely large. In this paper, we proposed a rice panicle detection and counting method based on deep learning which was especially designed for detecting rice panicles in rice field images with large image size. Different object detectors were compared and YOLOv5 was selected with MAPE of 3.44% and accuracy of 92.77%. Specifically, we proposed a new method for removing repeated detections and proved that the method outperformed the existing NMS methods. The proposed method was proved to be robust and accurate for counting panicles in field rice images of different illumination, rice accessions, and image input size. Also, the proposed method performed well on UAV images. In addition, an open-access and user-friendly web portal was developed for rice researchers to use the proposed method conveniently.

Fine-tuning OsCPK18/OsCPK4 activity via genome editing of phosphorylation motif improves rice yield and immunity.

Plants have evolved complex signalling networks to regulate growth and defence responses under an ever-changing environment. However, the molecular mechanisms underlying the growth-defence tradeoff are largely unclear. We previously reported that rice CALCIUM-DEPENDENT PROTEIN KINASE 18 (OsCPK18) and MITOGEN-ACTIVATED PROTEIN KINASE 5 (OsMPK5) mutually phosphorylate each other and that OsCPK18 phosphorylates and positively regulates OsMPK5 to suppress rice immunity. In this study, we found that OsCPK18 and its paralog OsCPK4 positively regulate plant height and yield-related traits. Further analysis reveals that OsCPK18 and OsMPK5 synergistically regulate defence-related genes but differentially regulate development-related genes. In vitro and in vivo kinase assays demonstrated that OsMPK5 phosphorylates C-terminal threonine (T505) and serine (S512) residues of OsCPK18 and OsCPK4, respectively. The kinase activity of OsCPK18T505D , in which T505 was replaced by aspartic acid to mimic T505 phosphorylation, displayed less calcium sensitivity than that of wild-type OsCPK18. Interestingly, editing the MAPK phosphorylation motif in OsCPK18 and its paralog OsCPK4, which deprives OsMPK5-mediated phosphorylation but retains calcium-dependent activation of kinase activity, simultaneously increases rice yields and immunity. This editing event also changed the last seven amino acid residues of OsCPK18 and attenuated its binding with OsMPK5. This study presents a new regulatory circuit that fine tunes the growth-defence tradeoff by modulating OsCPK18/4 activity and suggests that CRISPR/Cas9-mediated engineering phosphorylation pathways could simultaneously improve crop yield and immunity.

A MITE variation-associated heat-inducible isoform of a heat-shock factor confers heat tolerance through regulation of JASMONATE ZIM-DOMAIN genes in rice.

High temperatures cause huge yield losses in rice. Heat-shock factors (Hsfs) are key transcription factors which regulate the expression of heat stress-responsive genes, but natural variation in and functional characterization of Hsfs have seldom been reported. A significant heat response locus was detected via a genome-wide association study (GWAS) using green leaf area as an indicative trait. A miniature inverted-repeat transposable element (MITE) in the promoter of a candidate gene, HTG3 (heat-tolerance gene on chromosome 3), was found to be significantly associated with heat-induced expression of HTG3 and heat tolerance (HT). The MITE-absent variant has been selected in heat-prone rice-growing regions. HTG3a is an alternatively spliced isoform encoding a functional Hsf, and experiments using overexpression and knockout rice lines showed that HTG3a positively regulates HT at both vegetative and reproductive stages. The HTG3-regulated genes were enriched for heat shock proteins and jasmonic acid signaling. Two heat-responsive JASMONATE ZIM-DOMAIN (JAZ) genes were confirmed to be directly upregulated by HTG3a, and one of them, OsJAZ9, positively regulates HT. We conclude that HTG3 plays an important role in HT through the regulation of JAZs and other heat-responsive genes. The MITE-absent allele may be valuable for HT breeding in rice.

A long transcript mutant of the rubisco activase gene RCA upregulated by the transcription factor Ghd2 enhances drought tolerance in rice.

The transcription factor Ghd2 increases rice yield potential under normal conditions and accelerates leaf senescence under drought stress. However, its mechanism on the regulation of leaf senescence under drought stress remains unclear. In the present study, to unveil the mechanism, one target of Ghd2, the Rubisco activase gene RCA, was identified through the combined analysis of Ghd2-CRISPR transcriptome data and Ghd2-overexpression microarray data. Ghd2 binds to the 'CACA' motif in the RCA promoter by its CCT domain and upregulates RCA expression. RCA has alternative transcripts, RCAS and RCAL, which are predominantly expressed under normal conditions and drought stress, respectively. Similar to Ghd2-overexpressing plants, RCAL-overexpressing plants were more sensitive to drought stress than the wild-type. However, the plants overexpressing RCAS showed a weak drought-sensitive phenotype. Moreover, RCAL knockdown and knockout plants did not show yield loss under normal conditions, but exhibited enhanced drought tolerance and delayed leaf senescence. The chlorophyll content, the free amino acid content and the expression of senescence-related genes in the RCAL mutant were lower than those in the wild-type plants under drought stress. In summary, Ghd2 induces leaf senescence by upregulating RCAL expression under drought stress, and the RCAL mutant has important values in breeding drought-tolerant varieties.

A FLASH pipeline for arrayed CRISPR library construction and the gene function discovery of rice receptor-like kinases.

Clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated gene editing is revolutionizing plant research and crop breeding. Here, we present an effective and streamlined pipeline for arrayed CRISPR library construction and demonstrate it is suitable for small- to large-scale genome editing in plants. This pipeline introduces artificial PCR fragment-length markers for distinguishing guide RNAs (gRNAs) (FLASH), and a group of 12 constructs harboring different FLASH tags are co-transformed into plants each time. The identities of gRNAs in Agrobacterium mixtures and transgenic plants can therefore be read out by detecting the FLASH tags, a process that requires only conventional PCR and gel electrophoresis rather than sequencing. We generated an arrayed CRISPR library targeting all 1,072 members of the receptor-like kinase (RLK) family in rice. One-shot transformation generated a mutant population that covers gRNAs targeting 955 RLKs, and 74.3% (710/955) of the target genes had three or more independent T0 lines. Our results indicate that the FLASH tags act as bona fide surrogates for the gRNAs and are tightly (92.1%) associated with frameshift mutations in the target genes. In addition, the FLASH pipeline allows for rapid identification of unintended editing events without corresponding T-DNA integrations and generates high-order mutants of closely related RLK genes. Furthermore, we showed that the RLK mutant library enables rapid discovery of defense-related RLK genes. This study introduces an effective pipeline for arrayed CRISPR library construction and provides genome-wide rice RLK mutant resources for functional genomics.