RiceXPro: a platform for monitoring gene expression in japonica rice grown under natural field conditions

Analysis of Germin-like Protein Genes (OsGLPs) Family in Rice Using Various In silico Approaches

Background:

Germin-like Proteins (GLPs) play an important role in various stresses. Rice contains 43 GLPs, among which many remain functionally unexplored. The computational analysis will provide significant insight into their function.

Objective:

To find various structural properties, functional importance, phylogeny and expression pattern of all OsGLPs using various bioinformatics tools.

Methods:

Physiochemical properties, sub-cellular localization, domain composition, Nglycosylation and Phosphorylation sites, and 3D structural models of the OsGLPs were predicted using various bioinformatics tools. Functional analysis was carried out with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Blast2GO servers. The expression profile of the OsGLPs was predicted by retrieving the data for expression values from tissuespecific and hormonal stressed array libraries of RiceXPro. Their phylogenetic relationship was computed using Molecular and Evolutionary Genetic Analysis (MEGA6) tool.

Results:

Most of the OsGLPs are stable in the cellular environment with a prominent expression in the extracellular region (57%) and plasma membrane (33%). Besides, 3 basic cupin domains, 7 more were reported, among which NTTNKVGSNVTLINV, FLLAALLALASWQAI, and MASSSF were common to 99% of the sequences, related to bacterial pathogenicity, peroxidase activity, and peptide signal activity, respectively. Structurally, OsGLPs are similar but functionally they are diverse with novel enzymatic activities of oxalate decarboxylase, lyase, peroxidase, and oxidoreductase. Expression analysis revealed prominent activities in the root, endosperm, and leaves. OsGLPs were strongly expressed by abscisic acid, auxin, gibberellin, cytokinin, and brassinosteroid. Phylogenetically they showed polyphyletic origin with a narrow genetic background of 0.05%. OsGLPs of chromosome 3, 8, and 12 are functionally more important due to their defensive role against various stresses through co-expression strategy.

Conclusion:

The analysis will help to utilize OsGLPs in future food programs.

Evolution of Abscisic Acid Signaling for Stress Responses to Toxic Metals and Metalloids

Toxic heavy metals and metalloids in agricultural ecosystems are crucial factors that limit global crop productivity and food safety. Industrial toxic heavy metals and metalloids such as cadmium, lead, and arsenic have contaminated large areas of arable land in the world and their accumulation in the edible parts of crops is causing serious health risks to humans and animals. Plants have co-evolved with various concentrations of these toxic metals and metalloids in soil and water. Some green plant species have significant innovations in key genes for the adaptation of abiotic stress tolerance pathways that are able to tolerate heavy metals and metalloids. Increasing evidence has demonstrated that phytohormone abscisic acid (ABA) plays a vital role in the alleviation of heavy metal and metalloid stresses in plants. Here, we trace the evolutionary origins of the key gene families connecting ABA signaling with tolerance to heavy metals and metalloids in green plants. We also summarize the molecular and physiological aspects of ABA in the uptake, root-to-shoot translocation, chelation, sequestration, reutilization, and accumulation of key heavy metals and metalloids in plants. The molecular evolution and interaction between the ABA signaling pathway and mechanisms for heavy metal and metalloid tolerance are highlighted in this review. Therefore, we propose that it is promising to manipulate ABA signaling in plant tissues to reduce the uptake and accumulation of toxic heavy metals and metalloids in crops through the application of ABA-producing bacteria or ABA analogues. This may lead to improvements in tolerance of major crops to heavy metals and metalloids.

Public-transcriptome-database-assisted selection and validation of reliable reference genes for qRT-PCR in rice

Accurate quantitative reverse transcription PCR (qRT-PCR) requires reliable reference genes whose expression does not vary in different tissues and developmental stages. However, few reliable reference genes are available for qRT-PCR in rice (Oryza sativa). Here, we established an effective strategy for identifying novel reference genes (NRGs) for reliable normalization of qRT-PCR data in various rice organs and developmental stages. We selected candidate NRGs using the Information Commons for Rice Database and confirmed their expression in Rice Expression Profile Database (RiceXPro) data. Genes with low variation (<2.5 cycle quantification) across tissues and developmental stages, and little fluctuation in expression in heatmaps from RiceXPro data were considered stable NRGs. To validate this strategy, we selected 11 candidate NRGs and calculated their expression stability in different spatio-temporal conditions using five programs, and compared these genes with five established reference genes (ERGs). Only one of the ERGs (UBQ5) was reliable and 10 of the candidate NRGs were more stable than the four remaining ERGs. Therefore, public transcriptomic databases are useful for identifying NRGs. We selected two NRGs, UFC1 (Homolog of UFM1-Conjugating Enzyme 1) and FhaB (Homolog of Adhesin FhaB) for qRT-PCR analysis in rice; their homologs might be suitable for other monocot plants.

The overexpression of rice ACYL-CoA-BINDING PROTEIN2 increases grain size and bran oil content in transgenic rice

As Oryza sativa (rice) seeds represent food for over three billion people worldwide, the identification of genes that enhance grain size and composition is much desired. Past reports have indicated that Arabidopsis thaliana acyl-CoA-binding proteins (ACBPs) are important in seed development but did not affect seed size. Herein, rice OsACBP2 was demonstrated not only to play a role in seed development and germination, but also to influence grain size. OsACBP2 mRNA accumulated in embryos and endosperm of germinating seeds in qRT-PCR analysis, while β-glucuronidase (GUS) assays on OsACBP2pro::GUS rice transformants showed GUS expression in embryos, as well as the scutellum and aleurone layer of germinating seeds. Deletion analysis of the OsACBP2 5'-flanking region revealed five copies of the seed cis-element, Skn-I-like motif (-1486/-1482, -956/-952, -939/-935, -826/-822, and -766/-762), and the removal of any adversely affected expression in seeds, thereby providing a molecular basis for OsACBP2 expression in seeds. When OsACBP2 function was investigated using osacbp2 mutants and transgenic rice overexpressing OsACBP2 (OsACBP2-OE), osacbp2 was retarded in germination, while OsACBP2-OEs performed better than the wild-type and vector-transformed controls, in germination, seedling growth, grain size and grain weight. Transmission electron microscopy of OsACBP2-OE mature seeds revealed an accumulation of oil bodies in the scutellum cells, while confocal laser scanning microscopy indicated oil accumulation in OsACBP2-OE aleurone tissues. Correspondingly, OsACBP2-OE seeds showed gain in triacylglycerols and long-chain fatty acids over the vector-transformed control. As dietary rice bran contains beneficial bioactive components, OsACBP2 appears to be a promising candidate for enriching seed nutritional value.

OsbHLH058 and OsbHLH059 transcription factors positively regulate iron deficiency responses in rice

Subgroup IVc basic helix-loop-helix transcription factors OsbHLH058 and OsbHLH059 positively regulate major iron deficiency responses in rice in a similar but distinct manner, putatively under partial control by OsHRZs. Under low iron availability, plants transcriptionally induce the expression of genes involved in iron uptake and translocation. OsHRZ1 and OsHRZ2 ubiquitin ligases negatively regulate this iron deficiency response in rice. The basic helix-loop-helix (bHLH) transcription factor OsbHLH060 interacts with OsHRZ1, and positively regulates iron deficiency-inducible genes. However, the functions of three other subgroup IVc bHLH transcription factors in rice, OsbHLH057, OsbHLH058, and OsbHLH059, have not yet been characterized. In the present study, we investigated the functions of OsbHLH058 and OsbHLH059 related to iron deficiency response. OsbHLH058 expression was repressed under iron deficiency, whereas the expression of OsbHLH057 and OsbHLH060 was moderately induced. Yeast two-hybrid analysis indicated that OsbHLH058 interacts with OsHRZ1 and OsHRZ2 more strongly than OsbHLH060, whereas OsbHLH059 showed no interaction. An in vitro ubiquitination assay detected no OsbHLH058 and OsbHLH060 ubiquitination by OsHRZ1 and OsHRZ2. Transgenic rice lines overexpressing OsbHLH058 showed tolerance for iron deficiency and higher iron concentration in seeds. These lines also showed enhanced expression of many iron deficiency-inducible genes involved in iron uptake and translocation under iron-sufficient conditions. Conversely, OsbHLH058 knockdown lines showed susceptibility to iron deficiency and reduced expression of many iron deficiency-inducible genes. OsbHLH059 knockdown lines were also susceptible to iron deficiency, and formed characteristic brownish regions in iron-deficient new leaves. OsbHLH059 knockdown lines also showed reduced expression of many iron deficiency-inducible genes. These results indicate that OsbHLH058 and OsbHLH059 positively regulate major iron deficiency responses in a similar but distinct manner, and that this function may be partially controlled by OsHRZs.

Suppression of Melatonin 2-Hydroxylase Increases Melatonin Production Leading to the Enhanced Abiotic Stress Tolerance against Cadmium, Senescence, Salt, and Tunicamycin in Rice Plants

Melatonin 2-hydroxylase (M2H) catalyzes the conversion of melatonin into 2hydroxymelatonin (2OHM), which is present in plants at a higher concentration than melatonin. Although M2H has been cloned, the in vivo function of its product is unknown. Here, we generated stable T₂ homozygous transgenic rice plants in which expression of endogenous M2H was suppressed (RNAi lines). However, we failed to generate M2H overexpression transgenic rice due to failure of somatic embryogenesis. The M2H transcript level showed a diurnal rhythm with a peak at night concomitantly with the peak concentration of 2OHM. RNAi rice showed a reduced M2H mRNA level and 2OHM and melatonin concentrations. The unexpected decrease in the melatonin concentration was caused by redirection of melatonin into cyclic 3hydroxymelatonin via a detour catabolic pathway. Thus, the decrease in the melatonin concentration in M2H RNAi rice led to slowed seedling growth and delayed germination. By contrast, the transient increase in the melatonin concentration was of greater magnitude in the M2H RNAi than the wild-type rice upon cadmium treatment due to possible suppression of melatonin degradation. Due to its higher concentration of melatonin, the M2H RNAi rice displayed tolerance to senescence, salt, and tunicamycin stresses. Therefore, the increase in the melatonin concentration caused by suppression of melatonin degradation or by overexpression of melatonin biosynthetic genes enhances stress tolerance in rice.

The Rice Actin-Binding Protein RMD Regulates Light-Dependent Shoot Gravitropism

Light and gravity are two key determinants in orientating plant stems for proper growth and development. The organization and dynamics of the actin cytoskeleton are essential for cell biology and critically regulated by actin-binding proteins. However, the role of actin cytoskeleton in shoot negative gravitropism remains controversial. In this work, we report that the actin-binding protein Rice Morphology Determinant (RMD) promotes reorganization of the actin cytoskeleton in rice (Oryza sativa) shoots. The changes in actin organization are associated with the ability of the rice shoots to respond to negative gravitropism. Here, light-grown rmd mutant shoots exhibited agravitropic phenotypes. By contrast, etiolated rmd shoots displayed normal negative shoot gravitropism. Furthermore, we show that RMD maintains an actin configuration that promotes statolith mobility in gravisensing endodermal cells, and for proper auxin distribution in light-grown, but not dark-grown, shoots. RMD gene expression is diurnally controlled and directly repressed by the phytochrome-interacting factor-like protein OsPIL16. Consequently, overexpression of OsPIL16 led to gravisensing and actin patterning defects that phenocopied the rmd mutant. Our findings outline a mechanism that links light signaling and gravity perception for straight shoot growth in rice.

Genome-wide analysis of the rice J-protein family: identification, genomic organization, and expression profiles under multiple stresses

J-proteins which function as molecular chaperone played critical roles in plant growth, development, and response to various environment stresses, but little was reported on this gene family in rice. Here, we identified 115 putative rice J-proteins and classified them into nine major clades (I–IX) according to their phylogenetic relationships. Gene-structure analysis revealed that each member of the same clade has same or similar exon–intron structure, and most rice J-protein genes of clade VII were intronless. Chromosomes mapping suggested that tandem duplication was occurred in evolution. Expression profile showed that the 61 rice J-protein genes were expressed in at least one tissue. The result implied that they could be involved in the process of rice growth and development. The RNA-sequencing data identified 96 differentially expressed genes, 59.38% (57/96), 67.71% (65/96), and 62.50% (60/96) genes were induced by heat stress, drought stress, and salt stress, respectively. The results indicated that J-protein genes could participated in rice response to different stresses. The findings in this study would provide a foundation for further analyzing the function of J-proteins in rice.

Genetic dissection of developmental responses of agro-morphological traits under different doses of nutrient fertilizers using high-density SNP markers

The production and productivity of rice (Oryza sativa L.) are primarily influenced by the application of the critical nutrients nitrogen (N), phosphorus (P), and potassium (K). However, excessive application of these fertilizers is detrimental to the environment and increases the cost of production. Hence, there is a need to develop varieties that simultaneously increase yields under both optimal and suboptimal rates of fertilizer application by maximizing nutrient use efficiency (NuUE). To unravel the hidden genetic variation and understand the molecular and physiological mechanisms of NuUE, three different mapping populations (MPs; BC1F5) derived from three donors (Haoannong, Cheng-Hui 448, and Zhong 413) and recipient Weed Tolerant Rice 1 were developed. A total of three favorable agronomic traits (FATs) were considered as the measure of NuUE. Analysis of variance and descriptive statistics indicated the existence of genetic variation for NuUE and quantitative inheritance of FATs. The genotypic data from single-nucleotide polymorphism (SNP) markers from Tunable Genotyping-By-Sequencing (tGBS) and phenotypic values were used for locating the genomic regions conferring NuUE. A total of 19 quantitative trait loci (QTLs) were detected, out of which 11 QTLs were putative on eight chromosomes, which individually explained 17.02% to 34.85% of the phenotypic variation. Notably, qLC-II_1 and qLC-II_11 detected at zero fertilizer application showed higher performance for LC under zero percentage of NPK fertilizer. The remarkable findings of the present study are that the detected QTLs were associated in building tolerance to low/no nutrient application and six candidate genes on chromosomes 2 and 5 within these putative QTLs were found associated with low nutrient tolerance and related to several physiological and metabolic pathways involved in abiotic stress tolerance. The identified superior introgressed lines (ILs) and trait-associated genetic regions can be effectively used in marker-assisted selection (MAS) for NuUE breeding programs.

Suppression of phospholipase D genes improves chalky grain production by high temperature during the grain-filling stage in rice

High temperature (HT) during the grain developing stage causes deleterious effects on rice quality resulting in mature grains with a chalky appearance. Phospholipase D (PLD) plays an important role in plants, including responses to environmental stresses. OsPLDα1, α3 and β2-knockdown (KD) plants showed decreased production of chalky grains at HT. HT ripening increased H2O2 accumulated in the developing grains. However, the increase was canceled by the knockdown of OsPLDβ2. Expression levels of OsCATA which is one of three rice catalase genes, in developing grains of OsPLDβ2-KD plants at 10 DAF were increased compared with that in vector-controls in HT growth conditions. Overexpression of OsCATA markedly suppressed the production of chalky grains in HT growth conditions. These results suggested that OsPLDβ2 functions as a negative regulator of the induction of OsCATA and is involved in the production of chalky grains in HT growth conditions.

Development of a visible marker trait based on leaf sheath-specific anthocyanin pigmentation applicable to various genotypes in rice

To overcome a limitation to the breeding of autogamous crops, recurrent selection using transgenic male sterility (RSUTMS) has been proposed. In this system, negatively or positively selectable marker traits are required along with dominant transgenic male sterility. Anthocyanin pigmentation is an excellent marker trait. Two regulatory genes for MYB and bHLH and a structural gene for DFR are required for anthocyanin pigmentation in rice. Therefore, to apply anthocyanin pigmentation as a marker trait in various rice genotypes, coordinated expression of the three genes is required. In this study, we developed a leaf sheath-specific promoter and introduced three genes-DFR and C1/Myb, driven by the 35S promoter, and OsB2/bHLH, driven by the leaf sheath-specific promoter-into the rice genome. Leaf sheath-specific pigmentation was confirmed in all seven genotypes tested, which included japonica and indica cultivars. Analysis of genome sequence data from 25 cultivars showed that the strategy of conferring leaf sheath-specific anthocyanin pigmentation by introduction of these three genes would be effective for a wide range of genotypes and will be applicable to RSUTMS.

Recent advances in the characterization of plant transcriptomes in response to drought, salinity, heat, and cold stress

Despite recent advancements in plant molecular biology and biotechnology, providing food security for an increasing world population remains a challenge. Drought (water scarcity), salinity, heat, and cold stress are considered major limiting factors that affect crop production both qualitatively and quantitatively. Therefore, the development of cost-effective and environmentally friendly strategies will be needed to resolve these agricultural problems. This will require a comprehensive understanding of transcriptomic alterations that occur in plants in response to varying levels of environmental stresses, singly and in combination. Here, we briefly discuss the current status and future challenges in plant research related to understanding transcriptional changes that occur in response to drought, salinity, heat, and cold stress.

Heterologous expression of rice RNA-binding glycine-rich (RBG) gene OsRBGD3 in transgenic Arabidopsis thaliana confers cold stress tolerance

Imparting cold stress tolerance to crops is a major challenge in subtropical agriculture. New genes conferring cold tolerance needs to be identified and characterised for sustainable crop production in low-temperature stress affected areas. Here we report functional characterisation of OsRBGD3, classified previously as a class D glycine-rich RNA recognition motif (RRM) containing proteins from a drought-tolerant Indica rice cultivar N22. The gene was isolated by screening yeast one-hybrid library using the minimal promoter region of the OsMYB38 that is necessary for cold stress-responsive expression. OsRBGD3 exhibited cold, drought and salt stress inductive expression in a drought tolerant N22 rice cultivar as compared with susceptible variety IR64. OsRBGD3 was found to be localised to both nuclear and cytoplasmic subcellular destinations. Constitutive overexpression of the OsRBGD3 in transgenic Arabidopsis conferred tolerance to cold stress. ABA sensitivity was also observed in transgenic lines suggesting the regulatory role of this gene in the ABA signalling pathway. OsRBGD3 overexpression also attributed to significant root development and early flowering in transgenics. Hence, OsRBGD3 could be an important target for developing cold tolerant early flowering rice and other crops' genotypes for increasing production in low temperature affected areas.

Detection of stable QTLs for grain protein content in rice (Oryza sativa L.) employing high throughput phenotyping and genotyping platforms

Lack of appropriate donors, non-utilization of high throughput phenotyping and genotyping platforms with high genotype × environment interaction restrained identification of robust QTLs for grain protein content (GPC) in rice. In the present investigation a  BC₃F₄ mapping population was developed using grain protein donor, ARC10075 and high-yielding cultivar Naveen and 190 lines were genotyped using 40 K Affimetrix custom SNP array with the objective to identify stable QTLs for protein content. Three of the identified QTLs, one for GPC (qGPC1.1) and the other two for single grain protein content (qSGPC2.1, qSGPC7.1) were stable over the environments explaining  13%, 14% and 7.8% of the phenotypic variances, respectively. Stability and repeatability of these additive QTLs were supported by the synergistic additive effects of multi-environmental-QTLs. One epistatic-QTL, independent of  the  main effect QTL was detected over the environment for SGPC. A few functional genes governing seed storage protein were hypothesised inside these identified QTLs. The qGPC1.1 was validated by NIR Spectroscopy-based high throughput phenotyping in BC₃F₅ population. Higher glutelin content was estimated in high-protein lines with the introgression of qGPC1.1 in telomeric region of short arm of chromosome 1. This was supported by the postulation of probable candidate gene inside this QTL region encoding glutelin family proteins.

Transcriptome monitoring visualizes growth stage-dependent nutrient status dynamics in rice under field conditions

Crop plants undergo morpho-physiological changes throughout the growth process in response to both the internal and the external environment, and that eventually determine the yield. The system-level adjustment of the morpho-physiological changes has remained largely unclear, however, especially in field conditions. Here, we reveal changes in nutrient status associated with tiller development and soil conditions based on the leaf transcriptome profile of rice (Oryza sativa) throughout the entire period of growth. We performed gene co-expression network analysis and identified three gene sets as indicators for monitoring the internal nitrogen and phosphorus status. Expression profiling reveals that the phosphorus starvation response is expressed during the tillering stage and is then switched off with the transition to nitrogen deficiency. Coincident with phosphorus status dynamics, the level of phosphate in the leaf is demonstrated to be low during the tillering stage and subsequently increases drastically. The phosphorus dynamics are genetically validated by analysing mutants with a defect in phosphorus homeostasis. Notably, we show that nitrogen limitation directly suppresses the phosphorus starvation response. Finally, the phosphorus starvation response is demonstrated to be activated in soil with a high phosphate retention capacity, without the visible phenotypes associated with phosphorus starvation. Our results reveal a growth stage- and soil condition-dependent reaction that requires phosphorus, which is expressed to promote the phosphorus uptake required for developing tillers and is directly adjusted by nitrogen status. A molecular framework for elucidating nutrient status dynamics under field conditions would provide insights into improving crop productivity.

The Rice Alpha-Amylase, Conserved Regulator of Seed Maturation and Germination

Alpha-amylase, the major form of amylase with secondary carbohydrate binding sites, is a crucial enzyme throughout the growth period and life cycle of angiosperm. In rice, alpha-amylase isozymes are critical for the formation of the storage starch granule during seed maturation and motivate the stored starch to nourish the developing seedling during seed germination which will directly affect the plant growth and field yield. Alpha-amylase has not yet been studied intensely to understand its classification, structure, expression trait, and expression regulation in rice and other crops. Among the 10-rice alpha-amylases, most were exclusively expressed in the developing seed embryo and induced in the seed germination process. During rice seed germination, the expression of alpha-amylase genes is known to be regulated negatively by sugar in embryos, however positively by gibberellin (GA) in endosperm through competitively binding to the specific promoter domain; besides, it is also controlled by a series of other abiotic or biotic factors, such as salinity. In this review, we overviewed the research progress of alpha-amylase with focus on seed germination and reflected on how in-depth work might elucidate its regulation and facilitate crop breeding as an efficient biomarker.

A web-based tool for the prediction of rice transcription factor function

Transcription factors (TFs) are an important class of regulatory molecules. Despite their importance, only a small number of genes encoding TFs have been characterized in Oryza sativa (rice), often because gene duplication and functional redundancy complicate their analysis. To address this challenge, we developed a web-based tool called the Rice Transcription Factor Phylogenomics Database (RTFDB) and demonstrate its application for predicting TF function. The RTFDB hosts transcriptome and co-expression analyses. Sources include high-throughput data from oligonucleotide microarray (Affymetrix and Agilent) as well as RNA-Seq-based expression profiles. We used the RTFDB to identify tissue-specific and stress-related gene expression. Subsequently, 273 genes preferentially expressed in specific tissues or organs, 455 genes showing a differential expression pattern in response to 4 abiotic stresses, 179 genes responsive to infection of various pathogens and 512 genes showing differential accumulation in response to various hormone treatments were identified through the meta-expression analysis. Pairwise Pearson correlation coefficient analysis between paralogous genes in a phylogenetic tree was used to assess their expression collinearity and thereby provides a hint on their genetic redundancy. Integrating transcriptome with the gene evolutionary information reveals the possible functional redundancy or dominance played by paralog genes in a highly duplicated genome such as rice. With this method, we estimated a predominant role for 83.3% (65/78) of the TF or transcriptional regulator genes that had been characterized via loss-of-function studies. In this regard, the proposed method is applicable for functional studies of other plant species with annotated genome.

A rice PHD-finger protein OsTITANIA, is a growth regulator that functions through elevating expression of transporter genes for multiple metals

Essential metal absorption for plant growth is mediated predominantly by metal-specific transporters, with expression that responds to the environmental or cellular conditions of specific metals. Differing from metal-specific regulation, we describe a constitutively expressed transcription factor that regulates the transport of several metals in rice. We characterized the rice mutant LOW CADMIUM 5 (LC5), which exhibited reduced growth and accumulation of essential metals (e.g., copper [Cu], zinc [Zn] and manganese [Mn]) in shoots. LC5 was dwarf and developed less tillers than the wild type, but the structure of vasculature was apparently normal. Molecular genetic analysis revealed that the causal gene of LC5 is an ortholog of the transcriptional regulator Arabidopsis thaliana TITANIA (TTA), known as a transcriptional regulator. Expression analyses demonstrated that the OsTTA gene encodes a nucleus-localized protein containing a plant homeodomain-finger (PHD-finger) domain and is expressed ubiquitously in rice plants. RNA sequencing and quantitative PCR analyses revealed that the mRNA accumulation of transporter genes for essential metals, including iron (Fe), Zn, or Mn, were substantially lower in LC5 roots than in the wild type. Unlike known transcription factors of metal transport regulation, OsTTA transcript accumulation was not affected by metal availability. In addition, the growth defect of LC5 was partially rescued by Fe, Zn, or Mn supplementation, respectively. Taken together, OsTTA is a constitutively expressed regulator of multiple metal transporter genes responsible for essential metals delivery to shoots for their normal growth.

Production of high oleic/low linoleic rice by genome editing

Rice bran oil (RBO) contains many valuable healthy constituents, including oleic acid. Improvement of the fatty acid composition in RBO, including an increase in the content of oleic acid, which helps suppress lifestyle disease, would increase health benefits. The enzyme fatty acid desaturase 2 (FAD2) catalyzes the conversion of oleic acid to linoleic acid in plants, and FAD2 mutants exhibit altered oleic and linoleic acid content in many crops. There are three functional FAD2 genes in the genome of rice (Oryza sativa L.), and, of these, expression of the OsFAD2-1 gene is highest in rice seeds. In order to produce high oleic/low linoleic RBO, we attempted to disrupt the OsFAD2-1 gene by CRISPR/Cas9-mediated targeted mutagenesis. We succeeded in the production of homozygous OsFAD2-1 knockout rice plants. The content of oleic acid increased to more than twice that of wild type, and, surprisingly, linoleic acid, a catabolite of oleic acid by FAD2, decreased dramatically to undetectable levels in fad2-1 mutant brown rice seeds. In this study, by genome editing based on genome information, we succeeded in the production of rice whose fatty acid composition is greatly improved. We suggest that CRISPR/Cas9-mediated mutagenesis of a major gene that shows dominant expression in the target tissue could be a powerful tool to improve target traits in a tissue-specific manner.

Development of transgenic male-sterile rice by using anther-specific promoters identified by comprehensive screening of the gene expression profile database 'RiceXPro'

Because genomic selection is designed for the population breeding of allogamous species, a successive outcrossing system is required for efficient use of genomic selection in autogamous crops, such as Oryza sativa L. (rice). Transgenic and dominant male-sterility is a suitable tool for efficient outcrossing of autogamous crops. Though there have been some reports of dominant male-sterile rice developed using transgenic technology, the flowering habit was substandard. Here, to isolate promoters that, when linked to a lethal gene, induce dominant male-sterility while retaining a good flowering habit, we identified 38 candidate genes with anther-specific expression by using the 'RiceXPro' database. We then evaluated the abilities of the near-upstream regions of these genes to induce male-sterility when linked to the lethal gene barnase and introduced into the rice cultivar 'Nipponbare'. Seven of the 38 promoters induced clear dominant male-sterility; promoters expressed in the later stage of anther development induced male-sterility while retaining better flowering habits when compared to ones expressed in the early stage. These seven promoters could potentially be used to facilitate development of an efficient outcross-based breeding system in rice.

Consensus Coexpression Network Analysis Identifies Key Regulators of Flower and Fruit Development in Wild Strawberry

The diploid strawberry, Fragaria vesca, is a developing model system for the economically important Rosaceae family. Strawberry fleshy fruit develops from the floral receptacle and its ripening is nonclimacteric. The external seed configuration of strawberry fruit facilitates the study of seed-to-fruit cross tissue communication, particularly phytohormone biosynthesis and transport. To investigate strawberry fruit development, we previously generated spatial and temporal transcriptome data profiling F. vesca flower and fruit development pre- and postfertilization. In this study, we combined 46 of our existing RNA-seq libraries to generate coexpression networks using the Weighted Gene Co-Expression Network Analysis package in R. We then applied a post-hoc consensus clustering approach and used bootstrapping to demonstrate consensus clustering's ability to produce robust and reproducible clusters. Further, we experimentally tested hypotheses based on the networks, including increased iron transport from the receptacle to the seed postfertilization and characterized a F. vesca floral mutant and its candidate gene. To increase their utility, the networks are presented in a web interface (www.fv.rosaceaefruits.org) for easy exploration and identification of coexpressed genes. Together, the work reported here illustrates ways to generate robust networks optimized for the mining of large transcriptome data sets, thereby providing a useful resource for hypothesis generation and experimental design in strawberry and related Rosaceae fruit crops.

Melatonin is involved in skotomorphogenesis by regulating brassinosteroid biosynthesis in rice plants

Serotonin N-acetyltransferase (SNAT) is the penultimate enzyme in melatonin biosynthesis catalyzing the conversion of serotonin into N-acetylserotonin. In plants, SNAT is encoded by 2 isogenes of which SNAT1 is constitutively expressed and its overexpression confers increased yield in rice. However, the role of SNAT2 remains to be clarified. In contrast to SNAT1, the diurnal rhythm of SNAT2 mRNA expression peaks at night. In this study, transgenic rice plants in which SNAT2 expression were suppressed by RNAi technology showed a decrease in melatonin and a dwarf phenotype with erect leaves, reminiscent of brassinosteroids (BR)-deficient mutants. Of note, the dwarf phenotype was dependent on the presence of dark, suggesting that melatonin is involved in dark growth (skotomorphogenesis). In support of this suggestion, SNAT2 RNAi lines exhibited photomorphogenic phenotypes such as inhibition of internodes and increased expression of light-inducible CAB genes in the dark. The causative gene for the melatonin-mediated BR biosynthetic gene was DWARF4, a rate-limiting BR biosynthetic gene. Exogenous melatonin treatment induced several BR biosynthetic genes, including DWARF4, D11, and RAVL1. As expected from the erect leaves, the SNAT2 RNAi lines produced less BR than the wild type. Our results show for the first time that melatonin is a positive regulator of dark growth or shade outgrowth by regulating BR biosynthesis in plants.

Global analysis of threonine metabolism genes unravel key players in rice to improve the abiotic stress tolerance

The diversity in plant metabolites with improved phytonutrients is essential to achieve global food security and sustainable crop yield. Our study using computational metabolomics genome wide association study (cmGWAS) reports on a comprehensive profiling of threonine (Thr) metabolite in rice. Sixteen abiotic stress responsive (AbSR) – Thr metabolite producing genes (ThrMPG), modulate metabolite levels and play a significant role determining both physiological and nutritional importance of rice. These AbSR-ThrMPG were computationally analysed for their protein properties using OryzaCyc through plant metabolic network analyser. A total of 1373 and 1028 SNPs were involved in complex traits and genomic variations. Comparative mapping of AbSR-ThrMPG revealed the chromosomal colinearity with C4 grass species. Further, computational expression pattern of these genes predicted a differential expression profiling in diverse developmental tissues. Protein interaction of protein coding gene sequences revealed that the abiotic stresses (AbS) are multigenic in nature. In silico expression of AbSR-ThrMPG determined the putative involvement in response to individual AbS. This is the first comprehensive genome wide study reporting on AbSR –ThrMPG analysis in rice. The results of this study provide a pivotal resource for further functional investigation of these key genes in the vital areas of manipulating AbS signaling in rice improvement.

ARMOUR - A Rice miRNA: mRNA Interaction Resource

ARMOUR was developed as A Rice miRNA:mRNA interaction resource. This informative and interactive database includes the experimentally validated expression profiles of miRNAs under different developmental and abiotic stress conditions across seven Indian rice cultivars. This comprehensive database covers 689 known and 1664 predicted novel miRNAs and their expression profiles in more than 38 different tissues or conditions along with their predicted/known target transcripts. The understanding of miRNA:mRNA interactome in regulation of functional cellular machinery is supported by the sequence information of the mature and hairpin structures. ARMOUR provides flexibility to users in querying the database using multiple ways like known gene identifiers, gene ontology identifiers, KEGG identifiers and also allows on the fly fold change analysis and sequence search query with inbuilt BLAST algorithm. ARMOUR database provides a cohesive platform for novel and mature miRNAs and their expression in different experimental conditions and allows searching for their interacting mRNA targets, GO annotation and their involvement in various biological pathways. The ARMOUR database includes a provision for adding more experimental data from users, with an aim to develop it as a platform for sharing and comparing experimental data contributed by research groups working on rice.

Shortened Basal Internodes Encodes a Gibberellin 2-Oxidase and Contributes to Lodging Resistance in Rice

Breeding semi-dwarf varieties to improve lodging resistance has been proven to be enormously successful in increasing grain yield since the advent of the "green revolution." However, the breeding of the majority of semi-dwarf rice varieties in Asia has been dependent mainly on genetic introduction of the mutant alleles of SD1, which encodes a gibberellin (GA) 20-oxidase, OsGA20ox2, for catalyzing GA biosynthesis. Here, we report a new rice lodging-resistance gene, Shortened Basal Internodes (SBI), which encodes a gibberellin 2-oxidase and specifically controls the elongation of culm basal internodes through deactivating GA activity. SBI is predominantly expressed in culm basal internodes. Genetic analyses indicate that SBI is a semi-dominant gene affecting rice height and lodging resistance. SBI allelic variants display different activities and are associated with the height of rice varieties. Breeding with higher activity of the SBI allele generates new rice varieties with improved lodging resistance and increased yield. The discovery of the SBI provides a desirable gene resource for producing semi-dwarf rice phenotypes and offers an effective strategy for breeding rice varieties with enhanced lodging resistance and high yield.

Heterologous expression of rice 9-cis-epoxycarotenoid dioxygenase 4 (OsNCED4) in Arabidopsis confers sugar oversensitivity and drought tolerance

BACKGROUND

The 9-cis-epoxycarotenoid dioxygenases OsNCED4 was cloned from rice in conjunction with OsNCED 1-3 and 5, of which 3 has been shown to function in ABA biosynthesis and alteration of leaf morphology. In higher plants, NCEDs have been shown to be key enzymes controlling ABA biosynthesis and belong to a differentially expressed gene family. Aside from OsNCED3, it remains largely unknown if other OsNCED genes are involved in ABA biosynthesis in rice. Thus, transgenic Arabidopsis plants overexpressing OsNCED4 were generated in the 129B08/nced3 mutant background to explore OsNCED4 function in ABA biosynthesis.

RESULTS

Heterologous expression of OsNCED4 in Arabidopsis increased ABA levels and altered plant size and leaf shape, delayed seed germination, caused sugar oversensitivity in post-germination growth, and enhanced tolerance to drought. The native OsNCED3 and OsNCED4 promoters were expressed in an overlapping pattern in rice seeds and young seedlings, suggesting possible functional redundancy between OsNCED3 and OsNCED4. At the one-leaf stage, similar regulation of OsNCED3 and OsNCED4 gene expression in roots or leaves in response to moderate salt stress (150 mM NaCl) was observed.

CONCLUSION

Like OsNCED3, OsNCED4 is functionally active in ABA biosynthesis in rice. OsNCED3 and OsNCED4 might play redundant roles in controlling ABA biosynthesis in rice, as suggested by GUS staining assay, but this should be further analyzed through complementation of rice NCED knockout mutants.

Meta-expression analysis of unannotated genes in rice and approaches for network construction to suggest the probable roles

This work suggests 2020 potential candidates in rice for the functional annotation of unannotated genes using meta-analysis of anatomical samples derived from microarray and RNA-seq technologies and this information will be useful to identify novel morphological agronomic traits. Although the genome of rice (Oryza sativa) has been sequenced, 14,365 genes are considered unannotated because they lack putative annotation information. According to the Rice Genome Annotation Project Database ( http://rice.plantbiology.msu.edu/ ), the proportion of functionally characterized unannotated genes (0.35%) is quite limited when compared with the approximately 3.9% of annotated genes with assigned putative functions. Researchers require additional information to help them investigate the molecular mechanisms associated with those unannotated genes. To determine which of them might regulate morphological or physiological traits in the rice genome, we conducted a meta-analysis of expression data that covered a wide range of tissue/organ samples. Overall, 2020 genes showed cultivar-, tissue-, or organ-preferential patterns of expression. Representative candidates from featured groups were validated by RT-PCR, and the GUS reporter system was used to validate the expression of genes that were clustered according to their leaf or root preference. Taking a molecular and genetics approach, we examined meta-expression data and found that 127 genes were differentially expressed between japonica and indica rice cultivars. This is potentially significant for future agronomic applications. We also used a T-DNA insertional mutant and performed a co-expression network analysis of Sword shape dwarf1 (SSD1), a gene that regulates cell division. This network was refined via RT-PCR analysis. Our results suggested that SSD1 represses the expression of four genes related to the processes of DNA replication or cell division and provides insight into possible molecular mechanisms. Together, these strategies present a valuable tool for in-depth characterization of currently unannotated genes.

Integrated omics analysis of root-preferred genes across diverse rice varieties including Japonica and indica cultivars

Plant root systems play essential roles in developmental processes, such as the absorption of water and inorganic nutrients, and structural support. Gene expression is affected by growth conditions and the genetic background of plants. To identify highly conserved root-preferred genes in rice across diverse growth conditions and varieties, we used two independent meta-anatomical expression profiles based on a large collection of Affymetrix and Agilent 44K microarray data sets available for public use. We then identified 684 loci with root-preferred expression, which were validated with in silico analysis using both meta-expression profiles. The expression patterns of four candidate genes were confirmed in vivo by monitoring expression of β-glucuronidase under control of the candidate-gene promoters, providing new tools to manipulate agronomic traits associated with roots. We also utilized real-time PCR to examine the root-preferential expression of 14 genes across four rice varieties, including japonica and indica cultivars. Using a database of rice genes with known functions, we identified the reported functions of 39 out of the 684 candidate genes. Sixteen genes are directly involved in root development, while the remaining are involved in processes indirectly related to root development (i.e., soil-stress tolerance or growth retardation). This indicates the importance of our candidate genes for studies on root development and function. Gene ontology enrichment analysis in the 'biological processes' category revealed that root-preferred genes in rice are closely associated with nutrient transport-related genes, indicating that the primary role of roots is the uptake of nutrients from soil. In addition, predicted protein-protein interaction analysis suggested a molecular network for root development composed of 215 interactions associated with 44 root-preferred or root development-related genes. Taken together, our data provide an important foundation for future research on root development in rice.

Pathway Editing Targets for Thiamine Biofortification in Rice Grains

Thiamine deficiency is common in populations consuming polished rice as a major source of carbohydrates. Thiamine is required to synthesize thiamine pyrophosphate (TPP), an essential cofactor of enzymes of central metabolism. Its biosynthesis pathway has been partially elucidated and the effect of overexpression of a few genes such as thi1 and thiC, on thiamine accumulation in rice has been reported. Based on current knowledge, this review focuses on the potential of gene editing in metabolic engineering of thiamine biosynthesis pathway to improve thiamine in rice grains. Candidate genes, suitable for modification of the structural part to evolve more efficient versions of enzymes in the pathway, are discussed. For example, adjacent cysteine residues may be introduced in the catalytic domain of thi4 to improve the turn over activity of thiamine thiazole synthase 2. Motif specific editing to modify promoter regulatory regions of genes is discussed to modulate gene expression. Editing cis acting regulatory elements in promoter region can shift the expression of transporters and thiamine binding proteins to endosperm. This can enhance dietary availability of thiamine from rice grains. Differential transcriptomics on rice varieties with contrasting grain thiamine and functional genomic studies will identify more strategic targets for editing in future. Developing functionally enhanced foods by biofortification is a sustainable approach to make diets wholesome.

Sequencing of bulks of segregants allows dissection of genetic control of amylose content in rice

Amylose content (AC) is a key quality trait in rice. A cross between Oryza glaberrima (African rice) and Oryza sativa (Asian rice) segregating for AC was analysed by sequencing bulks of individuals with high and low AC. SNP associated with the granule bound starch synthase (GBSS1) locus on chromosome 6 were polymorphic between the bulks. In particular, a G/A SNP that would result in an Asp to Asn mutation was identified. This amino acid substitution may be responsible for differences in GBSS activity as it is adjacent to a disulphide linkage conserved in all grass GBSS proteins. Other polymorphisms in genomic regions closely surrounding this variation may be the result of linkage drag. In addition to the variant in the starch biosynthesis gene, SNP on chromosomes 1 and 11 linked to AC was also identified. SNP was found in the genes encoding the NAC and CCAAT-HAP5 transcription factors that have previously been linked to starch biosynthesis. This study has demonstrated that the approach of sequencing bulks was able to identify genes on different chromosomes associated with this complex trait.

RiceMetaSys for salt and drought stress responsive genes in rice: a web interface for crop improvement

Background

Genome-wide microarray has enabled development of robust databases for functional genomics studies in rice. However, such databases do not directly cater to the needs of breeders. Here, we have attempted to develop a web interface which combines the information from functional genomic studies across different genetic backgrounds with DNA markers so that they can be readily deployed in crop improvement. In the current version of the database, we have included drought and salinity stress studies since these two are the major abiotic stresses in rice.

Results

RiceMetaSys, a user-friendly and freely available web interface provides comprehensive information on salt responsive genes (SRGs) and drought responsive genes (DRGs) across genotypes, crop development stages and tissues, identified from multiple microarray datasets. ‘Physical position search’ is an attractive tool for those using QTL based approach for dissecting tolerance to salt and drought stress since it can provide the list of SRGs and DRGs in any physical interval. To identify robust candidate genes for use in crop improvement, the ‘common genes across varieties’ search tool is useful. Graphical visualization of expression profiles across genes and rice genotypes has been enabled to facilitate the user and to make the comparisons more impactful. Simple Sequence Repeat (SSR) search in the SRGs and DRGs is a valuable tool for fine mapping and marker assisted selection since it provides primers for survey of polymorphism. An external link to intron specific markers is also provided for this purpose. Bulk retrieval of data without any limit has been enabled in case of locus and SSR search.

Conclusions

The aim of this database is to facilitate users with a simple and straight-forward search options for identification of robust candidate genes from among thousands of SRGs and DRGs so as to facilitate linking variation in expression profiles to variation in phenotype.

Database URL: http://14.139.229.201

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1846-y) contains supplementary material, which is available to authorized users.

Construction and Optimization of a Large Gene Coexpression Network in Maize Using RNA-Seq Data

With the emergence of massively parallel sequencing, genomewide expression data production has reached an unprecedented level. This abundance of data has greatly facilitated maize research, but may not be amenable to traditional analysis techniques that were optimized for other data types. Using publicly available data, a gene coexpression network (GCN) can be constructed and used for gene function prediction, candidate gene selection, and improving understanding of regulatory pathways. Several GCN studies have been done in maize (Zea mays), mostly using microarray datasets. To build an optimal GCN from plant materials RNA-Seq data, parameters for expression data normalization and network inference were evaluated. A comprehensive evaluation of these two parameters and a ranked aggregation strategy on network performance, using libraries from 1266 maize samples, were conducted. Three normalization methods and 10 inference methods, including six correlation and four mutual information methods, were tested. The three normalization methods had very similar performance. For network inference, correlation methods performed better than mutual information methods at some genes. Increasing sample size also had a positive effect on GCN. Aggregating single networks together resulted in improved performance compared to single networks.

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

Flag leaves (FL) and second leaves (SL) in rice show differential aging patterns during monocarpic senescence. Coordination of aging programs between FL and SL is important for grain yield and quality. However, the molecular bases for differential aging programs between FL and SL have not been systematically explored in rice. Here, we performed mRNA-sequencing of FL and SL at six time points during grain-filling and identified four molecular bases for differential aging programs between FL and SL: phenylpropanoid biosynthesis, photosynthesis, amino acid (AA) transport, and hormone response. Of them, photosynthesis (carbon assimilation) and AA transport (nitrogen remobilization) predominantly occurred in FL and SL, respectively, during grain-filling. Unlike other molecular bases, AA transport showed consistent differential expression patterns between FL and SL in independent samples. Moreover, long-distance AA transporters showed invariant differential expression patterns between FL and SL after panicle removal, which was consistent to invariant differential nitrogen contents between FL and SL after panicle removal. Therefore, our results suggest that the supplies of carbon and nitrogen to seeds is functionally segregated between FL and SL and that long-distance AA transport is an invariant core program for high nitrogen remobilization in SL.

An integrated hyperspectral imaging and genome-wide association analysis platform provides spectral and genetic insights into the natural variation in rice

With progress of genetic sequencing technology, plant genomics has experienced rapid development and subsequently triggered the progress of plant phenomics. In this study, a high-throughput hyperspectral imaging system (HHIS) was developed to obtain 1,540 hyperspectral indices at whole-plant level during tillering, heading, and ripening stages. These indices were used to quantify traditional agronomic traits and to explore genetic variation. We performed genome-wide association study (GWAS) of these indices and traditional agronomic traits in a global rice collection of 529 accessions. With the genome-level suggestive P-value threshold, 989 loci were identified. Of the 1,540 indices, we detected 502 significant indices (designated as hyper-traits) that exhibited phenotypic and genetic relationship with traditional agronomic traits and had high heritability. Many hyper-trait-associated loci could not be detected using traditional agronomic traits. For example, we identified a candidate gene controlling chlorophyll content (Chl). This gene, which was not identified based on Chl, was significantly associated with a chlorophyll-related hyper-trait in GWAS and was demonstrated to control Chl. Moreover, our study demonstrates that red edge (680-760 nm) is vital for rice research for phenotypic and genetic insights. Thus, combination of HHIS and GWAS provides a novel platform for dissection of complex traits and for crop breeding.

Transcriptional regulatory network of WOX11 is involved in the control of crown root development, cytokinin signals, and redox in rice

The rice root system is mainly composed of post-embryonic shoot-borne roots called crown roots. WOX11, encoding a WUSCHEL-related homeobox domain transcription factor, is a key regulator of crown root growth and development in rice (Oryza sativa. L). In addition to specifically activating crown root development, WOX11 is also involved in lateral root initiation, root hair formation, and abiotic stresses. However, the gene regulatory network downstream of WOX11 remains largely unknown. Here, we studied the transcriptome of wox11 root tips by RNA-Seq and determined direct WOX11-binding targets by bioinformatic and biochemical analysis. The transcriptomic analysis revealed 664 differentially expressed genes, which covered a wide range of functions related to root development, cytokinin homeostasis/signaling, stress response, and redox metabolic processes. Bioinformatic analysis also revealed that the WOX11-binding motif was distributed over 41% (273/664) of the differentially expressed genes, and was mostly enriched in the promoter and intron regions. We used qRT-PCR and/or in situ hybridization to confirm co-expression of some of the WOX11-regulated genes in crown root development. We also used electrophoretic mobility shift assay and chromatin immunoprecipitation with anti-WOX11 antibody to validate direct regulation of these genes by WOX11. The analysis identified several genes that acted downstream of WOX11 in controlling crown root formation, cytokinin signaling, stress response, and redox metabolism. This work built a hierarchical regulatory model of WOX11 in rice crown root development.

Rice Expression Database (RED): An integrated RNA-Seq-derived gene expression database for rice

Rice is one of the most important stable food as well as a monocotyledonous model organism for the plant research community. Here, we present RED (Rice Expression Database; http://expression.ic4r.org), an integrated database of rice gene expression profiles derived entirely from RNA-Seq data. RED features a comprehensive collection of 284 high-quality RNA-Seq experiments, integrates a large number of gene expression profiles and covers a wide range of rice growth stages as well as various treatments. Based on massive expression profiles, RED provides a list of housekeeping and tissue-specific genes and dynamically constructs co-expression networks for gene(s) of interest. Besides, it provides user-friendly web interfaces for querying, browsing and visualizing expression profiles of concerned genes. Together, as a core resource in BIG Data Center, RED bears great utility for characterizing the function of rice genes and better understanding important biological processes and mechanisms underlying complex agronomic traits in rice.

High Temperature-Induced Expression of Rice α-Amylases in Developing Endosperm Produces Chalky Grains

Global warming impairs grain filling in rice and reduces starch accumulation in the endosperm, leading to chalky-appearing grains, which damages their market value. We found previously that high temperature-induced expression of starch-lytic α-amylases during ripening is crucial for grain chalkiness. Because the rice genome carries at least eight functional α-amylase genes, identification of the α-amylase(s) that contribute most strongly to the production of chalky grains could accelerate efficient breeding. To identify α-amylase genes responsible for the production of chalky grains, we characterized the histological expression pattern of eight α-amylase genes and the influences of their overexpression on grain appearance and carbohydrate components through a series of experiments with transgenic rice plants. The promoter activity of most α-amylase genes was elevated to various extents at high temperature. Among them, the expression of Amy1A and Amy3C was induced in the internal, especially basal to dorsal, region of developing endosperm, whereas that of Amy3D was confined near the ventral aleurone. These regions coincided with the site of occurrence of chalkiness, which was in clear contrast to conventionally known expression patterns of the enzyme in the scutellum and aleurone during seed germination. Furthermore, overexpression of α-amylase genes, except for Amy3E, in developing endosperm produced various degrees of chalky grains without heat exposure, whereas that of Amy3E yielded normal translucent grains, as was the case in the vector control, even though Amy3E-overexpressing grains contained enhanced α-amylase activities. The weight of the chalky grains was decreased due to reduced amounts of starch, and microscopic observation of the chalky part of these grains revealed that their endosperm consisted of loosely packed round starch granules that had numerous pits on their surface, confirming the hydrolysis of the starch reserve by α-amylases. Moreover, the chalky grains contained increased amounts of soluble sugars including maltooligosaccharides at the expense of starch. The integrated analyses proposed that expression of Amy1A, Amy3C, and Amy3D at the specific regions of the developing endosperm could generate the chalkiness. This finding provides the fundamental knowledge to narrow down the targets for the development of high temperature-tolerant premium rice.

Assessment of Five Chilling Tolerance Traits and GWAS Mapping in Rice Using the USDA Mini-Core Collection

Rice (Oryza sativa L.) is often exposed to cool temperatures during spring planting in temperate climates. A better understanding of genetic pathways regulating chilling tolerance will enable breeders to develop varieties with improved tolerance during germination and young seedling stages. To dissect chilling tolerance, five assays were developed; one assay for the germination stage, one assay for the germination and seedling stage, and three for the seedling stage. Based on these assays, five chilling tolerance indices were calculated and assessed using 202 O. sativa accessions from the Rice Mini-Core (RMC) collection. Significant differences between RMC accessions made the five indices suitable for genome-wide association study (GWAS) based quantitative trait loci (QTL) mapping. For young seedling stage indices, japonica and indica subspecies clustered into chilling tolerant and chilling sensitive accessions, respectively, while both subspecies had similar low temperature germinability distributions. Indica subspecies were shown to have chilling acclimation potential. GWAS mapping uncovered 48 QTL at 39 chromosome regions distributed across all 12 rice chromosomes. Interestingly, there was no overlap between the germination and seedling stage QTL. Also, 18 QTL and 32 QTL were in regions discovered in previously reported bi-parental and GWAS based QTL mapping studies, respectively. Two novel low temperature seedling survivability (LTSS)-QTL, qLTSS3-4 and qLTSS4-1, were not in a previously reported QTL region. QTL with strong effect alleles identified in this study will be useful for marker assisted breeding efforts to improve chilling tolerance in rice cultivars and enhance gene discovery for chilling tolerance.

Jasmonate regulates juvenile-to-adult phase transition in rice

Juvenile-to-adult phase transition is an important shift for the acquisition of adult vegetative characteristics and subsequent reproductive competence. We identified a recessive precocious (pre) mutant exhibiting a long leaf phenotype in rice. The long leaf phenotype is conspicuous in the second to the fourth leaves, which are juvenile and juvenile-to-adult transition leaves. We found that morphological and physiological traits, such as midrib formation, shoot meristem size, photosynthetic rate and plastochron, in juvenile and juvenile-to-adult transition stages of the pre mutant have precociously acquired adult characteristics. In agreement with these results, expression patterns of miR156 and miR172, which are microRNAs regulating phase change, support the accelerated juvenile-to-adult phase change in the pre mutant. The mutated gene encodes an allene oxide synthase (OsAOS1), which is a key enzyme for the biosynthesis of jasmonic acid (JA). The pre mutant showed a low level of JA and enhanced sensitivity to gibberellic acid, which promotes the phase change in some plant species. We also show that prolonged plastochron in the pre mutant is caused by accelerated PLASTOCHRON1 (PLA1) function. The present study reveals a substantial role of JA as a negative regulator of vegetative phase change.

Combining mapping of physiological quantitative trait loci and transcriptome for cold tolerance for counteracting male sterility induced by low temperatures during reproductive stage in rice

Male sterility induced by low temperatures (LTs) during the reproductive stage is a major constraint for temperate zone rice. To detect physiological quantitative trait loci (QTLs), we modeled genotypic variation in the physiological processes involved in low temperature spikelet sterility on the basis of anther length (AL), a proxy for microspore and pollen grain number per anther. The model accounted for 83% of the genotypic variation in potential AL at normal temperature and the ability to maintain AL at LT. We tested the model on 208 recombinant inbred lines of cold-tolerant 'Tohoku-PL3' (PL3) × cold-sensitive 'Akihikari' (AH) for 2 years. QTLs for spikelet fertility (FRT) at LT were detected on chromosomes 5 (QTL for Cold Tolerance at Reproductive stage, qCTR5) and 12 (qCTR12). qCTR12 was annotated with the ability to maintain AL under LTs. qCTR5 was in a region shared with QTLs for culm length and heading date. Genome-wide expression analysis showed 798 genes differentially expressed in the spikelets between the parents at LTs. Of these, 12 were near qCTR5 and 23 were near qCTR12. Gene expression analysis confirmed two candidate genes for qCTR5 (O-methyltransferase ZRP4, Os05g0515600; beta-1,3-glucanase-like protein, Os05g0535100) and one for qCTR12 (conserved hypothetical protein, Os12g0550600). Nucleotide polymorphisms (21 deletions, 2 insertions and 10 single nucleotide polymorphisms) in PL3 were found near the candidate conserved hypothetical protein (Os12g0550600) and upstream in PL3, but not in AH. Haplotype analysis revealed that this gene came from 'Kuchum'. The combination of mapping physiological QTLs with gene expression analysis can be extended to identify other genes for abiotic stress response in cereals.

Rice OsVAMP714, a membrane-trafficking protein localized to the chloroplast and vacuolar membrane, is involved in resistance to rice blast disease

Membrane trafficking plays pivotal roles in many cellular processes including plant immunity. Here, we report the characterization of OsVAMP714, an intracellular SNARE protein, focusing on its role in resistance to rice blast disease caused by the fungal pathogen Magnaporthe oryzae. Disease resistance tests using OsVAMP714 knockdown and overexpressing rice plants demonstrated the involvement of OsVAMP714 in blast resistance. The overexpression of OsVAMP7111, whose product is highly homologous to OsVAMP714, did not enhance blast resistance to rice, implying a potential specificity of OsVAMP714 to blast resistance. OsVAMP714 was localized to the chloroplast in mesophyll cells and to the cellular periphery in epidermal cells of transgenic rice plant leaves. We showed that chloroplast localization is critical for the normal OsVAMP714 functioning in blast resistance by analyzing the rice plants overexpressing OsVAMP714 mutants whose products did not localize in the chloroplast. We also found that OsVAMP714 was located in the vacuolar membrane surrounding the invasive hyphae of M. oryzae. Furthermore, we showed that OsVAMP714 overexpression promotes leaf sheath elongation and that the first 19 amino acids, which are highly conserved between animal and plant VAMP7 proteins, are crucial for normal rice plant growths. Our studies imply that the OsVAMP714-mediated trafficking pathway plays an important role in rice blast resistance as well as in the vegetative growth of rice.

Are Symbiotic Methanotrophs Key Microbes for N Acquisition in Paddy Rice Root?

The relationships between biogeochemical processes and microbial functions in rice (Oryza sativa) paddies have been the focus of a large number of studies. A mechanistic understanding of methane–nitrogen (CH₄–N) cycle interactions is a key unresolved issue in research on rice paddies. This minireview is an opinion paper for highlighting the mechanisms underlying the interactions between biogeochemical processes and plant-associated microbes based on recent metagenomic, metaproteomic, and isotope analyses. A rice symbiotic gene, relevant to rhizobial nodulation and mycorrhization in plants, likely accommodates diazotrophic methanotrophs or the associated bacterial community in root tissues under low-N fertilizer management, which may permit rice plants to acquire N via N₂ fixation. The amount of N fixed in rice roots was previously estimated to be approximately 12% of plant N based on measurements of ¹⁵N natural abundance in a paddy field experiment. Community analyses also indicate that methanotroph populations in rice roots are susceptible to environmental conditions such as the microclimate of rice paddies. Therefore, CH₄ oxidation by methanotrophs is a driving force in shaping bacterial communities in rice roots grown in CH₄-rich environments. Based on these findings, we propose a hypothesis with unanswered questions to describe the interplay between rice plants, root microbiomes, and their biogeochemical functions (CH₄ oxidation and N₂ fixation).

Amyloplast Membrane Protein SUBSTANDARD STARCH GRAIN6 Controls Starch Grain Size in Rice Endosperm

Starch is a biologically and commercially important polymer of glucose. Starch is organized into starch grains (SGs) inside amyloplasts. The SG size differs depending on the plant species and is one of the most important factors for industrial applications of starch. There is limited information on genetic factors regulating SG sizes. In this study, we report the rice (Oryza sativa) mutant substandard starch grain6 (ssg6), which develops enlarged SGs in endosperm. Enlarged SGs are observed starting at 3 d after flowering. During endosperm development, a number of smaller SGs appear and coexist with enlarged SGs in the same cells. The ssg6 mutation also affects SG morphologies in pollen. The SSG6 gene was identified by map-based cloning and microarray analysis. SSG6 encodes a protein homologous to aminotransferase. SSG6 differs from other rice homologs in that it has a transmembrane domain. SSG6-green fluorescent protein is localized in the amyloplast membrane surrounding SGs in rice endosperm, pollen, and pericarp. The results of this study suggest that SSG6 is a novel protein that controls SG size. SSG6 will be a useful molecular tool for future starch breeding and applications.

Detection of QTLs for cold tolerance of rice cultivar 'Kuchum' and effect of QTL pyramiding

A QTL for cold tolerance at the booting stage of rice cultivar 'Kuchum' was detected and delimited into a 1.36 Mb region, and a cold-tolerant line was developed by QTL pyramiding. Low temperature in summer causes pollen sterility in rice, resulting in a serious loss of yield. The second most widely grown rice cultivar in Japan, 'Hitomebore', has been developed as a cultivar highly tolerant to low temperature at the booting stage. However, even 'Hitomebore' exhibits sterility at a temperature lower than 18.5 °C. Further improvement of cold tolerance of rice is required. In the present study, QTLs for cold tolerance in a Bhutanese rice variety, 'Kuchum', were analyzed using backcrossed progenies and a major QTL, named qCT-4, was detected on chromosome 4. Evaluating cold tolerance of seven types of near isogenic lines having 'Kuchum' alleles around qCT-4 with a 'Hitomebore' genetic background, qCT-4 was delimited to a region of ca. 1.36 Mb between DNA markers 9_1 and 10_13. Homozygous 'Kuchum' alleles at qCT-4 showed an effect of increasing seed fertility by ca. 10 % under cold-water treatment. Near isogenic lines of 'Hitomebore' having 'Silewah' alleles of Ctb1 and Ctb2 and a 'Hokkai PL9' allele of qCTB8 did not exhibit higher cold tolerance than that of 'Hitomebore'. On the other hand, a qLTB3 allele derived from a Chinese cultivar 'Lijiangxintuanheigu' increased cold tolerance of 'Hitomebore', and pyramiding of the qCT-4 allele and the qLTB3 allele further increased seed fertility under cold-water treatment. Since NILs of 'Hitomebore' with the 'Kuchum' allele of qCT-4 were highly similar to 'Hitomebore' in other agronomic traits, the qCT-4 allele is considered to be useful for developing a cold-tolerant cultivar.

Genome-wide analysis of spatiotemporal gene expression patterns during early embryogenesis in rice

Embryogenesis in rice is different from that of most dicotolydonous plants in that it shows a non-stereotypic cell division pattern, formation of dorsal-ventral polarity, and endogenous initiation of the radicle. To reveal the transcriptional features associated with developmental events during rice early embryogenesis, we used microarray analysis coupled with laser microdissection to obtain both spatial and temporal transcription profiles. Our results allowed us to determine spatial expression foci for each expressed gene in the globular embryo, which revealed the importance of phytohormone-related genes and a suite of transcription factors to early embryogenesis. Our analysis showed the polarized expression of a small number of genes along the apical-basal and dorsal-ventral axes in the globular embryo, which tended to fluctuate in later developmental stages. We also analyzed gene expression patterns in the early globular embryo and how this relates to expression in embryonic organs at later stages. We confirmed the accuracy of the expression patterns found by microarray analysis of embryo subdomains using in situ hybridization. Our study identified homologous genes from Arabidopsis thaliana with known functions in embryogenesis in addition to unique and uncharacterized genes that show polarized expression patterns during embryogenesis. The results of this study are presented in a database to provide a framework for spatiotemporal gene expression during rice embryogenesis, to serve as a resource for future functional analysis of genes, and as a basis for comparative studies of plant embryogenesis.

Nucleolar DEAD-Box RNA Helicase TOGR1 Regulates Thermotolerant Growth as a Pre-rRNA Chaperone in Rice

Plants have evolved a considerable number of intrinsic tolerance strategies to acclimate to ambient temperature increase. However, their molecular mechanisms remain largely obscure. Here we report a DEAD-box RNA helicase, TOGR1 (Thermotolerant Growth Required1), prerequisite for rice growth themotolerance. Regulated by both temperature and the circadian clock, its expression is tightly coupled to daily temperature fluctuations and its helicase activities directly promoted by temperature increase. Located in the nucleolus and associated with the small subunit (SSU) pre-rRNA processome, TOGR1 maintains a normal rRNA homeostasis at high temperature. Natural variation in its transcript level is positively correlated with plant height and its overexpression significantly improves rice growth under hot conditions. Our findings reveal a novel molecular mechanism of RNA helicase as a key chaperone for rRNA homeostasis required for rice thermotolerant growth and provide a potential strategy to breed heat-tolerant crops by modulating the expression of TOGR1 and its orthologs.

Global warming is increasingly posing negative impacts on crop productivity. In this study, we report a nucleolar-located RNA helicase TOGR1 for thermotolerant growth in rice. TOGR1 maintains pre-rRNA homeostasis under high temperature by securing a proper pre-rRNA structure via elevating its helicase activity. Its expression is high temperature inducible with an afternoon peak expression, consistent with a high temperature anticipation of the circadian clock. Transcriptome analysis revealed that TOGR1 is essential in coordinating primary metabolisms to support thermotolerant growth. Importantly, an enhanced expression of TOGR1 significantly increased biomass of rice. Our findings reveal a novel role of a RNA helicase in thermotolerance and provide a potential strategy to breed heat-tolerant rice cultivars and possibly other heat-tolerant crops.

Regulatory role of FZP in the determination of panicle branching and spikelet formation in rice

FRIZZLE PANICLE (FZP) and RFL/ABERRANT PANICLE ORGANIZATION 2 (APO2) play important roles in regulating the ABCDE floral organ identity genes. However, the relationships among FZP and these floral identity genes in the regulation of panicle formation remain unclear. Here, we used the novel mutant fzp-11, wild-type and FZP-overexpressing plants to compare the expression of these genes during panicle development by real-time PCR and in situ hybridization. The results indicate that FZP is a major negative regulator of RFL/APO2 and determines the transition from panicle branching to spikelet formation. Moreover, overexpression of FZP severely represses axillary meristem formation in both the vegetative and reproductive phases and the outgrowth of secondary branches in panicle. FZP overexpression positively regulates the expression of a subset of the class B genes, AGL6 genes (OsMADS6 and OsMADS17) as well as class E genes (OsMADS1, OsMADS7 and OsMADS8) in floral meristem (FM). Thus, it suggested that FZP could specify floral organ identity by regulating the related OsMADS-box genes.

Comparative Transcriptional Profiling of Melatonin Synthesis and Catabolic Genes Indicates the Possible Role of Melatonin in Developmental and Stress Responses in Rice

As a well-known animal hormone, melatonin (N-acetyl-5-methoxytryptamine) is also involved in multiple plant biological processes, especially in various stress responses. Rice is one of the most important crops, and melatonin is taken in by many people everyday from rice. However, the transcriptional profiling of melatonin-related genes in rice is largely unknown. In this study, the expression patterns of 11 melatonin related genes in rice in different periods, tissues, in response to different treatments were synthetically analyzed using published microarray data. These results suggest that the melatonin-related genes may play important and dual roles in rice developmental stages. We highlight the commonly regulation of rice melatonin-related genes by abscisic acid (ABA), jasmonic acid (JA), various abiotic stresses and pathogen infection, indicating the possible role of these genes in multiple stress responses and underlying crosstalks of plant hormones, especially ABA and JA. Taken together, this study may provide insight into the association among melatonin biosynthesis and catabolic pathway, plant development and stress responses in rice. The profile analysis identified candidate genes for further functional characterization in circadian rhythm and specific stress responses.