A microarray analysis of the rice transcriptome and its comparison to Arabidopsis

Transcriptomic responses of oil palm (Elaeis guineensis) stem to waterlogging at plantation in relation to precipitation seasonality

Global warming-induced climate change causes significant agricultural problems by increasing the incidence of drought and flooding events. Waterlogging is an inevitable consequence of these changes but its effects on oil palms have received little attention and are poorly understood. Recent waterlogging studies have focused on oil palm seedlings, with particular emphasis on phenology. However, the transcriptomic waterlogging response of mature oil palms remains elusive in real environments. We therefore investigated transcriptomic changes over time in adult oil palms at plantations over a two-year period with pronounced seasonal variation in precipitation. A significant transcriptional waterlogging response was observed in the oil palm stem core but not in leaf samples when gene expression was correlated with cumulative precipitation over two-day periods. Pathways and processes upregulated or enriched in the stem core response included hypoxia, ethylene signaling, and carbon metabolism. Post-waterlogging recovery in oil palms was found to be associated with responses to heat stress and carotenoid biosynthesis. Nineteen transcription factors (TFs) potentially involved in the waterlogging response of mature oil palms were also identified. These data provide new insights into the transcriptomic responses of planted oil palms to waterlogging and offer valuable guidance on the sensitivity of oil palm plantations to future climate changes.

Rice RS2-9, which is bound by transcription factor OSH1, blocks enhancer-promoter interactions in plants

Insulators characterized in Drosophila and mammals have been shown to play a key role in the restriction of promiscuous enhancer-promoter interactions, as well as reshaping the topological landscape of chromosomes. Yet the role of insulators in plants remains poorly understood, in large part because of a lack of well-characterized insulators and binding factor(s). In this study, we isolated a 1.2-kb RS2-9 insulator from the Oryza sativa (rice) genome that can, when interposed between an enhancer and promoter, efficiently block the activation function of both constitutive and floral organ-specific enhancers in transgenic Arabidopsis and Nicotiana tabacum (tobacco). In the rice genome, the genes flanking RS2-9 exhibit an absence of mutual transcriptional interactions, as well as a lack of histone modification spread. We further determined that O. sativa Homeobox 1 (OSH1) bound two regions of RS2-9, as well as over 50 000 additional sites in the rice genome, the majority of which resided in intergenic regions. Mutation of one of the two OSH1-binding sites in RS2-9 impaired insulation activity by up to 60%, whereas the mutation of both binding sites virtually abolished insulator function. We also demonstrated that OSH1 binding sites were associated with 72% of the boundaries of topologically associated domains (TADs) identified in the rice genome, which is comparable to the 77% of TAD boundaries bound by the insulator CCCTC-binding factor (CTCF) in mammals. Taken together, our findings indicate that OSH1-RS2-9 acts as a true insulator in plants, and highlight a potential role for OSH1 in gene insulation and topological organization in plant genomes.

Functional Omics Identifies Serine Hydrolases That Mobilize Storage Lipids during Rice Seed Germination

Elucidating proteolipidome dynamics is crucial for understanding the roles of these molecules in plant physiology and disease. Sequence-based functional annotation of the protein is inadequate, since protein activities depend on posttranslational modification. In this study, we applied a gel-free activity-based protein profiling approach to unravel the active lipases, including other Serine hydrolases (SHs), expressed during seed germination in rice (Oryza sativa). We successfully mapped the active sites of 43 active SHs encompassing lipases/esterases, GDSL lipases, proteases, Ser carboxypeptidases, ABHD protein, pectin acetylesterase, and other SHs. The mRNA expression levels of those genes encoding the identified SHs were monitored using microarray analysis. The lipidome analysis revealed distinct patterns of molecular species distribution in individual lipid classes and displayed the metabolic connections between lipid mobilization and rice seedling growth. Changes in the mobilization of storage lipids and their molecular species remodeling were correlated with the expression of the identified lipases and their lipase activity in a time-dependent manner. The physiological significance of the identified SHs was explored during biotic stress with Fusarium verticillioides infection. The fungal infection significantly reduced lipase activity and lipid mobilization, thus impairing the rice seedling. Collectively, our data demonstrate application of the functional proteome strategy along with the shotgun lipidome approach for the identification of active SHs, and thus for deciphering the role of lipid homeostasis during rice seed germination.

Mechanisms of Cryptochrome-Mediated Photoresponses in Plants

Cryptochromes are blue-light receptors that mediate photoresponses in plants. The genomes of most land plants encode two clades of cryptochromes, CRY1 and CRY2, which mediate distinct and overlapping photoresponses within the same species and between different plant species. Photoresponsive protein-protein interaction is the primary mode of signal transduction of cryptochromes. Cryptochromes exist as physiologically inactive monomers in the dark; the absorption of photons leads to conformational change and cryptochrome homooligomerization, which alters the affinity of cryptochromes interacting with cryptochrome-interacting proteins to form various cryptochrome complexes. These cryptochrome complexes, collectively referred to as the cryptochrome complexome, regulate transcription or stability of photoresponsive proteins to modulate plant growth and development. The activity of cryptochromes is regulated by photooligomerization; dark monomerization; cryptochrome regulatory proteins; and cryptochrome phosphorylation, ubiquitination, and degradation. Most of the more than 30 presently known cryptochrome-interacting proteins are either regulated by other photoreceptors or physically interactingwith the protein complexes of other photoreceptors. Some cryptochrome-interacting proteins are also hormonal signaling or regulatory proteins. These two mechanisms enable cryptochromes to integrate blue-light signals with other internal and external signals to optimize plant growth and development.

Integrated Genome-Scale Analysis Identifies Novel Genes and Networks Underlying Senescence in Maize

Premature senescence in annual crops reduces yield, while delayed senescence, termed stay-green, imposes positive and negative impacts on yield and nutrition quality. Despite its importance, scant information is available on the genetic architecture of senescence in maize (Zea mays) and other cereals. We combined a systematic characterization of natural diversity for senescence in maize and coexpression networks derived from transcriptome analysis of normally senescing and stay-green lines. Sixty-four candidate genes were identified by genome-wide association study (GWAS), and 14 of these genes are supported by additional evidence for involvement in senescence-related processes including proteolysis, sugar transport and signaling, and sink activity. Eight of the GWAS candidates, independently supported by a coexpression network underlying stay-green, include a trehalose-6-phosphate synthase, a NAC transcription factor, and two xylan biosynthetic enzymes. Source-sink communication and the activity of cell walls as a secondary sink emerge as key determinants of stay-green. Mutant analysis supports the role of a candidate encoding Cys protease in stay-green in Arabidopsis (Arabidopsis thaliana), and analysis of natural alleles suggests a similar role in maize. This study provides a foundation for enhanced understanding and manipulation of senescence for increasing carbon yield, nutritional quality, and stress tolerance of maize and other cereals.

Investigation of Drought and Salinity Tolerance Related Genes and their Regulatory Mechanisms in Arabidopsis (Arabidopsis thaliana)

Background:

The development of genome microarrays of the model plant;Arabidopsis thaliana, with increasing repositories of publicly available data and high-throughput data analysis tools, has opened new avenues to genome-wide systemic analysis of plant responses to environmental stresses.

Objective:

To identify differentially expressed genes and their regulatory networks inArabidopsis thalianaunder harsh environmental condition.

Methods:

Two replications of eight microarray data sets were derived from two different tissues (root and shoot) and two different time courses (control and 24 hours after the beginning of stress occurrence) for comparative data analysis through various bioinformatics tools.

Results:

Under drought stress, 2558 gene accessions in root and 3691 in shoot tissues had significantly differential expression with respect to control condition. Likewise, under salinity stress 9078 gene accessions in root and 5785 in shoot tissues were discriminated between stressed and non-stressed conditions. Furthermore, the transcription regulatory activity of differentially expressed genes was mainly due to hormone, light, circadian and stress responsivecis-acting regulatory elements among which ABRE, ERE, P-box, TATC-box, CGTCA-motif, GARE-motif, TGACG-motif, GAG-motif, GA-motif, GATA- motif, TCT-motif, GT1-motif, Box 4, G-Box, I-box, LAMP-element, Sp1, MBS, TC-rich repeats, TCA-element and HSE were the most important elements in the identified up-regulated genes.

Conclusion:

The results of the high-throughput comparative analyses in this study provide more options for plant breeders and give an insight into genes andcis-acting regulatory elements involved in plant response to drought and salinity stresses in strategic crops such as cereals.

Analysis of Genes with Alternatively Spliced Transcripts in the Leaf, Root, Panicle and Seed of Rice Using a Long Oligomer Microarray and RNA-Seq

Pre-mRNA splicing further increases protein diversity acquired through evolution. The underlying driving forces for this phenomenon are unknown, especially in terms of gene expression. A rice alternatively spliced transcript detection microarray (ASDM) and RNA sequencing (RNA-Seq) were applied to differentiate the transcriptome of 4 representative organs of Oryza sativa L. cv. Ilmi: leaves, roots, 1-cm-stage panicles and young seeds at 21 days after pollination. Comparison of data obtained by microarray and RNA-Seq showed a bell-shaped distribution and a co-lineation for highly expressed genes. Transcripts were classified according to the degree of organ enrichment using a coefficient value (CV, the ratio of the standard deviation to the mean values): highly variable (CVI), variable (CVII), and constitutive (CVIII) groups. A higher index of the portion of loci with alternatively spliced transcripts in a group (IAST) value was observed for the constitutive group. Genes of the highly variable group showed the characteristics of the examined organs, and alternatively spliced transcripts tended to exhibit the same organ specificity or less organ preferences, with avoidance of 'organ distinctness'. In addition, within a locus, a tendency of higher expression was found for transcripts with a longer coding sequence (CDS), and a spliced intron was the most commonly found type of alternative splicing for an extended CDS. Thus, pre-mRNA splicing might have evolved to retain maximum functionality in terms of organ preference and multiplicity.

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

Background

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

Results

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

Conclusions

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

Database URL: http://14.139.229.201

Electronic supplementary material

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

Bioinformatic Identification of Conserved Cis-Sequences in Coregulated Genes

Bioinformatics tools can be employed to identify conserved cis-sequences in sets of coregulated plant genes because more and more gene expression and genomic sequence data become available. Knowledge on the specific cis-sequences, their enrichment and arrangement within promoters, facilitates the design of functional synthetic plant promoters that are responsive to specific stresses. The present chapter illustrates an example for the bioinformatic identification of conserved Arabidopsis thaliana cis-sequences enriched in drought stress-responsive genes. This workflow can be applied for the identification of cis-sequences in any sets of coregulated genes. The workflow includes detailed protocols to determine sets of coregulated genes, to extract the corresponding promoter sequences, and how to install and run a software package to identify overrepresented motifs. Further bioinformatic analyses that can be performed with the results are discussed.

Deep transcriptome sequencing provides new insights into the structural and functional organization of the wheat genome

Background

Because of its size, allohexaploid nature, and high repeat content, the bread wheat genome is a good model to study the impact of the genome structure on gene organization, function, and regulation. However, because of the lack of a reference genome sequence, such studies have long been hampered and our knowledge of the wheat gene space is still limited. The access to the reference sequence of the wheat chromosome 3B provided us with an opportunity to study the wheat transcriptome and its relationships to genome and gene structure at a level that has never been reached before.

Results

By combining this sequence with RNA-seq data, we construct a fine transcriptome map of the chromosome 3B. More than 8,800 transcription sites are identified, that are distributed throughout the entire chromosome. Expression level, expression breadth, alternative splicing as well as several structural features of genes, including transcript length, number of exons, and cumulative intron length are investigated. Our analysis reveals a non-monotonic relationship between gene expression and structure and leads to the hypothesis that gene structure is determined by its function, whereas gene expression is subject to energetic cost. Moreover, we observe a recombination-based partitioning at the gene structure and function level.

Conclusions

Our analysis provides new insights into the relationships between gene and genome structure and function. It reveals mechanisms conserved with other plant species as well as superimposed evolutionary forces that shaped the wheat gene space, likely participating in wheat adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0601-9) contains supplementary material, which is available to authorized users.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii)

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

Comparative transcriptional profiling of three super-hybrid rice combinations

Utilization of heterosis has significantly increased rice yields. However, its mechanism remains unclear. In this study, comparative transcriptional profiles of three super-hybrid rice combinations, LY2163, LY2186 and LYP9, at the flowering and filling stages, were created using rice whole-genome oligonucleotide microarray. The LY2163, LY2186 and LYP9 hybrids yielded 1193, 1630 and 1046 differentially expressed genes (DGs), accounting for 3.2%, 4.4% and 2.8% of the total number of genes (36,926), respectively, after using the z-test (p < 0.01). Functional category analysis showed that the DGs in each hybrid combination were mainly classified into the carbohydrate metabolism and energy metabolism categories. Further analysis of the metabolic pathways showed that DGs were significantly enriched in the carbon fixation pathway (p < 0.01) for all three combinations. Over 80% of the DGs were located in rice quantitative trait loci (QTLs) of the Gramene database, of which more than 90% were located in the yield related QTLs in all three combinations, which suggested that there was a correlation between DGs and rice heterosis. Pathway Studio analysis showed the presence of DGs in the circadian regulatory network of all three hybrid combinations, which suggested that the circadian clock had a role in rice heterosis. Our results provide information that can help to elucidate the molecular mechanism underlying rice heterosis.

Transcriptome profiling analysis using rice oligonucleotide microarrays

Transcriptome analysis using oligonucleotide microarrays is a powerful tool for detecting changes in genome-wide transcripts under a given biological condition. Although the rice genome sequence is available, the number of functionally characterized genes in rice is still very limited. Genome-wide transcriptome analysis is a useful tool for elucidating the functions of rice genes that have not yet been determined. Currently, more than 3,000 arrays are publicly available. Here, we introduce methods for genome-wide transcriptome analysis in rice.

Comparative co-expression analysis in plant biology

The analysis of gene expression data generated by high-throughput microarray transcript profiling experiments has shown that transcriptionally coordinated genes are often functionally related. Based on large-scale expression compendia grouping multiple experiments, this guilt-by-association principle has been applied to study modular gene programmes, identify cis-regulatory elements or predict functions for unknown genes in different model plants. Recently, several studies have demonstrated how, through the integration of gene homology and expression information, correlated gene expression patterns can be compared between species. The incorporation of detailed functional annotations as well as experimental data describing protein-protein interactions, phenotypes or tissue specific expression, provides an invaluable source of information to identify conserved gene modules and translate biological knowledge from model organisms to crops. In this review, we describe the different steps required to systematically compare expression data across species. Apart from the technical challenges to compute and display expression networks from multiple species, some future applications of plant comparative transcriptomics are highlighted.

Identification of potential antisense transcripts in rice using conventional microarray

Natural antisense transcripts (NATs) are endogenous transcripts that contain reverse complementary sequences to other RNAs (usually called sense transcripts). NATs regulate the expression of sense transcripts in a wide range of species. The identification and analysis of NATs are the prerequisite to elucidate their functions. Microarray is a genome-wide method to detect gene expression. However, conventional microarrays do not contain the specific probes of NATs; thus, they cannot be utilized to detect NATs. In this article, we developed a novel method to identify potential NATs with the conventional microarrays. In this method of our study, we labeled the first strand cDNA from one sample with Cy5 and labeled the second strand cDNA from another sample with Cy3, and then hybridized these labeled samples with oligonucleotide microarray. Using this method, we identified 920 potential NATs in rice variety Nipponbare. Among these potential NATs, 88 of them were confirmed by either full-length cDNA or orientated ESTs (expressed sequence tags). This is the first time that a conventional oligonucleotide microarray was employed to identify NATs in rice.

Evaluation of alternative RNA labeling protocols for transcript profiling with Arabidopsis AGRONOMICS1 tiling arrays

Microarrays are routine tools for transcript profiling, and genomic tiling arrays such as the Arabidopsis AGRONOMICS1 arrays have been found to be highly suitable for such experiments because changes in genome annotation can be easily integrated at the data analysis level. In a transcript profiling experiment, RNA labeling is a critical step, most often initiated by oligo-dT-primed reverse transcription. Although this has been found to be a robust and reliable method, very long transcripts or non-polyadenylated transcripts might be labeled inefficiently. In this study, we first provide data handling methods to analyze AGRONOMICS1 tiling microarrays based on the TAIR10 genome annotation. Second, we describe methods to easily quantify antisense transcripts on such tiling arrays. Third, we test a random-primed RNA labeling method, and find that on AGRONOMICS1 arrays this method has similar general performance as the conventional oligo-dT-primed method. In contrast to the latter, however, the former works considerably better for long transcripts and for non-polyadenylated transcripts such as found in mitochondria and plastids. We propose that researchers interested in organelle function use the random-primed method to unleash the full potential of genomic tiling arrays.

Tomato plants overexpressing cryptochrome 2 reveal altered expression of energy and stress-related gene products in response to diurnal cues

In order to sense and respond to the fluctuating light conditions, higher plants possess several families of photoreceptors, such as phytochromes (PHYs), cryptochromes (CRYs) and phototropins. CRYs are responsible for photomorphogenesis and play a role in circadian, developmental and adaptive growth regulation of plants. In tomato (Solanum lycopersicum), CRY2 controls vegetative development, flowering time, fruit antioxidant content as well as the diurnal transcription of several other photoreceptor genes. We applied large-scale molecular approaches to identify altered transcripts and proteins in tomato wild-type (WT) versus a CRY2 overexpressing transgenic genotype, under a diurnal rhythm. Our results showed that tomato CRY2 profoundly affects both gene and protein expression in response to daily light cycle. Particularly altered molecular pathways are related to biotic/abiotic stress, photosynthesis, including components of the light and dark reactions and of starch and sucrose biosynthesis, as well as to secondary metabolism, such as phenylpropanoid, phenolic and flavonoid/anthocyanin biosynthesis pathways. One of the most interesting results is the coordinated up-regulation, in the transgenic genotype, of a consistent number of transcripts and proteins involved in photorespiration and photosynthesis. It is conceivable that light modulates the energetic metabolism of tomato through a fine CRY2-mediated transcriptional control.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

Expression dynamics of metabolic and regulatory components across stages of panicle and seed development in indica rice

Carefully analyzed expression profiles can serve as a valuable reference for deciphering gene functions. We exploited the potential of whole genome microarrays to measure the spatial and temporal expression profiles of rice genes in 19 stages of vegetative and reproductive development. We could verify expression of 22,980 genes in at least one of the tissues. Differential expression analysis with respect to five vegetative tissues and preceding stages of development revealed reproductive stage-preferential/-specific genes. By using subtractive logic, we identified 354 and 456 genes expressing specifically during panicle and seed development, respectively. The metabolic/hormonal pathways and transcription factor families playing key role in reproductive development were elucidated after overlaying the expression data on the public databases and manually curated list of transcription factors, respectively. During floral meristem differentiation (P1) and male meiosis (P3), the genes involved in jasmonic acid and phenylpropanoid biosynthesis were significantly upregulated. P6 stage of panicle, containing mature gametophytes, exhibited enrichment of transcripts involved in homogalacturonon degradation. Genes regulating auxin biosynthesis were induced during early seed development. We validated the stage-specificity of regulatory regions of three panicle-specific genes, OsAGO3, OsSub42, and RTS, and an early seed-specific gene, XYH, in transgenic rice. The data generated here provides a snapshot of the underlying complexity of the gene networks regulating rice reproductive development.

Transcriptome analysis of rice mature root tissue and root tips in early development by massive parallel sequencing

Despite the major physiological dissimilarities between mature root regions and their tips, differences in their gene expression profiles remain largely unexplored. In this research, the transcriptome of rice (Oryza sativa L. subsp. japonica) mature root tissue and root tips was monitored using mRNA-Seq at two time points. Almost 50 million 76 bp reads were mapped onto the rice genome sequence, expression patterns for different tissues and time points were investigated, and at least 1106 novel transcriptionally active regions (nTARs) expressed in rice root tissue were detected. More than 30 000 genes were found to be expressed in rice roots, among which were 1761 root-enriched and 306 tip-enriched transcripts. Mature root tissue appears to respond more strongly to external stimuli than tips, showing a higher expression of, for instance, auxin-responsive and abscisic acid-responsive genes, as well as the phenylpropanoid pathway and photosynthesis upon light. The root tip-enriched transcripts are mainly involved in mitochondrial electron transport, organelle development, secondary metabolism, DNA replication and metabolism, translation, and cellular component organization. During root maturation, genes involved in cell wall biosynthesis and modification, response to oxidative stress, and secondary metabolism were activated. For some nTARs, a potential role in root development can be put forward based on homology to genes involved in CLAVATA signalling, cell cycle regulators, and hormone signalling. A subset of differentially expressed genes and novel transcripts was confirmed using (quantitative) reverse transcription-PCR. These results uncover previously unrecognized tissue-specific expression profiles and provide an interesting starting point to study the different regulation of transcribed regions of these tissues.

Rice functional genomics research: progress and implications for crop genetic improvement

Rice is a staple food crop and has become a reference of monocot plant for functional genomic research. With the availability of high quality rice genome sequence, there has been rapid accumulation of functional genomic resources, including: large mutant libraries by T-DNA insertion, transposon tagging, and chemical mutagenesis; global expression profiles of the genes in the entire life cycle of rice growth and development; full-length cDNAs for both indica and japonica rice; sequences from resequencing large numbers of diverse germplasm accessions. Such resource development has greatly accelerated gene cloning. By the end of 2010, over 600 genes had been cloned using various methods. Many of the genes control agriculturally useful traits such as yield, grain quality, resistances to biotic and abiotic stresses, and nutrient-use efficiency, thus have potential utility in crop genetic improvement. This review was aimed to provide a comprehensive summary of such progress. We also presented our perspective for future studies.

Comparative network analysis reveals that tissue specificity and gene function are important factors influencing the mode of expression evolution in Arabidopsis and rice

Microarray experiments have yielded massive amounts of expression information measured under various conditions for the model species Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Expression compendia grouping multiple experiments make it possible to define correlated gene expression patterns within one species and to study how expression has evolved between species. We developed a robust framework to measure expression context conservation (ECC) and found, by analyzing 4,630 pairs of orthologous Arabidopsis and rice genes, that 77% showed conserved coexpression. Examples of nonconserved ECC categories suggested a link between regulatory evolution and environmental adaptations and included genes involved in signal transduction, response to different abiotic stresses, and hormone stimuli. To identify genomic features that influence expression evolution, we analyzed the relationship between ECC, tissue specificity, and protein evolution. Tissue-specific genes showed higher expression conservation compared with broadly expressed genes but were fast evolving at the protein level. No significant correlation was found between protein and expression evolution, implying that both modes of gene evolution are not strongly coupled in plants. By integration of cis-regulatory elements, many ECC conserved genes were significantly enriched for shared DNA motifs, hinting at the conservation of ancestral regulatory interactions in both model species. Surprisingly, for several tissue-specific genes, patterns of concerted network evolution were observed, unveiling conserved coexpression in the absence of conservation of tissue specificity. These findings demonstrate that orthologs inferred through sequence similarity in many cases do not share similar biological functions and highlight the importance of incorporating expression information when comparing genes across species.

Evolutionary transients in the rice transcriptome

In the canonical version of evolution by gene duplication, one copy is kept unaltered while the other is free to evolve. This process of evolutionary experimentation can persist for millions of years. Since it is so short lived in comparison to the lifetime of the core genes that make up the majority of most genomes, a substantial fraction of the genome and the transcriptome may—in principle—be attributable to what we will refer to as “evolutionary transients”, referring here to both the process and the genes that have gone or are undergoing this process. Using the rice gene set as a test case, we argue that this phenomenon goes a long way towards explaining why there are so many more rice genes than Arabidopsis genes, and why most excess rice genes show low similarity to eudicots.

Transcriptional characteristics of Xa21-mediated defense responses in rice

Bacterial blight, caused by Xanthomonas oryzae pv. oryzae (Xoo), is the most destructive bacterial disease of rice. The cloned rice gene Xa21 confers resistance to a broad spectrum of Xoo races. To identify genes involved in Xa21-mediated immunity, a whole-genome oligonucleotide microarray of rice was used to profile the expression of rice genes between incompatible interactions and mock treatments at 0, 4, 8, 24, 72 and 120 h post inoculation (hpi) or between incompatible and compatible interactions at 4 hpi, respectively. A total of 441 differentially expressed genes, designated as XDGs (Xa21 mediated differentially expressed genes), were identified. Based on their functional annotations, the XDGs were assigned to 14 categories, including defense-related, signaling, transcriptional regulators. Most of the defense-related genes belonged to the pathogenesis-related gene family, which was induced dramatically at 72 and 120 hpi. Interestingly, most signaling and transcriptional regulator genes were downregulated at 4 and 8 hpi, suggesting that negative regulation of cellular signaling may play a role in the Xa21-mediated defense response. Comparison of expression profiles between Xa21- and other R gene-mediated defense systems revealed interesting common responses. Representative XDGs with supporting evidences were also discussed.

Identification of rhizome-specific genes by genome-wide differential expression analysis in Oryza longistaminata

BACKGROUND

Rhizomatousness is a key component of perenniality of many grasses that contribute to competitiveness and invasiveness of many noxious grass weeds, but can potentially be used to develop perennial cereal crops for sustainable farmers in hilly areas of tropical Asia. Oryza longistaminata, a perennial wild rice with strong rhizomes, has been used as the model species for genetic and molecular dissection of rhizome development and in breeding efforts to transfer rhizome-related traits into annual rice species. In this study, an effort was taken to get insights into the genes and molecular mechanisms underlying the rhizomatous trait in O. longistaminata by comparative analysis of the genome-wide tissue-specific gene expression patterns of five different tissues of O. longistaminata using the Affymetrix GeneChip Rice Genome Array.

RESULTS

A total of 2,566 tissue-specific genes were identified in five different tissues of O. longistaminata, including 58 and 61 unique genes that were specifically expressed in the rhizome tips (RT) and internodes (RI), respectively. In addition, 162 genes were up-regulated and 261 genes were down-regulated in RT compared to the shoot tips. Six distinct cis-regulatory elements (CGACG, GCCGCC, GAGAC, AACGG, CATGCA, and TAAAG) were found to be significantly more abundant in the promoter regions of genes differentially expressed in RT than in the promoter regions of genes uniformly expressed in all other tissues. Many of the RT and/or RI specifically or differentially expressed genes were located in the QTL regions associated with rhizome expression, rhizome abundance and rhizome growth-related traits in O. longistaminata and thus are good candidate genes for these QTLs.

CONCLUSION

The initiation and development of the rhizomatous trait in O. longistaminata are controlled by very complex gene networks involving several plant hormones and regulatory genes, different members of gene families showing tissue specificity and their regulated pathways. Auxin/IAA appears to act as a negative regulator in rhizome development, while GA acts as the activator in rhizome development. Co-localization of the genes specifically expressed in rhizome tips and rhizome internodes with the QTLs for rhizome traits identified a large set of candidate genes for rhizome initiation and development in rice for further confirmation.

Common and distinct organ and stress responsive transcriptomic patterns in Oryza sativa and Arabidopsis thaliana

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

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

Elucidating the function of all predicted genes in rice remains as the ultimate goal in cereal genomics in order to ensure the development of improved varieties that will sustain an expanding world population. We constructed a gene expression database (RiceXPro, URL: http://ricexpro.dna.affrc.go.jp/) to provide an overview of the transcriptional changes throughout the growth of the rice plant in the field. RiceXPro contains two data sets corresponding to spatiotemporal gene expression profiles of various organs and tissues, and continuous gene expression profiles of leaf from transplanting to harvesting. A user-friendly web interface enables the extraction of specific gene expression profiles by keyword and chromosome search, and basic data analysis, thereby providing useful information as to the organ/tissue and developmental stage specificity of expression of a particular gene. Analysis tools such as t-test, calculation of fold change and degree of correlation facilitate the comparison of expression profiles between two random samples and the prediction of function of uncharacterized genes. As a repository of expression data encompassing growth in the field, this database can provide baseline information of genes that underlie various agronomically important traits in rice.

Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-seq

The functional complexity of the rice transcriptome is not yet fully elucidated, despite many studies having reported the use of DNA microarrays. Next-generation DNA sequencing technologies provide a powerful approach for mapping and quantifying the transcriptome, termed RNA sequencing (RNA-seq). In this study, we applied RNA-seq to globally sample transcripts of the cultivated rice Oryza sativa indica and japonica subspecies for resolving the whole-genome transcription profiles. We identified 15,708 novel transcriptional active regions (nTARs), of which 51.7% have no homolog to public protein data and >63% are putative single-exon transcripts, which are highly different from protein-coding genes (<20%). We found that approximately 48% of rice genes show alternative splicing patterns, a percentage considerably higher than previous estimations. On the basis of the available rice gene models, 83.1% (46,472 genes) of the current rice gene models were validated by RNA-seq, and 6228 genes were identified to be extended at the 5' and/or 3' ends by at least 50 bp. Comparative transcriptome analysis demonstrated that 3464 genes exhibited differential expression patterns. The ratio of SNPs with nonsynonymous/synonymous mutations was nearly 1:1.06. In total, we interrogated and compared transcriptomes of the two rice subspecies to reveal the overall transcriptional landscape at maximal resolution.

Genome-scale analysis and comparison of gene expression profiles in developing and germinated pollen in Oryza sativa

Background

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.

Results

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.

Conclusions

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.

Gene expression profiles in rice roots under low phosphorus stress

Phosphorus (P), an important plant macronutrient, is a component of key molecules such as nucleic acids, phospholipids and ATP. P is often the limiting nutrient for crop yield potential because of the low concentration of soluble P that can be absorbed directly by plant. Plants have evolved a series of molecular and morphological adaptations to cope with P limitation. However, the molecular bases of these responses to P deficiency have not been thoroughly elucidated. In this report, the gene expression profiles of low-P-tolerant rice Zhongzao 18 (Oryza sativa ssp. Indica) and not-low-P-tolerant rice Lagrue (Oryza sativa ssp. Indica) roots at 6 h, 24 h and 72 h under low P stress were investigated and compared with a control (normal P conditions) profile, using a DNA chip of 60,000 oligos (70 mer) that represented all putative genes of the rice genome. A total of 1,518 and 2,358 genes exhibited alterations in expression in response to low P stress in at least one of the three time points in rice Zhongzao 18 and rice Lagrue, respectively. The differentially expressed genes included those involved in phosphate (Pi) transportation, transportations except for Pi transportation, phosphatase, enzymes other than phosphatase, primary metabolism, secondary metabolism and so on. Several genes involved in glycolysis and TCA cycle were up-regulated during the early stages of low P treatment in rice Zhongzao 18 roots, but not in rice Lagrue roots. The results may provide useful information to further studies of the molecular mechanism of plant adaptation to low P and thus facilitate research in improving P utilization in crop species.

A dynamic gene expression atlas covering the entire life cycle of rice

Growth and development of a plant are controlled by programmed expression of suits of genes at the appropriate time, tissue and abundance. Although genomic resources have been developed rapidly in recent years in rice, a model plant for cereal genome research, data of gene expression profiling are still insufficient to relate the developmental processes to transcriptomes, leaving a large gap between the genome sequence and phenotype. In this study, we generated genome-wide expression data by hybridizing 190 Affymetrix GeneChip Rice Genome Arrays with RNA from 39 tissues collected throughout the life cycle of the rice plant from two varieties, Zhenshan 97 and Minghui 63. Analyses of the global transcriptomes revealed many interesting features of dynamic patterns of gene expression across the tissues and stages. In total, 38 793 probe sets were detected as expressed and 69% of the expressed transcripts showed significantly variable expression levels among tissues/organs. We found that similarity of transcriptomes among organs corresponded well to their developmental relatedness. About 5.2% of the expressed transcripts showed tissue-specific expression in one or both varieties and 22.7% of the transcripts exhibited constitutive expression including 19 genes with high and stable expression in all the tissues. This dataset provided a versatile resource for plant genomic research, which can be used for associating the transcriptomes to the developmental processes, understanding the regulatory network of these processes, tracing the expression profile of individual genes and identifying reference genes for quantitative expression analyses.

AGRONOMICS1: a new resource for Arabidopsis transcriptome profiling

Transcriptome profiling has become a routine tool in biology. For Arabidopsis (Arabidopsis thaliana), the Affymetrix ATH1 expression array is most commonly used, but it lacks about one-third of all annotated genes present in the reference strain. An alternative are tiling arrays, but previous designs have not allowed the simultaneous analysis of both strands on a single array. We introduce AGRONOMICS1, a new Affymetrix Arabidopsis microarray that contains the complete paths of both genome strands, with on average one 25mer probe per 35-bp genome sequence window. In addition, the new AGRONOMICS1 array contains all perfect match probes from the original ATH1 array, allowing for seamless integration of the very large existing ATH1 knowledge base. The AGRONOMICS1 array can be used for diverse functional genomics applications such as reliable expression profiling of more than 30,000 genes, detection of alternative splicing, and chromatin immunoprecipitation coupled to microarrays (ChIP-chip). Here, we describe the design of the array and compare its performance with that of the ATH1 array. We find results from both microarrays to be of similar quality, but AGRONOMICS1 arrays yield robust expression information for many more genes, as expected. Analysis of the ATH1 probes on AGRONOMICS1 arrays produces results that closely mirror those of ATH1 arrays. Finally, the AGRONOMICS1 array is shown to be useful for ChIP-chip experiments. We show that heterochromatic H3K9me2 is strongly confined to the gene body of target genes in euchromatic chromosome regions, suggesting that spreading of heterochromatin is limited outside of pericentromeric regions.

Metabolomic screening applied to rice FOX Arabidopsis lines leads to the identification of a gene-changing nitrogen metabolism

Plant metabolomics developed as a powerful tool to examine gene functions and to gain deeper insight into the physiology of the plant cell. In this study, we screened Arabidopsis lines overexpressing rice full-length (FL) cDNAs (rice FOX Arabidopsis lines) using a gas chromatography-time-of-flight mass spectrometry (GC-TOF/MS)-based technique to identify rice genes that caused metabolic changes. This screening system allows fast and reliable identification of candidate lines showing altered metabolite profiles. We performed metabolomic and transcriptomic analysis of a rice FOX Arabidopsis line that harbored the FL cDNA of the rice ortholog of the Lateral Organ Boundaries (LOB) Domain (LBD)/Asymmetric Leaves2-like (ASL) gene of Arabidopsis, At-LBD37/ASL39. The investigated rice FOX Arabidopsis line showed prominent changes in the levels of metabolites related to nitrogen metabolism. The transcriptomic data as well as the results from the metabolite analysis of the Arabidopsis At-LBD37/ASL39-overexpressor plants were consistent with these findings. Furthermore, the metabolomic and transcriptomic analysis of the Os-LBD37/ASL39-overexpressing rice plants indicated that Os-LBD37/ASL39 is associated with processes related to nitrogen metabolism in rice. Thus, the combination of a metabolomics-based screening method and a gain-of-function approach is useful for rapid characterization of novel genes in both Arabidopsis and rice.

Phytochromes are the sole photoreceptors for perceiving red/far-red light in rice

Phytochromes are believed to be solely responsible for red and far-red light perception, but this has never been definitively tested. To directly address this hypothesis, a phytochrome triple mutant (phyAphyBphyC) was generated in rice (Oryza sativa L. cv. Nipponbare) and its responses to red and far-red light were monitored. Since rice only has three phytochrome genes (PHYA, PHYB and PHYC), this mutant is completely lacking any phytochrome. Rice seedlings grown in the dark develop long coleoptiles while undergoing regular circumnutation. The phytochrome triple mutants also show this characteristic skotomorphogenesis, even under continuous red or far-red light. The morphology of the triple mutant seedlings grown under red or far-red light appears completely the same as etiolated seedlings, and they show no expression of the light-induced genes. This is direct evidence demonstrating that phytochromes are the sole photoreceptors for perceiving red and far-red light, at least during rice seedling establishment. Furthermore, the shape of the triple mutant plants was dramatically altered. Most remarkably, triple mutants extend their internodes even during the vegetative growth stage, which is a time during which wild-type rice plants never elongate their internodes. The triple mutants also flowered very early under long day conditions and set very few seeds due to incomplete male sterility. These data indicate that phytochromes play an important role in maximizing photosynthetic abilities during the vegetative growth stage in rice.

Comparing genomic expression patterns across plant species reveals highly diverged transcriptional dynamics in response to salt stress

Background

Rice and barley are both members of Poaceae (grass family) but have a marked difference in salt tolerance. The molecular mechanism underlying this difference was previously unexplored. This study employs a comparative genomics approach to identify analogous and contrasting gene expression patterns between rice and barley.

Results

A hierarchical clustering approach identified several interesting expression trajectories among rice and barley genotypes. There were no major conserved expression patterns between the two species in response to salt stress. A wheat salt-stress dataset was queried for comparison with rice and barley. Roughly one-third of the salt-stress responses of barley were conserved with wheat while overlap between wheat and rice was minimal. These results demonstrate that, at transcriptome level, rice is strikingly different compared to the more closely related barley and wheat. This apparent lack of analogous transcriptional programs in response to salt stress is further highlighted through close examination of genes associated with root growth and development.

Conclusion

The analysis provides support for the hypothesis that conservation of transcriptional signatures in response to environmental cues depends on the genetic similarity among the genotypes within a species, and on the phylogenetic distance between the species.

Expansion mechanisms and functional annotations of hypothetical genes in the rice genome

In each completely sequenced genome, 30% to 50% of genes are annotated as uncharacterized hypothetical genes. In the rice (Oryza sativa) genome, 10,918 hypothetical genes were annotated in the latest version (release 6) of the Michigan State University rice genome annotation. We have implemented an integrative approach to analyze their duplication/expansion and function. The analyses show that tandem/segmental duplication and transposition/retrotransposition have significantly contributed to the expansion of hypothetical genes despite their different contribution rates. A total of 3,769 hypothetical genes have been detected from retrogene, tandem, segmental, Pack-MULE, or long terminated direct repeat-related duplication/expansion. The nonsynonymous substitutions per site and synonymous substitutions per site analyses showed that 21.65% of them were still functional, accounting for 7.47% of total hypothetical genes. Global expression analyses have identified 1,672 expressed hypothetical genes. Among them, 415 genes might function in a developmental stage-specific manner. Antisense strand expression and small RNA analyses have demonstrated that a high percentage of these hypothetical genes might play important roles in negatively regulating gene expression. Homologous searches against Arabidopsis (Arabidopsis thaliana), maize (Zea mays), sorghum (Sorghum bicolor), and indica rice genomes suggest that most of the hypothetical genes could be annotated from recently evolved genomic sequences. These data advance the understanding of rice hypothetical genes as being involved in lineage-specific expansion and that they function in a specific developmental stage. Our analyses also provide a valuable means to facilitate the characterization and functional annotation of hypothetical genes in other organisms.

Comparative transcriptomics in the Triticeae

Background

Barley and particularly wheat are two grass species of immense agricultural importance. In spite of polyploidization events within the latter, studies have shown that genotypically and phenotypically these species are very closely related and, indeed, fertile hybrids can be created by interbreeding. The advent of two genome-scale Affymetrix GeneChips now allows studies of the comparison of their transcriptomes.

Results

We have used the Wheat GeneChip to create a "gene expression atlas" for the wheat transcriptome (cv. Chinese Spring). For this, we chose mRNA from a range of tissues and developmental stages closely mirroring a comparable study carried out for barley (cv. Morex) using the Barley1 GeneChip. This, together with large-scale clustering of the probesets from the two GeneChips into "homologous groups", has allowed us to perform a genomic-scale comparative study of expression patterns in these two species. We explore the influence of the polyploidy of wheat on the results obtained with the Wheat GeneChip and quantify the correlation between conservation in gene sequence and gene expression in wheat and barley. In addition, we show how the conservation of expression patterns can be used to elucidate, probeset by probeset, the reliability of the Wheat GeneChip.

Conclusion

While there are many differences in expression on the level of individual genes and tissues, we demonstrate that the wheat and barley transcriptomes appear highly correlated. This finding is significant not only because given small evolutionary distance between the two species it is widely expected, but also because it demonstrates that it is possible to use the two GeneChips for comparative studies. This is the case even though their probeset composition reflects rather different design principles as well as, of course, the present incomplete knowledge of the gene content of the two species. We also show that, in general, the Wheat GeneChip is not able to distinguish contributions from individual homoeologs. Furthermore, the comparison between the two species leads us to conclude that the conservation of both gene sequence as well as gene expression is positively correlated with absolute expression levels, presumably reflecting increased selection pressure on genes coding for proteins present at high levels. In addition, the results indicate the presence of a correlation between sequence and expression conservation within the Triticeae.

Discovery of cis-elements between sorghum and rice using co-expression and evolutionary conservation

BACKGROUND

The spatiotemporal regulation of gene expression largely depends on the presence and absence of cis-regulatory sites in the promoter. In the economically highly important grass family, our knowledge of transcription factor binding sites and transcriptional networks is still very limited. With the completion of the sorghum genome and the available rice genome sequence, comparative promoter analyses now allow genome-scale detection of conserved cis-elements.

RESULTS

In this study, we identified thousands of phylogenetic footprints conserved between orthologous rice and sorghum upstream regions that are supported by co-expression information derived from three different rice expression data sets. In a complementary approach, cis-motifs were discovered by their highly conserved co-occurrence in syntenic promoter pairs. Sequence conservation and matches to known plant motifs support our findings. Expression similarities of gene pairs positively correlate with the number of motifs that are shared by gene pairs and corroborate the importance of similar promoter architectures for concerted regulation. This strongly suggests that these motifs function in the regulation of transcript levels in rice and, presumably also in sorghum.

CONCLUSION

Our work provides the first large-scale collection of cis-elements for rice and sorghum and can serve as a paradigm for cis-element analysis through comparative genomics in grasses in general.

A transcriptomic analysis of superhybrid rice LYP9 and its parents

By using a whole-genome oligonucleotide microarray, designed based on known and predicted indica rice genes, we investigated transcriptome profiles in developing leaves and panicles of superhybrid rice LYP9 and its parental cultivars 93-11 and PA64s. We detected 22,266 expressed genes out of 36,926 total genes set collectively from 7 tissues, including leaves at seedling and tillering stages, flag leaves at booting, heading, flowering, and filling stages, and panicles at filling stage. Clustering results showed that the F1 hybrid's expression profiles resembled those of its parental lines more than that which lies between the 2 parental lines. Out of the total gene set, 7,078 genes are shared by all sampled tissues and 3,926 genes (10.6% of the total gene set) are differentially expressed genes (DG). As we divided DG into those between the parents (DG(PP)) and between the hybrid and its parents (DG(HP)), the comparative results showed that genes in the categories of energy metabolism and transport are enriched in DG(HP) rather than in DG(PP). In addition, we correlated the concurrence of DG and yield-related quantitative trait loci, providing a potential group of heterosis-related genes.

The dynamics of photosynthesis

Despite recent elucidation of the three-dimensional structure of major photosynthetic complexes, our understanding of light energy conversion in plant chloroplasts and microalgae under physiological conditions requires exploring the dynamics of photosynthesis. The photosynthetic apparatus is a flexible molecular machine that can acclimate to metabolic and light fluctuations in a matter of seconds and minutes. On a longer time scale, changes in environmental cues trigger acclimation responses that elicit intracellular signaling between the nucleo-cytosol and chloroplast resulting in modification of the biogenesis of the photosynthetic machinery. Here we attempt to integrate well-established knowledge on the functional flexibility of light-harvesting and electron transfer processes, which has greatly benefited from genetic approaches, with data derived from the wealth of recent transcriptomic and proteomic studies of acclimation responses in photosynthetic eukaroytes.

A transcriptome atlas of rice cell types uncovers cellular, functional and developmental hierarchies

The functions of the plant body rely on interactions among distinct and nonequivalent cell types. The comparison of transcriptomes from different cell types should expose the transcriptional networks that underlie cellular attributes and contributions. Using laser microdissection and microarray profiling, we have produced a cell type transcriptome atlas that includes 40 cell types from rice (Oryza sativa) shoot, root and germinating seed at several developmental stages, providing patterns of cell specificity for individual genes and gene classes. Cell type comparisons uncovered previously unrecognized properties, including cell-specific promoter motifs and coexpressed cognate binding factor candidates, interaction partner candidates and hormone response centers. We inferred developmental regulatory hierarchies of gene expression in specific cell types by comparison of several stages within root, shoot and embryo.

[A place in tears]

Alternative splicing creates a diversity of gene products in higher eukaryotes. Twenty-five percent (1,583/6,371) of predicted alternatively spliced transcripts can be detected using the NSF45K rice whole-genome oligonucleotide array. We used the NSF45K array to assess differential expression patterns of 507 loci showing at least a twofold change in expression between light- and dark-grown seedlings. At least 42% of these loci show evidence of alternative splicing in aerial seedling tissue of Oryza sativa ssp. japonica cv. Nipponbare. Most alternative splice forms display the same pattern of regulation as the primary, or most highly expressed, transcript; however, splice forms for ten loci, represented by 35 oligos, display opposite expression patterns in the light vs. dark. We found similar evidence of alternative splicing events in Affymetrix microarray data for Nipponbare rice treated with the causative agent of fungal rice blast, Magnaporthe grisea. This strategy for analyzing alternative splicing in microarray data will enable delineation of the diversity of splicing in rice.

Characterization of the genome expression trends in the heading-stage panicle of six rice lineages

To study how changes in gene regulation shape phenotypic variations in rice, we performed a comparative analysis of genome expression in the heading-stage panicle from six lineages of cultivated and wild rice, including Oryza sativa subsp. indica, japonica and javanica, O. nivara , O. rufipogon and O. glaberrima. While nearly three-quarters of the genes are expressed at a constant level in all six lineages, a large portion of the genome, ranging from 1767 to 4489 genes, exhibited differential expression between Asian domesticated and wild rice with repression or down-regulation of genome expression in Asian cultivated rice as the dominant trend. Importantly, we found this repression was achieved to a large extent by the differential expression of a single member of paralogous gene families. Functional analysis of the differentially expressed genes revealed that genes related to catabolism are repressed while genes related to anabolism up-regulated. Finally, we observed that distinct evolutionary forces may have acted on gene expression and the coding sequences in the examined rice lineages.

Refinement of light-responsive transcript lists using rice oligonucleotide arrays: evaluation of gene-redundancy

Studies of gene function are often hampered by gene-redundancy, especially in organisms with large genomes such as rice (Oryza sativa). We present an approach for using transcriptomics data to focus functional studies and address redundancy. To this end, we have constructed and validated an inexpensive and publicly available rice oligonucleotide near-whole genome array, called the rice NSF45K array. We generated expression profiles for light- vs. dark-grown rice leaf tissue and validated the biological significance of the data by analyzing sources of variation and confirming expression trends with reverse transcription polymerase chain reaction. We examined trends in the data by evaluating enrichment of gene ontology terms at multiple false discovery rate thresholds. To compare data generated with the NSF45K array with published results, we developed publicly available, web-based tools (www.ricearray.org). The Oligo and EST Anatomy Viewer enables visualization of EST-based expression profiling data for all genes on the array. The Rice Multi-platform Microarray Search Tool facilitates comparison of gene expression profiles across multiple rice microarray platforms. Finally, we incorporated gene expression and biochemical pathway data to reduce the number of candidate gene products putatively participating in the eight steps of the photorespiration pathway from 52 to 10, based on expression levels of putatively functionally redundant genes. We confirmed the efficacy of this method to cope with redundancy by correctly predicting participation in photorespiration of a gene with five paralogs. Applying these methods will accelerate rice functional genomics.

A genome-wide transcription analysis reveals a close correlation of promoter INDEL polymorphism and heterotic gene expression in rice hybrids

Heterosis, or hybrid vigor, refers to the phenomenon in which hybrid progeny of two inbred varieties exhibits enhanced growth or agronomic performance. Although a century-long history of research has generated several hypotheses regarding the genetic basis of heterosis, the molecular mechanisms underlying heterosis and heterotic gene expression remain elusive. Here, we report a genome-wide gene expression analysis of two heterotic crosses in rice, taking advantage of its fully sequenced genomes. Approximately 7-9% of the genes were differentially expressed in the seedling shoots from two sets of heterotic crosses, including many transcription factor genes, and exhibited multiple modes of gene action. Comparison of the putative promoter regions of the ortholog genes between inbred parents revealed extensive sequence variation, particularly small insertions/deletions (INDELs), many of which result in the formation/disruption of putative cis-regulatory elements. Together, these results suggest that a combinatorial interplay between expression of transcription factors and polymorphic promoter cis-regulatory elements in the hybrids is one plausible molecular mechanism underlying heterotic gene action and thus heterosis in rice.

Identification and functional analysis of light-responsive unique genes and gene family members in rice

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.