Population-level exploration of alternative splicing and its unique role in controlling agronomic traits of rice.

Alternative splicing (AS) plays crucial roles in regulating various biological processes in plants. However, the genetic mechanisms underlying AS and its role in controlling important agronomic traits in rice (Oryza sativa) remain poorly understood. In this study, we explored AS in rice leaves and panicles using the rice minicore collection. Our analysis revealed a high level of transcript isoform diversity, with approximately one fifth of potential isoforms acting as major transcripts in both tissues. Regarding the genetic mechanism of AS, we found that the splicing of 833 genes in the leaf and 1,230 genes in the panicle was affected by cis-genetic variation. Twenty-one percent of these AS events could only be explained by large structural variations. Approximately 77.5% of genes with significant splicing quantitative trait loci (sGenes) exhibited tissue-specific regulation, and AS can cause 26.9% (leaf) and 23.6% (panicle) of sGenes to have altered, lost or gained functional domains. Additionally, through splicing-phenotype association analysis, we identified phosphate-starvation induced RING-type E3 ligase (OsPIE1; LOC_Os01g72480), whose splicing ratio was significantly associated with plant height. In summary, this study provides an understanding of AS in rice and its contribution to the regulation of important agronomic traits.

mRNA cleavage by 21-nucleotide phasiRNAs determines temperature-sensitive male sterility in rice.

Temperature-sensitive male sterility is one of the core components for hybrid rice (Oryza sativa) breeding based on the 2-line system. We previously found that knockout of ARGONAUTE 1d (AGO1d) causes temperature-sensitive male sterility in rice by influencing phased small interfering RNA (phasiRNA) biogenesis and function. However, the specific phasiRNAs and their targets underlying the temperature-sensitive male sterility in the ago1d mutant remain unknown. Here, we demonstrate that the ago1d mutant displays normal female fertility but complete male sterility at low temperature. Through a multiomics analysis of small RNA (sRNA), degradome, and transcriptome, we found that 21-nt phasiRNAs account for the greatest proportion of the 21-nt sRNA species in rice anthers and are sensitive to low temperature and markedly downregulated in the ago1d mutant. Moreover, we found that 21-nt phasiRNAs are essential for the mRNA cleavage of a set of fertility- and cold tolerance-associated genes, such as Earlier Degraded Tapetum 1 (EDT1), Tapetum Degeneration Retardation (TDR), OsPCF5, and OsTCP21, directly or indirectly determined by AGO1d-mediated gene silencing. The loss of function of 21-nt phasiRNAs can result in upregulation of their targets and causes varying degrees of defects in male fertility and grain setting. Our results highlight the essential functions of 21-nt phasiRNAs in temperature-sensitive male sterility in rice and suggest their promising application in 2-line hybrid rice breeding in the future.

A rice variation map derived from 10 548 rice accessions reveals the importance of rare variants.

Detailed knowledge of the genetic variations in diverse crop populations forms the basis for genetic crop improvement and gene functional studies. In the present study, we analyzed a large rice population with a total of 10 548 accessions to construct a rice super-population variation map (RSPVM), consisting of 54 378 986 single nucleotide polymorphisms, 11 119 947 insertion/deletion mutations and 184 736 presence/absence variations. Assessment of variation detection efficiency for different population sizes revealed a sharp increase of all types of variation as the population size increased and a gradual saturation of that after the population size reached 10 000. Variant frequency analysis indicated that ∼90% of the obtained variants were rare, and would therefore likely be difficult to detect in a relatively small population. Among the rare variants, only 2.7% were predicted to be deleterious. Population structure, genetic diversity and gene functional polymorphism of this large population were evaluated based on different subsets of RSPVM, demonstrating the great potential of RSPVM for use in downstream applications. Our study provides both a rich genetic basis for understanding natural rice variations and a powerful tool for exploiting great potential of rare variants in future rice research, including population genetics and functional genomics.

Programmable RNA N6-methyladenosine editing with CRISPR/dCas13a in plants.

N6-methyladenonsine (m6A) is the most prevalent internal modification of messenger RNA (mRNA) and plays critical roles in mRNA processing and metabolism. However, perturbation of individual m6A modification to reveal its function and the phenotypic effects is still lacking in plants. Here, we describe the construction and characterization of programmable m6A editing tools by fusing the m6A writers, the core catalytic domain of the MTA and MTB complex, and the AlkB homologue 5 (ALKBH5) eraser, to catalytically dead Cas13a (dCas13a) to edit individual m6A sites on mRNAs. We demonstrated that our m6A editors could efficiently and specifically deposit and remove m6A modifications on specific RNA transcripts in both Nicotiana benthamiana and Arabidopsis thaliana. Moreover, we found that targeting SHORT-ROOT (SHR) transcripts with a methylation editor could significantly increase its m6A levels with limited off-target effects and promote its degradation. This leads to a boost in plant growth with enlarged leaves and roots, increased plant height, plant biomass, and total grain weight in Arabidopsis. Collectively, these findings suggest that our programmable m6A editing tools can be applied to study the functions of individual m6A modifications in plants, and may also have potential applications for future crop improvement.

A centromere map based on super pan-genome highlights the structure and function of rice centromeres.

Rice (Oryza sativa) is a significant crop worldwide with a genome shaped by various evolutionary factors. Rice centromeres are crucial for chromosome segregation, and contain some unreported genes. Due to the diverse and complex centromere region, a comprehensive understanding of rice centromere structure and function at the population level is needed. We constructed a high-quality centromere map based on the rice super pan-genome consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. We showed that rice centromeres have diverse satellite repeat CentO, which vary across chromosomes and subpopulations, reflecting their distinct evolutionary patterns. We also revealed that long terminal repeats (LTRs), especially young Gypsy-type LTRs, are abundant in the peripheral CentO-enriched regions and drive rice centromere expansion and evolution. Furthermore, high-quality genome assembly and complete telomere-to-telomere (T2T) reference genome enable us to obtain more centromeric genome information despite mapping and cloning of centromere genes being challenging. We investigated the association between structural variations and gene expression in the rice centromere. A centromere gene, OsMAB, which positively regulates rice tiller number, was further confirmed by expression quantitative trait loci, haplotype analysis and clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 methods. By revealing the new insights into the evolutionary patterns and biological roles of rice centromeres, our finding will facilitate future research on centromere biology and crop improvement.

Identification of natural allelic variation in TTL1 controlling thermotolerance and grain size by a rice super pan-genome.

Continuously increasing global temperatures present great challenges to food security. Grain size, one of the critical components determining grain yield in rice (Oryza sativa L.), is a prime target for genetic breeding. Thus, there is an immediate need for genetic improvement in rice to maintain grain yield under heat stress. However, quantitative trait loci (QTLs) endowing heat stress tolerance and grain size in rice are extremely rare. Here, we identified a novel negative regulator with pleiotropic effects, Thermo-Tolerance and grain Length 1 (TTL1), from the super pan-genomic and transcriptomic data. Loss-of-function mutations in TTL1 enhanced heat tolerance, and caused an increase in grain size by coordinating cell expansion and proliferation. TTL1 was shown to function as a transcriptional regulator and localized to the nucleus and cell membrane. Furthermore, haplotype analysis showed that hapL and hapS of TTL1 were obviously correlated with variations of thermotolerance and grain size in a core collection of cultivars. Genome evolution analysis of available rice germplasms suggested that TTL1 was selected during domestication of the indica and japonica rice subspecies, but still had much breeding potential for increasing grain length and thermotolerance. These findings provide insights into TTL1 as a novel potential target for the development of high-yield and thermotolerant rice varieties.

Temperature-sensitive male sterility in rice determined by the roles of AGO1d in reproductive phasiRNA biogenesis and function.

Phased secondary siRNAs (phasiRNAs) are broadly present in the reproductive tissues of flowering plants, with spatial-temporal specificity. However, the ARGONAUTE (AGO) proteins associated with phasiRNAs and their miRNA triggers remain elusive. Here, through histological and high-throughput sequencing analyses, we show that rice AGO1d, which is specifically expressed in anther wall cells before and during meiosis, associates with both miR2118 and miR2275 to mediate phasiRNA biogenesis. AGO1d preferentially binds to miR2118-triggered 21-nucleotide (nt) phasiRNAs with a 5'-terminal uridine, suggesting a dual role in phasiRNA biogenesis and function. Depletion of AGO1d causes a reduction of 21- and 24-nt phasiRNAs and temperature-sensitive male sterility. At lower temperatures, anthers of the ago1d mutant predominantly show excessive tapetal cells with little starch accumulation during pollen formation, possibly caused by the dysregulation of cell metabolism. These results uncover an essential role of AGO1d in rice anther development at lower temperatures and demonstrate coordinative roles of AGO proteins during reproductive phasiRNA biogenesis and function.

A quantitative trait locus GW6 controls rice grain size and yield through the gibberellin pathway.

Grain size is one of the essential components determining rice yield and is a target for both domestication and artificial breeding. Gibberellins (GAs) are diterpenoid phytohormones that influence diverse aspects of plant growth and development. Several quantitative trait loci (QTLs) have been identified that control grain size through phytohormone regulation. However, little is known about the role of GAs in the control of grain size. Here we report the cloning and characterization of a QTL, GW6 (GRAIN WIDTH 6), which encodes a GA-regulated GAST family protein and positively regulates grain width and weight. GW6 is highly expressed in the young panicle and increases grain width by promoting cell expansion in the spikelet hull. Knockout of GW6 exhibits reduced grain size and weight, whereas overexpression of GW6 results in increased grain size and weight. GW6 is induced by GA and its knockout downregulates the expression of GA biosynthesis genes and decreases GA content in the young panicle. We found that a natural variation in the cis element CAAT-box in the promoter of GW6 is associated with its expression level and grain width and weight. Furthermore, introduction of GW6 to Oryza indica variety HJX74 can lead to a 10.44% increase in rice grain yield, indicating that GW6 has great potential to improve grain yield in rice.

GRAIN SIZE AND NUMBER1 Negatively Regulates the OsMKKK10-OsMKK4-OsMPK6 Cascade to Coordinate the Trade-off between Grain Number per Panicle and Grain Size in Rice.

Grain number and size are interactive agronomic traits that determine grain yield. However, the molecular mechanisms responsible for coordinating the trade-off between these traits remain elusive. Here, we characterized the rice (Oryza sativa) grain size and number1 (gsn1) mutant, which has larger grains but sparser panicles than the wild type due to disordered localized cell differentiation and proliferation. GSN1 encodes the mitogen-activated protein kinase phosphatase OsMKP1, a dual-specificity phosphatase of unknown function. Reduced expression of GSN1 resulted in larger and fewer grains, whereas increased expression resulted in more grains but reduced grain size. GSN1 directly interacts with and inactivates the mitogen-activated protein kinase OsMPK6 via dephosphorylation. Consistent with this finding, the suppression of mitogen-activated protein kinase genes OsMPK6, OsMKK4, and OsMKKK10 separately resulted in denser panicles and smaller grains, which rescued the mutant gsn1 phenotypes. Therefore, OsMKKK10-OsMKK4-OsMPK6 participates in panicle morphogenesis and acts on a common pathway in rice. We confirmed that GSN1 is a negative regulator of the OsMKKK10-OsMKK4-OsMPK6 cascade that determines panicle architecture. The GSN1-MAPK module coordinates the trade-off between grain number and grain size by integrating localized cell differentiation and proliferation. These findings provide important insights into the developmental plasticity of the panicle and a potential means to improve crop yields.

NAL8 encodes a prohibitin that contributes to leaf and spikelet development by regulating mitochondria and chloroplasts stability in rice.

BACKGROUND: Leaf morphology and spikelet number are two important traits associated with grain yield. To understand how genes coordinating with sink and sources of cereal crops is important for grain yield improvement guidance. Although many researches focus on leaf morphology or grain number in rice, the regulating molecular mechanisms are still unclear. RESULTS: In this study, we identified a prohibitin complex 2α subunit, NAL8, that contributes to multiple developmental process and is required for normal leaf width and spikelet number at the reproductive stage in rice. These results were consistent with the ubiquitous expression pattern of NAL8 gene. We used genetic complementation, CRISPR/Cas9 gene editing system, RNAi gene silenced system and overexpressing system to generate transgenic plants for confirming the fuctions of NAL8. Mutation of NAL8 causes a reduction in the number of plastoglobules and shrunken thylakoids in chloroplasts, resulting in reduced cell division. In addition, the auxin levels in nal8 mutants are higher than in TQ, while the cytokinin levels are lower than in TQ. Moreover, RNA-sequencing and proteomics analysis shows that NAL8 is involved in multiple hormone signaling pathways as well as photosynthesis in chloroplasts and respiration in mitochondria. CONCLUSIONS: Our findings provide new insights into the way that NAL8 functions as a molecular chaperone in regulating plant leaf morphology and spikelet number through its effects on mitochondria and chloroplasts associated with cell division.

The MRE11-ATM-SOG1 DNA damage signaling pathway confers rice immunity to Xanthomonas oryzae.

Plants are constantly exposed to microbial pathogens in the environment. One branch of innate plant immunity is mediated by cell-membrane-localized receptors, but less is known about associations between DNA damage and plant immune responses. Here, we show that rice (Oryza sativa) mesophyll cells are prone to DNA double-stranded breaks (DSBs) in response to ZJ173, a strain of Xanthomonas oryzae pv. oryzae (Xoo). The DSB signal transducer ataxia telangiectasia mutated (ATM), but not the ATM and Rad3-related branch, confers resistance against Xoo. Mechanistically, the MRE11-ATM module phosphorylates suppressor of gamma response 1 (SOG1), which activates several phenylpropanoid pathway genes and prompts downstream phytoalexin biosynthesis during Xoo infection. Intriguingly, overexpression of the topoisomerase gene TOP6A3 causes a switch from the classic non-homologous end joining (NHEJ) pathway to the alternative NHEJ and homologous recombination pathways at Xoo-induced DSBs. The enhanced ATM signaling of the alternative NHEJ pathway strengthens the SOG1-regulated phenylpropanoid pathway and thereby boosts Xoo-induced phytoalexin biosynthesis in TOP6A3-OE1 overexpression lines. Overall, the MRE11-ATM-SOG1 pathway serves as a prime example of plant-pathogen interactions that occur via host non-specific recognition. The function of TOP6-facilitated ATM signaling in the defense response makes it a promising target for breeding of rice germplasm that exhibits resistance to bacterial blight disease without a growth penalty.

A pan-TE map highlights transposable elements underlying domestication and agronomic traits in Asian rice.

Transposable elements (TEs) are ubiquitous genomic components and hard to study due to being highly repetitive. Here we assembled 232 chromosome-level genomes based on long-read sequencing data. Coupling the 232 genomes with 15 existing assemblies, we developed a pan-TE map comprising both cultivated and wild Asian rice. We detected 177 084 high-quality TE variations and inferred their derived state using outgroups. We found TEs were one source of phenotypic variation during rice domestication and differentiation. We identified 1246 genes whose expression variation was associated with TEs but not single-nucleotide polymorphisms (SNPs), such as OsRbohB, and validated OsRbohB's relative expression activity using a dual-Luciferase (LUC) reporter assays system. Our pan-TE map allowed us to detect multiple novel loci associated with agronomic traits. Collectively, our findings highlight the contributions of TEs to domestication, differentiation and agronomic traits in rice, and there is massive potential for gene cloning and molecular breeding by the high-quality Asian pan-TE map we generated.

Uncovering key salt-tolerant regulators through a combined eQTL and GWAS analysis using the super pan-genome in rice.

For sessile plants, gene expression plays a pivotal role in responding to salinity stress by activating or suppressing specific genes. However, our knowledge of genetic variations governing gene expression in response to salt stress remains limited in natural germplasm. Through transcriptome analysis of the Global Mini-Core Rice Collection consisting of a panel of 202 accessions, we identified 22 345 and 27 610 expression quantitative trait loci associated with the expression of 7787 and 9361 eGenes under normal and salt-stress conditions, respectively, leveraging the super pan-genome map. Notably, combined with genome-wide association studies, we swiftly pinpointed the potential candidate gene STG5-a major salt-tolerant locus known as qSTS5. Intriguingly, STG5 is required for maintaining Na+/K+ homeostasis by directly regulating the transcription of multiple members of the OsHKT gene family. Our study sheds light on how genetic variants influence the dynamic changes in gene expression responding to salinity stress and provides a valuable resource for the mining of salt-tolerant genes in the future.

TurboID-mediated proximity labeling for screening interacting proteins of FIP37 in Arabidopsis.

Proximity labeling was recently developed to detect protein-protein interactions and members of subcellular multiprotein structures in living cells. Proximity labeling is conducted by fusing an engineered enzyme with catalytic activity, such as biotin ligase, to a protein of interest (bait protein) to biotinylate adjacent proteins. The biotinylated protein can be purified by streptavidin beads, and identified by mass spectrometry (MS). TurboID is an engineered biotin ligase with high catalytic efficiency, which is used for proximity labeling. Although TurboID-based proximity labeling technology has been successfully established in mammals, its application in plant systems is limited. Here, we report the usage of TurboID for proximity labeling of FIP37, a core member of m6A methyltransferase complex, to identify FIP37 interacting proteins in Arabidopsis thaliana. By analyzing the MS data, we found 214 proteins biotinylated by GFP-TurboID-FIP37 fusion, including five components of m6A methyltransferase complex that have been previously confirmed. Therefore, the identified proteins may include potential proteins directly involved in the m6A pathway or functionally related to m6A-coupled mRNA processing due to spatial proximity. Moreover, we demonstrated the feasibility of proximity labeling technology in plant epitranscriptomics study, thereby expanding the application of this technology to more subjects of plant research.

Synthesis and anticancer activity in vitro and in vivo evaluation of iridium(III) complexes on mouse melanoma B16 cells.

Combining the ligand NPIP (2-(2-nitrophenyl)-1H-imidazo[4,5-f][1,10]phenanthroline) with piq (1-phenylisoquinoline) and bzq (benzo[h]quinolone) gave [Ir(piq)2(NPIP)](PF6) (Ir1), and [Ir(bzq)2(NPIP)](PF6) (Ir2). The newly synthesized complexes were characterized by high-resolution mass spectrometry (HRMS), 1H NMR and 13C NMR. The complexes showed high antiproliferative activity against B16 cells. Three-dimensional (3D) cell model in vitro was used to evaluate the inhibitory effect of iridium (III) complex on B16 cells. The cellular uptake, mitochondrial localization, and intracellular distribution of the drugs confirmed that the iridium (III) complexes targeted the mitochondria, and the complexes can lead to the loss of mitochondrial membrane potential (MMP), increases the intracellular ROS content, further induces apoptosis. We also found that Ir1 and Ir2 can trigger the release of damage-associated molecular patterns (DAMPs) (cell surface calreticulin (CRT), heat-shock protein 70 (HSP70) and high mobility group box 1 (HMGB1)). In addition, Ir1 and Ir2 inhibited glutathione (GSH) synthesis and thus induced oxidative stress, Ir1 and Ir2 promoted malondialdehyde (MDA) production which is the stable metabolite of lipid peroxidation products. Finally, mice xenograft assay was performed to demonstrate that the complex shows higher antitumor activity in vivo than cisplatin. The inhibitory rates for cisplatin and Ir1 are 38.95% and 69.67%, respectively.

UDP-glucosyltransferase regulates grain size and abiotic stress tolerance associated with metabolic flux redirection in rice.

Grain size is an important component trait of grain yield, which is frequently threatened by abiotic stress. However, little is known about how grain yield and abiotic stress tolerance are regulated. Here, we characterize GSA1, a quantitative trait locus (QTL) regulating grain size and abiotic stress tolerance associated with metabolic flux redirection. GSA1 encodes a UDP-glucosyltransferase, which exhibits glucosyltransferase activity toward flavonoids and monolignols. GSA1 regulates grain size by modulating cell proliferation and expansion, which are regulated by flavonoid-mediated auxin levels and related gene expression. GSA1 is required for the redirection of metabolic flux from lignin biosynthesis to flavonoid biosynthesis under abiotic stress and the accumulation of flavonoid glycosides, which protect rice against abiotic stress. GSA1 overexpression results in larger grains and enhanced abiotic stress tolerance. Our findings provide insights into the regulation of grain size and abiotic stress tolerance associated with metabolic flux redirection and a potential means to improve crops.

Translational Regulation of Plant Response to High Temperature by a Dual-Function tRNAHis Guanylyltransferase in Rice.

As sessile organisms, plants have evolved numerous strategies to acclimate to changes in environmental temperature. However, the molecular basis of this acclimation remains largely unclear. In this study we identified a tRNAHis guanylyltransferase, AET1, which contributes to the modification of pre-tRNAHis and is required for normal growth under high-temperature conditions in rice. Interestingly, AET1 possibly interacts with both RACK1A and eIF3h in the endoplasmic reticulum. Notably, AET1 can directly bind to OsARF mRNAs including the uORFs of OsARF19 and OsARF23, indicating that AET1 is associated with translation regulation. Furthermore, polysome profiling assays suggest that the translational status remains unaffected in the aet1 mutant, but that the translational efficiency of OsARF19 and OsARF23 is reduced; moreover, OsARF23 protein levels are obviously decreased in the aet1 mutant under high temperature, implying that AET1 regulates auxin signaling in response to high temperature. Our findings provide new insights into the molecular mechanisms whereby AET1 regulates the environmental temperature response in rice by playing a dual role in tRNA modification and translational control.