Molecular bases for differential aging programs between flag and second leaves during grain-filling in rice

Integrated transcriptomics and miRNAomics provide insights into the complex multi-tiered regulatory networks associated with coleoptile senescence in rice

Coleoptile is the small conical, short-lived, sheath-like organ that safeguards the first leaf and shoot apex in cereals. It is also the first leaf-like organ to senesce that provides nutrition to the developing shoot and is, therefore, believed to play a crucial role in seedling establishment in rice and other grasses. Though histochemical studies have helped in understanding the pattern of cell death in senescing rice coleoptiles, genome-wide expression changes during coleoptile senescence have not yet been explored. With an aim to investigate the gene regulation underlying the coleoptile senescence (CS), we performed a combinatorial whole genome expression analysis by sequencing transcriptome and miRNAome of senescing coleoptiles. Transcriptome analysis revealed extensive reprogramming of 3439 genes belonging to several categories, the most prominent of which encoded for transporters, transcription factors (TFs), signaling components, cell wall organization enzymes, redox homeostasis, stress response and hormone metabolism. Small RNA sequencing identified 41 known and 21 novel miRNAs that were differentially expressed during CS. Comparison of gene expression and miRNA profiles generated for CS with publicly available leaf senescence (LS) datasets revealed that the two aging programs are remarkably distinct at molecular level in rice. Integration of expression data of transcriptome and miRNAome identified high confidence 140 miRNA-mRNA pairs forming 42 modules, thereby demonstrating multi-tiered regulation of CS. The present study has generated a comprehensive resource of the molecular networks that enrich our understanding of the fundamental pathways regulating coleoptile senescence in rice.

Genome-wide H3K9 acetylation level increases with age-dependent senescence of flag leaf in rice

Flag leaf senescence is an important biological process that drives the remobilization of nutrients to the growing organs of rice. Leaf senescence is controlled by genetic information via gene expression and histone modification, but the precise mechanism is as yet unclear. Here, we analysed genome-wide acetylated lysine residue 9 of histone H3 (H3K9ac) enrichment by chromatin immunoprecipitation-sequencing (ChIP-seq), and examined its association with transcriptomes by RNA-seq during flag leaf aging in rice (Oryza sativa). We found that genome-wide H3K9 acetylation levels increased with age-dependent senescence in rice flag leaf, and there was a positive correlation between the density and breadth of H3K9ac with gene expression and transcript elongation. During flag leaf aging, we observed 1249 up-regulated differentially expressed genes (DEGs) and 996 down-regulated DEGs, showing a strong relationship between temporal changes in gene expression and gain/loss of H3K9ac. We produced a landscape of H3K9 acetylation-modified gene expression targets that include known senescence-associated genes, metabolism-related genes, as well as miRNA biosynthesis-related genes. Our findings reveal a complex regulatory network of metabolism- and senescence-related pathways mediated by H3K9ac, and elucidate patterns of H3K9ac-mediated regulation of gene expression during flag leaf aging in rice.

A Comparative Transcriptomic Meta-Analysis Revealed Conserved Key Genes and Regulatory Networks Involved in Drought Tolerance in Cereal Crops

Drought affects plant growth and development, causing severe yield losses, especially in cereal crops. The identification of genes involved in drought tolerance is crucial for the development of drought-tolerant crops. The aim of this study was to identify genes that are conserved key players for conferring drought tolerance in cereals. By comparing the transcriptomic changes between tolerant and susceptible genotypes in four Gramineae species, we identified 69 conserved drought tolerant-related (CDT) genes that are potentially involved in the drought tolerance of all of the analysed species. The CDT genes are principally involved in stress response, photosynthesis, chlorophyll biogenesis, secondary metabolism, jasmonic acid signalling, and cellular transport. Twenty CDT genes are not yet characterized and can be novel candidates for drought tolerance. The k-means clustering analysis of expression data highlighted the prominent roles of photosynthesis and leaf senescence-related mechanisms in differentiating the drought response between tolerant and sensitive genotypes. In addition, we identified specific transcription factors that could regulate the expression of photosynthesis and leaf senescence-related genes. Our analysis suggests that the balance between the induction of leaf senescence and maintenance of photosynthesis during drought plays a major role in tolerance. Fine-tuning of CDT gene expression modulation by specific transcription factors can be the key to improving drought tolerance in cereals.

Association mapping of drought tolerance and agronomic traits in rice (Oryza sativa L.) landraces

BACKGROUND

Asian cultivars were predominantly represented in global rice panel selected for sequencing and to identify novel alleles for drought tolerance. Diverse genetic resources adapted to Indian subcontinent were not represented much in spite harboring useful alleles that could improve agronomic traits, stress resilience and productivity. These rice accessions are valuable genetic resource in developing rice varieties suited to different rice ecosystem that experiences varying drought stress level, and at different crop stages. A core collection of rice germplasm adapted to Southwestern Indian peninsular genotyped using SSR markers and characterized by contrasting water regimes to associate genomic regions for physiological, root traits and yield related traits. Genotyping-By-Sequencing of selected accessions within the diverse panel revealed haplotype variation in genic content within genomic regions mapped for physiological, morphological and root traits.

RESULTS

Diverse rice panel (99 accessions) were evaluated in field and measurements on plant physiological, root traits and yield related traits were made over five different seasons experiencing varying drought stress intensity at different crop stages. Traits like chlorophyll stability index, leaf rolling, days to 50% flowering, chlorophyll content, root volume and root biomass were identified as best predictors of grain yield under stress. Association mapping revealed genetic variation among accessions and revealed 14 genomic targets associated with different physiological, root and plant production traits. Certain accessions were found to have beneficial allele to improve traits, plant height, root length and spikelet fertility, that contribute to the grain yield under stress. Genomic characterization of eleven accessions revealed haplotype variation within key genomic targets on chromosomes 1, 4, 6 and 11 for potential use as molecular markers to combine drought avoidance and tolerance traits. Genes mined within the genomic QTL intervals identified were prioritized based on tissue specific expression level in publicly available rice transcriptome data.

CONCLUSION

The genetic and genomic resources identified will enable combining traits with agronomic value to optimize yield under stress and hasten trait introgression into elite cultivars. Alleles associated with plant height, specific leaf area, root length from PTB8 and spikelet fertility and grain weight from PTB26 can be harnessed in future rice breeding program.

Current Understanding of Leaf Senescence in Rice

Leaf senescence, which is the last developmental phase of plant growth, is controlled by multiple genetic and environmental factors. Leaf yellowing is a visual indicator of senescence due to the loss of the green pigment chlorophyll. During senescence, the methodical disassembly of macromolecules occurs, facilitating nutrient recycling and translocation from the sink to the source organs, which is critical for plant fitness and productivity. Leaf senescence is a complex and tightly regulated process, with coordinated actions of multiple pathways, responding to a sophisticated integration of leaf age and various environmental signals. Many studies have been carried out to understand the leaf senescence-associated molecular mechanisms including the chlorophyll breakdown, phytohormonal and transcriptional regulation, interaction with environmental signals, and associated metabolic changes. The metabolic reprogramming and nutrient recycling occurring during leaf senescence highlight the fundamental role of this developmental stage for the nutrient economy at the whole plant level. The strong impact of the senescence-associated nutrient remobilization on cereal productivity and grain quality is of interest in many breeding programs. This review summarizes our current knowledge in rice on (i) the actors of chlorophyll degradation, (ii) the identification of stay-green genotypes, (iii) the identification of transcription factors involved in the regulation of leaf senescence, (iv) the roles of leaf-senescence-associated nitrogen enzymes on plant performance, and (v) stress-induced senescence. Compiling the different advances obtained on rice leaf senescence will provide a framework for future rice breeding strategies to improve grain yield.

Variation between rice accessions in photosynthetic induction in flag leaves and underlying mechanisms

Several breeding initiatives have sought to improve flag leaf performance as its health and physiology are closely correlated to rice yield. Previous studies have described natural variation of photosynthesis for flag leaves; however, none has examined their performance under the non-steady-state conditions that prevail in crop fields. Photosynthetic induction is the transient response of photosynthesis to a change from low to high light. Rice flag leaf photosynthesis was measured in both steady- and non-steady-state conditions to characterize natural variation. Between the lowest and highest performing accession, there was a 152% difference for average CO2 assimilation during induction (Ā300), a 77% difference for average intrinsic water use efficiency during induction (iWUEavg), and a 185% difference for the speed of induction (IT50), indicating plentiful variation. No significant correlation was found between steady- and non-steady-state photosynthetic traits. Additionally, measures of neither steady-state nor non-steady-state photosynthesis of flag leaves correlated with the same measures of leaves in the vegetative growth stage, with the exception of iWUEavg. Photosynthetic induction was measured at six [CO2], to determine biochemical and diffusive limitations to photosynthesis in vivo. Photosynthetic induction in rice flag leaves was limited primarily by biochemistry.

ACCELERATED CELL DEATH 6 Acts on Natural Leaf Senescence and Nitrogen Fluxes in Arabidopsis

As the last step of leaf development, senescence is a molecular process involving cell death mechanism. Leaf senescence is trigged by both internal age-dependent factors and environmental stresses. It must be tightly regulated for the plant to adopt a proper response to environmental variation and to allow the plant to recycle nutrients stored in senescing organs. However, little is known about factors that regulate both nutrients fluxes and plant senescence. Taking advantage of variation for natural leaf senescence between Arabidopsis thaliana accessions, Col-0 and Ct-1, we did a fine mapping of a quantitative trait loci for leaf senescence and identified ACCELERATED CELL DEATH 6 (ACD6) as the causal gene. Using two near-isogeneic lines, differing solely around the ACD6 locus, we showed that ACD6 regulates rosette growth, leaf chlorophyll content, as well as leaf nitrogen and carbon percentages. To unravel the role of ACD6 in N remobilization, the two isogenic lines and acd6 mutant were grown and labeled with ¹⁵N at the vegetative stage in order to determine ¹⁵N partitioning between plant organs at harvest. Results showed that N remobilization efficiency was significantly lower in all the genotypes with lower ACD6 activity irrespective of plant growth and productivity. Measurement of N uptake at vegetative and reproductive stages revealed that ACD6 did not modify N uptake efficiency but enhanced nitrogen translocation from root to silique. In this study, we have evidenced a new role of ACD6 in regulating both sequential and monocarpic senescences and disrupting the balance between N remobilization and N uptake that is required for a good seed filling.

Effect of exogenously-applied abscisic acid, putrescine and hydrogen peroxide on drought tolerance of barley

The objective of this study was to identify the effect of abscisic acid (ABA), putrescine (Put) and hydrogen peroxide (H2O2) foliar pre-treatment on drought tolerance of barley. Despite water limitation, ABA-sprayed plants preserved increased water content, photosynthetic efficiency of PSII (ΦPSII) and CO2 assimilation rate (Pn) compared to untreated stressed plants. The ABA-treated plants presented also the lowest rate of lipid peroxidation (MDA), lowered the rate of PSII primary acceptor reduction (1 – qP) and increased the yield of regulated energy dissipation (NPQ) with higher accumulation of PGRL1 (PROTON GRADIENT REGULATION LIKE1) protein. These plants preserved a similar level of photochemical efficiency and the rate of electron transport of PSII (ETRII) to the well-watered samples. The significantly less pronounced response was observed in Put-sprayed samples under drought. Additionally, the combined effects of drought and H2O2 application increased the 1 – qP and quantum yield of non-regulated energy dissipation in PSII (ΦNO) and reduced the accumulation of Rubisco activase (RCA). In conclusion, ABA foliar application allowed to balance water retention and preserve antioxidant capacity resulting in efficient photosynthesis and the restricted risk of oxidative damage under drought. Neither hydrogen peroxide nor putrescine has been able to ameliorate drought stress as effectively as ABA.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

Natural variations at the Stay-Green gene promoter control lifespan and yield in rice cultivars

Increased grain yield will be critical to meet the growing demand for food, and could be achieved by delaying crop senescence. Here, via quantitative trait locus (QTL) mapping, we uncover the genetic basis underlying distinct life cycles and senescence patterns of two rice subspecies, indica and japonica. Promoter variations in the Stay-Green (OsSGR) gene encoding the chlorophyll-degrading Mg⁺⁺-dechelatase were found to trigger higher and earlier induction of OsSGR in indica, which accelerated senescence of indica rice cultivars. The indica-type promoter is present in a progenitor subspecies O. nivara and thus was acquired early during the evolution of rapid cycling trait in rice subspecies. Japonica OsSGR alleles introgressed into indica-type cultivars in Korean rice fields lead to delayed senescence, with increased grain yield and enhanced photosynthetic competence. Taken together, these data establish that naturally occurring OsSGR promoter and related lifespan variations can be exploited in breeding programs to augment rice yield.

A RING-Type E3 Ubiquitin Ligase, OsGW2, Controls Chlorophyll Content and Dark-Induced Senescence in Rice

Leaf senescence is the final stage of plant development. Many internal and external factors affect the senescence process in rice (Oryza sativa L.). In this study, we identified qCC2, a major quantitative trait locus (QTL) for chlorophyll content using a population derived from an interspecific cross between O. sativa (cv. Hwaseong) and Oryza grandiglumis. The O. grandiglumis allele at qCC2 increased chlorophyll content and delayed senescence. GW2 encoding E3 ubiquitin ligase in the qCC2 region was selected as a candidate for qCC2. To determine if GW2 is allelic to qCC2, a gw2-knockout mutant (gw2-ko) was examined using a dark-induced senescence assay. gw2-ko showed delayed leaf senescence in the dark with down-regulated expression of senescence-associated genes (SAGs) and chlorophyll degradation genes (CDGs). The association of the GW2 genotype with the delayed senescence phenotype was confirmed in an F₂ population. RNA-seq analysis was conducted to investigate 30-day-old leaf transcriptome dynamics in Hwaseong and a backcross inbred line-CR2002-under dark treatment. This resulted in the identification of genes involved in phytohormone signaling and associated with senescence. These results suggested that transcriptional regulation was associated with delayed senescence in CR2002, and RING-type E3 ubiquitin ligase GW2 was a positive regulator of leaf senescence in rice.

Differential Expression of Genes at Panicle Initiation and Grain Filling Stages Implied in Heterosis of Rice Hybrids

RNA-Seq technology was used to analyze the transcriptome of two rice hybrids, Ajay (based on wild-abortive (WA)-cytoplasm) and Rajalaxmi (based on Kalinga-cytoplasm), and their respective parents at the panicle initiation (PI) and grain filling (GF) stages. Around 293 and 302 million high quality paired-end reads of Ajay and Rajalaxmi, respectively, were generated and aligned against the Nipponbare reference genome. Transcriptome profiling of Ajay revealed 2814 and 4819 differentially expressed genes (DEGs) at the PI and GF stages, respectively, as compared to its parents. In the case of Rajalaxmi, 660 and 5264 DEGs were identified at PI and GF stages, respectively. Functionally relevant DEGs were selected for validation through qRT-PCR, which were found to be co-related with the expression patterns to RNA-seq. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated significant DEGs enriched for energy metabolism pathways, such as photosynthesis, oxidative phosphorylation, and carbon fixation, at the PI stage, while carbohydrate metabolism-related pathways, such as glycolysis and starch and sucrose metabolism, were significantly involved at the GF stage. Many genes involved in energy metabolism exhibited upregulation at the PI stage, whereas the genes involved in carbohydrate biosynthesis had higher expression at the GF stage. The majority of the DEGs were successfully mapped to know yield related rice quantitative trait loci (QTLs). A set of important transcription factors (TFs) was found to be encoded by the identified DEGs. Our results indicated that a complex interplay of several genes in different pathways contributes to higher yield and vigor in rice hybrids.

Identification of Transcription Factors Regulating Senescence in Wheat through Gene Regulatory Network Modelling

Senescence is a tightly regulated developmental program coordinated by transcription factors. Identifying these transcription factors in crops will provide opportunities to tailor the senescence process to different environmental conditions and regulate the balance between yield and grain nutrient content. Here, we use ten time points of gene expression data along with gene network modeling to identify transcription factors regulating senescence in polyploid wheat (Triticum aestivum). We observe two main phases of transcriptional changes during senescence: early down-regulation of housekeeping functions and metabolic processes followed by up-regulation of transport and hormone-related genes. These two phases are largely conserved with Arabidopsis (Arabidopsis thaliana), although the individual genes underlying these changes are often not orthologous. We have identified transcription factor families associated with these early and later waves of differential expression. Using gene regulatory network modeling, we identified candidate transcription factors that may control senescence. Using independent, publicly available datasets, we found that the most highly ranked candidate genes in the network were enriched for senescence-related functions compared with all genes in the network. We validated the function of one of these candidate transcription factors in senescence using wheat chemically induced mutants. This study lays the groundwork to understand the transcription factors that regulate senescence in polyploid wheat and exemplifies the integration of time-series data with publicly available expression atlases and networks to identify candidate regulatory genes.

Rice Overexpressing OsNUC1-S Reveals Differential Gene Expression Leading to Yield Loss Reduction after Salt Stress at the Booting Stage

Rice nucleolin (OsNUC1), consisting of two isoforms, OsNUC1-L and OsNUC1-S, is a multifunctional protein involved in salt-stress tolerance. Here, OsNUC1-S’s function was investigated using transgenic rice lines overexpressing OsNUC1-S. Under non-stress conditions, the transgenic lines showed a lower yield, but higher net photosynthesis rates, stomatal conductance, and transpiration rates than wild type only in the second leaves, while in the flag leaves, these parameters were similar among the lines. However, under salt-stress conditions at the booting stage, the higher yields in transgenic lines were detected. Moreover, the gas exchange parameters of the transgenic lines were higher in both flag and second leaves, suggesting a role for OsNUC1-S overexpression in photosynthesis adaptation under salt-stress conditions. Moreover, the overexpression lines could maintain light-saturation points under salt-stress conditions, while a decrease in the light-saturation point owing to salt stress was found in wild type. Based on a transcriptome comparison between wild type and a transgenic line, after 3 and 9 days of salt stress, the significantly differentially expressed genes were enriched in the metabolic process of nucleic acid and macromolecule, photosynthesis, water transport, and cellular homeostasis processes, leading to the better performance of photosynthetic processes under salt-stress conditions at the booting stage.