OsGSTU34, a Bz2-like anthocyanin-related glutathione transferase transporter, is essential for rice (Oryza sativa L.) organs coloration.

Anthocyanins are a flavonoid compound known as one of the most important chromogenic substances. They play several functions, including health promotion and sustaining plants during adverse conditions. They are synthesized at the endoplasmic reticulum and sequestered in the vacuole. In this work, we generated knock-out lines of OsGSTU34, a glutathione transporter's tau gene family, with no transgene line and off-target through CRISPR/Cas9 mutagenesis and highlighted the loss of pigmentation in rice flowers, leaves, stems, shoots, and caryopsis. The anthocyanin quantification in the wild-type BLWT and mutant line BLG34-8 caryopsis showed that cyanidin-3-O-glucoside (C3G) and peonidin-3-O-glucoside (P3G) were almost undetectable in the mutant line. A tandem mass tag (TMT) labeling proteomic analysis was conducted to elucidate the proteomic changes in the BLWT and BLG34-8. The result revealed that 1175 proteins were altered, including 408 that were down-regulated and 767 that were upregulated. The accumulation of the OsGSTU34-related protein (Q8L576), along with several anthocyanin-related proteins, was down-regulated. The enrichment analysis showed that the down-regulated proteins were enriched in different pathways, among which the phenylpropanoid biosynthesis pathway, flavonoid biosynthesis metabolites, and anthocyanin biosynthesis pathway. Protein interaction network prediction revealed that glutathione-S-transferase (Q8L576) was connected to the proteins involved in the flavonoid and anthocyanin biosynthesis pathways, such as flavanone 3-dioxygenase 1 (Q7XM21), leucoanthocyanidin dioxygenase 1 (Q93VC3), 4-coumarate-CoA ligase 2 (Q42982), phenylalanine ammonia-lyase (P14717), chalcone synthase 1 (Q2R3A1), and 4-coumarate-CoA ligase 5 (Q6ZAC1). However, the expression of the most important anthocyanin biosynthesis gene was not altered, suggesting that only the transport mechanism was affected. Our findings highlight new insight into the anthocyanin pigmentation in black rice and provide new perspectives for future research.

Integrative HPLC profiling and transcriptome analysis revealed insights into anthocyanin accumulation and key genes at three developmental stages of black rice (Oryza sativa. L) caryopsis.

Introduction: Anthocyanins are plants' secondary metabolites belonging to the flavonoid class with potential health-promoting properties. They are greatly employed in the food industry as natural alternative food colorants for dairy and ready-to-eat desserts and pH indicators. These tremendous advantages make them economically important with increasing market trends. Black rice is a rich source of anthocyanin that can be used to ensure food and nutritional security around the world. However, research on anthocyanin accumulation and gene expression during rice caryopsis development is lacking. Methods: In this study, we combined high-performance liquid chromatography (HPLC) and transcriptome analysis to profile the changes in anthocyanin content and gene expression dynamics at three developmental stages (milky, doughy, and mature). Results: Our results showed that anthocyanin accumulation started to be visible seven days after flowering (DAF), increased rapidly from milky (11 DAF) to dough stage, then started decreasing after the peak was attained at 18 DAF. RNA-seq showed that 519 out of 14889, 477 out of 17914, and 1614 out of 18810 genes were uniquely expressed in the milky, doughy, and mature stages, respectively. We performed three pairwise comparisons: milky vs. dough, milky vs. mature, and dough vs. mature, and identified 6753, 9540, and 2531 DEGs, respectively. The DEGs' abundance was higher in milky vs. mature, with 5527 up-regulated genes and 4013 down-regulated genes, while it was smaller in dough vs. mature, with 1419 up-regulated genes and 1112 down-regulated DEGs. This result was consistent with the changes in anthocyanin profiling, and the expression of structural, regulatory, and transporter genes involved in anthocyanin biosynthesis showed their highest expression at the dough stage. Through the gene expression profile and protein interaction network, we deciphered six main contributors of the anthocyanin peak observed at dough stage, including OsANS, OsDFR, OsGSTU34, OsMYB3, OsbHLH015, and OsWD40-50. Discussion: This study is the first to report the investigation of anthocyanin and gene expression at three developmental stages of black rice caryopsis. The findings of this study could aid in predicting the best harvesting time to achieve maximum anthocyanin content and the best time to collect samples for various gene expression analysis, laying the groundwork for future research into the molecular mechanisms underlying rice caryopsis coloration.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.).

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F2 population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways' biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.