Pollen development from the microspore involves a series of coordinated cellular events, and the resulting mature pollen has a specialized function to quickly germinate, produce a polar-growth pollen tube derived from the vegetative cell, and deliver two sperm cells into the embryo sac for double fertilization. The gene expression profiles of developing and germinated pollen have been characterised by use of the eudicot model plant Arabidopsis. Rice, one of the most important cereal crops, has been used as an excellent monocot model. A comprehensive analysis of transcriptome profiles of developing and germinated pollen in rice is important to understand the conserved and diverse mechanism underlying pollen development and germination in eudicots and monocots.
We used Affymetrix GeneChip® Rice Genome Array to comprehensively analyzed the dynamic changes in the transcriptomes of rice pollen at five sequential developmental stages from microspores to germinated pollen. Among the 51,279 transcripts on the array, we found 25,062 pollen-preferential transcripts, among which 2,203 were development stage-enriched. The diversity of transcripts decreased greatly from microspores to mature and germinated pollen, whereas the number of stage-enriched transcripts displayed a "U-type" change, with the lowest at the bicellular pollen stage; and a transition of overrepresented stage-enriched transcript groups associated with different functional categories, which indicates a shift in gene expression program at the bicellular pollen stage. About 54% of the now-annotated rice F-box protein genes were expressed preferentially in pollen. The transcriptome profile of germinated pollen was significantly and positively correlated with that of mature pollen. Analysis of expression profiles and coexpressed features of the pollen-preferential transcripts related to cell cycle, transcription, the ubiquitin/26S proteasome system, phytohormone signalling, the kinase system and defense/stress response revealed five expression patterns, which are compatible with changes in major cellular events during pollen development and germination. A comparison of pollen transcriptomes between rice and Arabidopsis revealed that 56.6% of the rice pollen preferential genes had homologs in Arabidopsis genome, but 63.4% of these homologs were expressed, with a small proportion being expressed preferentially, in Arabidopsis pollen. Rice and Arabidopsis pollen had non-conservative transcription factors each.
Our results demonstrated that rice pollen expressed a set of reduced but specific transcripts in comparison with vegetative tissues, and the number of stage-enriched transcripts displayed a "U-type" change during pollen development, with the lowest at the bicellular pollen stage. These features are conserved in rice and Arabidopsis. The shift in gene expression program at the bicellular pollen stage may be important to the transition from earlier cell division to later pollen maturity. Pollen at maturity pre-synthesized transcripts needed for germination and early pollen tube growth. The transcription regulation associated with pollen development would have divergence between the two species. Our results also provide novel insights into the molecular program and key components of the regulatory network regulating pollen development and germination.
Duplications are very common in the evolution of plant genomes, explaining the high number of members in plant gene families. New genes born after duplication can undergo pseudogenization, neofunctionalization or subfunctionalization. Rice is a model for functional genomics research, an important crop for human nutrition and a target for biofortification. Increased zinc and iron content in the rice grain could be achieved by manipulation of metal transporters. Here, we describe the ZINC-INDUCED FACILITATOR-LIKE (ZIFL) gene family in plants, and characterize the genomic structure and expression of rice paralogs, which are highly affected by segmental duplication.
Sequences of sixty-eight ZIFL genes, from nine plant species, were comparatively analyzed. Although related to MSF_1 proteins, ZIFL protein sequences consistently grouped separately. Specific ZIFL sequence signatures were identified. Monocots harbor a larger number of ZIFL genes in their genomes than dicots, probably a result of a lineage-specific expansion. The rice ZIFL paralogs were named OsZIFL1 to OsZIFL13 and characterized. The genomic organization of the rice ZIFL genes seems to be highly influenced by segmental and tandem duplications and concerted evolution, as rice genome contains five highly similar ZIFL gene pairs. Most rice ZIFL promoters are enriched for the core sequence of the Fe-deficiency-related box IDE1. Gene expression analyses of different plant organs, growth stages and treatments, both from our qPCR data and from microarray databases, revealed that the duplicated ZIFL gene pairs are mostly co-expressed. Transcripts of OsZIFL4, OsZIFL5, OsZIFL7, and OsZIFL12 accumulate in response to Zn-excess and Fe-deficiency in roots, two stresses with partially overlapping responses.
We suggest that ZIFL genes have different evolutionary histories in monocot and dicot lineages. In rice, concerted evolution affected ZIFL duplicated genes, possibly maintaining similar expression patterns between pairs. The enrichment for IDE1 boxes in rice ZIFL gene promoters suggests a role in Zn-excess and Fe-deficiency up-regulation of ZIFL transcripts. Moreover, this is the first description of the ZIFL gene family in plants and the basis for functional studies on this family, which may play important roles in Zn and Fe homeostasis in plants.
Plants face a variety of environmental stresses and have evolved molecular mechanisms to survive these challenges. One of these stresses is low oxygen conditions, which can occur under flooding conditions. Rice (Oryza sativa) is somewhat unique for its ability to tolerate and even germinate under low to no oxygen conditions. In this study, we examined global transcriptomic responses over the course of germination and in response to low oxygen and other abiotic stress in rice and Arabidopsis (Arabidopsis thaliana). Over 150 microarray datasets were analyzed in parallel to determine just how unique the low oxygen response is in rice. Comparison of aerobic germination in rice and Arabidopsis, with anaerobic germination in rice revealed conserved transcriptomic responses that are not only conserved across both species but also occur in the absence of oxygen in rice. Thus, these genes may represent functions necessary for the developmental progression of germination, whether or not oxygen is present in rice. Analysis of genes that responded differently in rice compared to Arabidopsis revealed responses specific to anaerobic germination in rice, including the down-regulation of genes encoding redox functions and up-regulation of receptor kinases. Comparison of a range of hypoxia/anoxia studies within and across Arabidopsis and rice revealed both conserved and species specific changes in gene expression (e.g., Arabidopsis specific up-regulation of WRKYs and rice specific down-regulation of heme), unveiling unique transcriptomic signatures of the low oxygen response. Lastly, a comparison of the low oxygen response with cold, salt, drought and heat stress revealed some similarity with the response to heat stress in Arabidopsis, which was not seen in rice. Comparison of these heat-responsive, abiotic stress marker genes in Arabidopsis with their rice orthologs revealed that while low oxygen may be perceived as an abiotic stress in Arabidopsis, this is not the case in rice.
low oxygen; transcriptomes; rice; Arabidopsis; microarrays
Plant nutrition is one of the important areas for improving the yield and quality in crops as well as non-crop plants. Potassium is an essential plant nutrient and is required in abundance for their proper growth and development. Potassium deficiency directly affects the plant growth and hence crop yield and production. Recently, potassium-dependent transcriptomic analysis has been performed in the model plant Arabidopsis, however in cereals and crop plants; such a transcriptome analysis has not been undertaken till date. In rice, the molecular mechanism for the regulation of potassium starvation responses has not been investigated in detail. Here, we present a combined physiological and whole genome transcriptomic study of rice seedlings exposed to a brief period of potassium deficiency then replenished with potassium. Our results reveal that the expressions of a diverse set of genes annotated with many distinct functions were altered under potassium deprivation. Our findings highlight altered expression patterns of potassium-responsive genes majorly involved in metabolic processes, stress responses, signaling pathways, transcriptional regulation, and transport of multiple molecules including K+. Interestingly, several genes responsive to low-potassium conditions show a reversal in expression upon resupply of potassium. The results of this study indicate that potassium deprivation leads to activation of multiple genes and gene networks, which may be acting in concert to sense the external potassium and mediate uptake, distribution and ultimately adaptation to low potassium conditions. The interplay of both upregulated and downregulated genes globally in response to potassium deprivation determines how plants cope with the stress of nutrient deficiency at different physiological as well as developmental stages of plants.
Arabidopsis thaliana is clearly established as the model plant species. Given the ever-growing demand for food, there is a need to translate the knowledge learned in Arabidopsis to agronomically important species, such as rice (Oryza sativa). To gain a comparative insight into the similarities and differences into how organs are built and how plants respond to stress, the transcriptomes of Arabidopsis and rice were compared at the level of gene orthology and functional categorisation.
Organ specific transcripts in rice and Arabidopsis display less overlap in terms of gene orthology compared to the orthology observed between both genomes. Although greater overlap in terms of functional classification was observed between root specific transcripts in rice and Arabidopsis, this did not extend to flower, leaf or seed specific transcripts. In contrast, the overall abiotic stress response transcriptome displayed a significantly greater overlap in terms of gene orthology compared to the orthology observed between both genomes. However, ~50% or less of these orthologues responded in a similar manner in both species. In fact, under cold and heat treatments as many or more orthologous genes responded in an opposite manner or were unchanged in one species compared to the other. Examples of transcripts that responded oppositely include several genes encoding proteins involved in stress and redox responses and non-symbiotic hemoglobins that play central roles in stress signalling pathways. The differences observed in the abiotic transcriptomes were mirrored in the presence of cis-acting regulatory elements in the promoter regions of stress responsive genes and the transcription factors that potentially bind these regulatory elements. Thus, both the abiotic transcriptome and its regulation differ between rice and Arabidopsis.
These results reveal significant divergence between Arabidopsis and rice, in terms of the abiotic stress response and its regulation. Both plants are shown to employ unique combinations of genes to achieve growth and stress responses. Comparison of these networks provides a more rational approach to translational studies that is based on the response observed in these two diverse plant models.
Functional redundancy limits detailed analysis of genes in many organisms. Here, we report a method to efficiently overcome this obstacle by combining gene expression data with analysis of gene-indexed mutants. Using a rice NSF45K oligo-microarray to compare 2-week-old light- and dark-grown rice leaf tissue, we identified 365 genes that showed significant 8-fold or greater induction in the light relative to dark conditions. We then screened collections of rice T-DNA insertional mutants to identify rice lines with mutations in the strongly light-induced genes. From this analysis, we identified 74 different lines comprising two independent mutant lines for each of 37 light-induced genes. This list was further refined by mining gene expression data to exclude genes that had potential functional redundancy due to co-expressed family members (12 genes) and genes that had inconsistent light responses across other publicly available microarray datasets (five genes). We next characterized the phenotypes of rice lines carrying mutations in ten of the remaining candidate genes and then carried out co-expression analysis associated with these genes. This analysis effectively provided candidate functions for two genes of previously unknown function and for one gene not directly linked to the tested biochemical pathways. These data demonstrate the efficiency of combining gene family-based expression profiles with analyses of insertional mutants to identify novel genes and their functions, even among members of multi-gene families.
Rice, a model monocot, is the first crop plant to have its entire genome sequenced. Although genome-wide transcriptome analysis tools and genome-wide, gene-indexed mutant collections have been generated for rice, the functions of only a handful of rice genes have been revealed thus far. Functional genomics approaches to studying crop plants like rice are much more labor-intensive and difficult in terms of maintaining the plants than when studying Arabidopsis, a model dicot. Here, we describe an efficient method for dissecting gene function in rice and other crop plants. We identified light response-related phenotypes for ten genes, the functions for which were previously unknown in rice. We also carried out co-expression analysis of 72 genes involved in specific biochemical pathways connected in lines carrying mutations in these ten genes. This analysis led to the identification of a novel set of genes likely involved in these pathways. The rapid progress of functional genomics in crops will significantly contribute to overcoming a food crisis in the near future.
Co-expression networks systematically constructed from large-scale transcriptome data reflect the interactions and functions of genes with similar expression patterns and are a powerful tool for the comprehensive understanding of biological events and mining of novel genes. In Arabidopsis (a model dicot plant), high-resolution co-expression networks have been constructed from very large microarray datasets and these are publicly available as online information resources. However, the available transcriptome data of rice (a model monocot plant) have been limited so far, making it difficult for rice researchers to achieve reliable co-expression analysis. In this study, we performed co-expression network analysis by using combined 44 K agilent microarray datasets of rice, which consisted of 33 laser microdissection (LM)-microarray datasets of anthers, and 143 spatiotemporal transcriptome datasets deposited in RicexPro. The entire data of the rice co-expression network, which was generated from the 176 microarray datasets by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off, contained 24,258 genes and 60,441 genes pairs. Using these datasets, we constructed high-resolution co-expression subnetworks of two specific biological events in the anther, “meiosis” and “pollen wall synthesis”. The meiosis network contained many known or putative meiotic genes, including genes related to meiosis initiation and recombination. In the pollen wall synthesis network, several candidate genes involved in the sporopollenin biosynthesis pathway were efficiently identified. Hence, these two subnetworks are important demonstrations of the efficiency of co-expression network analysis in rice. Our co-expression analysis included the separated transcriptomes of pollen and tapetum cells in the anther, which are able to provide precise information on transcriptional regulation during male gametophyte development in rice. The co-expression network data presented here is a useful resource for rice researchers to elucidate important and complex biological events.
Circadian clocks provide an adaptive advantage through anticipation of daily and seasonal environmental changes. In plants, the central clock oscillator is regulated by several interlocking feedback loops. It was shown that a substantial proportion of the Arabidopsis genome cycles with phases of peak expression covering the entire day. Synchronized transcriptome cycling is driven through an extensive network of diurnal and clock-regulated transcription factors and their target cis-regulatory elements. Study of the cycling transcriptome in other plant species could thus help elucidate the similarities and differences and identify hubs of regulation common to monocot and dicot plants.
Using a combination of oligonucleotide microarrays and data mining pipelines, we examined daily rhythms in gene expression in one monocotyledonous and one dicotyledonous plant, rice and poplar, respectively. Cycling transcriptomes were interrogated under different diurnal (driven) and circadian (free running) light and temperature conditions. Collectively, photocycles and thermocycles regulated about 60% of the expressed nuclear genes in rice and poplar. Depending on the condition tested, up to one third of oscillating Arabidopsis-poplar-rice orthologs were phased within three hours of each other suggesting a high degree of conservation in terms of rhythmic gene expression. We identified clusters of rhythmically co-expressed genes and searched their promoter sequences to identify phase-specific cis-elements, including elements that were conserved in the promoters of Arabidopsis, poplar, and rice.
Our results show that the cycling patterns of many circadian clock genes are highly conserved across poplar, rice, and Arabidopsis. The expression of many orthologous genes in key metabolic and regulatory pathways is diurnal and/or circadian regulated and phased to similar times of day. Our results confirm previous findings in Arabidopsis of three major classes of cis-regulatory modules within the plant circadian network: the morning (ME, GBOX), evening (EE, GATA), and midnight (PBX/TBX/SBX) modules. Identification of identical overrepresented motifs in the promoters of cycling genes from different species suggests that the core diurnal/circadian cis-regulatory network is deeply conserved between mono- and dicotyledonous species.
The response regulators represent the elements of bacterial two-component system and have been characterized from dicot plants like Arabidopsis but little information is available on the monocots, including the cereal crops. The aim of this study was to characterize type-A response regulator genes from rice, and to investigate their expression in various organs as well as in response to different hormones, including cytokinin, and environmental stimuli.
By analysis of the whole genome sequence of rice, we have identified ten genes encoding type-A response regulators based upon their high sequence identity within the receiver domain. The exon-intron organization, intron-phasing as well as chromosomal location of all the RT-PCR amplified rice (Oryza sativa) response regulator (OsRR) genes have been analyzed. The transcripts of OsRR genes could be detected by real-time PCR in all organs of the light- and dark-grown rice seedlings/plants, although there were quantitative differences. The steady-state transcript levels of most of the OsRR genes increased rapidly (within 15 min) on exogenous cytokinin application even in the presence of cycloheximide. Moreover, the expression of the OsRR6 gene was enhanced in rice seedlings exposed to salinity, dehydration and low temperature stress.
Ten type-A response regulator genes identified in rice, the model monocot plant, show overlapping/differential expression patterns in various organs and in response to light. The induction of OsRR genes by cytokinin even in the absence of de novo protein synthesis qualifies them to be primary cytokinin response genes. The induction of OsRR6 in response to different environmental stimuli indicates its role in cross-talk between abiotic stress and cytokinin signaling. These results provide a foundation for further investigations on specific as well as overlapping cellular functions of type-A response regulators in rice.
The bZIP class Abscisic acid Responsive Element (ABRE)-binding factor, OSBZ8 (38.5 kD) has been considered to regulate ABA-mediated transcription in the suspension cultured cells of japonica rice. Still, nothing is known about the expression of OSBZ8 at protein level in vegetative tissue of salt sensitive and salt tolerant rice plants. In our previous study, Electrophoretic Mobility Shift Assay (EMSA) of [32P]ABRE-DNA and nuclear extracts prepared from the lamina of Pokkali rice plants has detected the presence of an ABRE-binding factor. Northern analysis has also detected salinity stress induced accumulation of transcripts for bZIP class of factor. Therefore, OSBZ8 was considered to play an important role in the regulation of transcription in the vegetative tissue of rice. The aim of this study is to find out whether OSBZ8 has any role in regulating the NaCl-stress induced gene expression in vegetative tissue and whether the expression of OSBZ8 factor directly correlates with the stress tolerance of different varieties of indica type rice.
Northern analysis of total RNA from roots and lamina of salt-sensitive M-I-48 and salt-tolerant Nonabokra, when probed with the N-terminal unique region of OSBZ8 (OSBZ8p, without the highly conserved basic region), a transcript of 1.3 kb hybridized and its level was much higher in tolerant cultivar. EMSA with Em1a, the strongest ABA Responsive Element till reported from the upstream of EmBP1, and the nuclear extracts from laminar tissue of untreated and salt-treated seedlings of three salt sensitive, one moderately sensitive and two salt tolerant indica rice cultivars showed specific binding of nuclear factor to ABRE element. Intensity of binding was low and inducible in salt sensitive rice cultivars while high and constitutive in salt tolerant cultivars. EMSA with 300 bp 5'upstream region of Rab16A gene, a well known salt stress and ABA-inducible gene of rice, showed formation of two complexes, again very weak in salt sensitive and strong in salt tolerant rice cultivar.
The bZIP factor OSBZ8 was found to be present in the ABRE-DNA: protein complex as shown by the supershift of the complex by the purified antiserum raised against OSBZ8p. Treatment of the seedlings with NaCl was found to enhance the complex formation, suggesting the regulation of OSBZ8 gene at both transcriptional and post-translational steps. Comparative EMSA with different varieties of rice suggests a positive correlation with the expression pattern of OSBZ8 and salt tolerance in rice cultivars.
MADS-box transcription factors, besides being involved in floral organ specification, have also been implicated in several aspects of plant growth and development. In recent years, there have been reports on genomic localization, protein motif structure, phylogenetic relationships, gene structure and expression of the entire MADS-box family in the model plant system, Arabidopsis. Though there have been some studies in rice as well, an analysis of the complete MADS-box family along with a comprehensive expression profiling was still awaited after the completion of rice genome sequencing. Furthermore, owing to the role of MADS-box family in flower development, an analysis involving structure, expression and functional aspects of MADS-box genes in rice and Arabidopsis was required to understand the role of this gene family in reproductive development.
A genome-wide molecular characterization and microarray-based expression profiling of the genes encoding MADS-box transcription factor family in rice is presented. Using a thorough annotation exercise, 75 MADS-box genes have been identified in rice and categorized into MIKCc, MIKC*, Mα, Mβ and Mγ groups based on phylogeny. Chromosomal localization of these genes reveals that 16 MADS-box genes, mostly MIKCc-type, are located within the duplicated segments of the rice genome, whereas most of the M-type genes, 20 in all, seem to have resulted from tandem duplications. Nine members belonging to the Mβ group, which was considered absent in monocots, have also been identified. The expression profiles of all the MADS-box genes have been analyzed under 11 temporal stages of panicle and seed development, three abiotic stress conditions, along with three stages of vegetative development. Transcripts for 31 genes accumulate preferentially in the reproductive phase, of which, 12 genes are specifically expressed in seeds, and six genes show expression specific to panicle development. Differential expression of seven genes under stress conditions is also evident. An attempt has been made to gain insight into plausible functions of rice MADS-box genes by collating the expression data of functionally validated genes in rice and Arabidopsis.
Only a limited number of MADS genes have been functionally validated in rice. A comprehensive annotation and transcriptome profiling undertaken in this investigation adds to our understanding of the involvement of MADS-box family genes during reproductive development and stress in rice and also provides the basis for selection of candidate genes for functional validation studies.
The protein phosphatase 2Cs (PP2Cs) from various organisms have been implicated to act as negative modulators of protein kinase pathways involved in diverse environmental stress responses and developmental processes. A genome-wide overview of the PP2C gene family in plants is not yet available.
A comprehensive computational analysis identified 80 and 78 PP2C genes in Arabidopsis thaliana (AtPP2Cs) and Oryza sativa (OsPP2Cs), respectively, which denotes the PP2C gene family as one of the largest families identified in plants. Phylogenic analysis divided PP2Cs in Arabidopsis and rice into 13 and 11 subfamilies, respectively, which are supported by the analyses of gene structures and protein motifs. Comparative analysis between the PP2C genes in Arabidopsis and rice identified common and lineage-specific subfamilies and potential 'gene birth-and-death' events. Gene duplication analysis reveals that whole genome and chromosomal segment duplications mainly contributed to the expansion of both OsPP2Cs and AtPP2Cs, but tandem or local duplication occurred less frequently in Arabidopsis than rice. Some protein motifs are widespread among the PP2C proteins, whereas some other motifs are specific to only one or two subfamilies. Expression pattern analysis suggests that 1) most PP2C genes play functional roles in multiple tissues in both species, 2) the induced expression of most genes in subfamily A by diverse stimuli indicates their primary role in stress tolerance, especially ABA response, and 3) the expression pattern of subfamily D members suggests that they may constitute positive regulators in ABA-mediated signaling pathways. The analyses of putative upstream regulatory elements by two approaches further support the functions of subfamily A in ABA signaling, and provide insights into the shared and different transcriptional regulation machineries in dicots and monocots.
This comparative genome-wide overview of the PP2C family in Arabidopsis and rice provides insights into the functions and regulatory mechanisms, as well as the evolution and divergence of the PP2C genes in dicots and monocots. Bioinformatics analyses suggest that plant PP2C proteins from different subfamilies participate in distinct signaling pathways. Our results have established a solid foundation for future studies on the functional divergence in different PP2C subfamilies.
Alternative splicing (AS) contributes significantly to protein diversity, by selectively using different combinations of exons of the same gene under certain circumstances. One particular type of AS is the use of alternative first exons (AFEs), which can have consequences far beyond the fine-tuning of protein functions. For example, AFEs may change the N-termini of proteins and thereby direct them to different cellular compartments. When alternative first exons are distant, they are usually associated with alternative promoters, thereby conferring an extra level of gene expression regulation. However, only few studies have examined the patterns of AFEs, and these analyses were mainly focused on mammalian genomes. Recent studies have shown that AFEs exist in the rice genome, and are regulated in a tissue-specific manner. Our current understanding of AFEs in plants is still limited, including important issues such as their regulation, contribution to protein diversity, and evolutionary conservation.
We systematically identified 1,378 and 645 AFE-containing clusters in rice and Arabidopsis, respectively. From our data sets, we identified two types of AFEs according to their genomic organisation. In genes with type I AFEs, the first exons are mutually exclusive, while most of the downstream exons are shared among alternative transcripts. Conversely, in genes with type II AFEs, the first exon of one gene structure is an internal exon of an alternative gene structure. The functionality analysis indicated about half and ~19% of the AFEs in Arabidopsis and rice could alter N-terminal protein sequences, and ~5% of the functional alteration in type II AFEs involved protein domain addition/deletion in both genomes. Expression analysis indicated that 20~66% of rice AFE clusters were tissue- and/or development- specifically transcribed, which is consistent with previous observations; however, a much smaller percentage of Arabidopsis AFEs was regulated in this manner, which suggests different regulation mechanisms of AFEs between rice and Arabidopsis. Statistical analysis of some features of AFE clusters, such as splice-site strength and secondary structure formation further revealed differences between these two species. Orthologous search of AFE-containing gene pairs detected only 19 gene pairs conserved between rice and Arabidopsis, accounting only for a few percent of AFE-containing clusters.
Our analysis of AFE-containing genes in rice and Arabidopsis indicates that AFEs have multiple functions, from regulating gene expression to generating protein diversity. Comparisons of AFE clusters revealed different features in the two plant species, which indicates that AFEs may have evolved independently after the separation of rice (a model monocot) and Arabidopsis (a model dicot).
Heat shock proteins (Hsps) constitute an important component in the heat shock response of all living systems. Among the various plant Hsps (i.e. Hsp100, Hsp90, Hsp70 and Hsp20), Hsp20 or small Hsps (sHsps) are expressed in maximal amounts under high temperature stress. The characteristic feature of the sHsps is the presence of α-crystallin domain (ACD) at the C-terminus. sHsps cooperate with Hsp100/Hsp70 and co-chaperones in ATP-dependent manner in preventing aggregation of cellular proteins and in their subsequent refolding. Database search was performed to investigate the sHsp gene family across rice genome sequence followed by comprehensive expression analysis of these genes.
We identified 40 α-crystallin domain containing genes in rice. Phylogenetic analysis showed that 23 out of these 40 genes constitute sHsps. The additional 17 genes containing ACD clustered with Acd proteins of Arabidopsis. Detailed scrutiny of 23 sHsp sequences enabled us to categorize these proteins in a revised scheme of classification constituting of 16 cytoplasmic/nuclear, 2 ER, 3 mitochondrial, 1 plastid and 1 peroxisomal genes. In the new classification proposed herein nucleo-cytoplasmic class of sHsps with 9 subfamilies is more complex in rice than in Arabidopsis. Strikingly, 17 of 23 rice sHsp genes were noted to be intronless. Expression analysis based on microarray and RT-PCR showed that 19 sHsp genes were upregulated by high temperature stress. Besides heat stress, expression of sHsp genes was up or downregulated by other abiotic and biotic stresses. In addition to stress regulation, various sHsp genes were differentially upregulated at different developmental stages of the rice plant. Majority of sHsp genes were expressed in seed.
We identified twenty three sHsp genes and seventeen Acd genes in rice. Three nucleocytoplasmic sHsp genes were found only in monocots. Analysis of expression profiling of sHsp genes revealed that these genes are differentially expressed under stress and at different stages in the life cycle of rice plant.
Gibberellin (GA) 2-oxidases play an important role in the GA catabolic pathway through 2β-hydroxylation. There are two classes of GA2oxs, i.e., a larger class of C19-GA2oxs and a smaller class of C20-GA2oxs. In this study, the gene encoding a GA 2-oxidase of rice, Oryza sativa GA 2-oxidase 5 (OsGA2ox5), was cloned and characterized. BLASTP analysis showed that OsGA2ox5 belongs to the C20-GA2oxs subfamily, a subfamily of GA2oxs acting on C20-GAs (GA12, GA53). Subcellular localization of OsGA2ox5-YFP in transiently transformed onion epidermal cells revealed the presence of this protein in both of the nucleus and cytoplasm. Real-time PCR analysis, along with GUS staining, revealed that OsGA2ox5 is expressed in the roots, culms, leaves, sheaths and panicles of rice. Rice plants overexpressing OsGA2ox5 exhibited dominant dwarf and GA-deficient phenotypes, with shorter stems and later development of reproductive organs than the wild type. The dwarfism phenotype was partially rescued by the application of exogenous GA3 at a concentration of 10 µM. Ectopic expression of OsGA2ox5 cDNA in Arabidopsis resulted in a similar phenotype. Real-time PCR assays revealed that both GA synthesis-related genes and GA signaling genes were expressed at higher levels in transgenic rice plants than in wild-type rice; OsGA3ox1, which encodes a key enzyme in the last step of the bioactive GAs synthesis pathway, was highly expressed in transgenic rice. The roots of OsGA2ox5-ox plants exhibited increased starch granule accumulation and gravity responses, revealing a role for GA in root starch granule development and gravity responses. Furthermore, rice and Arabidopsis plants overexpressing OsGA2ox5 were more resistant to high-salinity stress than wild-type plants. These results suggest that OsGA2ox5 plays important roles in GAs homeostasis, development, gravity responses and stress tolerance in rice.
The development of rice (Oryza sativa) seed is closely associated with assimilates storage and plant yield, and is fine controlled by complex regulatory networks. Exhaustive transcriptome analysis of developing rice embryo and endosperm will help to characterize the genes possibly involved in the regulation of seed development and provide clues of yield and quality improvement.
Our analysis showed that genes involved in metabolism regulation, hormone response and cellular organization processes are predominantly expressed during rice development. Interestingly, 191 transcription factor (TF)-encoding genes are predominantly expressed in seed and 59 TFs are regulated during seed development, some of which are homologs of seed-specific TFs or regulators of Arabidopsis seed development. Gene co-expression network analysis showed these TFs associated with multiple cellular and metabolism pathways, indicating a complex regulation of rice seed development. Further, by employing a cold-resistant cultivar Hanfeng (HF), genome-wide analyses of seed transcriptome at normal and low temperature reveal that rice seed is sensitive to low temperature at early stage and many genes associated with seed development are down-regulated by low temperature, indicating that the delayed development of rice seed by low temperature is mainly caused by the inhibition of the development-related genes. The transcriptional response of seed and seedling to low temperature is different, and the differential expressions of genes in signaling and metabolism pathways may contribute to the chilling tolerance of HF during seed development.
These results provide informative clues and will significantly improve the understanding of rice seed development regulation and the mechanism of cold response in rice seed.
In plants, complex regulatory mechanisms are at the core of physiological and developmental processes. The phytohormone abscisic acid (ABA) is involved in the regulation of various such processes, including stomatal closure, seed and bud dormancy, and physiological responses to cold, drought and salinity stress. The underlying tissue or plant-wide control circuits often include combinatorial gene regulatory mechanisms and networks that we are only beginning to unravel with the help of new molecular tools. The increasing availability of genomic sequences and gene expression data enables us to dissect ABA regulatory mechanisms at the individual gene expression level. In this paper we used an in-silico-based approach directed towards genome-wide prediction and identification of specific features of ABA-responsive elements. In particular we analysed the genome-wide occurrence and positional arrangements of two well-described ABA-responsive cis-regulatory elements (CREs), ABRE and CE3, in thale cress (Arabidopsis thaliana) and rice (Oryza sativa).
Our results show that Arabidopsis and rice use the ABA-responsive elements ABRE and CE3 distinctively. Earlier reports for various monocots have identified CE3 as a coupling element (CE) associated with ABRE. Surprisingly, we found that while ABRE is equally abundant in both species, CE3 is practically absent in Arabidopsis. ABRE-ABRE pairs are common in both genomes, suggesting that these can form functional ABA-responsive complexes (ABRCs) in Arabidopsis and rice. Furthermore, we detected distinct combinations, orientation patterns and DNA strand preferences of ABRE and CE3 motifs in rice gene promoters.
Our computational analyses revealed distinct recruitment patterns of ABA-responsive CREs in upstream sequences of Arabidopsis and rice. The apparent absence of CE3s in Arabidopsis suggests that another CE pairs with ABRE to establish a functional ABRC capable of interacting with transcription factors. Further studies will be needed to test whether the observed differences are extrapolatable to monocots and dicots in general, and to understand how they contribute to the fine-tuning of the hormonal response. The outcome of our investigation can now be used to direct future experimentation designed to further dissect the ABA-dependent regulatory networks.
Cereal endosperm represents 60% of the calories consumed by human beings worldwide. In addition, cereals also serve as the primary feedstock for livestock. However, the regulatory mechanism of cereal endosperm and seed development is largely unknown. Polycomb complex has been shown to play a key role in the regulation of endosperm development in Arabidopsis, but its role in cereal endosperm development remains obscure. Additionally, the enzyme activities of the polycomb complexes have not been demonstrated in plants. Here we purified the rice OsFIE2-polycomb complex using tandem affinity purification and demonstrated its specific H3 methyltransferase activity. We found that the OsFIE2 gene product was responsible for H3K27me3 production specifically in vivo. Genetic studies showed that a reduction of OsFIE2 expression led to smaller seeds, partially filled seeds, and partial loss of seed dormancy. Gene expression and proteomics analyses found that the starch synthesis rate limiting step enzyme and multiple storage proteins are down-regulated in OsFIE2 reduction lines. Genome wide ChIP–Seq data analysis shows that H3K27me3 is associated with many genes in the young seeds. The H3K27me3 modification and gene expression in a key helix-loop-helix transcription factor is shown to be regulated by OsFIE2. Our results suggest that OsFIE2-polycomb complex positively regulates rice endosperm development and grain filling via a mechanism highly different from that in Arabidopsis.
Rice is the staple food for over half of the world's population and an important feedstock for livestock. The rice grain is mainly endosperm tissue. The regulatory mechanism of rice endosperm development is still largely unknown thus far. Understanding the underlying mechanism will lead to crop yield and quality improvement in the long term, besides gaining new knowledge. Polycomb complex is a protein complex with a potential role in endosperm development according to prior publications. In this manuscript, we purified the rice OsFIE2-polycomb protein complex and demonstrated the enzyme activity of the complex. Genetic studies showed that a reduction of polycomb group gene OsFIE2 expression led to smaller seeds, partially filled seeds, and seed germination before seed maturation. Gene expression and proteomics analyses found that the starch synthesis rate limiting step enzyme and multiple storage proteins are down-regulated while a key transcription factor is up-regulated in OsFIE2 reduction lines. In addition, we identified many loci in the rice genome whose histone proteins are modified by the polycomb complex enzyme via a method called ChIP–Seq. Our results demonstrate that OsFIE2-polycomb complex positively regulates rice grain development via a mechanism distinct from that in Arabidopsis and provide new insight into the regulation of rice grain development.
Plant diurnal oscillation is a 24-hour period based variation. The correlation between diurnal genes and biological pathways was widely revealed by microarray analysis in different species. Rice (Oryza sativa) is the major food staple for about half of the world's population. The rice flag leaf is essential in providing photosynthates to the grain filling. However, there is still no comprehensive view about the diurnal transcriptome for rice leaves. In this study, we applied rice microarray to monitor the rhythmically expressed genes in rice seedling and flag leaves. We developed a new computational analysis approach and identified 6,266 (10.96%) diurnal probe sets in seedling leaves, 13,773 (24.08%) diurnal probe sets in flag leaves. About 65% of overall transcription factors were identified as flag leaf preferred. In seedling leaves, the peak of phase distribution was from 2:00am to 4:00am, whereas in flag leaves, the peak was from 8:00pm to 2:00am. The diurnal phase distribution analysis of gene ontology (GO) and cis-element enrichment indicated that, some important processes were waken by the light, such as photosynthesis and abiotic stimulus, while some genes related to the nuclear and ribosome involved processes were active mostly during the switch time of light to dark. The starch and sucrose metabolism pathway genes also showed diurnal phase. We conducted comparison analysis between Arabidopsis and rice leaf transcriptome throughout the diurnal cycle. In summary, our analysis approach is feasible for relatively unbiased identification of diurnal transcripts, efficiently detecting some special periodic patterns with non-sinusoidal periodic patterns. Compared to the rice flag leaves, the gene transcription levels of seedling leaves were relatively limited to the diurnal rhythm. Our comprehensive microarray analysis of seedling and flag leaves of rice provided an overview of the rice diurnal transcriptome and indicated some diurnal regulated biological processes and key functional pathways in rice.
Plants respond to low temperature through an intricately coordinated transcriptional network. The CBF/DREB-regulated network of genes has been shown to play a prominent role in freeze-tolerance of Arabidopsis through the process of cold acclimation (CA). Recent evidence also showed that the CBF/DREB regulon is not unique to CA but evolutionarily conserved between chilling-insensitive (temperate) and chilling-sensitive (warm-season) plants. In this study, the wide contrast in chilling sensitivity between indica and japonica rice was used as model to identify other regulatory clusters by integrative analysis of promoter architecture (ab initio) and gene expression profiles.
Transcriptome analysis in chilling tolerant japonica rice identified a subset of 121 'early response' genes that were upregulated during the initial 24 hours at 10°C. Among this group were four transcription factors including ROS-bZIP1 and another larger sub-group with a common feature of having as1/ocs-like elements in their promoters. Cold-induction of ROS-bZIP1 preceded the induction of as1/ocs-like element-containing genes and they were also induced by exogenous H2O2 at ambient temperature. Coordinated expression patterns and similar promoter architectures among the 'early response' genes suggest that they belong to a potential regulon (ROS-bZIP – as1/ocs regulatory module) that responds to elevated levels of ROS during chilling stress. Cultivar-specific expression signatures of the candidate genes indicate a positive correlation between the activity of the putative regulon and genotypic variation in chilling tolerance.
A hypothetical model of an ROS-mediated regulon (ROS-bZIP – as1/ocs) triggered by chilling stress was assembled in rice. Based on the current results, it appears that this regulon is independent of ABA and CBF/DREB, and that its activation has an important contribution in configuring the rapid responses of rice seedlings to chilling stress.
Cis-natural antisense transcripts (cis-NATs) are RNAs transcribed from the antisense strand of a gene locus, and are complementary to the RNA transcribed from the sense strand. Common techniques including microarray approach and analysis of transcriptome databases are the major ways to globally identify cis-NATs in various eukaryotic organisms. Genome-wide in silico analysis has identified a large number of cis-NATs that may generate endogenous short interfering RNAs (nat-siRNAs), which participate in important biogenesis mechanisms for transcriptional and post-transcriptional regulation in rice. However, the transcriptomes are yet to be deeply sequenced to comprehensively investigate cis-NATs.
We applied high-throughput strand-specific complementary DNA sequencing technology (ssRNA-seq) to deeply sequence mRNA for assessing sense and antisense transcripts that were derived under salt, drought and cold stresses, and normal conditions, in the model plant rice (Oryza sativa). Combined with RAP-DB genome annotation (the Rice Annotation Project Database build-5 data set), 76,013 transcripts corresponding to 45,844 unique gene loci were assembled, in which 4873 gene loci were newly identified. Of 3819 putative rice cis-NATs, 2292 were detected as expressed and giving rise to small RNAs from their overlapping regions through integrated analysis of ssRNA-seq data and small RNA data. Among them, 503 cis-NATs seemed to be associated with specific conditions. The deep sequence data from isolated epidermal cells of rice seedlings further showed that 54.0% of cis-NATs were expressed simultaneously in a population of homogenous cells. Nearly 9.7% of rice transcripts were involved in one-to-one or many-to-many cis-NATs formation. Furthermore, only 17.4-34.7% of 223 many-to-many cis-NAT groups were all expressed and generated nat-siRNAs, indicating that only some cis-NAT groups may be involved in complex regulatory networks.
Our study profiles an abundance of cis-NATs and nat-siRNAs in rice. These data are valuable for gaining insight into the complex function of the rice transcriptome.
Oryza sativa; Cis-NATs; Nat-siRNAs; SsRNA-seq; Transcriptome
Glycosyl hydrolase family 1 (GH1) β-glucosidases have been implicated in physiologically important processes in plants, such as response to biotic and abiotic stresses, defense against herbivores, activation of phytohormones, lignification, and cell wall remodeling. Plant GH1 β-glucosidases are encoded by a multigene family, so we predicted the structures of the genes and the properties of their protein products, and characterized their phylogenetic relationship to other plant GH1 members, their expression and the activity of one of them, to begin to decipher their roles in rice.
Forty GH1 genes could be identified in rice databases, including 2 possible endophyte genes, 2 likely pseudogenes, 2 gene fragments, and 34 apparently competent rice glycosidase genes. Phylogenetic analysis revealed that GH1 members with closely related sequences have similar gene structures and are often clustered together on the same chromosome. Most of the genes appear to have been derived from duplications that occurred after the divergence of rice and Arabidopsis thaliana lineages from their common ancestor, and the two plants share only 8 common gene lineages. At least 31 GH1 genes are expressed in a range of organs and stages of rice, based on the cDNA and EST sequences in public databases. The cDNA of the Os4bglu12 gene, which encodes a protein identical at 40 of 44 amino acid residues with the N-terminal sequence of a cell wall-bound enzyme previously purified from germinating rice, was isolated by RT-PCR from rice seedlings. A thioredoxin-Os4bglu12 fusion protein expressed in Escherichia coli efficiently hydrolyzed β-(1,4)-linked oligosaccharides of 3–6 glucose residues and laminaribiose.
Careful analysis of the database sequences produced more reliable rice GH1 gene structure and protein product predictions. Since most of these genes diverged after the divergence of the ancestors of rice and Arabidopsis thaliana, only a few of their functions could be implied from those of GH1 enzymes from Arabidopsis and other dicots. This implies that analysis of GH1 enzymes in monocots is necessary to understand their function in the major grain crops. To begin this analysis, Os4bglu12 β-glucosidase was characterized and found to have high exoglucanase activity, consistent with a role in cell wall metabolism.
Iron (Fe) is the most limiting micronutrient element for crop production in alkaline soils. A number of transcription factors involved in regulating Fe uptake from soil and transport in plants have been identified. Analysis of transcriptome data from Oryza sativa grown under limiting Fe conditions reveals that transcript abundances of several genes encoding transcription factors are altered by Fe availability. These transcription factors are putative regulators of Fe deficiency responses.
Transcript abundance of one nuclear located basic helix-loop-helix family transcription factor, OsIRO3, is up-regulated from 25- to 90-fold under Fe deficiency in both root and shoot respectively. The expression of OsIRO3 is specifically induced by Fe deficiency, and not by other micronutrient deficiencies. Transgenic rice plants over-expressing OsIRO3 were hypersensitive to Fe deficiency, indicating that the Fe deficiency response was compromised. Furthermore, the Fe concentration in shoots of transgenic rice plants over-expressing OsIRO3 was less than that in wild-type plants. Analysis of the transcript abundances of genes normally induced by Fe deficiency in OsIRO3 over-expressing plants indicated their induction was markedly suppressed.
A novel Fe regulated bHLH transcription factor (OsIRO3) that plays an important role for Fe homeostasis in rice was identified. The inhibitory effect of OsIRO3 over-expression on Fe deficiency response gene expression combined with hypersensitivity of OsIRO3 over-expression lines to low Fe suggest that OsIRO3 is a negative regulator of the Fe deficiency response in rice.
Rice blast disease caused by Magnaporthe oryzae is one of the most serious diseases of cultivated rice (Oryza sativa L.) in most rice-growing regions of the world. In order to investigate early response genes in rice, we utilized the transcriptome analysis approach using a 300 K tilling microarray to rice leaves infected with compatible and incompatible M. oryzae strains. Prior to the microarray experiment, total RNA was validated by measuring the differential expression of rice defense-related marker genes (chitinase 2, barwin, PBZ1, and PR-10) by RT-PCR, and phytoalexins (sakuranetin and momilactone A) with HPLC. Microarray analysis revealed that 231 genes were up-regulated (>2 fold change, p < 0.05) in the incompatible interaction compared to the compatible one. Highly expressed genes were functionally characterized into metabolic processes and oxidation-reduction categories. The oxidative stress response was induced in both early and later infection stages. Biotic stress overview from MapMan analysis revealed that the phytohormone ethylene as well as signaling molecules jasmonic acid and salicylic acid is important for defense gene regulation. WRKY and Myb transcription factors were also involved in signal transduction processes. Additionally, receptor-like kinases were more likely associated with the defense response, and their expression patterns were validated by RT-PCR. Our results suggest that candidate genes, including receptor-like protein kinases, may play a key role in disease resistance against M. oryzae attack.
defense response; Magnaporthe oryzae; MapMan analysis; rice; transcriptomics
Plants are simultaneously exposed to multiple stresses resulting in enormous changes in the molecular landscape within the cell. Identification and characterization of the synergistic and antagonistic components of stress response mechanisms contributing to the cross talk between stresses is of high priority to explore and enhance multiple stress responses. To this end, we performed meta-analysis of drought (abiotic), bacterial (biotic) stress response in rice and Arabidopsis by analyzing a total of 386 microarray samples belonging to 20 microarray studies and identified approximately 3100 and 900 DEGs in rice and Arabidopsis, respectively. About 38.5% (1214) and 28.7% (272) DEGs were common to drought and bacterial stresses in rice and Arabidopsis, respectively. A majority of these common DEGs showed conserved expression status in both stresses. Gene ontology enrichment analysis clearly demarcated the response and regulation of various plant hormones and related biological processes. Fatty acid metabolism and biosynthesis of alkaloids were upregulated and, nitrogen metabolism and photosynthesis was downregulated in both stress conditions. WRKY transcription family genes were highly enriched in all upregulated gene sets while ‘CO-like’ TF family showed inverse relationship of expression between drought and bacterial stresses. Weighted gene co-expression network analysis divided DEG sets into multiple modules that show high co-expression and identified stress specific hub genes with high connectivity. Detection of consensus modules based on DEGs common to drought and bacterial stress revealed 9 and 4 modules in rice and Arabidopsis, respectively, with conserved and reversed co-expression patterns.