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

Exploration of the pearl millet phospholipase gene family to identify potential candidates for grain quality traits

BACKGROUND

Phospholipases constitute a diverse category of enzymes responsible for the breakdown of phospholipids. Their involvement in signal transduction with a pivotal role in plant development and stress responses is well documented.

RESULTS

In the present investigation, a thorough genome-wide analysis revealed that the pearl millet genome contains at least 44 phospholipase genes distributed across its 7 chromosomes, with chromosome one harbouring the highest number of these genes. The synteny analysis suggested a close genetic relationship of pearl millet phospholipases with that of foxtail millet and sorghum. All identified genes were examined to unravel their gene structures, protein attributes, cis-regulatory elements, and expression patterns in two pearl millet genotypes contrasting for rancidity. All the phospholipases have a high alpha-helix content and distorted regions within the predicted secondary structures. Moreover, many of these enzymes possess binding sites for both metal and non-metal ligands. Additionally, the putative promoter regions associated with these genes exhibit multiple copies of cis-elements specifically responsive to biotic and abiotic stress factors and signaling molecules. The transcriptional profiling of 44 phospholipase genes in two genotypes contrasting for rancidity across six key tissues during pearl millet growth revealed a predominant expression in grains, followed by seed coat and endosperm. Specifically, the genes PgPLD-alpha1-1, PgPLD-alpha1-5, PgPLD-delta1-7a, PgPLA1-II-1a, and PgPLD-delta1-2a exhibited notable expression in grains of both the genotypes while showing negligible expression in the other five tissues. The sequence alignment of putative promoters revealed several variations including SNPs and InDels. These variations resulted in modifications to the corresponding cis-acting elements, forming distinct transcription factor binding sites suggesting the transcriptional-level regulation for these five genes in pearl millet.

CONCLUSIONS

The current study utilized a genome-wide computational analysis to characterize the phospholipase gene family in pearl millet. A comprehensive expression profile of 44 phospholipases led to the identification of five grain-specific candidates. This underscores a potential role for at least these five genes in grain quality traits including the regulation of rancidity in pearl millet. Therefore, this study marks the first exploration highlighting the possible impact of phospholipases towards enhancing agronomic traits in pearl millet.

OsbZIP1 regulates phosphorus uptake and nitrogen utilization, contributing to improved yield

Increasing nutrient uptake and use efficiency in plants can contribute to improved crop yields and reduce the demand for fertilizers in crop production. In this study, we characterized a rice mutant, 88n which showed long roots under low nitrogen (N) or phosphorus (P) conditions. Low expression levels of N transporter genes were observed in 88n root, and total N concentration in 88n shoots were decreased, however, C concentrations and shoot dry weight in 88n were comparable to that in WT. Therefore, 88n showed high nitrogen utilization efficiency (NUtE). mRNA accumulation of Pi transporter genes was higher in 88n roots, and Pi concentration and uptake activity were higher in 88n than in WT. Therefore, 88n also showed high phosphorus uptake efficiency (PUpE). Molecular genetic analysis revealed that the causal gene of 88n phenotypes was OsbZIP1, a monocot-specific ortholog of the A. thaliana bZIP transcription factor HY5. Similar to the hy5 mutant, chlorophyll content in roots was decreased and root angle was shallower in 88n than in WT. Finally, we tested the yield of 88n in paddy fields over 3 years because 88n mutant plants showed higher PUpE and NUtE activity and different root architecture at the seedling stage. 88n showed large panicles and increased panicle weight/plant. Taken together, a mutation in OsbZIP1 could contribute to improved crop yields.

Utilization of the Winkler scale of plants using big data temperature presented by the Korea Meteorological Administration

Introduction

Rice is an important food source that can provide a stable supply of calories for most people around the world. However, owing to the recent rapid temperature rise, we are facing social issues related to the increase in the Winkler scale. In this study, a strategy for screening potential candidate genes related to the yield according to the Winkler scale is presented, and the possibility of using a candidate gene identified through sequence haplotype and homology analysis as a breeding source is suggested.

Methods

QTL for the Winkler scale was identified using a population of 120 double haploids derived from a cross between Cheongchoneg, Indica, and Nagdong, Japonica.

Results and discussion

A total of 79 candidate genes were detected in the identified QTL region, and OsHAq8 was finally screened. Through haplotype analysis, OsHAq8 was derived from the Indica group and orthologous to Graminae's activator of Hsp90 ATPase, suggesting that it is a candidate gene involved in yield according to temperature during the growing period. The expression level of OsHAq8 increased as the Winkler scale increased. The findings of this study can serve as a crucial indicator for predicting harvest time and grain quality while achieving a stable yield through marker selection and adaptation to climate change. Climate change occurs more frequently. In these situations, it is very important to predict harvest time and apply relevant candidate genes to breeding. The candidate genes presented in this study can be effectively applied to rice breeding in preparation for climate change.

Analysis of Germin-like protein genes family in Vitis vinifera (VvGLPs) using various in silico approaches

Abstract Germin-like proteins (GLPs) play an important role against various stresses. Vitis vinifera L. genome contains 7 GLPs; many of them are functionally unexplored. However, the computational analysis may provide important new insight into their function. Currently, physicochemical properties, subcellular localization, domain architectures, 3D structures, N-glycosylation & phosphorylation sites, and phylogeney of the VvGLPs were investigated using the latest computational tools. Their functions were predicted using the Search tool for the retrieval of interacting genes/proteins (STRING) and Blast2Go servers. Most of the VvGLPs were extracellular (43%) in nature but also showed periplasmic (29%), plasma membrane (14%), and mitochondrial- or chloroplast-specific (14%) expression. The functional analysis predicted unique enzymatic activities for these proteins including terpene synthase, isoprenoid synthase, lipoxygenase, phosphate permease, receptor kinase, and hydrolases generally mediated by Mn+ cation. VvGLPs showed similarity in the overall structure, shape, and position of the cupin domain. Functionally, VvGLPs control and regulate the production of secondary metabolites to cope with various stresses. Phylogenetically VvGLP1, -3, -4, -5, and VvGLP7 showed greater similarity due to duplication while VvGLP2 and VvGLP6 revealed a distant relationship. Promoter analysis revealed the presence of diverse cis-regulatory elements among which CAAT box, MYB, MYC, unnamed-4 were common to all of them. The analysis will help to utilize VvGLPs and their promoters in future food programs by developing resistant cultivars against various biotic (Erysiphe necator and in Powdery Mildew etc.) and abiotic (Salt, drought, heat, dehydration, etc.) stresses.

Meta-QTL s and haplotypes for efficient zinc biofortification of rice

Biofortification of rice with improved grain zinc (Zn) content is the most sustainable and cost-effective approach to address Zn malnutrition in Asia. Genomics-assisted breeding using precise and consistent Zn quantitative trait loci (QTLs), genes, and haplotypes can fast-track the development of Zn biofortified rice varieties. We conducted the meta-analysis of 155 Zn QTLs reported from 26 different studies. Results revealed 57 meta-QTLs with a significant reduction of 63.2% and 80% in the number and confidence interval of the Zn QTLs, respectively. Meta-quantitative trait loci (MQTLs) regions were found to be enriched with diverse metal homeostasis genes; at least 11 MQTLs were colocated with 20 known major genes involved in the production of root exudates, metal uptake, transport, partitioning, and loading into grains in rice. These genes were differentially expressed in vegetative and reproductive tissues, and a complex web of interactions were observed among them. We identified superior haplotypes and their combinations for nine candidate genes (CGs), and the frequency and allelic effects of superior haplotypes varied in different subgroups. The precise MQTLs with high phenotypic variance, CGs, and superior haplotypes identified in our study are useful for an efficient Zn biofortification of rice and to ensure Zn as an essential component of all the future rice varieties through mainstreaming of Zn breeding.

Integrating omics databases for enhanced crop breeding

Crop plant breeding involves selecting and developing new plant varieties with desirable traits such as increased yield, improved disease resistance, and enhanced nutritional value. With the development of high-throughput technologies, such as genomics, transcriptomics, and metabolomics, crop breeding has entered a new era. However, to effectively use these technologies, integration of multi-omics data from different databases is required. Integration of omics data provides a comprehensive understanding of the biological processes underlying plant traits and their interactions. This review highlights the importance of integrating omics databases in crop plant breeding, discusses available omics data and databases, describes integration challenges, and highlights recent developments and potential benefits. Taken together, the integration of omics databases is a critical step towards enhancing crop plant breeding and improving global food security.

A cell wall-localized cytokinin/purine riboside nucleosidase is involved in apoplastic cytokinin metabolism in Oryza sativa

In the final step of cytokinin biosynthesis, the main pathway is the elimination of a ribose-phosphate moiety from the cytokinin nucleotide precursor by phosphoribohydrolase, an enzyme encoded by a gene named LONELY GUY (LOG). This reaction accounts for most of the cytokinin supply needed for regulating plant growth and development. In contrast, the LOG-independent pathway, in which dephosphorylation and deribosylation sequentially occur, is also thought to play a role in cytokinin biosynthesis, but the gene entity and physiological contribution have been elusive. In this study, we profiled the phytohormone content of chromosome segment substitution lines of Oryza sativa and searched for genes affecting the endogenous levels of cytokinin ribosides by quantitative trait loci analysis. Our approach identified a gene encoding an enzyme that catalyzes the deribosylation of cytokinin nucleoside precursors and other purine nucleosides. The cytokinin/purine riboside nucleosidase 1 (CPN1) we identified is a cell wall-localized protein. Loss-of-function mutations (cpn1) were created by inserting a Tos17-retrotransposon that altered the cytokinin composition in seedling shoots and leaf apoplastic fluid. The cpn1 mutation also abolished cytokinin riboside nucleosidase activity in leaf extracts and attenuated the trans-zeatin riboside-responsive expression of cytokinin marker genes. Grain yield of the mutants declined due to altered panicle morphology under field-grown conditions. These results suggest that the cell wall-localized LOG-independent cytokinin activating pathway catalyzed by CPN1 plays a role in cytokinin control of rice growth. Our finding broadens our spatial perspective of the cytokinin metabolic system.

Function identification and characterization of Oryza sativa ZRT and IRT-like proteins computationally for nutrition and biofortification in rice

Zinc plays a very critical role and function in all organisms. Its deficiency can cause a serious issue. In Oryza sativa, the ZRT/IRT transporter-like proteins play a role in the zinc metal uptake and transport. Few OsZIPs genes have been validated and characterized for their biological functions and most of OsZIPs are not well physiologically, biochemically and phenotypically characterized. In the current study, they analyzed for their function through subcellular localization, phylogenetic analysis, homology modeling, expression analysis, protein-protein interaction (PPI) network prediction, and prediction of their binding sites. Hierarchical clustering of OsZIP genes based on different anatomical parts and developmental stages also orthologs prediction was identified. The presence of SNPs, SSRs, ESTs, FSTs, MPSS, and SAGE tags were analyzed for useful development of markers. SNPs were identified in all OsZIPs genes and each gene was further classified based on their number and position in the 3'UTR and 5'UTR regions of the gene-specific sequences. Binding clusters and their location on the protein sequences were predicted. We found Changing in residues number and position which were due to partial overlapping and sequence alignment, but they share the same mechanism of binding and transporting Zinc. A wide range of CRISPR Cas9 gRNAs was designed based on single nucleotide polymorphism (SNP) for each OsZIP transporter gene for well-function identification and characterization with genome-wide association studies. Hence this study would provide useful information, understanding, and predicting molecular insights for the future studies that will help for improvement of nutritional quality of rice varieties.Communicated by Ramaswamy H. Sarma.

Genome-Wide Association Study Reveals the Genetic Basis of Total Flavonoid Content in Brown Rice

Flavonoids have anti-inflammatory, antioxidative, and anticarcinogenic effects. Breeding rice varieties rich in flavonoids can prevent chronic diseases such as cancer and cardio-cerebrovascular diseases. However, most of the genes reported are known to regulate flavonoid content in leaves or seedlings. To further elucidate the genetic basis of flavonoid content in rice grains and identify germplasm rich in flavonoids in grains, a set of rice core collections containing 633 accessions from 32 countries was used to determine total flavonoid content (TFC) in brown rice. We identified ten excellent germplasms with TFC exceeding 300 mg/100 g. Using a compressed mixed linear model, a total of 53 quantitative trait loci (QTLs) were detected through a genome-wide association study (GWAS). By combining linkage disequilibrium (LD) analysis, location of significant single nucleotide polymorphisms (SNPs), gene expression, and haplotype analysis, eight candidate genes were identified from two important QTLs (qTFC1-6 and qTFC9-7), among which LOC_Os01g59440 and LOC_Os09g24260 are the most likely candidate genes. We also analyzed the geographic distribution and breeding utilization of favorable haplotypes of the two genes. Our findings provide insights into the genetic basis of TFC in brown rice and could facilitate the breeding of flavonoid-rich varieties, which may be a prevention and adjuvant treatment for cancer and cardio-cerebrovascular diseases.

In-silico analysis of heat shock transcription factor (OsHSF) gene family in rice (Oryza sativa L.)

BACKGROUND

One of the most important cash crops worldwide is rice (Oryza sativa L.). Under varying climatic conditions, however, its yield is negatively affected. In order to create rice varieties that are resilient to abiotic stress, it is essential to explore the factors that control rice growth, development, and are source of resistance. HSFs (heat shock transcription factors) control a variety of plant biological processes and responses to environmental stress. The in-silico analysis offers a platform for thorough genome-wide identification of OsHSF genes in the rice genome.

RESULTS

In this study, 25 randomly dispersed HSF genes with significant DNA binding domains (DBD) were found in the rice genome. According to a gene structural analysis, all members of the OsHSF family share Gly-66, Phe-67, Lys-69, Trp-75, Glu-76, Phe-77, Ala-78, Phe-82, Ile-93, and Arg-96. Rice HSF family genes are widely distributed in the vegetative organs, first in the roots and then in the leaf and stem; in contrast, in reproductive tissues, the embryo and lemma exhibit the highest levels of gene expression. According to chromosomal localization, tandem duplication and repetition may have aided in the development of novel genes in the rice genome. OsHSFs have a significant role in the regulation of gene expression, regulation in primary metabolism and tolerance to environmental stress, according to gene networking analyses.

CONCLUSION

Six genes viz; Os01g39020, Os01g53220, Os03g25080, Os01g54550, Os02g13800 and Os10g28340 were annotated as promising genes. This study provides novel insights for functional studies on the OsHSFs in rice breeding programs. With the ultimate goal of enhancing crops, the data collected in this survey will be valuable for performing genomic research to pinpoint the specific function of the HSF gene during stress responses.

Fantastic genes: where and how to find them? Exploiting rice genetic resources for the improvement of yield, tolerance, and resistance to a wide array of stresses in rice

Rice production is a critical component of global food security. To date, rice is grown in over 100 countries and is the primary source of food for more than 3 billion people. Despite its importance, rice production is facing numerous challenges that threaten its future viability. One of the primary problems is the advent of climate change. The changing climatic conditions greatly affect the growth and productivity of rice crop and the quality of rice yield. Similarly, biotic stresses brought about by pathogen and pest infestations are greatly affecting the productivity of rice. To address these issues, the utilization of rice genetic resources is necessary to map, identify, and understand the genetics of important agronomic traits. This review paper highlights the role of rice genetic resources for developing high-yielding and stress-tolerant rice varieties. The integration of genetic, genomic, and phenomic tools in rice breeding programs has led to the development of high-yielding and stress-tolerant rice varieties. The collaboration of multidisciplinary teams of experts, sustainable farming practices, and extension services for farmers is essential for accelerating the development of high-yielding and stress-tolerant rice varieties.

Serotonin: A frontline player in plant growth and stress responses

Serotonin is a well-studied pineal hormone that functions as a neurotransmitter in mammals and is found in varying amounts in diverse plant species. By modulating gene and phytohormonal crosstalk, serotonin has a significant role in plant growth and stress response, including root, shoot, flowering, morphogenesis, and adaptability responses to numerous environmental signals. Despite its prevalence and importance in plant growth and development, its molecular action, regulation and signalling processes remain unknown. Here, we highlight the current knowledge of the role of serotonin-mediated regulation of plant growth and stress response. We focus on serotonin and its regulatory connections with phytohormonal crosstalk and address their possible functions in coordinating diverse phytohormonal responses during distinct developmental phases, correlating with melatonin. Additionally, we have also discussed the possible role of microRNAs (miRNAs) in the regulation of serotonin biosynthesis. In summary, serotonin may act as a node molecule to coordinate the balance between plant growth and stress response, which may shed light on finding its key regulatory pathways for uncovering its mysterious molecular network.

Natural variation in Tiller Number 1 affects its interaction with TIF1 to regulate tillering in rice

Tiller number per plant-a cardinal component of ideal plant architecture-affects grain yield potential. Thus, alleles positively affecting tillering must be mined to promote genetic improvement. Here, we report a Tiller Number 1 (TN1) protein harbouring a bromo-adjacent homology domain and RNA recognition motifs, identified through genome-wide association study of tiller numbers. Natural variation in TN1 affects its interaction with TIF1 (TN1 interaction factor 1) to affect DWARF14 expression and negatively regulate tiller number in rice. Further analysis of variations in TN1 among indica genotypes according to geographical distribution revealed that low-tillering varieties with TN1-hap^L are concentrated in Southeast Asia and East Asia, whereas high-tillering varieties with TN1-hap^H are concentrated in South Asia. Taken together, these results indicate that TN1 is a tillering regulatory factor whose alleles present apparent preferential utilization across geographical regions. Our findings advance the molecular understanding of tiller development.

Fine Mapping of the Affecting Tillering and Plant Height Gene CHA-1 in Rice

The plant architecture of rice is an important factor affecting yield. Strigolactones (SLs) are newly discovered carotenoid-derived plant hormones that play an important role in rice plant architecture. In this study, a high-tillering dwarf mutant, CHA-1, was identified by spatial mutagenesis. CHA-1 was located in the region of 31.52-31.55 MB on chromosome 1 by map-based cloning. Compared with the wild-type THZ, the CHA-1 mutant showed that ACCAC replaced TGGT in the coding region of the candidate gene LOC_Os01g54810, leading to premature termination of expression. Genetic complementation experiments proved that LOC_Os01g54810 was CHA-1, which encodes a putative member of Class III lipase. Expression analysis showed that CHA-1 was constitutively expressed in various organs of rice. Compared with those in THZ, the expression levels of the D17 and D10 genes were significantly downregulated in the CHA-1 mutant. In addition, the concentrations of ent-2'-epi-5-deoxystrigol (epi-5DS) in the root exudates of the CHA-1 mutant was significantly reduced compared with that of THZ, and exogenous application of GR24 inhibited the tillering of the CHA-1 mutant. These results suggest that CHA-1 influences rice architecture by affecting SL biosynthesis.

Time of Day Analysis over a Field Grown Developmental Time Course in Rice

Plants integrate time of day (TOD) information over an entire season to ensure optimal growth, flowering time, and grain fill. However, most TOD expression studies have focused on a limited number of combinations of daylength and temperature under laboratory conditions. Here, an Oryza sativa (rice) expression study that followed TOD expression in the field over an entire growing season was re-analyzed. Similar to Arabidopsis thaliana, almost all rice genes have a TOD-specific expression over the developmental time course. As has been suggested in other grasses, thermocycles were a stronger cue for TOD expression than the photocycles over the growing season. All the core circadian clock genes display consistent TOD expression over the season with the interesting exception that the two grass paralogs of EARLY FLOWERING 3 (ELF3) display a distinct phasing based on the interaction between thermo- and photo-cycles. The dataset also revealed how specific pathways are modulated to distinct TOD over the season consistent with the changing biology. The data presented here provide a resource for researchers to study how TOD expression changes under natural conditions over a developmental time course, which will guide approaches to engineer more resilient and prolific crops.

Natural variation of RGN1a regulates grain number per panicle in japonica rice

The grain number per panicle (GNP) is an important yield component. Identifying naturally favorable variations in GNP will benefit high-yield rice breeding. Here, we performed a genome-wide association study using a mini-core collection of 266 cultivated rice accessions with deep sequencing data and investigated the phenotype for three years. Three genes, i.e., TOTOU1 (TUT1), Grain height date 7 (Ghd7), and Days to heading 7/Grain height date 7.1/Pseudo-Response Regulator37 (DTH7/Ghd7.1/OsPRR37), which regulate GNP, were found in the quantitative trait loci (QTL) identified in this study. A stable QTL, qGNP1.3, which showed a strong correlation with variations in GNP, was repeatedly detected. After functional and transgenic phenotype analysis, we identified a novel gene, regulator of grain number 1a (RGN1a), which codes for protein kinase, controlling GNP in rice. The RGN1a mutation caused 37.2%, 27.8%, 51.2%, and 25.5% decreases in grain number, primary branch number per panicle, secondary branch number per panicle, and panicle length, respectively. Furthermore, breeding utilization analysis revealed that the additive effects of the dominant allelic variants of RGN1a and DTH7 played a significant role in increasing the grain number per panicle in japonica rice. Our findings enrich the gene pool and provide an effective strategy for the genetic improvement of grain numbers.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

Effect of Exogenous Glycine Betaine on the Germination of Tomato Seeds under Cold Stress

Cold stress is known to influence tomato growth, development, and yield. In this study, we analyzed the germination of tomato seeds treated with exogenous glycine betaine (GB) at a low temperature (14 °C). The results showed that cold stress inhibited tomato seed germination, and pretreatment with exogenous GB reduced this inhibition and enhanced the germination rate (GR), germination index (GI), and viability of tomato seeds at low temperatures. Analysis of gene expression and metabolism revealed that GB positively regulated endogenous hormone gibberellin (GA) content and negatively regulated abscisic acid (ABA) content, while GB reduced the starch content in the seeds by up-regulating the amylase gene expression. Gene expression analysis showed that the key genes (SlSOD, SlPOD, and SlchlAPX) involved in reactive oxygen species (ROS) scavenging systems were up-regulated in GB-pretreated tomato seeds compared with the control. At the same time, levels of malondialdehyde and hydrogen peroxide were significantly lower, while the proline content and peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) levels were elevated compared with those in the control. These results demonstrate that exogenous GB as a positive regulator effectively alleviated the inhibition of tomato seed germination under cold stress by different signal pathways.

CRISPR/Cas9 suppression of OsAT10, a rice BAHD acyltransferase, reduces p-coumaric acid incorporation into arabinoxylan without increasing saccharification

Ester-linked hydroxycinnamic acids ferulic acid (FA) and para-coumaric acid (p-CA) play important roles in crosslinking within cell wall arabinoxylans (AX) and between AX and lignin in grass cell walls. The addition of hydroxycinnamates to AX, is mediated by the Mitchell clade of BAHD acyl-coenzyme A-utilizing transferases. Overexpression of OsAT10 (a Mitchell clade BAHD acyl transferase) in rice, has previously been shown to increase p-CA content in AX in leaves and stems, leading to increased cell wall digestibility, potentially associated with a concomitant decrease in FA content. To investigate the physiological role of OsAT10 we established CRISPR/Cas9 rice knock-out mutants devoid of OsAT10. Our analysis of hydroxycinnamic acid content in wild type plants revealed that AX associated p-CA is found almost exclusively in rice husks, with very little found in other tissues. Mutant plants were essentially devoid of ester-linked p-CA associated with AX, indicating that OsAT10 represents the major enzyme responsible for the addition of p-CA to arabinoxylan in rice plants. We found no change in the digestibility of rice husk lacking AX-associated p-CA, suggesting that the changes in digestibility seen in OsAT10 overexpressing plants were solely due to compensatory decreases in AX-associated FA.

Golgi-Localized OsFPN1 is Involved in Co and Ni Transport and Their Detoxification in Rice

Cobalt (Co) and nickel (Ni) are beneficial and essential elements for plants, respectively, with the latter required for urease activity, which hydrolyzes urea into ammonium in plants. However, excess Co and Ni are toxic to plants and their transport mechanisms in rice are unclear. Here, we analyzed an ethyl methanesulfonate (EMS)-mutagenized rice mutant, 1187_n, with increased Co and Ni contents in its brown rice and shoots. 1187_n has a mutation in OsFPN1, which was correlated with a high Co and Ni phenotype in F₂ crosses between the parental line and mutant. In addition, CRISPR/Cas9 mutants exhibited a phenotype similar to that of 1187_n, demonstrating that OsFPN1 is the causal gene of the mutant. In addition to the high Co and Ni in brown rice and shoots, the mutant also exhibited high Co and Ni concentrations in the xylem sap, but low concentrations in the roots, suggesting that OsFPN1 is involved in the root-to-shoot translocation of Co and Ni. The growth of 1187_n and CRISPR/Cas9 lines were suppressed under high Co and Ni condition, indicating OsFPN1 is required for the normal growth under high Co and Ni. An OsFPN1-green fluorescent protein (GFP) fusion protein was localized to the Golgi apparatus. Yeast carrying GFP-OsFPN1 increased sensitivity to high Co contents and decreased Co and Ni accumulation. These results suggest that OsFPN1 can transport Co and Ni and is vital detoxification in rice.

Global Integrated Genomic and Transcriptomic Analyses of MYB Transcription Factor Superfamily in C3 Model Plant Oryza sativa (L.) Unravel Potential Candidates Involved in Abiotic Stress Signaling

Plant transcription factors (TFs) are significant players in transcriptional regulations, signal transduction, and constitute an integral part of signaling networks. MYB TFs are major TF superfamilies that play pivotal roles in regulation of transcriptional reprogramming, physiological processes, and abiotic stress (AbS) responses. To explore the understanding of MYB TFs, genome and transcriptome-wide identification was performed in the C3 model plant, Oryza sativa (OsMYB). This study retrieved 114 OsMYB TFs that were computationally analyzed for their expression profiling, gene organization, cis-acting elements, and physicochemical properties. Based on the microarray datasets, six OsMYB genes which were sorted out and identified by a differential expression pattern were noted in various tissues. Systematic expression profiling of OsMYB TFs showed their meta-differential expression of different AbS treatments, spatio-temporal gene expression of various tissues and their growth in the field, and gene expression profiling in responses to phytohormones. In addition, the circular ideogram of OsMYB genes in related C4 grass plants conferred the gene synteny. Protein–protein interactions of these genes revealed the molecular crosstalk of OsMYB TFs. Transcriptional analysis (qPCR) of six OsMYB players in response to drought and salinity stress suggested the involvement in individual and combined AbS responses. To decipher how these OsMYB play functional roles in AbS dynamics, further research is a prerequisite.

The basic leucine zipper transcription factor OsbZIP83 and the glutaredoxins OsGRX6 and OsGRX9 facilitate rice iron utilization under the control of OsHRZ ubiquitin ligases

Under low iron availability, plants induce the expression of various genes for iron uptake and translocation. The rice (Oryza sativa) ubiquitin ligases OsHRZ1 and OsHRZ2 cause overall repression of these iron-related genes at the transcript level, but their protein-level regulation is unclear. We conducted a proteome analysis to identify key regulators whose abundance was regulated by OsHRZs at the protein level. In response to iron deficiency or OsHRZ knockdown, many genes showed differential regulation between the transcript and protein levels, including the TGA-type basic leucine zipper transcription factor OsbZIP83. We also identified two glutaredoxins, OsGRX6 and OsGRX9, as OsHRZ-interacting proteins in yeast and plant cells. OsGRX6 also interacted with OsbZIP83. Our in vitro degradation assay suggested that OsbZIP83, OsGRX6 and OsGRX9 proteins are subjected to 26S proteasome- and OsHRZ-dependent degradation. Proteome analysis and our in vitro degradation assay also suggested that OsbZIP83 protein was preferentially degraded under iron-deficient conditions in rice roots. Transgenic rice lines overexpressing OsGRX9 and OsbZIP83 showed improved tolerance to iron deficiency. Expression of iron-related genes was affected in the OsGRX9 and OsGRX6 knockdown lines, suggesting disturbed iron utilization and signaling. OsbZIP83 overexpression lines showed enhanced expression of OsYSL2 and OsNAS3, which are involved in internal iron translocation, in addition to OsGRX9 and genes related to phytoalexin biosynthesis and the salicylic acid pathway. The results suggest that OsbZIP83, OsGRX6 and OsGRX9 facilitate iron utilization downstream of the OsHRZ pathway.

Identification of a unique allele in the quantitative trait locus for crown root number in japonica rice from Japan using genome-wide association studies

To explore the genetic resources that could be utilized to help improve root system architecture phenotypes in rice (Oryza sativa), we have conducted genome-wide association studies to investigate maximum root length and crown root number in 135 10-day-old Japanese rice accessions grown hydroponically. We identified a quantitative trait locus for crown root number at approximately 32.7 Mbp on chromosome 4 and designated it qNCR1 (quantitative trait locus for Number of Crown Root 1). A linkage disequilibrium map around qNCR1 suggested that three candidate genes are involved in crown root number: a cullin (LOC_Os04g55030), a gibberellin 20 oxidase 8 (LOC_Os04g55070), and a cyclic nucleotide-gated ion channel (LOC_Os04g55080). The combination of haplotypes for each gene was designated as a haploblock, and haploblocks 1, 2, and 3 were defined. Compared to haploblock 1, the accessions with haploblocks 2 and 3 had fewer crown roots; approximately 5% and 10% reductions in 10-day-old plants and 15% and 25% reductions in 42-day-old plants, respectively. A Japanese leading variety Koshihikari and its progenies harbored haploblock 3. Their crown root number could potentially be improved using haploblocks 1 and 2.

OsJAZ11 regulates spikelet and seed development in rice

Seed size is one of the major determinants of seed weight and eventually, crop yield. As the global population is increasing beyond the capacity of current food production, enhancing seed size is a key target for crop breeders. Despite the identification of several genes and QTLs, current understanding about the molecular regulation of seed size/weight remains fragmentary. In the present study, we report novel role of a jasmonic acid (JA) signaling repressor, OsJAZ11 controlling rice seed width and weight. Transgenic rice lines overexpressing OsJAZ11 exhibited up to a 14% increase in seed width and ~30% increase in seed weight compared to wild type (WT). Constitutive expression of OsJAZ11 dramatically influenced spikelet morphogenesis leading to extra glume-like structures, open hull, and abnormal numbers of floral organs. Furthermore, overexpression lines accumulated higher JA levels in spikelets and developing seeds. Expression studies uncovered altered expression of JA biosynthesis/signaling and MADS box genes in overexpression lines compared to WT. Yeast two-hybrid and pull-down assays revealed that OsJAZ11 interacts with OsMADS29 and OsMADS68. Remarkably, expression of OsGW7, a key negative regulator of grain size, was significantly reduced in overexpression lines. We propose that OsJAZ11 participates in the regulation of seed size and spikelet development by coordinating the expression of JA-related, OsGW7 and MADS genes.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

The submergence tolerance regulator SUB1A differentially coordinates molecular adaptation to submergence in mature and growing leaves of rice (Oryza sativa L.)

A typical adaptive response to submergence regulated by SUB1A, the ethylene-responsive transcription factor gene, is the restricted elongation of the uppermost leaves. However, the molecular and physiological functions of SUB1A have been characterized using entire shoot tissues, most of which are mature leaves that do not elongate under submergence. We aimed to identify leaf-type-specific and overlapping adaptations coordinated in SUB1A-dependent and -independent manners. To this end, we compared the transcriptomic and hormonal responses to submergence between mature and growing leaves using rice genotypes with and without SUB1A. Monosaccharide, branched-chain amino acid, and nucleoside metabolism, associated with ATP synthesis, were commonly activated in both leaf types regardless of genotype. In both leaf types, pathways involved in carbohydrate and nitrogen metabolism were suppressed by SUB1A, with more severe restriction in growing leaves that have a greater energy demand if SUB1A is absent. In growing leaves, accumulation of and responsiveness to growth-regulating hormones were properly modulated by SUB1A, which correlated with restricted elongation. In mature leaves, submergence-induced auxin accumulation was suppressed by SUB1A. This study demonstrates that different sets of hormonal pathways, both of which are modulated by SUB1A, contribute to distinct adaptive responses to submergence in mature and growing rice leaves.

NARROW AND DWARF LEAF 1, the Ortholog of Arabidopsis ENHANCER OF SHOOT REGENERATION1/DORNRÖSCHEN, Mediates Leaf Development and Maintenance of the Shoot Apical Meristem in Oryza sativa L

The molecular basis for leaf development, a major focus in developmental biology, remains unclear in the monocotyledonous grass, rice (Oryza sativa). Here, we performed a mutant screen in rice and identified an AP2-type transcription factor family protein, NARROW AND DWARF LEAF1 (NDL1). NDL1 is the ortholog of Arabidopsis thaliana (subsequently called Arabidopsis) ENHANCER OF SHOOT REGENERATION1 (ESR1)/DORNRÖSCHEN (DRN) and mediates leaf development and maintenance of the shoot apical meristem (SAM). Loss of function of NDL1 results in bladeless leaves and SAMs that are flat, rather than dome-shaped, and lack cell proliferation activity. This loss of function also causes reduced auxin signaling. Moreover, as is the case with Arabidopsis ESR1/DRN, NDL1 plays crucial roles in shoot regeneration. Importantly, we found that NDL1 is not expressed in the SAM but is expressed in leaf primordia. We propose that NDL1 cell autonomously regulates leaf development, but non-cell autonomously regulates SAM maintenance in rice.

Abscisic acid is required for exodermal suberization to form a barrier to radial oxygen loss in the adventitious roots of rice (Oryza sativa)

To acclimate to waterlogged conditions, wetland plants form a barrier to radial oxygen loss (ROL) that can enhance oxygen transport to the root apex. We hypothesized that one or more hormones are involved in the induction of the barrier and searched for such hormones in rice. We previously identified 98 genes that were tissue-specifically upregulated during ROL barrier formation in rice. The RiceXPro database showed that most of these genes were highly enhanced by exogenous abscisic acid (ABA). We then examined the effect of ABA on ROL barrier formation by using an ABA biosynthesis inhibitor (fluridone, FLU), by applying exogenous ABA and by examining a mutant with a defective ABA biosynthesis gene (osaba1). FLU suppressed barrier formation in a stagnant solution that mimics waterlogged soil. Under aerobic conditions, rice does not naturally form a barrier, but 24 h of ABA treatment induced barrier formation. osaba1 did not form a barrier under stagnant conditions, but the application of ABA rescued the barrier. In parallel with ROL barrier formation, suberin lamellae formed in the exodermis. These findings strongly suggest that ABA is an inducer of suberin lamellae formation in the exodermis, resulting in an ROL barrier formation in rice.

Localized expression of the Dwarf14-like2a gene in rice roots on infection of arbuscular mycorrhizal fungus and hydrolysis of rac-GR24 by the encoded protein

Strigolactones (SLs) are plant hormones that control diverse aspects of the shoot and root growth and are exuded into the soil as recruitment signals for arbuscular mycorrhizal (AM) fungi. SL signaling in plants is transduced via the α/β-hydrolase receptor Dwarf14 (D14). The D14 family consists of D14, Dwarf14-like (D14L), and Dwarf14-like 2 (D14L2) clades in rice. The D14L receptor is known to condition pre-symbiotic perception of AM fungi. In this study, it was found that the Dwarf14-like2a (D14L2a) gene expression was significantly induced by AM fungal colonization. The transcript of D14L2a appeared not only in mature arbuscule-containing cells but also in epidermal/cortical cells at an early colonization stage and near the elongating intercellular hyphae. D14L2a transcript was detected normally in mycorrhizal roots of str1-2 mutant that form stunted arbuscules, suggesting that the gene expression is independent of arbuscule development. Moreover, the recombinant D14L2a protein exhibited hydrolase activity of synthetic SL, rac-GR24. Based on these results, we discussed the role of D14L2 in the establishment of AM symbiosis.

Genotype-dependent and heat-induced grain chalkiness in rice correlates with the expression patterns of starch biosynthesis genes

Starch biosynthesis is a complex process underlying grain chalkiness in rice in a genotype-dependent manner. Coordinated expression of starch biosynthesis genes is important for producing translucent rice grains, while disruption in this process leads to opaque or chalky grains. To better understand the dynamics of starch biosynthesis genes in grain chalkiness, six rice genotypes showing variable chalk levels were subjected to gene expression analysis during reproductive stages. In the chalky genotypes, peak expression of the large subunit genes of ADP-glucose pyrophosphorylase (AGPase), encoding the first key step in starch biosynthesis, occurred in the stages before grain filling commenced, creating a gap with the upregulation of starch synthase genes, granule bound starch synthase I (GBSSI) and starch synthase IIA (SSIIA). Whereas, in low-chalk genotypes, AGPase large subunit genes expressed at later stages, generally following the expression patterns of GBSSI and SSIIA. However, heat treatment altered the expression in a genotype-dependent manner that was accompanied by transformed grain morphology and increased chalkiness. The suppression of AGPase subunit genes during early grain filling stages was observed in the chalky genotypes or upon heat treatment, which could result in a limited pool of ADP-Glucose for synthesizing amylose and amylopectin, the major components of the starch. This suboptimal starch biosynthesis process could subsequently lead to inefficient grain filling and air pockets that contribute to chalkiness. In summary, this study suggests a mechanism of grain chalkiness based on the expression patterns of the starch biosynthesis genes in rice.

Genome-Wide Association Study of the Genetic Basis of Effective Tiller Number in Rice

BACKGROUND

Effective tiller number (ETN) has a pivotal role in determination of rice (Oryza sativa L.) grain yield. ETN is a complex quantitative trait regulated by both genetic and environmental factors. Despite multiple tillering-related genes have been cloned previously, few of them have been utilized in practical breeding programs.

RESULTS

In this study, we conducted a genome-wide association study (GWAS) for ETN using a panel of 490 rice accessions derived from the 3 K rice genomes project. Thirty eight ETN-associated QTLs were identified, interestingly, four of which colocalized with the OsAAP1, DWL2, NAL1, and OsWRKY74 gene previously reported to be involved in rice tillering regulation. Haplotype (Hap) analysis revealed that Hap5 of OsAAP1, Hap3 and 6 of DWL2, Hap2 of NAL1, and Hap3 and 4 of OsWRKY74 are favorable alleles for ETN. Pyramiding favorable alleles of all these four genes had more enhancement in ETN than accessions harboring the favorable allele of only one gene. Moreover, we identified 25 novel candidate genes which might also affect ETN, and the positive association between expression levels of the OsPILS6b gene and ETN was validated by RT-qPCR. Furthermore, transcriptome analysis on data released on public database revealed that most ETN-associated genes showed a relatively high expression from 21 days after transplanting (DAT) to 49 DAT and decreased since then. This unique expression pattern of ETN-associated genes may contribute to the transition from vegetative to reproductive growth of tillers.

CONCLUSIONS

Our results revealed that GWAS is a feasible way to mine ETN-associated genes. The candidate genes and favorable alleles identified in this study have the potential application value in rice molecular breeding for high ETN and grain yield.

Function of the pseudo phosphotransfer proteins has diverged between rice and Arabidopsis

The phytohormone cytokinin plays a significant role in nearly all aspects of plant growth and development. Cytokinin signaling has primarily been studied in the dicot model Arabidopsis, with relatively little work done in monocots, which include rice (Oryza sativa) and other cereals of agronomic importance. The cytokinin signaling pathway is a phosphorelay comprised of the histidine kinase receptors, the authentic histidine phosphotransfer proteins (AHPs) and type-B response regulators (RRs). Two negative regulators of cytokinin signaling have been identified: the type-A RRs, which are cytokinin primary response genes, and the pseudo histidine phosphotransfer proteins (PHPs), which lack the His residue required for phosphorelay. Here, we describe the role of the rice PHP genes. Phylogenic analysis indicates that the PHPs are generally first found in the genomes of gymnosperms and that they arose independently in monocots and dicots. Consistent with this, the three rice PHPs fail to complement an Arabidopsis php mutant (aphp1/ahp6). Disruption of the three rice PHPs results in a molecular phenotype consistent with these elements acting as negative regulators of cytokinin signaling, including the induction of a number of type-A RR and cytokinin oxidase genes. The triple php mutant affects multiple aspects of rice growth and development, including shoot morphology, panicle architecture, and seed fill. In contrast to Arabidopsis, disruption of the rice PHPs does not affect root vascular patterning, suggesting that while many aspects of key signaling networks are conserved between monocots and dicots, the roles of at least some cytokinin signaling elements are distinct.

Iron deficiency-inducible peptide-coding genes OsIMA1 and OsIMA2 positively regulate a major pathway of iron uptake and translocation in rice

Under low iron (Fe) availability, plants transcriptionally induce various genes responsible for Fe uptake and translocation to obtain adequate amounts of Fe. Although transcription factors and ubiquitin ligases involved in these Fe deficiency responses have been identified, the mechanisms coordinating these pathways have not been clarified in rice. Recently identified Fe-deficiency-inducible IRON MAN (IMA)/FE UPTAKE-INDUCING PEPTIDE (FEP) positively regulates many Fe-deficiency-inducible genes for Fe uptake in Arabidopsis. Here, we report that the expression of two IMA/FEP genes in rice, OsIMA1 and OsIMA2, is strongly induced under Fe deficiency, positively regulated by the transcription factors IDEF1, OsbHLH058, and OsbHLH059, as well as OsIMA1 and OsIMA2 themselves, and negatively regulated by HRZ ubiquitin ligases. Overexpression of OsIMA1 or OsIMA2 in rice conferred tolerance to Fe deficiency and accumulation of Fe in leaves and seeds. These OsIMA-overexpressing rice exhibited enhanced expression of all of the known Fe-deficiency-inducible genes involved in Fe uptake and translocation, except for OsYSL2, a Fe-nicotianamine transporter gene, in roots but not in leaves. Knockdown of OsIMA1 or OsIMA2 caused minor effects, including repression of some Fe uptake- and translocation-related genes in OsIMA1 knockdown roots. These results indicate that OsIMA1 and OsIMA2 play key roles in enhancing the major pathway of the Fe deficiency response in rice.

Genome-wide analysis of spatiotemporal expression patterns during rice leaf development

BACKGROUND

Rice leaves consist of three distinct regions along a proximal-distal axis, namely the leaf blade, sheath, and blade-sheath boundary region. Each region has a unique morphology and function, but the genetic programs underlying the development of each region are poorly understood. To fully elucidate rice leaf development and discover genes with unique functions in rice and grasses, it is crucial to explore genome-wide transcriptional profiles during the development of the three regions.

RESULTS

In this study, we performed microarray analysis to profile the spatial and temporal patterns of gene expression in the rice leaf using dissected parts of leaves sampled in broad developmental stages. The dynamics in each region revealed that the transcriptomes changed dramatically throughout the progress of tissue differentiation, and those of the leaf blade and sheath differed greatly at the mature stage. Cluster analysis of expression patterns among leaf parts revealed groups of genes that may be involved in specific biological processes related to rice leaf development. Moreover, we found novel genes potentially involved in rice leaf development using a combination of transcriptome data and in situ hybridization, and analyzed their spatial expression patterns at high resolution. We successfully identified multiple genes that exhibit localized expression in tissues characteristic of rice or grass leaves.

CONCLUSIONS

Although the genetic mechanisms of leaf development have been elucidated in several eudicots, direct application of that information to rice and grasses is not appropriate due to the morphological and developmental differences between them. Our analysis provides not only insights into the development of rice leaves but also expression profiles that serve as a valuable resource for gene discovery. The genes and gene clusters identified in this study may facilitate future research on the unique developmental mechanisms of rice leaves.

Genomics Analyses Reveal Unique Classification, Population Structure and Novel Allele of Neo-Tetraploid Rice

BACKGROUND

Neo-tetraploid rice (NTR) is a useful new germplasm that developed from the descendants of the autotetraploid rice (ATR) hybrids. NTR showed improved fertility and yield potential, and produced high yield heterosis when crossed with indica ATR for commercial utilization. However, their classification, population structure and genomic feature remain elusive.

RESULTS

Here, high-depth genome resequencing data of 15 NTRs and 18 ATRs, together with 38 publicly available data of diploid rice accessions, were analyzed to conduct classification, population structure and haplotype analyses. Five subpopulations were detected and NTRs were clustered into one independent group that was adjacent to japonica subspecies, which maybe the reason for high heterosis when NTRs crossed with indica ATRs. Haplotype patterns of 717 key genes that associated with yield and other agronomic traits were revealed in these NTRs. Moreover, a novel specific SNP variation was detected in the first exon of HSP101, a known heat-inducible gene, which was conserved in all NTRs but absent in ATRs, 3KRG and RiceVarMap2 databases. The novel allele was named as HSP101-1, which was confirmed to be a heat response factor by qRT-PCR, and knockout of HSP101-1 significantly decreased the thermotolerance capacity of NTR. Interestingly, HSP101-1 was also specifically expressed in the anthers of NTR at pre-meiotic and meiosis stages under optimal environment without heat stress, and its loss-of-function mutant showed significant decrease in fertility of NTR.

CONCLUSION

The construction of first genomic variation repository and the revelation of population structure provide invaluable information for optimizing the designs of tetraploid rice breeding. The detection of specific genomic variations offered useful genomic markers and new directions to resolve high fertility mechanism of NTR.

Transcriptome-based approaches for clarification of nutritional responses and improvement of crop production

Genome-wide transcriptome profiling is a powerful tool for identifying key genes and pathways involved in plant development and physiological processes. This review summarizes studies that have used transcriptome profiling mainly in rice to focus on responses to macronutrients such as nitrogen, phosphorus and potassium, and spatio-temporal root profiling in relation to the regulation of root system architecture as well as nutrient uptake and transport. We also discuss strategies based on meta- and co-expression analyses with different attributed transcriptome data, which can be used for investigating the regulatory mechanisms and dynamics of nutritional responses and adaptation, and speculate on further advances in transcriptome profiling that could have potential application to crop breeding and cultivation.

Cytological Observations and Bulked-Segregant Analysis Coupled Global Genome Sequencing Reveal Two Genes Associated with Pollen Fertility in Tetraploid Rice

Neo-tetraploid rice with high fertility is a useful germplasm for polyploid rice breeding, which was developed from the crossing of different autotetraploid rice lines. However, little information is available on the molecular mechanism underlying the fertility of neo-tetraploid rice. Here, two contrasting populations of tetraploid rice, including one with high fertility (hereafter referred to as JG) and another with low fertility (hereafter referred to as JD), were generated by crossing Huaduo 3 (H3), a high fertility neo-tetraploid rice that was developed by crossing Jackson-4x with 96025-4x, and Huajingxian74-4x (T452), a low fertility autotetraploid rice parent. Cytological, global genome sequencing-based bulked-segregant (BSA-seq) and CRISPR/Cas9 technology were employed to study the genes associated with pollen fertility in neo-tetraploid rice. The embryo sacs of JG and JD lines were normal; however, pollen fertility was low in JD, which led to scarce fertilization and low seed setting. Cytological observations displayed low pollen fertility (25.1%) and approximately 31.3 and 27.2% chromosome lagging at metaphase I and II, and 28.8 and 24.8% chromosome straggling at anaphase I and II in JD, respectively. BSA-seq of F_2–3 generations and RNA-seq of F₄ generation detected a common fragment, i.e., 18,915,234–19,500,000, at chromosome 7, which was comprised of 78 genes associated with fertility. Among 78 genes, 9 genes had been known to be involved in meiosis and pollen development. Two mutants ny1 (LOC_Os07g32406) and ny2 (LOC_Os07g32040) were generated by CRISPR/Cas9 knockout in neo-tetraploid rice, and which exhibited low pollen fertility and abnormal chromosome behavior. Our study revealed that two unknown genes, LOC_Os07g32406 (NY1) and LOC_Os07g32040 (NY2) play an important role in pollen development of neo-tetraploid rice and provides a new perspective about the genetic mechanisms of fertility in polyploid rice.

Targeted Transgene Expression in Rice Using a Callus Strong Promoter for Selectable Marker Gene Control

Precise expression of a transgene in the desired manner is important for plant genetic engineering and gene function deciphering, but it is a challenge to obtain specific transgene expression free from the interference of the constitutive promoters used to express the selectable marker gene, such as the Cauliflower mosaic virus (CaMV) 35S promoter. So, the solutions to avoid these inappropriate regulations are largely demanded. In this study, we report the characterization of a callus strong promoter (CSP1) in rice and its application for accurate transgene expression. Our results indicate that the high expression of the CSP1 promoter in the callus enables efficient selection of hygromycin equivalent to that provided by the CaMV 35S promoter, whereas its expression in other tissues is low. To evaluate possible leaky effects, the expression of a β-glucuronidase reporter driven by six specific promoters involving hormone signaling, pathogen response, cell fate determination, and proliferation was observed in transgenic rice plants generated by CSP1-mediated selection. Distinct β-glucuronidase expression was found consistently in most of the transgenic lines obtained for each promoter. In addition, we applied these specific marker lines to investigate the root cellular responses to exogenous cytokinin and auxin treatment. The results reveal that the root growth inhibition by cytokinin was differently regulated at high and low concentrations. In summary, we have established the feasibility of using callus-specific promoter-dependent selection to mitigate the transgene misexpression in rice. By enabling efficient transformation, rice plants with reliable transgene expression will be easily acquired for broad applications.

A fast, efficient and high-throughput procedure involving laser microdissection and RT droplet digital PCR for tissue-specific expression profiling of rice roots

Background

In rice, the cortex and outer tissues play a key role in submergence tolerance. The cortex differentiates into aerenchyma, which are air-containing cavities that allow the flow of oxygen from shoots to roots, whereas exodermis suberification and sclerenchyma lignification limit oxygen loss from the mature parts of roots by forming a barrier to root oxygen loss (ROL). The genes and their networks involved in the cellular identity and differentiation of these tissues remain poorly understood. Identification and characterization of key regulators of aerenchyma and ROL barrier formation require determination of the specific expression profiles of these tissues.

Results

We optimized an approach combining laser microdissection (LM) and droplet digital RT-PCR (ddRT-PCR) for high-throughput identification of tissue-specific expression profiles. The developed protocol enables rapid (within 3 days) extraction of high-quality RNA from root tissues with a low contamination rate. We also demonstrated the possibility of extracting RNAs from paraffin blocks stored at 4 °C without any loss of quality. We included a detailed troubleshooting guide that should allow future users to adapt the proposed protocol to other tissues and/or species. We demonstrated that our protocol, which combines LM with ddRT-PCR, can be used as a complementary tool to in situ hybridization for tissue-specific characterization of gene expression even with a low RNA concentration input. We illustrated the efficiency of the proposed approach by validating three of four potential tissue-specific candidate genes detailed in the RiceXpro database.

Conclusion

The detailed protocol and the critical steps required to optimize its use for other species will democratize tissue-specific transcriptome approaches combining LM with ddRT-PCR for analyses of plants.

Defects in the rice aconitase-encoding OsACO1 gene alter iron homeostasis

Rice aconitase gene OsACO1 is involved in the iron deficiency-signaling pathway for the expression of iron deficiency-inducible genes, either thorough enzyme activity or possible specific RNA binding for post-transcriptional regulation. Iron (Fe) is an essential element for virtually all living organisms. When plants are deficient in Fe, Fe acquisition systems are activated to maintain Fe homeostasis, and this regulation is mainly executed at the gene transcription level. Many molecules responsible for Fe uptake, translocation, and storage in plants have been identified and characterized. However, how plants sense Fe status within cells and then induce a transcriptional response is still unclear. In the present study, we found that knockdown of the OsACO1 gene, which encodes an aconitase in rice, leads to the down-regulation of selected Fe deficiency-inducible genes involved in Fe uptake and translocation in roots, and a decrease in Fe concentration in leaves, even when grown under Fe-sufficient conditions. OsACO1 knockdown plants showed a delayed transcriptional response to Fe deficiency compared to wild-type plants. In contrast, overexpression of OsACO1 resulted in the opposite effects. These results suggest that OsACO1 is situated upstream of the Fe deficiency-signaling pathway. Furthermore, we found that the OsACO1 protein potentially has RNA-binding activity. In vitro screening of RNA interactions with OsACO1 revealed that RNA potentially forms a unique stem-loop structure that interacts with OsACO1 via a conserved GGUGG motif within the loop structure. These results suggest that OsACO1 regulate Fe deficiency response either thorough enzyme activity catalyzing isomerization of citrate, or specific RNA binding for post-transcriptional regulation.

Genetic architecture to cause dynamic change in tiller and panicle numbers revealed by genome-wide association study and transcriptome profile in rice

Panicle number (PN) is one of the three yield components in rice. As one of the most unstable traits, the dynamic change in tiller number (DCTN) may determine the final PN. However, the genetic basis of DCTN and its relationship with PN remain unclear. Here, 377 deeply re-sequenced rice accessions were used to perform genome-wide association studies (GWAS) for tiller/PN. It was found that the DCTN pattern rather than maximum tiller number or effective tiller ratio is the determinant factor of high PN. The DCTN pattern that affords more panicles exhibits a period of stable tillering peak between 30 and 45 days after transplant (called DT30 and DT45, respectively), which was believed as an ideal pattern contributing to the steady transition from tiller development to panicle development (ST-TtP). Consistently, quantitative trait loci (QTL) expressed near DT30-DT45 were especially critical to the rice DCTN and in supporting the ST-TtP. The spatio-temporal expression analysis showed that the expression pattern of keeping relatively high expression in root at 24:00 (R24-P2) from about DT30 to DT45 is a typical expression pattern of cloned tiller genes, and the candidate genes with R24-P2 can facilitate the prediction of PN. Moreover, gene OsSAUR27 was identified by an integrated approach combining GWAS, bi-parental QTL mapping and transcription. These findings related to the genetic basis underlying the DCTN will provide the genetic theory in making appropriate decisions on field management, and in developing new varieties with high PN and ideal dynamic plant architecture.

Transcriptome profiling of laser-captured crown root primordia reveals new pathways activated during early stages of crown root formation in rice

Crown roots constitute the main part of the rice root system. Several key genes involved in crown root initiation and development have been identified by functional genomics approaches. Nevertheless, these approaches are impaired by functional redundancy and mutant lethality. To overcome these limitations, organ targeted transcriptome analysis can help to identify genes involved in crown root formation and early development. In this study, we generated an atlas of genes expressed in developing crown root primordia in comparison with adjacent stem cortical tissue at three different developmental stages before emergence, using laser capture microdissection. We identified 3975 genes differentially expressed in crown root primordia. About 30% of them were expressed at the three developmental stages, whereas 10.5%, 19.5% and 12.8% were specifically expressed at the early, intermediate and late stages, respectively. Sorting them by functional ontology highlighted an active transcriptional switch during the process of crown root primordia formation. Cross-analysis with other rice root development-related datasets revealed genes encoding transcription factors, chromatin remodeling factors, peptide growth factors, and cell wall remodeling enzymes that are likely to play a key role during crown root primordia formation. This atlas constitutes an open primary data resource for further studies on the regulation of crown root initiation and development.

Transcriptome and Gene Editing Analyses Reveal MOF1a Defect Alters the Expression of Genes Associated with Tapetum Development and Chromosome Behavior at Meiosis Stage Resulting in Low Pollen Fertility of Tetraploid Rice

Autotetraploid rice is a useful rice germplasm for polyploid rice breeding. However, low fertility limits its commercial production. A neo-tetraploid rice with high fertility was developed from the progenies of crossing between autotetraploid lines by our research group. Our previous study showed that a myeloblastosis (MYB) transcription factor, MOF1, might be associated with the pollen development in tetraploid rice. However, little information is available about its role in pollen development in tetraploid rice. Here, we identified a new haplotype of MOF1 from neo-tetraploid rice and marked it as MOF1a. Transcriptome and qRT-PCR analysis demonstrated that MOF1a highly expressed in anthers, and displayed differential expression in neo-tetraploid rice compared to tetraploid rice line with low pollen fertility. The mutant (mof1a) of MOF1a, which was generated by CRISPR/Cas9 knockout in neo-tetraploid rice, showed low pollen fertility, and also exhibited abnormal tapetum and middle layer development, and defective chromosome behaviors during meiosis. A total of 13 tapetal related genes were found to be up-regulated in meiotic anthers of MOF1a compared with wild type plants by RNA-seq analysis, including CYP703A3, PTC1, and OsABCG26, which had been demonstrated to affect tapetal development. Moreover, 335 meiosis-related genes displayed differential expression patterns at same stage, including nine important meiosis-related genes, such as metallothionein OsMT1a. These results demonstrated that MOF1a plays an important role in pollen development and provides a foundation for understanding the molecular mechanism underlying MOF1a in reproduction of tetraploid rice.

Reproductive tissues-specific meta-QTLs and candidate genes for development of heat-tolerant rice cultivars

By integrating genetics and genomics data, reproductive tissues-specific and heat stress responsive 35 meta-QTLs and 45 candidate genes were identified, which could be exploited through marker-assisted breeding for fast-track development of heat-tolerant rice cultivars. Rice holds the key to future food security. In rice-growing areas, temperature has already reached an optimum level for growth, hence, any further increase due to global climate change could significantly reduce rice yield. Several mapping studies have identified a plethora of reproductive tissue-specific and heat stress associated inconsistent quantitative trait loci (QTL), which could be exploited for improvement of heat tolerance. In this study, we performed a meta-analysis on previously reported QTLs and identified 35 most consistent meta-QTLs (MQTLs) across diverse genetic backgrounds and environments. Genetic and physical intervals of nearly 66% MQTLs were narrower than 5 cM and 2 Mb respectively, indicating hotspot genomic regions for heat tolerance. Comparative analyses of MQTLs underlying genes with microarray and RNA-seq based transcriptomic data sets revealed a core set of 45 heat-responsive genes, among which 24 were reproductive tissue-specific and have not been studied in detail before. Remarkably, all these genes corresponded to various stress associated functions, ranging from abiotic stress sensing to regulating plant stress responses, and included heat-shock genes (OsBiP2, OsMed37_1), transcription factors (OsNAS3, OsTEF1, OsWRKY10, OsWRKY21), transmembrane transporters (OsAAP7A, OsAMT2;1), sugar metabolizing (OsSUS4, α-Gal III) and abiotic stress (OsRCI2-7, SRWD1) genes. Functional data evidences from Arabidopsis heat-shock genes also suggest that OsBIP2 may be associated with thermotolerance of pollen tubes under heat stress conditions. Furthermore, promoters of identified genes were enriched with heat, dehydration, pollen and sugar responsive cis-acting regulatory elements, proposing a common regulatory mechanism might exist in rice for mitigating reproductive stage heat stress. These findings strongly support our results and provide new candidate genes for fast-track development of heat-tolerant rice cultivars.

Gene regulatory network and its constituent transcription factors that control nitrogen-deficiency responses in rice

Increase in the nitrogen (N)-use efficiency and optimization of N response in crop species are urgently needed. Although transcription factor-based genetic engineering is a promising approach for achieving these goals, transcription factors that play key roles in the response to N deficiency have not been studied extensively. Here, we performed RNA-seq analysis of root samples of 20 Asian rice (Oryza sativa) accessions with differential nutrient uptake. Data obtained from plants exposed to N-replete and N-deficient conditions were subjected to coexpression analysis and machine learning-based pathway inference to dissect the gene regulatory network required for the response to N deficiency. Four transcription factors, including members of the G2-like and bZIP families, were predicted to function as key regulators of gene transcription within the network in response to N deficiency. Cotransfection assays validated inferred novel regulatory pathways, and further analyses using genome-edited knockout lines suggested that these transcription factors are important for N-deficiency responses in planta. Many of the N deficiency-responsive genes, including those encoding key regulators within the network, were coordinately regulated by transcription factors belonging to different families. Transcription factors identified in this study could be valuable for the modification of N response and metabolism.

Production Assessment and Genome Comparison Revealed High Yield Potential and Novel Specific Alleles Associated with Fertility and Yield in Neo-Tetraploid Rice

BACKGROUND

Neo-tetraploid rice (NTR) is a new tetraploid rice germplasm that developed from the crossing and directional selection of different autotetraploid rice lines, which showed high fertility and promising yield potential. However, systematic yield assessment, genome composition and functional variations associated with fertility and yield remain elusive.

RESULTS

Two season's field trials of 15 NTRs and 27 autotetraploid rice (ATR) lines revealed that the improvement of YPP (yield per plant, 4.45 g increase) were significantly associated with the increase of SS (seed setting, 29.44% increase), and yield and seed setting of NTRs improved significantly compared to parental lines. Whole genome resequencing of 13 NTR sister lines and their parents at about 48.63 depth were conducted and genome compositions were illustrated using inherited chromosomal blocks. Interestingly, 222 non-parental genes were detected between NTRs and their low fertility parental lines, which were conserved in 13 NTRs. These genes were overlapped with yield and fertility QTLs, and RNA-Seq analysis revealed that 81 of them were enriched in reproductive tissues. CRISPR/Cas9 gene knockout was conducted for 9 non-parental genes to validate their function. Knockout mutants showed on an average 25.63% and 4.88 g decrease in SS and YPP, respectively. Notably, some mutants showed interesting phenotypes, e.g., kin7l (kinesin motor gene) and kin14m (kinesin motor gene), bzr3 (BES1/BZR1 homolog) and nrfg4 (neo-tetraploid rice fertility related gene) exhibited 44.65%, 24.30%, 24.42% and 28.33% decrease in SS and 8.81 g, 4.71 g, 5.90 g, 6.22 g reduction in YPP, respectively.

CONCLUSION

Comparative genomics provides insights into genome composition of neo-tetraploid rice and the genes associated with fertility and yield will play important role to reveal molecular mechanisms for the improvement of tetraploid rice.

Novel allelic variation in the Phospholipase D alpha1 gene (OsPLDα1) of wild Oryza species implies to its low expression in rice bran

Rice bran, a by-product after milling, is a rich source of phytonutrients like oryzanols, tocopherols, tocotrienols, phytosterols, and dietary fibers. Moreover, exceptional properties of the rice bran oil make it unparalleled to other vegetable oils. However, a lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) causes rancidity and ‘stale flavor’ in the oil, and thus limits the rice bran usage for human consumption. To improve the rice bran quality, sequence based allele mining at OsPLDα1 locus (3.6 Kb) was performed across 48 accessions representing 11 wild Oryza species, 8 accessions of African cultivated rice, and 7 Oryza sativa cultivars. From comparative sequence analysis, 216 SNPs and 30 InDels were detected at the OsPLDα1 locus. Phylogenetic analysis revealed 20 OsPLDα1 cDNA variants which further translated into 12 protein variants. The O. officinalis protein variant, when compared to Nipponbare, showed maximum variability comprising 22 amino acid substitutions and absence of two peptides and two β-sheets. Further, expression profiling indicated significant differences in transcript abundance within as well as between the OsPLDα1 variants. Also, a new OsPLDα1 transcript variant having third exon missing in it, Os01t0172400-06, has been revealed. An O. officinalis accession (IRGC101152) had lowest gene expression which suggests the presence of novel allele, named as OsPLDα1-1a (GenBank accession no. MF966931). The identified novel allele could be further deployed in the breeding programs to overcome rice bran rancidity in elite cultivars.

Transcriptomic Insights into Mechanisms of Early Seed Maturation in the Garden Pea (Pisum sativum L.)

The garden pea (Pisum sativum L.) is a legume crop of immense economic value. Extensive breeding has led to the emergence of numerous pea varieties, of which some are distinguished by accelerated development in various stages of ontogenesis. One such trait is rapid seed maturation, which, despite novel insights into the genetic control of seed development in legumes, remains poorly studied. This article presents an attempt to dissect mechanisms of early maturation in the pea line Sprint-2 by means of whole transcriptome RNA sequencing in two developmental stages. By using a de novo assembly approach, we have obtained a reference transcriptome of 25,756 non-redundant entries expressed in pea seeds at either 10 or 20 days after pollination. Differential expression in Sprint-2 seeds has affected 13,056 transcripts. A comparison of the two pea lines with a common maturation rate demonstrates that while at 10 days after pollination, Sprint-2 seeds show development retardation linked to intensive photosynthesis, morphogenesis, and cell division, and those at 20 days show a rapid onset of desiccation marked by the cessation of translation and cell anabolism and accumulation of dehydration-protective and -storage moieties. Further inspection of certain transcript functional categories, including the chromatin constituent, transcription regulation, protein turnover, and hormonal regulation, has revealed transcriptomic trends unique to specific stages and cultivars. Among other remarkable features, Sprint-2 demonstrated an enhanced expression of transposable element-associated open reading frames and an altered expression of major maturation regulators and DNA methyltransferase genes. To the best of our knowledge, this is the first comparative transcriptomic study in which the issue of the seed maturation rate is addressed.

Evaluation of Whole-Genome Sequence, Genetic Diversity, and Agronomic Traits of Basmati Rice (Oryza sativa L.)

Basmati is considered a unique varietal group of rice (Oryza sativa L.) because of its aroma and superior grain quality. Previous genetic analyses of rice showed that most of the Basmati varieties are classified into the aromatic group. Despite various efforts, genomic relationship of Basmati rice with other varietal groups and genomic variation in Basmati rice are yet to be understood. In the present study, we resequenced the whole genome of three traditional Basmati varieties at a coverage of more than 25X using Illumina HiSeq2500 and mapped the obtained sequences to the reference genome sequences of Nipponbare (japonica rice), Kasalath (aus rice), and Zhenshan 97 (indica rice). Comparison of these sequences revealed common single nucleotide polymorphisms (SNPs) in the genic regions of three Basmati varieties. Analysis of these SNPs revealed that Basmati varieties showed fewer sequence variations compared with the aus group than with the japonica and indica groups. Gene ontology (GO) enrichment analysis indicated that SNPs were present in genes with various biological, molecular, and cellular functions. Additionally, functional annotation of the Basmati mutated gene cluster shared by Nipponbare, Kasalath, and Zhenshan 97 was found to be associated with the metabolic process involved in the cellular aromatic compound, suggesting that aroma is an important specific genomic feature of Basmati varieties. Furthermore, 30 traditional Basmati varieties were classified into three different groups, aromatic (22 varieties), aus (four varieties), and indica (four varieties), based on genome-wide SNPs. All 22 aromatic Basmati varieties harbored the fragrant-inducing Badh2 allele. We also performed comparative analysis of 13 key agronomic and grain quality traits of Basmati rice and other rice varieties. Three traits including length-to-width ratio of grain (L/W ratio), panicle length (PL), and amylose content (AC) showed significant (P < 0.05 and P < 0.01) differences between the aromatic and indica/aus groups. Comparative analysis of genome structure, based on genome sequence variation and GO analysis, revealed that the Basmati genome was derived mostly from the aus and japonica groups. Overall, whole-genome sequence data and genetic diversity information obtained in this study will serve as an important resource for molecular breeding and genetic analysis of Basmati varieties.

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.