Identification and Characterization of Differentially Expressed Genes in Inferior and Superior Spikelets of Rice Cultivars with Contrasting Panicle-Compactness and Grain-Filling Properties

Genes determining panicle morphology and grain quality in rice (Oryza sativa)

The world's increase in rice (Oryza sativa L.) production is not keeping up with the increase in its population. To boost the introduction of new high-yielding cultivars, knowledge is being gained on the genes and quantitative trait loci (QTLs) determining the panicle phenotype. The important are those determining yield of the crop, such as grain numbers per panicle and size and weight of the grains. Biochemical and molecular functions of many of them are understood in some details. Among these, OsCKX2 and OsSPL14 have been shown to increase panicle branching and grain numbers when overexpressed. Furthermore, miRNAs appear to play an important role in determining the panicle morphology by regulating the expressions of the genes like OsSPL14 and GRF4 involved in panicle branching and grain numbers and length. Mutations also greatly influence the grain shape and size. However, the information gained so far on the genetic regulation of grain filling and panicle morphology has not been successfully put into commercial application. Furthermore, the identification of the gene(s)/QTLs regulating panicle compactness is still lacking, which may enable the researchers to convert a compact-panicle cultivar into a lax/open one, and thereby increasing the chances of enhancing the yield of a desired compact-panicle cultivar obtained by the breeding effort.

Effect of Panicle Morphology on Grain Filling and Rice Yield: Genetic Control and Molecular Regulation

The demand for rice is likely to increase approximately 1.5 times by the year 2050. In contrast, the rice production is stagnant since the past decade as the ongoing rice breeding program is unable to increase the production further, primarily because of the problem in grain filling. Investigations have revealed several reasons for poor filling of the grains in the inferior spikelets of the compact panicle, which are otherwise genetically competent to develop into well-filled grains. Among these, the important reasons are 1) poor activities of the starch biosynthesizing enzymes, 2) high ethylene production leading to inhibition in expressions of the starch biosynthesizing enzymes, 3) insufficient division of the endosperm cells and endoreduplication of their nuclei, 4) low accumulation of cytokinins and indole-3-acetic acid (IAA) that promote grain filling, and 5) altered expressions of the miRNAs unfavorable for grain filling. At the genetic level, several genes/QTLs linked to the yield traits have been identified, but the information so far has not been put into perspective toward increasing the rice production. Keeping in view the genetic competency of the inferior spikelets to develop into well-filled grains and based on the findings from the recent research studies, improving grain filling in these spikelets seems plausible through the following biotechnological interventions: 1) spikelet-specific knockdown of the genes involved in ethylene synthesis and overexpression of β-CAS (β-cyanoalanine) for enhanced scavenging of CN⁻ formed as a byproduct of ethylene biosynthesis; 2) designing molecular means for increased accumulation of cytokinins, abscisic acid (ABA), and IAA in the caryopses; 3) manipulation of expression of the transcription factors like MYC and OsbZIP58 to drive the expression of the starch biosynthesizing enzymes; 4) spikelet-specific overexpression of the cyclins like CycB;1 and CycH;1 for promoting endosperm cell division; and 5) the targeted increase in accumulation of ABA in the straw during the grain filling stage for increased carbon resource remobilization to the grains. Identification of genes determining panicle compactness could also lead to an increase in rice yield through conversion of a compact-panicle into a lax/open one. These efforts have the ability to increase rice production by as much as 30%, which could be more than the set production target by the year 2050.

Biochemical and molecular processes contributing to grain filling and yield in rice

The increase in much required rice production through breeding programmes is on decline. The primary reason being poor filling of grains in the basal spikelets of the heavy and compact panicle rice developed. These spikelets are genetically competent to develop into well filled grains, but fail to do so because the carbohydrate assimilates available to them remain unutilized, reportedly due to poor activities of the starch biosynthesizing enzymes, high production of ethylene leading to enhanced synthesis of the downstream signaling component RSR1 protein that inhibits GBSS1 activity, poor endosperm cell division and endoreduplication of the endosperm nuclei, altered expression of the transcription factors influencing grain filling, enhanced expression and phosphorylation of 14-3-3 proteins, poor expression of the seed storage proteins, reduced synthesis of the hormones like cytokinins and IAA that promote grain filling, and altered expression of miRNAs preventing their normal role in grain filling. Since the basal spikelets are genetically competent to develop into well filled mature grains, biotechnological interventions in terms of spikelet-specific overexpression of the genes encoding enzymes involved in grain filling and/or knockdown/overexpression of the genes influencing the activities of the starch biosynthesizing enzymes, various cell cycle events and hormone biosynthesis could increase rice production by as much as 30%, much more than the set production target of 800 mmt. Application of these biotechnological interventions in the heavy and compact panicle cultivars producing grains of desired quality would also maintain the quality of the grains having demand in market besides increasing the rice production per se.

Transcriptomic identification of long noncoding RNAs and their hormone-associated nearby coding genes involved in the differential development of caryopses localized on different branches in rice

Long noncoding RNAs (lncRNAs) play important regulatory roles in caryopsis development and grain size in rice. However, whether there exist differences in lncRNA expression between caryopses located on primary branches (CPB) and caryopses located on secondary branches (CSB) that contribute to their differential development remains elusive. Here, we performed transcriptome-wide analysis to identify 2,273 lncRNAs expressed in CPB and CSB at 0, 5, 12, and 20 days after flowering (DAF). Although these lncRNAs were widely distributed, the majority were located in intergenic regions of the 12 rice chromosomes. Based on gene expression cluster analysis, lncRNAs expressed in CPB and CSB were clustered into two subtypes in a position-independent manner: one includes 0- and 5-DAF CPB and CSB, and 12-DAF CSB; the second includes 12-DAF CPB and 20-DAF CPB and CSB. Furthermore, according to the expression value of each lncRNA, K-means cluster analysis revealed 135 early-stage, 116 middle-stage, and 114 late-stage expression-delayed lncRNAs in CSB. Then, we analyzed the expression values of the expression-delayed lncRNAs and nearby coding genes (100 kb upstream and downstream of the lncRNAs), and found 631 lncRNA-mRNA pairs, including 258 lncRNAs and 571 nearby coding genes, some of which are related to hormone-regulated grain development. These results suggested that expression-delayed lncRNAs in CSB may regulate the development of CPB and CSB, providing insight into the mechanism underlying the developmental differences between CPB and CSB, and the differences in grain yield.

Identification of novel QTLs for grain fertility and associated traits to decipher poor grain filling of basal spikelets in dense panicle rice

High grain number is positively correlated with grain yield in rice, but it is compromised because of poor filling of basal spikelets in dense panicle bearing numerous spikelets. The phenomenon that turns the basal spikelets of compact panicle sterile in rice is largely unknown. In order to understand the factor(s) that possibly determines such spikelet sterility in compact panicle cultivars, QTLs and candidate genes were identified for spikelet fertility and associated traits like panicle compactness, and ethylene production that significantly influences the grain filling using recombinant inbred lines developed from a cross between indica rice cultivars, PDK Shriram (compact, high spikelet number) and Heera (lax, low spikelet number). Novel QTLs, qSFP1.1, qSFP3.1, and qSFP6.1 for spikelet fertility percentage; qIGS3.2 and qIGS4.1 for panicle compactness; and qETH1.2, qETH3.1, and qETH4.1 for ethylene production were consistently identified in both kharif seasons of 2017 and 2018. The comparative expression analysis of candidate genes like ERF3, AP2-like ethylene-responsive transcription factor, EREBP, GBSS1, E3 ubiquitin-protein ligase GW2, and LRR receptor-like serine/threonine-protein kinase ERL1 associated with identified QTLs revealed their role in poor grain filling of basal spikelets in a dense panicle. These candidate genes thus could be important for improving grain filling in compact-panicle rice cultivars through biotechnological interventions.

MicroRNAs modulate ethylene induced retrograde signal for rice endosperm starch biosynthesis by default expression of transcriptome

Control of stage specific spike in ethylene production at anthesis has been a vauable route to potentially enhance genetic ceiling for grain filling of rice spikelet. A number of genes controlling ethylene homeostasis and starch synthesis have been identified so long, but lack of credible information on master modulation of gene expression by miRNAs and their target genes associated with hormonal dynamics obfuscate mechanisms controlling genotype difference in quantum of grain filling. The confusion accounts for consequent shrinkage of options for yield manipulation. In a two by two factorial design, miRNA regulation of spikelet specific grain development in low against high sterile recombinant inbred lines of rice Oryza sativa L. namely CR 3856-62-11-3-1-1-1-1-1-1 (SR 157) and CR 3856-63-1-1-1-1-1-1-1 (SR 159) respectively, and inferior verses superior spikelets were compared during first 10 days after anthesis. Grain filling was poorer in SR159 than SR157 and inferior spikelets in the former were most vulnerable. Between the cultivars, overall expression of unique miRNAs with targets on ethylene pathway genes was higher in SR159 than SR157 and the situation was opposite for auxin pathway genes. Precision analysis in psTarget server database identified up-regulation of MIR2877 and MIR530-5p having Os11t0141000-02 and Os07t0239400-01 (PP2A regulatory subunit-like protein and ethylene-responsive small GTP-binding proteins) and MIR396h having Os01t0643300-02 (an auxin efflux carrier protein) and Os01t0643300-01 (a PIN1-like auxin transport protein), as targets with highest probability at anthesis and 5 days after anthesis respectively, in the inferior spikelet and the fold change values of DGE matched with pattern of gene expression (relative transcript level) in the qRT-PCR studies conducted for relevant miRNAs and protein factors for ethylene and auxin signalling. In conclusion, epigenetic regulation of both auxin and ethylene homeostasis control grain filling of rice spikelet was established, but evidences were more robust for the latter.

Cell cycle events and expression of cell cycle regulators are determining factors in differential grain filling in rice spikelets based on their spatial location on compact panicles

Rice being a staple crop for human, its production is required to be increased significantly, particularly keeping in view the expected world's population of 9.6 billion by the year 2050. In this context, although the rice breeding programs have been successful in increasing the number of spikelets per panicle, the basal spikelets remain poorly filled, undermining the yield potential. The present study also found the grain filling to bear negative correlation with the panicle grain density. The poorly filled basal spikelets of the compact-panicle cultivars showed a lower endosperm cell division rate and ploidy status of the endosperm nuclei coupled with no significant greater expression of CYCB;1 and CYCH;1 compared with the apical spikelets, unlike that observed in the lax-panicle cultivars, which might have prevented them from overcoming apical dominance. Significantly greater expression of CYCB2;2 in the basal spikelets than in the apical spikelets might also have prevented the former to enter into endoreduplication. Furthermore, expression studies of KRPs in the caryopses revealed that a higher expression of KRP;1 and KRP;4 in the basal spikelets than in the apical spikelets of the compact-panicle cultivars could also be detrimental to grain filling in the former, as KRPs form complex primarily with CDKA-CYCD that promotes S-phase activity and G1/S transition, and thus inhibits endosperm cell division. The study indicates that targeted manipulation of expression of CYCB1;1, CYCB2;2, CYCH1;1, KRP;1 and KRP4 in the basal spikelets of the compact-panicle cultivars may significantly improve their yield performance.

Study of expressions of miRNAs in the spikelets based on their spatial location on panicle in rice cultivars provided insight into their influence on grain development

Development of rice cultivars bearing numerous spikelets by breeding approach to increase the yearly production of rice to approximately 800 million metric tons to feed the ever increasing population of the world accompanies poor grain filling in the inferior spikelets preventing achievement of the yield potential. As the initial stages of caryopses development are of much importance for grain filling, spatio-temporal expressions of the miRNAs were studied during these periods in the spikelets of a compact-panicle rice cultivar, Oryza sativa cv. Mahalaxmi, bearing numerous spikelets per panicle to understand the reason of poor grain filling at the level of the initial biochemical events. Differential expression of several known miRNAs between the superior and inferior spikelets suggested great difference in metabolism related to grain filling in the spikelets based on their spatial location on compact panicle. Expressions of five known and four novel miRNAs were validated by Northern. Their targets included the enzymes directly involved in starch biosynthesis like sucrose synthase, starch synthase and pullulanase, besides others. Spatio-temporal expression studies of these miRNAs in the spikelets of Mahalaxmi revealed a pattern of mostly a greater expression in the inferior spikelets compared with the superior ones concomitant with an inverse expression of the target genes, which was not observed in the lax-panicle cultivar Upahar. The study thus revealed that the grain filling in rice is greatly regulated by miRNAs, and these miRNAs or their target genes could be considered for biotechnological interventions for improving grain filling in the rice cultivars of interest.

Transcriptomic analysis reveals the contribution of auxin on the differentially developed caryopses on primary and secondary branches in rice

Generally, the caryopses located on proximal secondary branches (CSB) have smaller grain size and slower and poorer filling rate than those on apical primary branches (CPB) in rice, which greatly limits the grain yield potential fulfillment. However, the key regulators determining the developmental differences between CPB and CSB remain elusive. Here, we have performed transcriptomic analysis in CPB and CSB at four developmental stages [0, 5, 12 and 20 days after fertilization (DAF)] using high-throughput RNA-sequencing technique. Based on gene expression cluster analysis, the genes expressed in CPB and CSB were clustered into two subtypes in a positional-independent manner: one includes 0- and 5-DAF CPB and CSB, and 12-DAF CSB; another includes 12-DAF CPB, 20-DAF CPB and CSB. Moreover, according to the expression value of each gene, K-mean cluster analysis showed that the K4 to K6 classifiers contain the genes highly expressed in 5-DAF CPB and 12-DAF CSB, which were enriched in DNA synthesis, protein synthesis and cell proliferation mainly responsible for grain size decision. Then, functional enrichment analysis in Gene Ontology database showed that auxin-related genes were relatively enriched, indicating that auxin might be the key determinant for gene expression in K4 to K6 classifiers. Finally, the application of exogenous IAA in CSB before fertilization promoted gene expression, caryopsis development and grain weight closer to that in CPB, providing a molecular framework to optimize CSB development and potential targets for increasing grain yield.

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

The 14-3-3 protein GF14f negatively affects grain filling of inferior spikelets of rice (Oryza sativa L.)

In rice (Oryza sativa L.), later flowering inferior spikelets (IS), which are located on proximal secondary branches, fill slowly and produce smaller and lighter grains than earlier flowering superior spikelets (SS). Many genes have been reported to be involved in poor grain filling of IS, however the underlying molecular mechanisms remain unclear. The present study determined that GF14f, a member of the 14-3-3 protein family, showed temporal and spatial differences in expression patterns between SS and IS. Using GF14f-RNAi plants, we observed that a reduction in GF14f expression in the endosperm resulted in a significant increase in both grain length and weight, which in turn improved grain yield. Furthermore, pull-down assays indicated that GF14f interacts with enzymes that are involved in sucrose breakdown, starch synthesis, tricarboxylic acid (TCA) cycle and glycolysis. At the same time, an increase in the activity of sucrose synthase (SuSase), adenosine diphosphate-glucose pyrophosphorylase (AGPase), and starch synthase (StSase) was observed in the GF14f-RNAi grains. Comprehensive analysis of the proteome and metabolite profiling revealed that the abundance of proteins related to the TCA cycle, and glycolysis increased in the GF14f-RNAi grains together with several carbohydrate intermediates. These results suggested that GF14f negatively affected grain development and filling, and the observed higher abundance of the GF14f protein in IS compared with SS may be responsible for poor IS grain filling. The study provides insights into the molecular mechanisms underlying poor grain filling of IS and suggests that GF14f could serve as a potential tool for improving rice grain filling.

SWATH-MS Quantitative Analysis of Proteins in the Rice Inferior and Superior Spikelets during Grain Filling

Modern rice cultivars have large panicle but their yield potential is often not fully achieved due to poor grain-filling of late-flowering inferior spikelets (IS). Our earlier work suggested a broad transcriptional reprogramming during grain filling and showed a difference in gene expression between IS and earlier-flowering superior spikelets (SS). However, the links between the abundances of transcripts and their corresponding proteins are unclear. In this study, a SWATH-MS (sequential window acquisition of all theoretical spectra-mass spectrometry) -based quantitative proteomic analysis has been applied to investigate SS and IS proteomes. A total of 304 proteins of widely differing functionality were observed to be differentially expressed between IS and SS. Detailed gene ontology analysis indicated that several biological processes including photosynthesis, protein metabolism, and energy metabolism are differentially regulated. Further correlation analysis revealed that abundances of most of the differentially expressed proteins are not correlated to the respective transcript levels, indicating that an extra layer of gene regulation which may exist during rice grain filling. Our findings raised an intriguing possibility that these candidate proteins may be crucial in determining the poor grain-filling of IS. Therefore, we hypothesize that the regulation of proteome changes not only occurs at the transcriptional, but also at the post-transcriptional level, during grain filling in rice.