A collection of 10,096 indica rice full-length cDNAs reveals highly expressed sequence divergence between Oryza sativa indica and japonica subspecies

Non-coding RNA expression analysis revealed the molecular mechanism of flag leaf heterosis in inter-subspecific hybrid rice

Heterosis has been used widespread in agriculture, but its molecular mechanism is inadequately understood. Plants have a large number of non-coding RNAs (ncRNAs), among them, functional ncRNAs that have been studied widely containing long non-coding RNA (lncRNA) and circular RNA (circRNA) that play a role in varied biological processes, as well as microRNA (miRNA), which can not only regulate the post-transcriptional expression of target genes, but also target lncRNA and circRNA then participate the competing endogenous RNA (ceRNA) regulatory network. However, the influence of these three ncRNAs and their regulatory relationships on heterosis is unknown in rice. In this study, the expression profile of ncRNAs and the ncRNA regulatory network related to heterosis were comprehensively analyzed in inter-subspecific hybrid rice. A total of 867 miRNAs, 3,278 lncRNAs and 2,521 circRNAs were identified in the hybrid and its parents. Analysis of the global profiles of these three types of ncRNAs indicated that significant differences existed in the distribution and sequence characteristics of the corresponding genes. The numbers of miRNA and lncRNA in hybrid were higher than those in its parents. A total of 784 ncRNAs (169 miRNAs, 573 lncRNAs and 42 circRNAs) showed differentially expressed in the hybrid, and their target/host genes were vital in stress tolerance, growth and development in rice. These discoveries suggested that the expression plasticity of ncRNA has an important role of inter-subspecific hybrid rice heterosis. It is worth mentioning that miRNAs exhibited substantially more variations between hybrid and parents compared with observed variation for lncRNA and circRNA. Non-additive expression ncRNAs and allele-specific expression genes-related ncRNAs in hybrid were provided in this study, and multiple sets of ncRNA regulatory networks closely related to heterosis were obtained. Meanwhile, heterosis-related regulatory networks of ceRNA (lncRNA and circRNA) and miRNA were also demonstrated.

Novel QTL Associated with Aerenchyma-Mediated Radial Oxygen Loss (ROL) in Rice (Oryza sativa L.) under Iron (II) Sulfide

In rice, high radial oxygen loss (ROL) has been associated with the reduction in the activity of methanogens, therefore reducing the formation of methane (CH₄) due to the abundance in application of nitrogen (N)-rich fertilizers. In this study, we evaluated the root growth behavior and ROL rate of a doubled haploid (DH) population (n = 117) and parental lines 93-11 (P1, indica) and Milyang352 (P2, japonica) in response to iron (II) sulfide (FeS). In addition, we performed a linkage mapping and quantitative trait locus (QTL) analysis on the same population for the target traits. The results of the phenotypic evaluation revealed that parental lines had distinctive root growth and ROL patterns, with 93-11 (indica) and Milyang352 (japonica) showing low and high ROL rates, respectively. This was also reflected in their derived population, indicating that 93.2% of the DH lines exhibited a high ROL rate and about 6.8% had a low ROL pattern. Furthermore, the QTL and linkage map analysis detected two QTLs associated with the control of ROL and root area on chromosomes 2 (qROL-2-1, 127 cM, logarithm of the odds (LOD) 3.04, phenotypic variation explained (PVE) 11.61%) and 8 (qRA-8-1, 97 cM, LOD 4.394, PVE 15.95%), respectively. The positive additive effect (2.532) of qROL-2-1 indicates that the allele from 93-11 contributed to the observed phenotypic variation for ROL. The breakthrough is that the qROL-2-1 harbors genes proposed to be involved in stress signaling, defense response mechanisms, and transcriptional regulation, among others. The qPCR results revealed that the majority of genes harbored by the qROL-2-1 recorded a higher transcript accumulation level in Milyang352 over time compared to 93-11. Another set of genes exhibited a high transcript abundance in P1 compared to P2, while a few were differentially regulated between both parents. Therefore, OsTCP7 and OsMYB21, OsARF8 genes encoding transcription factors (TFs), coupled with OsTRX, OsWBC8, and OsLRR2 are suggested to play important roles in the positive regulation of ROL in rice. However, the recorded differential expression of OsDEF7 and OsEXPA, and the decrease in OsNIP2, Oscb5, and OsPLIM2a TF expression between parental lines proposes them as being involved in the control of oxygen flux level in rice roots.

Early Domestication History of Asian Rice Revealed by Mutations and Genome-Wide Analysis of Gene Genealogies

BACKGROUND

Asian rice (Oryza sativa L.) has been a model plant but its cultivation history is inadequately understood, and its origin still under debate. Several enigmas remain, including how this annual crop shifted its growth habit from its perennial ancestor, O. rufipogon, why genetic divergence between indica and japonica appears older than the history of human domestication, and why some domestication genes do not show signals of introgression between subgroups. Addressing these issues may benefit both basic research and rice breeding.

RESULTS

Gene genealogy-based mutation (GGM) analysis shows that history of Asian rice is divided into two phases (Phase I and II) of about equal lengths. Mutations occurred earlier than the partition of indica and japonica to Os genome mark Phase-I period. We diagnosed 91 such mutations among 101 genes sampled across 12 chromosomes of Asian rice and its wild relatives. Positive selection, detected more at 5' regions than at coding regions of some of the genes, involved 22 loci (e.g., An-1, SH4, Rc, Hd3a, GL3.2, OsMYB3, OsDFR, and OsMYB15), which affected traits from easy harvesting, grain color, flowering time, productivity, to likely taste and tolerance. Phase-I mutations of OsMYB3, OsHd3a and OsDFR were experimentally tested and all caused enhanced functions of the genes in vivo. Phase-II period features separate cultivations, lineage-specific selection, and expanded domestication to more genes. Further genomic analysis, along with phenotypic comparisons, indicates that O. sativa is hybrid progeny of O. rufipogon and O. nivara, inherited slightly more genes of O. rufipogon. Congruently, modern alleles of the sampled genes are approximately 6% ancient, 38% uni-specific, 40% bi-specific (mixed), and 15% new after accumulating significant mutations. Results of sequencing surveys across modern cultivars/landraces indicate locus-specific usages of various alleles while confirming the associated mutations.

CONCLUSIONS

Asian rice was initially domesticated as one crop and later separate selection mediated by human resulted in its major subgroups. This history and the hybrid origin well explain previous puzzles. Positive selection, particularly in 5' regions, was the major force underlying trait domestication. Locus-specific domestication can be characterized and the result may facilitate breeders in developing better rice varieties in future.

Chilling tolerance in rice: Past and present

Rice is generally sensitive to chilling stress, which seriously affects growth and yield. Since early in the last century, considerable efforts have been made to understand the physiological and molecular mechanisms underlying the response to chilling stress and improve rice chilling tolerance. Here, we review the research trends and advances in this field. The phenotypic and biochemical changes caused by cold stress and the physiological explanations are briefly summarized. Using published data from the past 20 years, we reviewed the past progress and important techniques in the identification of quantitative trait loci (QTL), novel genes, and cellular pathways involved in rice chilling tolerance. The advent of novel technologies has significantly advanced studies of cold tolerance, and the characterization of QTLs, key genes, and molecular modules have sped up molecular design breeding for cold tolerance in rice varieties. In addition to gene function studies based on overexpression or artificially generated mutants, elucidating natural allelic variation in specific backgrounds is emerging as a novel approach for the study of cold tolerance in rice, and the superior alleles identified using this approach can directly facilitate breeding.

A New Way of Rice Breeding: Polyploid Rice Breeding

Polyploid rice, first discovered by Japanese scientist Eiiti Nakamori in 1933, has a history of nearly 90 years. In the following years, polyploid rice studies have mainly focused on innovations in breeding theory, induction technology and the creation of new germplasm, the analysis of agronomic traits and nutritional components, the study of gametophyte development and reproduction characteristics, DNA methylation modification and gene expression regulation, distant hybridization and utilization among subspecies, species and genomes. In recent years, PMeS lines and neo-tetraploid rice lines with stable high seed setting rate characteristics have been successively selected, breaking through the bottleneck of low seed setting rate of polyploid rice. Following, a series of theoretical and applied studies on high seed setting rate tetraploid rice were carried out. This has pushed research on polyploid rice to a new stage, opening new prospects for polyploid rice breeding.

Rapid and pervasive development- and tissue-specific homeolog expression partitioning in newly formed inter-subspecific rice segmental allotetraploids

Background

In diverse plant taxa, whole-genome duplication (WGD) events are major sources of phenotypic novelty. Studies of gene expression in synthetic polyploids have shown immediate expression and functional partitioning of duplicated genes among different tissues. Many studies of the tissue-specific homeolog expression partitioning have focused on allopolyploids that have very different parental genomes, while few studies have focused on autopolyploids or allopolyploids that have similar parental genomes.

Results

In this study, we used a set of reciprocal F1 hybrids and synthetic tetraploids constructed from subspecies (japonica and indica) of Asian rice (Oryza sativa L.) as a model to gain insights into the expression partitioning of homeologs among tissues in a developmental context. We assayed the tissue-specific silencing (TSS) of the parental homeologs of 30 key genes in the hybrids and tetraploids relative to the in vitro “hybrids” (parental mixes) using Sequenom MassARRAY. We found that the parental mix and synthetic tetraploids had higher frequencies of homeolog TSS than the F1, revealing an instantaneous role of WGD on homeolog expression partitioning.

Conclusions

Our observations contradicted those of previous studies in which newly formed allopolyploids had a low TSS frequency, similar to that of F1 hybrids, suggesting that the impact of WGD on homeolog expression requires a longer time to manifest. In addition, we found that the TSS frequency in the tetraploids varied at different growth stages and that roots had a much higher frequency of TSS than leaves, which indicated that developmental and metabolic traits may influence the expression states of duplicated genes in newly formed plant polyploids.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5150-7) contains supplementary material, which is available to authorized users.

Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice

The rich genetic diversity in Oryza sativa and Oryza rufipogon serves as the main sources in rice breeding. Large-scale resequencing has been undertaken to discover allelic variants in rice, but much of the information for genetic variation is often lost by direct mapping of short sequence reads onto the O. sativa japonica Nipponbare reference genome. Here we constructed a pan-genome dataset of the O. sativa-O. rufipogon species complex through deep sequencing and de novo assembly of 66 divergent accessions. Intergenomic comparisons identified 23 million sequence variants in the rice genome. This catalog of sequence variations includes many known quantitative trait nucleotides and will be helpful in pinpointing new causal variants that underlie complex traits. In particular, we systemically investigated the whole set of coding genes using this pan-genome data, which revealed extensive presence and absence of variation among rice accessions. This pan-genome resource will further promote evolutionary and functional studies in rice.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica

BACKGROUND

The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate.

RESULTS

Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the "pan-genome" of three divergent rice varieties and document several megabases of each genome absent in the other two.

CONCLUSIONS

Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

Analyses of Ca2+ dynamics using a ubiquitin-10 promoter-driven Yellow Cameleon 3.6 indicator reveal reliable transgene expression and differences in cytoplasmic Ca2+ responses in Arabidopsis and rice (Oryza sativa) roots

Ca(2+) signatures are central to developmental processes and adaptive responses in plants. However, high-resolution studies of Ca(2+) dynamics using genetically encoded Ca(2+) indicators (GECIs) such as Yellow Cameleon (YC) proteins have so far not been conducted in important model crops such as rice (Oryza sativa). We conducted a comparative study of 35S and ubiquitin-10 (UBQ10) promoter functionality in Arabidopsis thaliana and O. sativa plants expressing the Ca(2+) indicator Yellow Cameleon 3.6 (YC3.6) under control of the UBQ10 or 35S promoter. Ca(2+) signatures in roots of both species were analyzed during exposure to hyperpolarization/depolarization cycles or in response to application of the amino acid glutamate. We found a superior performance of the UBQ10 promoter with regard to expression pattern, levels and expression stabilities in both species. We observed remarkable differences between the two species in the spatiotemporal parameters of the observed Ca(2+) signatures. Rice appeared in general to respond with a lower maximal signal amplitude but greatly increased signal duration when compared with Arabidopsis. Our results identify important advantages to using the UBQ10 promoter in Arabidopsis and rice and in T-DNA mutant backgrounds. Moreover, the observed differences in Ca(2+) signaling in the two species underscore the need for comparative studies to achieve a comprehensive understanding of Ca(2+) signaling in plants.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica

Background

The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate.

Results

Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the ‘pan-genome’ of three divergent rice varieties and document several megabases of each genome absent in the other two.

Conclusions

Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

Electronic supplementary material

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

Identification and molecular characterization of tissue-preferred rice genes and their upstream regularly sequences on a genome-wide level

BACKGROUND

Gene upstream regularly sequences (URSs) can be used as one of the tools to annotate the biological functions of corresponding genes. In addition, tissue-preferred URSs are frequently used to drive the transgene expression exclusively in targeted tissues during plant transgenesis. Although many rice URSs have been molecularly characterized, it is still necessary and valuable to identify URSs that will benefit plant transformation and aid in analyzing gene function.

RESULTS

In this study, we identified and characterized root-, seed-, leaf-, and panicle-preferred genes on a genome-wide level in rice. Subsequently, their expression patterns were confirmed through quantitative real-time RT-PCR (qRT-PCR) by randomly selecting 9candidate tissue-preferred genes. In addition, 5 tissue-preferred URSs were characterized by investigating the URS::GUS transgenic plants. Of these URS::GUS analyses, the transgenic plants harboring LOC_Os03g11350 URS::GUS construct showed the GUS activity only in young pollen. In contrast, when LOC_Os10g22450 URS was used to drive the reporter GUS gene, the GUS activity was detected only in mature pollen. Interestingly, the LOC_Os10g34360 URS was found to be vascular bundle preferred and its activities were restricted only to vascular bundles of leaves, roots and florets. In addition, we have also identified two URSs from genes LOC_Os02G15090 and LOC_Os06g31070 expressed in a seed-preferred manner showing the highest expression levels of GUS activities in mature seeds.

CONCLUSION

By genome-wide analysis, we have identified tissue-preferred URSs, five of which were further characterized using transgenic plants harboring URS::GUS constructs. These data might provide some evidence for possible functions of the genes and be a valuable resource for tissue-preferred candidate URSs for plant transgenesis.

Genome-wide disruption of gene expression in allopolyploids but not hybrids of rice subspecies

Hybridization and polyploidization are prominent processes in plant evolution. Hybrids and allopolyploids typically exhibit radically altered gene expression patterns relative to their parents, a phenomenon termed “transcriptomic shock.” To distinguish the effects of hybridization from polyploidization on coregulation of divergent alleles, we analyzed expression of parental copies (homoeologs) of 11,608 genes using RNA-seq-based transcriptome profiling in reciprocal hybrids and tetraploids constructed from subspecies japonica and indica of Asian rice (Oryza sativa L.). The diploid hybrids and their derived allopolyploids differ dramatically in morphology, despite having the same suite of genes and genic proportions. Allelic and homoeolog-specific transcripts were unequivocally diagnosed in the hybrids and tetraploids based on parent-specific SNPs. Compared with the in silico hybrid (parental mix), the range of progenitor expression divergence was significantly reduced in both reciprocally generated F1 hybrids, presumably due to the ameliorating effects of a common trans environment on divergent cis-factors. In contrast, parental expression differences were greatly elaborated at the polyploid level, which we propose is a consequence of stoichiometric disruptions associated with the numerous chromosomal packaging and volumetric changes accompanying nascent polyploidy. We speculate that the emergent property of “whole genome doubling” has repercussions that reverberate throughout the transcriptome and downstream, ultimately generating altered phenotypes. This perspective may yield insight into the nature of adaptation and the origin of evolutionary novelty accompanying polyploidy.

The rice OsLTP6 gene promoter directs anther-specific expression by a combination of positive and negative regulatory elements

Characterization of tissue-specific plant gene promoters will benefit genetic improvement in crops. Here, we isolated a novel rice anther-specific plant lipid transfer protein (OsLTP6) gene through high through-put expressional profiling. The promoter of OsLTP6 was introduced to the upstream of the uidA gene, which encodes β-glucuronidase (GUS), and transformed into rice plants for functional analysis. Histochemical and fluorometric GUS assay showed that GUS was specifically expressed in the anthers and pollens in OsLTP6 promoter::uidA transgenic plants. Transverse section of the rice anther further indicated that the OsLTP6 promoter directed the reporter gene specifically expressed in anther tapetum. To identify regulatory elements within OsLTP6 promoter region, four progressive deletions of the OsLTP6 promoter were constructed. The results indicated that the OsLTP6 promoter achieved anther-specific expression through a combination of positive and negative regulatory elements. A 26-bp motif upstream of TATA box was a key transcriptional activator for OsLTP6 gene. CAAT box and GTGA box were the putative motifs to increase the transcription level to full expression. Two negative regulatory elements were also found in two distinct regions within this promoter. They repressed the expression in leaf and stem, respectively. These results revealed the regulating complexity of anther-specific expression.

A BAC physical map of aus rice cultivar 'Kasalath', and the map-based genomic sequence of 'Kasalath' chromosome 1

Comparative analysis using available genomic resources within closely related species is an effective way to investigate genomic sequence and structural diversity. Rice (Oryza sativa L.) has undergone significant physiological and morphological changes during its domestication and local adaptation. We present a complete bacterial artificial chromosome (BAC) physical map for the aus rice cultivar 'Kasalath', which covers 90% of the sequence of temperate japonica rice cultivar 'Nipponbare'. Examination of physical distances between computational and experimental measurements of 'Kasalath' BAC insert size revealed the presence of more than 500 genomic regions that appear to have significant chromosome structural changes between the two cultivars. In particular, a genomic region on the long arm of 'Kasalath' chromosome 11 carrying a disease-resistance gene cluster was greatly expanded relative to the 'Nipponbare' genome. We also decoded 41.37 Mb of high-quality genomic sequence from 'Kasalath' chromosome 1. Extensive comparisons of chromosome 1 between 'Kasalath' and 'Nipponbare' led to the discovery of 317,843 single-nucleotide polymorphisms (SNPs) and 66,331 insertion/deletion (indel) sites. Nearly two-thirds of the expressed genes on rice chromosome 1 carried natural variations involving SNPs and/or indels that resulted in substitutions, insertions or deletions of amino acids in one cultivar relative to the other. We also observed gain and loss of genes caused by large indels. This study provides an important framework and an invaluable dataset for further understanding of the molecular mechanisms underlying the evolution and functions of the rice genome.

Genomic resources for functional analyses of the rice genome

With the availability of the rice genome sequence, rice research communities are entering a new era of plant functional genomics. The last decade has seen rapid worldwide progress on establishing platforms for rice functional genomic research. These platforms offer practical toolkits and genomic resources for high-throughput identification of genes and pathways. In this review, we summarize available genomic resources for functional analyses of the rice genome. These genomic resources include high-quality bacterial artificial chromosome libraries, large-scale expression sequence tags, full-length cDNA collections, large amounts of data on global expression profiles, various mutant libraries and integrated bioinformatics databases. We not only present the current status of genomic resources but also discuss their usage in elucidating gene functions of the rice genome.

Mutant resources for the functional analysis of the rice genome

Rice is one of the most important crops worldwide, both as a staple food and as a model system for genomic research. In order to systematically assign functions to all predicted genes in the rice genome, a large number of rice mutant lines, including those created by T-DNA insertion, Ds/dSpm tagging, Tos17 tagging, and chemical/irradiation mutagenesis, have been generated by groups around the world. In this study, we have reviewed the current status of mutant resources for functional analysis of the rice genome. A total of 246 566 flanking sequence tags from rice mutant libraries with T-DNA, Ds/dSpm, or Tos17 insertion have been collected and analyzed. The results show that, among 211 470 unique hits, inserts located in the genic region account for 68.16%, and 60.49% of nuclear genes contain at least one insertion. Currently, 57% of non-transposable-element-related genes in rice have insertional tags. In addition, chemical/irradiation-induced rice mutant libraries have contributed a lot to both gene identification and new technology for the identification of mutant sites. In this review, we summarize how these tools have been used to generate a large collection of mutants. In addition, we discuss the merits of classic mutation strategies. In order to achieve saturation of mutagenesis in rice, DNA targeting, and new resources like RiceFox for gene functional identification are reviewed from a perspective of the future generation of rice mutant resources.

Visual analysis of transcriptome data in the context of anatomical structures and biological networks

The complexity and temporal as well as spatial resolution of transcriptome datasets is constantly increasing due to extensive technological developments. Here we present methods for advanced visualization and intuitive exploration of transcriptomics data as necessary prerequisites in order to facilitate the gain of biological knowledge. Color-coding of structural images based on the expression level enables a fast visual data analysis in the background of the examined biological system. The network-based exploration of these visualizations allows for comparative analysis of genes with specific transcript patterns and supports the extraction of functional relationships even from large datasets. In order to illustrate the presented methods, the tool HIVE was applied for visualization and exploration of database-retrieved expression data for master regulators of Arabidopsis thaliana flower and seed development in the context of corresponding tissue-specific regulatory networks.

Advances in omics and bioinformatics tools for systems analyses of plant functions

Omics and bioinformatics are essential to understanding the molecular systems that underlie various plant functions. Recent game-changing sequencing technologies have revitalized sequencing approaches in genomics and have produced opportunities for various emerging analytical applications. Driven by technological advances, several new omics layers such as the interactome, epigenome and hormonome have emerged. Furthermore, in several plant species, the development of omics resources has progressed to address particular biological properties of individual species. Integration of knowledge from omics-based research is an emerging issue as researchers seek to identify significance, gain biological insights and promote translational research. From these perspectives, we provide this review of the emerging aspects of plant systems research based on omics and bioinformatics analyses together with their associated resources and technological advances.

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.

Assigning biological functions to rice genes by genome annotation, expression analysis and mutagenesis

Rice is the first cereal genome to be completely sequenced. Since the completion of its genome sequencing, considerable progress has been made in multiple areas including the whole genome annotation, gene expression profiling, mutant collection, etc. Here, we summarize the current status of rice genome annotation and review the methodology of assigning biological functions to hundreds of thousands of rice genes as well as discuss the major limitations and the future perspective in rice functional genomics. Available data analysis shows that the rice genome encodes around 32,000 protein-coding genes. Expression analysis revealed at least 31,000 genes with expression evidence from full-length cDNA/EST collection or other transcript profiling. In addition, we have summarized various strategies to generate mutant population including natural, physical, chemical, T-DNA, transposon/retrotransposon or gene silencing based mutagenesis. Currently, more than 1 million of mutants have been generated and 27,551 of them have their flanking sequence tags. To assign biological functions to hundreds of thousands of rice genes, global co-operations are required, various genetic resources should be more easily accessible and diverse data from transcriptomics, proteomics, epigenetics, comparative genomics and bioinformatics should be integrated to better understand the functions of these genes and their regulatory mechanisms.

Efficient plant gene identification based on interspecies mapping of full-length cDNAs

We present an annotation pipeline that accurately predicts exon–intron structures and protein-coding sequences (CDSs) on the basis of full-length cDNAs (FLcDNAs). This annotation pipeline was used to identify genes in 10 plant genomes. In particular, we show that interspecies mapping of FLcDNAs to genomes is of great value in fully utilizing FLcDNA resources whose availability is limited to several species. Because low sequence conservation at 5′- and 3′-ends of FLcDNAs between different species tends to result in truncated CDSs, we developed an improved algorithm to identify complete CDSs by the extension of both ends of truncated CDSs. Interspecies mapping of 71 801 monocot FLcDNAs to the Oryza sativa genome led to the detection of 22 142 protein-coding regions. Moreover, in comparing two mapping programs and three ab initio prediction programs, we found that our pipeline was more capable of identifying complete CDSs. As demonstrated by monocot interspecies mapping, in which nucleotide identity between FLcDNAs and the genome was ∼80%, the resultant inferred CDSs were sufficiently accurate. Finally, we applied both inter- and intraspecies mapping to 10 monocot and dicot genomes and identified genes in 210 551 loci. Interspecies mapping of FLcDNAs is expected to effectively predict genes and CDSs in newly sequenced genomes.

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-wide characterization of the biggest grass, bamboo, based on 10,608 putative full-length cDNA sequences

BACKGROUND

With the availability of rice and sorghum genome sequences and ongoing efforts to sequence genomes of other cereal and energy crops, the grass family (Poaceae) has become a model system for comparative genomics and for better understanding gene and genome evolution that underlies phenotypic and ecological divergence of plants. While the genomic resources have accumulated rapidly for almost all major lineages of grasses, bamboo remains the only large subfamily of Poaceae with little genomic information available in databases, which seriously hampers our ability to take a full advantage of the wealth of grass genomic data for effective comparative studies.

RESULTS

Here we report the cloning and sequencing of 10,608 putative full length cDNAs (FL-cDNAs) primarily from Moso bamboo, Phyllostachys heterocycla cv. pubescens, a large woody bamboo with the highest ecological and economic values of all bamboos. This represents the third largest FL-cDNA collection to date of all plant species, and provides the first insight into the gene and genome structures of bamboos. We developed a Moso bamboo genomic resource database that so far contained the sequences of 10,608 putative FL-cDNAs and nearly 38,000 expressed sequence tags (ESTs) generated in this study.

CONCLUSION

Analysis of FL-cDNA sequences show that bamboo diverged from its close relatives such as rice, wheat, and barley through an adaptive radiation. A comparative analysis of the lignin biosynthesis pathway between bamboo and rice suggested that genes encoding caffeoyl-CoA O-methyltransferase may serve as targets for genetic manipulation of lignin content to reduce pollutants generated from bamboo pulping.

Natural and artificial mutants as valuable resources for functional genomics and molecular breeding

With the completion of rice genome sequencing, large collection of expression data and the great efforts in annotating rice genomes, the next challenge is to systematically assign functions to all predicted genes in the genome. The generations and collections of mutants at the genome-wide level form technological platform of functional genomics. In this study, we have reviewed currently employed tools to generate such mutant populations. These tools include natural, physical, chemical, tissue culture, T-DNA, transposon or gene silencing based mutagenesis. We also reviewed how these tools were used to generate a large collection of mutants and how these mutants can be screened and detected for functional analysis of a gene. The data suggested that the current population of mutants might be large enough to tag all predicted genes. However, the collection of flanking sequencing tags (FSTs) is limited due to the relatively higher cost. Thus, we have proposed a new strategy to generate gene-silencing mutants at the genome-wide level. Due to the large collection of insertion mutants, the next step to rice functional genomics should be focusing on functional characterization of tagged genes by detailed survey of corresponding mutants. Additionally, we also evaluated the utilization of these mutants as valuable resources for molecular breeding.

Genomics and bioinformatics resources for crop improvement

Recent remarkable innovations in platforms for omics-based research and application development provide crucial resources to promote research in model and applied plant species. A combinatorial approach using multiple omics platforms and integration of their outcomes is now an effective strategy for clarifying molecular systems integral to improving plant productivity. Furthermore, promotion of comparative genomics among model and applied plants allows us to grasp the biological properties of each species and to accelerate gene discovery and functional analyses of genes. Bioinformatics platforms and their associated databases are also essential for the effective design of approaches making the best use of genomic resources, including resource integration. We review recent advances in research platforms and resources in plant omics together with related databases and advances in technology.

Gene expression profiles deciphering rice phenotypic variation between Nipponbare (Japonica) and 93-11 (Indica) during oxidative stress

Rice is a very important food staple that feeds more than half the world's population. Two major Asian cultivated rice (Oryza sativa L.) subspecies, japonica and indica, show significant phenotypic variation in their stress responses. However, the molecular mechanisms underlying this phenotypic variation are still largely unknown. A common link among different stresses is that they produce an oxidative burst and result in an increase of reactive oxygen species (ROS). In this study, methyl viologen (MV) as a ROS agent was applied to investigate the rice oxidative stress response. We observed that 93-11 (indica) seedlings exhibited leaf senescence with severe lesions under MV treatment compared to Nipponbare (japonica). Whole-genome microarray experiments were conducted, and 1,062 probe sets were identified with gene expression level polymorphisms between the two rice cultivars in addition to differential expression under MV treatment, which were assigned as Core Intersectional Probesets (CIPs). These CIPs were analyzed by gene ontology (GO) and highlighted with enrichment GO terms related to toxin and oxidative stress responses as well as other responses. These GO term-enriched genes of the CIPs include glutathine S-transferases (GSTs), P450, plant defense genes, and secondary metabolism related genes such as chalcone synthase (CHS). Further insertion/deletion (InDel) and regulatory element analyses for these identified CIPs suggested that there may be some eQTL hotspots related to oxidative stress in the rice genome, such as GST genes encoded on chromosome 10. In addition, we identified a group of marker genes individuating the japonica and indica subspecies. In summary, we developed a new strategy combining biological experiments and data mining to study the possible molecular mechanism of phenotypic variation during oxidative stress between Nipponbare and 93-11. This study will aid in the analysis of the molecular basis of quantitative traits.

Overexpression of a NAC transcription factor enhances rice drought and salt tolerance

The plant-specific NAC (NAM, ATAF1/2, CUC2) transcription factors play diverse roles in plant development and stress responses. In this study, a rice NAC gene, ONAC045, was functionally characterized, especially with regard to its role in abiotic stress resistance. Expression analysis revealed that ONAC045 was induced by drought, high salt, and low temperature stresses, and abscisic acid (ABA) treatment in leaves and roots. Transcriptional activation assay in yeast indicated that ONAC045 functioned as a transcriptional activator. Transient expression of GFP-ONAC045 in onion epidermal cells revealed that ONAC045 protein was localized in the nucleus. Transgenic rice plants overexpressing ONAC045 showed enhanced tolerance to drought and salt treatments. Two stress-responsive genes were upregulated in transgenic rice. Together, these results suggest that ONAC045 encodes a novel stress-responsive NAC transcription factor and is potential useful for engineering drought and salt tolerant rice.

RICD: a rice indica cDNA database resource for rice functional genomics

BACKGROUND

The Oryza sativa L. indica subspecies is the most widely cultivated rice. During the last few years, we have collected over 20,000 putative full-length cDNAs and over 40,000 ESTs isolated from various cDNA libraries of two indica varieties Guangluai 4 and Minghui 63. A database of the rice indica cDNAs was therefore built to provide a comprehensive web data source for searching and retrieving the indica cDNA clones.

RESULTS

Rice Indica cDNA Database (RICD) is an online MySQL-PHP driven database with a user-friendly web interface. It allows investigators to query the cDNA clones by keyword, genome position, nucleotide or protein sequence, and putative function. It also provides a series of information, including sequences, protein domain annotations, similarity search results, SNPs and InDels information, and hyperlinks to gene annotation in both The Rice Annotation Project Database (RAP-DB) and The TIGR Rice Genome Annotation Resource, expression atlas in RiceGE and variation report in Gramene of each cDNA.

CONCLUSION

The online rice indica cDNA database provides cDNA resource with comprehensive information to researchers for functional analysis of indica subspecies and for comparative genomics. The RICD database is available through our website http://www.ncgr.ac.cn/ricd.

Collection and comparative analysis of 1888 full-length cDNAs from wild rice Oryza rufipogon Griff. W1943

A huge amount of cDNA and EST resources have been developed for cultivated rice species Oryza sativa; however, only few cDNA resources are available for wild rice species. In this study, we isolated and completely sequenced 1888 putative full-length cDNA (FLcDNA) clones from wild rice Oryza rufipogon Griff. W1943 for comparative analysis between wild and cultivated rice species. Two cDNA libraries were constructed from 3-week-old leaf samples under either normal or cold-treated conditions. Homology searching of these cDNA sequences revealed that >96.8% of the wild rice cDNAs were matched to the cultivated rice O. sativa ssp. japonica cv. Nipponbare genome sequence. However, <22% of them were fully matched to the cv. Nipponbare genome sequence. The comparative analysis showed that O. rufipogon W1943 had greater similarity to O. sativa ssp. japonica than to ssp. indica cultivars. In addition, 17 novel rice cDNAs were identified, and 41 putative tissue-specific expression genes were defined through searching the rice massively parallel signature-sequencing database. In conclusion, these FLcDNA clones are a resource for further function verification and could be broadly utilized in rice biological studies.

Genome-wide analysis of transposon insertion polymorphisms reveals intraspecific variation in cultivated rice

Insertions and precise eliminations of transposable elements generated numerous transposon insertion polymorphisms (TIPs) in rice (Oryza sativa). We observed that TIPs represent more than 50% of large insertions and deletions (>100 bp) in the rice genome. Using a comparative genomic approach, we identified 2,041 TIPs between the genomes of two cultivars, japonica Nipponbare and indica 93-11. We also identified 691 TIPs between Nipponbare and indica Guangluai 4 in the 23-Mb collinear regions of chromosome 4. Among them, retrotransposon-based insertion polymorphisms were used to reveal the evolutionary relationships of these three cultivars. Our conservative estimates suggest that the TIPs generated approximately 14% of the genomic DNA sequence differences between subspecies indica and japonica. It was also found that more than 10% of TIPs were located in expressed gene regions, representing an important source of genetic variation. Transcript evidence implies that these TIPs induced a series of genetic differences between two subspecies, including interrupting host genes, creating different expression forms, drastically changing intron length, and affecting expression levels of adjacent genes. These analyses provide genome-wide insights into evolutionary history and genetic variation of rice.

Rice 2020: a call for an international coordinated effort in rice functional genomics

We describe a call for an international coordinated effort in rice functional genomics in the form of a project named RICE2020. The mission of the project will be: to determine the function of every gene in the rice genome by the year 2020, to identify functional diversity of alleles for agriculturally useful genes from the primary gene pool of rice, and to apply the findings of functional genomics research to rice genetic improvement.

Characterization and comparison of intron structure and alternative splicing between Medicago truncatula, Populus trichocarpa, Arabidopsis and rice

Alignment of transcripts and genome sequences yielded a set of alternatively spliced transcripts in four angiosperm genomes: three dicotyledon species Medicago truncatula (Medicago), Populus trichocarpa (poplar) and Arabidopsis thaliana (Arabidopsis), and the monocotyledon Oryzae sativa (rice). Intron retention was the predominant mode of alternative splicing (AS) in each species, consistent with previous reports for Arabidopsis and rice. We analyzed the structure of 5'-splice junctions and observed commonalities between species. There was dependency of base composition between sites flanking the 5'-splice junction, with the potential to create a subset of splice sites that interact more weakly or strongly than average with U1 snRNA. Such altered nucleotide composition was correlated with splicing fidelity in all four species. For Medicago, poplar and Arabidopsis, but not in rice, alternative splicing was most prevalent for introns with decreased UA content, consistent with lower UA content for monocot introns and potentially reflecting evolved differences in splicing mechanisms. Similarly, the occurrence of AS between transcript Gene Ontology categories was positively correlated between Arabidopsis and Medicago, with no correlation between dicots and rice. Analysis of within-species paralogs and between-species reciprocal best-hit homologs yielded rare cases of potentially conserved AS events. Reverse transcriptase PCR and amplicon sequencing were used to confirm a subset of the in silico-predicted AS events within Medicago, as well as to characterize conserved AS events between Medicago and Arabidopsis.

Strategies for developing Green Super Rice

From a global viewpoint, a number of challenges need to be met for sustainable rice production: (i) increasingly severe occurrence of insects and diseases and indiscriminate pesticide applications; (ii) high pressure for yield increase and overuse of fertilizers; (iii) water shortage and increasingly frequent occurrence of drought; and (iv) extensive cultivation in marginal lands. A combination of approaches based on the recent advances in genomic research has been formulated to address these challenges, with the long-term goal to develop rice cultivars referred to as Green Super Rice. On the premise of continued yield increase and quality improvement, Green Super Rice should possess resistances to multiple insects and diseases, high nutrient efficiency, and drought resistance, promising to greatly reduce the consumption of pesticides, chemical fertilizers, and water. Large efforts have been focused on identifying germplasms and discovering genes for resistance to diseases and insects, N- and P-use efficiency, drought resistance, grain quality, and yield. The approaches adopted include screening of germplasm collections and mutant libraries, gene discovery and identification, microarray analysis of differentially regulated genes under stressed conditions, and functional test of candidate genes by transgenic analysis. Genes for almost all of the traits have now been isolated in a global perspective and are gradually incorporated into genetic backgrounds of elite cultivars by molecular marker-assisted selection or transformation. It is anticipated that such strategies and efforts would eventually lead to the development of Green Super Rice.