Deep sequencing reveals the complex and coordinated transcriptional regulation of genes related to grain quality in rice cultivars

Transcriptome and Metabolome Analyses Reveal New Insights into the Regulatory Mechanism of Head Milled Rice Rate

The head milled rice rate (HMRR) is the most important trait of milling quality, which affects the final yield and quality of rice. However, few genes related to HMRR have been identified and the regulatory mechanism of HMRR remains elusive. In this study, we performed a comparative analysis integrating the transcriptome sequencing of developing seeds at the grain-filling stage and a metabolome analysis of brown rice between two groups of accessions with contrasting performances in HMRR. A total of 768 differentially expressed genes (DEGs) were identified between the transcriptome profiles of low-HMRR and high-HMRR accessions. In comparison to the high-HMRR accessions, 655 DEGs were up-regulated in the low-HMRR accessions, which was 4.79 folds higher than the number of down-regulated genes. These up-regulated DEGs were enriched in various metabolic and biosynthetic processes, oxidation reduction, phosphorylation, ion transport and ATP-related processes. However, the 113 down-regulated DEGs in the low-HMRR accessions were concentrated in carbohydrate metabolic processes, cell-death-related processes and defense response. Among the 30 differential metabolites, 20 and 10 metabolites were down-/up-regulated, respectively, in the accessions with low HMRR. In addition, 10 differential metabolites, including five metabolites of the shikimate pathway and five metabolites of the pyruvate pathway, were integrated into two separate pathways, starting from sucrose. Our global analysis of HMRR provides an invaluable resource for a better understanding of the molecular mechanism underlying the genetic regulation of HMRR.

The Function of DNA Demethylase Gene ROS1a Null Mutant on Seed Development in Rice (Oryza Sativa) Using the CRISPR/CAS9 System

The endosperm is the main nutrient source in cereals for humans, as it is a highly specialized storage organ for starch, lipids, and proteins, and plays an essential role in seed growth and development. Active DNA demethylation regulates plant developmental processes and is ensured by cytosine methylation (5-meC) DNA glycosylase enzymes. To find out the role of OsROS1a in seed development, the null mutant of OsROS1a was generated using the CRISPR/Cas9 system. The null mutant of OsROS1a was stable and heritable, which affects the major agronomic traits, particularly in rice seeds. The null mutant of OsROS1a showed longer and narrower grains, and seeds were deformed containing an underdeveloped and less-starch-producing endosperm with slightly irregularly shaped embryos. In contrast to the transparent grains of the wild type, the grains of the null mutant of OsROS1a were slightly opaque and rounded starch granules, with uneven shapes, sizes, and surfaces. A total of 723 differential expression genes (DEGs) were detected in the null mutant of OsROS1a by RNA-Seq, of which 290 were downregulated and 433 were upregulated. The gene ontology (GO) terms with the top 20 enrichment factors were visualized for cellular components, biological processes, and molecular functions. The key genes that are enriched for these GO terms include starch synthesis genes (OsSSIIa and OsSSIIIa) and cellulose synthesis genes (CESA2, CESA3, CESA6, and CESA8). Genes encoding polysaccharides and glutelin were found to be downregulated in the mutant endosperm. The glutelins were further verified by SDS-PAGE, suggesting that glutelin genes could be involved in the null mutant of OsROS1a seed phenotype and OsROS1a could have the key role in the regulation of glutelins. Furthermore, 378 differentially alternative splicing (AS) genes were identified in the null mutant of OsROS1a, suggesting that the OsROS1a gene has an impact on AS events. Our findings indicated that the function on rice endosperm development in the null mutant of OsROS1a could be influenced through regulating gene expression and AS, which could provide the base to properly understand the molecular mechanism related to the OsROS1a gene in the regulation of rice seed development.

Study of expressions of miRNAs in the spikelets based on their spatial location on panicle in rice cultivars provided insight into their influence on grain development

Development of rice cultivars bearing numerous spikelets by breeding approach to increase the yearly production of rice to approximately 800 million metric tons to feed the ever increasing population of the world accompanies poor grain filling in the inferior spikelets preventing achievement of the yield potential. As the initial stages of caryopses development are of much importance for grain filling, spatio-temporal expressions of the miRNAs were studied during these periods in the spikelets of a compact-panicle rice cultivar, Oryza sativa cv. Mahalaxmi, bearing numerous spikelets per panicle to understand the reason of poor grain filling at the level of the initial biochemical events. Differential expression of several known miRNAs between the superior and inferior spikelets suggested great difference in metabolism related to grain filling in the spikelets based on their spatial location on compact panicle. Expressions of five known and four novel miRNAs were validated by Northern. Their targets included the enzymes directly involved in starch biosynthesis like sucrose synthase, starch synthase and pullulanase, besides others. Spatio-temporal expression studies of these miRNAs in the spikelets of Mahalaxmi revealed a pattern of mostly a greater expression in the inferior spikelets compared with the superior ones concomitant with an inverse expression of the target genes, which was not observed in the lax-panicle cultivar Upahar. The study thus revealed that the grain filling in rice is greatly regulated by miRNAs, and these miRNAs or their target genes could be considered for biotechnological interventions for improving grain filling in the rice cultivars of interest.

Exploring the Relationships Between Yield and Yield-Related Traits for Rice Varieties Released in China From 1978 to 2017

Despite evidence from previous case studies showing that agronomic traits partially determine the resulting yield of different rice (Oryza sativa L.) varieties, it remains unclear whether this is true at the ecotype level. Here, an extensive dataset of the traits of 7686 rice varieties, released in China from 1978 to 2017, was used to study the relationship between yield and other agronomic traits. We assessed the association between yield and other agronomic traits for four different rice ecotypes, i.e., indica inbred, indica hybrid, japonica inbred, and japonica hybrid. We found that associations between agronomic traits and yield were ecotype-dependent. For both the indica inbred and indica hybrid ecotypes, we found that greater values of certain traits, including the filled grain number per panicle, 1000-grain-weight, plant height, panicle length, grains per panicle, seed setting rate, long growth period, low panicle number per unit area, and low seed length/width ratio, have accounted for high grain yield. In the japonica inbred and japonica hybrid ecotypes, we found that only high panicle number per unit area and long growth period led to high grain yield. Indirectly, growth period consistently had a positive effect on yield in all ecotypes, and plant height had a positive effect on yield for the indicas and japonica inbred only. Plant height had a negative effect for the japonica hybrid. Altogether, our findings potentially have valuable implications for improving the breeds of rice ecotypes.

Rice grain nutritional traits and their enhancement using relevant genes and QTLs through advanced approaches

BACKGROUND

Rice breeding program needs to focus on development of nutrient dense rice for value addition and helping in reducing malnutrition. Mineral and vitamin deficiency related problems are common in the majority of the population and more specific to developing countries as their staple food is rice.

RESULTS

Genes and QTLs are recently known for the nutritional quality of rice. By comprehensive literature survey and public domain database, we provided a critical review on nutritional aspects like grain protein and amino acid content, vitamins and minerals, glycemic index value, phenolic and flavonoid compounds, phytic acid, zinc and iron content along with QTLs linked to these traits. In addition, achievements through transgenic and advanced genomic approaches have been discussed. The information available on genes and/or QTLs involved in enhancement of micronutrient element and amino acids are summarized with graphical representation.

CONCLUSION

Compatible QTLs/genes may be combined together to design a desirable genotype with superior in multiple grain quality traits. The comprehensive review will be helpful to develop nutrient dense rice cultivars by integrating molecular markers and transgenic assisted breeding approaches with classical breeding.

Comparative proteomics of the superior and inferior spikelets at the early grain filling stage in rice cultivars contrast for panicle compactness and ethylene evolution

The breeding programmes in rice aimed at increasing the number of spikelets per panicle have been accompanied by poor grain filling in the inferior spikelets of large panicle rice, leading to yield disadvantage. The present study attempted to understand the reason for differential grain filling in the inferior and superior spikelets by comparative proteomics considering a compact-panicle rice cultivar Mahalaxmi and a lax-panicle rice cultivar Upahar, which show poor and good grain filling, respectively. An initial study of two rice cultivars for panicle compactness and grain filling revealed an inverse correlation between the two parameters. It was further observed that the panicle compactness in Mahalaxmi was associated with a higher evolution of ethylene by the spikelets, both superior and inferior, compared with the lax-panicle Upahar. The proteomic studies revealed that the superior and inferior spikelets of Mahalaxmi differentially expressed 21 proteins that were also expressed in Upahar. However, in Upahar, only two of these proteins were differentially expressed between the superior and inferior spikelets, indicating that the metabolic activities of the spikelets occupying the superior and inferior positions on the panicle were very different in Mahalaxmi compared with those in Upahar. Among the proteins that were downregulated in the inferior spikelets compared with the superior ones in Mahalaxmi were importin-α, elongation factor 1-β and cell division control protein 48, which are essential for cell cycle progression and cell division. Low expression of these proteins might inhibit endosperm cell division in the inferior spikelets, limiting their sink capacity and leading to poor grain filling compared to that in the superior spikelets. The poor grain filling in Mahalaxmi may also be a result of the high evolution of ethylene in the inferior spikelets, which is reflected from the observation that these spikelets showed significantly higher expression of S-adenosylmethionine synthase and the gene encoding the enzyme than the superior spikelets in this cultivar, but not in Upahar; S-adenosynlmethionine synthase catalyses the synthesis of S-adenosylmethionine, the precursor of ethylene biosynthesis.

Rice false smut fungus hijacks the rice nutrients supply by blocking and mimicking the fertilization of rice ovary

Rice false smut disease is an increasing threat to rice production in the world. Despite of best efforts, research for the infection of the fungus has yielded equivocal and conflicting results about where and how the infection is initiated and developed. Here we show a stepwise infection pattern and sophisticated regulation during this process. Initial infection occurred on the filaments, which prevented the production of mature pollen thus blocked the pollination. In the following days, the pathogen invaded the stigmas and styles, occasionally the ovaries. Expression analysis indicated that the fungus mimicked a successful fertilization process and enabled the continuous supply of nutrients for fungus to produce false smut balls. The stepwise infection of flower organs and mimicry of ovary fertilization unveiled in this study guided the rice plant into supplying nutrients for false smut ball development and represents a new and unique biological process of host pathogen interactions.

Proteomic Dissection of Endosperm Starch Granule Associated Proteins Reveals a Network Coordinating Starch Biosynthesis and Amino Acid Metabolism and Glycolysis in Rice Endosperms

Starch biosynthesis and starch granule packaging in cereal endosperms involve a coordinated action of starch biosynthesis enzymes and coordination with other metabolisms. Because directly binding to starch granules, starch granule-associated proteins (SGAPs) are essential to understand the underlying mechanisms, however the information on SGAPs remains largely unknown. Here, we dissected developmentally changed SGAPs from developing rice endosperms from 10 to 20 days after flowering (DAF). Starch granule packaging was not completed at 10 DAF, and was finished in the central endosperm at 15 DAF and in the whole endosperm at 20 DAF. Proteomic analysis with two-dimensional differential in-gel electrophoresis and mass spectrometry revealed 115 developmentally changed SGAPs, representing 37 unique proteins. 65% of the unique proteins had isoforms. 39% of the identified SGAPs were involved in starch biosynthesis with main functions in polyglucan elongation and granule structure trimming. Almost all proteins involved in starch biosynthesis, amino acid biosynthesis, glycolysis, protein folding, and PPDK pathways increased abundance as the endosperm developed, and were predicted in an interaction network. The network represents an important mechanism to orchestrate carbon partitioning among starch biosynthesis, amino acid biosynthesis and glycolysis for efficient starch and protein storage. These results provide novel insights into mechanisms of starch biosynthesis and its coordination with amino acid metabolisms and glycolysis in cereal endosperms.

De novo analysis of transcriptome reveals genes associated with leaf abscission in sugarcane (Saccharum officinarum L.)

Background

Sugarcane (Saccharum officinarum L.) is an important sugar crop which belongs to the grass family and can be used for fuel ethanol production. The growing demands for sugar and biofuel is asking for breeding a sugarcane variety that can shed their leaves during the maturity time due to the increasing cost on sugarcane harvest.

Results

To determine leaf abscission related genes in sugarcane, we generated 524,328,950 paired reads with RNA-Seq and profiled the transcriptome of new born leaves of leaf abscission sugarcane varieties (Q1 and T) and leaf packaging sugarcane varieties (Q2 and B). Initially, 275,018 transcripts were assembled with N50 of 1,177 bp. Next, the transcriptome was annotated by mapping them to NR, UniProtKB/Swiss-Prot, Gene Ontology and KEGG pathway databases. Further, we used TransDecoder and Trinotate to obtain the likely proteins and annotate them in terms of known proteins, protein domains, signal peptides, transmembrane regions and rRNA transcripts. Different expression analysis showed 1,202 transcripts were up regulated in leaf abscission sugarcane varieties, relatively to the leaf packaging sugarcane varieties. Functional analysis told us 62, 38 and 10 upregulated transcripts were involved in plant-pathogen interaction, response to stress and abscisic acid associated pathways, respectively. The upregulation of transcripts encoding 4 disease resistance proteins (RPM1, RPP13, RGA2, and RGA4), 6 ABC transporter G family members and 16 transcription factors including WRK33 and heat stress transcription factors indicate they may be used as candidate genes for sugarcane breeding. The expression levels of transcripts were validated by qRT-PCR. In addition, we characterized 3,722 SNPs between leaf abscission and leaf packaging sugarcane plants.

Conclusion

Our results showed leaf abscission associated genes in sugarcane during the maturity period. The output of this study provides a valuable resource for future genetic and genomic studies in sugarcane.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2552-2) contains supplementary material, which is available to authorized users.

An RNA-Seq-based reference transcriptome for Citrus

Previous RNA-Seq studies in citrus have been focused on physiological processes relevant to fruit quality and productivity of the major species, especially sweet orange. Less attention has been paid to vegetative or reproductive tissues, while most Citrus species have never been analysed. In this work, we characterized the transcriptome of vegetative and reproductive tissues from 12 Citrus species from all main phylogenetic groups. Our aims were to acquire a complete view of the citrus transcriptome landscape, to improve previous functional annotations and to obtain genetic markers associated with genes of agronomic interest. 28 samples were used for RNA-Seq analysis, obtained from 12 Citrus species: C. medica, C. aurantifolia, C. limon, C. bergamia, C. clementina, C. deliciosa, C. reshni, C. maxima, C. paradisi, C. aurantium, C. sinensis and Poncirus trifoliata. Four different organs were analysed: root, phloem, leaf and flower. A total of 3421 million Illumina reads were produced and mapped against the reference C. clementina genome sequence. Transcript discovery pipeline revealed 3326 new genes, the number of genes with alternative splicing was increased to 19,739, and a total of 73,797 transcripts were identified. Differential expression studies between the four tissues showed that gene expression is overall related to the physiological function of the specific organs above any other variable. Variants discovery analysis revealed the presence of indels and SNPs in genes associated with fruit quality and productivity. Pivotal pathways in citrus such as those of flavonoids, flavonols, ethylene and auxin were also analysed in detail.

Transcriptomic and Metabolomic Studies Disclose Key Metabolism Pathways Contributing to Well-maintained Photosynthesis under the Drought and the Consequent Drought-Tolerance in Rice

In contrast to wild species, drought-tolerance in crops requires a fully functional metabolism during drought (particularly photosynthetic processes). However, the link between drought-tolerance, photosynthetic regulation during drought, and the associated transcript and metabolic foundation, remains largely unknown. For this study, we used two rice cultivars with contrasting drought-tolerance (the drought-intolerant cultivar IRAT109 and the drought-tolerant cultivar IAC1246) to explore transcript and metabolic responses to long-term drought. The drought-tolerant cultivar represented higher osmotic adjustment and antioxidant capacity, as well as higher relative photosynthesis rate under a progressive drought stress occurred in a modified field with shallow soil-layers. A total of 4059 and 2677 differentially expressed genes (DEGs) were identified in IRAT109 and IAC1246 between the drought and well-watered conditions, respectively. A total of 69 and 47 differential metabolites (DMs) were identified between the two treatments in IRAT109 and IAC1246, respectively. Compared to IRAT109, the DEGs of IAC1246 displayed enhanced regulatory amplitude during drought. We found significant correlations between DEGs and the osmolality and total antioxidant capacity (AOC) of both cultivars. During the early stages of drought, we detected up-regulation of DEGs in IAC1246 related to photosynthesis, in accordance with its higher relative photosynthesis rate. The contents of six differential metabolites were correlated with the osmotic potential and AOC. Moreover, they were differently regulated between the two cultivars. Particularly, up-regulations of 4-hydroxycinnamic acid and ferulic acid were consistent with the performance of photosynthesis-related DEGs at the early stages of drought in IAC1246. Therefore, 4-hydroxycinnamic acid and ferulic acid were considered as key metabolites for rice drought-tolerance. DEGs involved in pathways of these metabolites are expected to be good candidate genes to improve drought-tolerance. In conclusion, well-maintained photosynthesis under drought should contribute to improved drought-tolerance in rice. Metabolites play vital roles in protecting photosynthesis under dehydration via osmotic adjustments and/or antioxidant mechanisms. A metabolite-based method was thus an effective way to explore drought candidate genes. Metabolic accompanied by transcript responses to drought stress should be further studied to find more useful metabolites, pathways, and genes.

Identification and Characterization of Differentially Expressed Genes in Inferior and Superior Spikelets of Rice Cultivars with Contrasting Panicle-Compactness and Grain-Filling Properties

Breeding programs for increasing spikelet number in rice have resulted in compactness of the panicle, accompanied by poor grain filling in inferior spikelets. Although the inefficient utilization of assimilate has been indicated as responsible for this poor grain filling, the underlying cause remains elusive. The current study utilized the suppression subtractive hybridization technique to identify 57 and 79 genes that overexpressed in the superior and inferior spikelets (with respect to each other), respectively, of the compact-panicle rice cultivar Mahalaxmi. Functional categorization of these differentially expressed genes revealed a marked metabolic difference between the spikelets according to their spatial location on the panicle. The expression of genes encoding seed storage proteins was dominant in inferior spikelets, whereas genes encoding regulatory proteins, such as serine-threonine kinase, zinc finger protein and E3 ligase, were highly expressed in superior spikelets. The expression patterns of these genes in the inferior and superior spikelets of Mahalaxmi were similar to those observed in another compact-panicle cultivar, OR-1918, but differed from those obtained in two lax-panicle cultivars, Upahar and Lalat. The results first suggest that the regulatory proteins abundantly expressed in the superior spikelets of compact-panicle cultivars and in both the superior and inferior spikelets of lax-panicle cultivars but poorly expressed in the inferior spikelets of compact-panicle cultivars promote grain filling. Second, the high expression of seed-storage proteins observed in the inferior spikelets of compact-panicle cultivars appears to inhibit the grain filling process. Third, the low expression of enzymes of the Krebs cycle in inferior spikelets compared with superior spikelets of compact-panicle cultivars is bound to lead to poor ATP generation in the former and consequently limit starch biosynthesis, an ATP-consuming process, resulting in poor grain filling.

Deep sequencing transcriptional fingerprinting of rice kernels for dissecting grain quality traits

BACKGROUND

Rice represents one the most important foods all over the world. In Europe, Italy is the first rice producer and Italian production is driven by tradition and quality. All main rice grain quality traits, like cooking properties, texture, gelatinization temperature, chalkiness and yield, are related to the content and composition of starch and seed-storage proteins in the endosperm and to grain shape. In addition, a number of nutraceutical compounds and allergens are known to have a significant effect on grain quality determination. To investigate the genetic bases underlying the qualitative differences that characterize traditional Italian rice cultivars, a comparative RNA-Seq-based transcriptomic analysis of developing caryopsis was conducted at 14 days after flowering on six popular Italian varieties (Carnaroli, Arborio, Balilla, Vialone Nano, Gigante Vercelli and Volano) phenotypically differing for qualitative grain-related traits.

RESULTS

Co-regulation analyses of differentially expressed genes showing the same expression patterns in the six genotypes highlighted clusters of loci up or down-regulated in specific varieties, with respect to the others. Among them, we detected loci involved in cell wall biosynthesis, protein metabolism and redox homeostasis, classes of genes affecting in chalkiness determination. Moreover, loci encoding for seed-storage proteins, allergens or involved in the biosynthesis of specific nutraceutical compounds were also present and specifically regulated in the different clusters. A wider investigation of all the DEGs detected in pair-wise comparisons revealed transcriptional variation, among the six genotypes, for quality-related loci involved in starch biosynthesis (e.g. GBSSI, starch synthases and AGPase), genes encoding for transcription factors, additional seed storage proteins, allergens or belonging to additional nutraceutical compounds biosynthetic pathways and loci affecting grain size. Putative functional SNPs associated to amylose content in starch, gelatinization temperature and grain size were also identified.

CONCLUSIONS

The present work represents a more extended phenotypic characterization of a set of rice accessions that present a wider genetic variability than described nowadays in literature. The results provide the first transcriptional picture for several of the grain quality differences observed among the Italian rice varieties analyzed and reveal that each variety is characterized by the over-expression of a peculiar set of loci affecting grain appearance and quality. A list of candidates and SNPs affecting specific grain properties has been identified offering a starting point for further works aimed to characterize genes and molecular markers for breeding programs.

Identification of candidate genes associated with positive and negative heterosis in rice

To identify the genes responsible for yield related traits, and heterosis, massively parallel signature sequencing (MPSS) libraries were constructed from leaves, roots and meristem tissues from the two parents, 'Nipponbare' and '93-11', and their F1 hybrid. From the MPSS libraries, 1-3 million signatures were obtained. Using cluster analysis, commonly and specifically expressed genes in the parents and their F1 hybrid were identified. To understand heterosis in the F1 hybrid, the differentially expressed genes in the F1 hybrid were mapped to yield related quantitative trait loci (QTL) regions using a linkage map constructed from 131 polymorphic simple sequence repeat markers with 266 recombinant inbred lines derived from a cross between Nipponbare and 93-11. QTLs were identified for yield related traits including days to heading, plant height, plant type, number of tillers, main panicle length, number of primary branches per main panicle, number of kernels per main panicle, total kernel weight per main panicle, 1000 grain weight and total grain yield per plant. Seventy one QTLs for these traits were mapped, of which 3 QTLs were novel. Many highly expressed chromatin-related genes in the F1 hybrid encoding histone demethylases, histone deacetylases, argonaute-like proteins and polycomb proteins were located in these yield QTL regions. A total of 336 highly expressed transcription factor (TF) genes belonging to 50 TF families were identified in the yield QTL intervals. These findings provide the starting genomic materials to elucidate the molecular basis of yield related traits and heterosis in rice.

Improvement of marker-based predictability of Apparent Amylose Content in japonica rice through GBSSI allele mining

BACKGROUND

Apparent Amylose Content (AAC), regulated by the Waxy gene, represents the key determinant of rice cooking properties. In occidental countries high AAC rice represents the most requested market class but the availability of molecular markers allowing specific selection of high AAC varieties is limited.

RESULTS

In this study, the effectiveness of available molecular markers in predicting AAC was evaluated in a collection of 127 rice accessions (125 japonica ssp. and 2 indica ssp.) characterized by AAC values from glutinous to 26%. The analyses highlighted the presence of several different allelic patterns identifiable by a few molecular markers, and two of them, i.e., the SNPs at intron1 and exon 6, were able to explain a maximum of 79.5% of AAC variation. However, the available molecular markers haplotypes did not provide tools for predicting accessions with AAC higher than 24.5%. To identify additional polymorphisms, the re-sequencing of the Waxy gene and 1kbp of the putative upstream regulatory region was performed in 21 genotypes representing all the AAC classes identified. Several previously un-characterized SNPs were identified and four of them were used to develop dCAPS markers.

CONCLUSIONS

The addition of the SNPs newly identified slightly increased the AAC explained variation and allowed the identification of a haplotype almost unequivocally associated to AAC higher than 24.5%. Haplotypes at the waxy locus were also associated to grain length and length/width (L/W) ratio. In particular, the SNP at the first intron, which identifies the Wxa and Wxb alleles, was associated with differences in the width of the grain, the L/W ratio and the length of the kernel, most likely as a result of human selection.

QTL map meets population genomics: an application to rice

Genes involved in the transition from wild to cultivated crop species should be of great agronomic importance. Population genomic approaches utilizing genome resequencing data have been recently applied for this purpose, although it only reports a large list of candidate genes with no biological information. Here, by resequencing more than 30 genomes altogether of wild rice Oryza rufipogon and cultivated rice O. sativa, we identified a number of regions with clear footprints of selection during the domestication process. We then focused on identifying candidate domestication genes in these regions by utilizing the wealth of QTL information in rice. We were able to identify a number of interesting candidates such as transcription factors that should control key domestication traits such as shattering, awn length, and seed dormancy. Other candidates include those that might have been related to the improvement of grain quality and those that might have been involved in the local adaptation to dry conditions and colder environments. Our study shows that population genomic approaches and QTL mapping information can be used together to identify genes that might be of agronomic importance.

Genome-wide analysis of endosperm-specific genes in rice

The endosperm of the cereal crop is an important nutrient source for humans. It also acts as a critical integrator of plant seed growth and development. Despite its importance, the comprehensive understanding in regulating of endosperm development in rice remains elusive. Here, we performed a genomic survey comprising the identification and functional characterization of the endosperm-specific genes (OsEnS) in rice using Affymetrix microarray data and Gene Ontology (GO) analysis. A total of 151 endosperm-specific genes were identified, and the expression patterns of 13 selected genes were confirmed by qRT-PCR analysis. Promoter regions of the endosperm-specific expression genes were analyzed by PLACE Signal Scan Search. The results indicated that some motifs were involved in endosperm-specific expression regulation, and some cis-elements were responsible for hormone regulation. The bootstrap analysis indicated that the RY repeat (CATGCA box) was over-represented in promoter regions of endosperm-specific expression genes. GO analysis indicated that these genes could be classified into 12 groups, namely, transcription factor, stress/defense, seed storage protein (SSP), carbohydrate and energy metabolism, seed maturation, protein metabolism, lipid metabolism, transport, cell wall related, hormone related, signal transduction, and one unclassified group. Taken together, our results provide informative clues for further functional characterization of the endosperm-specific genes, which facilitate the understanding of the molecular mechanism in rice endosperm development.

WikiPathways for plants: a community pathway curation portal and a case study in rice and arabidopsis seed development networks

BACKGROUND

Next-generation sequencing and 'omics' platforms are used extensively in plant biology research to unravel new genomes and study their interactions with abiotic and biotic agents in the growth environment. Despite the availability of a large and growing number of genomic data sets, there are only limited resources providing highly-curated and up-to-date metabolic and regulatory networks for plant pathways.

RESULTS

Using PathVisio, a pathway editor tool associated with WikiPathways, we created a gene interaction network of 430 rice (Oryza sativa) genes involved in the seed development process by curating interactions reported in the published literature. We then applied an InParanoid-based homology search to these genes and used the resulting gene clusters to identify 351 Arabidopsis thaliana genes. Using this list of homologous genes, we constructed a seed development network in Arabidopsis by processing the gene list and the rice network through a Perl utility software called Pathway GeneSWAPPER developed by us. In order to demonstrate the utility of these networks in generating testable hypotheses and preliminary analysis prior to more in-depth downstream analysis, we used the expression viewer and statistical analysis features of PathVisio to analyze publicly-available and published microarray gene expression data sets on diurnal photoperiod response and the seed development time course to discover patterns of coexpressed genes found in the rice and Arabidopsis seed development networks. These seed development networks described herein, along with other plant pathways and networks, are freely available on the plant pathways portal at WikiPathways (http://plants.wikipathways.org).

CONCLUSION

In collaboration with the WikiPathways project we present a community curation and analysis platform for plant biologists where registered users can freely create, edit, share and monitor pathways supported by published literature. We describe the curation and annotation of a seed development network in rice, and the projection of a similar, gene homology-based network in Arabidopsis. We also demonstrate the utility of the Pathway GeneSWAPPER (PGS) application in saving valuable time and labor when a reference network in one species compiled in GPML format is used to project a similar network in another species based on gene homology.

Transcriptomic analysis of rice (Oryza sativa) endosperm using the RNA-Seq technique

The endosperm plays an important role in seed formation and germination, especially in rice (Oryza sativa). We used a high-throughput sequencing technique (RNA-Seq) to reveal the molecular mechanisms involved in rice endosperm development. Three cDNA libraries were taken from rice endosperm at 3, 6 and 10 days after pollination (DAP), which resulted in the detection of 21,596, 20,910 and 19,459 expressed gens, respectively. By ERANGE, we identified 10,371 differentially expressed genes (log(2)Ratio ≥1, FDR ≤0.001). The results were compared against three public databases (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes and MapMan) in order to annotate the gene descriptions, associate them with gene ontology terms and to assign each to pathways. A large number of genes related to ribosomes, the spliceosome and oxidative phosphorylation were found to be expressed in the early and middle stages. Plant hormone, galactose metabolism and carbon fixation related genes showed a significant increase in expression at the middle stage, whereas genes for defense against disease or response to stress as well as genes for starch/sucrose metabolism were strongly expressed during the later stages of endosperm development. Interestingly, most metabolic pathways were down-regulated between 3 and 10 DAP except for those involved in the accumulation of material, such as starch/sucrose and protein metabolism. We also identified the expression of 1,118 putative transcription factor genes in endosperm development. The RNA-Seq results provide further systematic understanding of rice endosperm development at a fine scale and a foundation for future studies.

Identification and characterization of in planta-expressed secreted effector proteins from Magnaporthe oryzae that induce cell death in rice

Interactions between rice and Magnaporthe oryzae involve the recognition of cellular components and the exchange of complex molecular signals from both partners. How these interactions occur in rice cells is still elusive. We employed robust-long serial analysis of gene expression, massively parallel signature sequencing, and sequencing by synthesis to examine transcriptome profiles of infected rice leaves. A total of 6,413 in planta-expressed fungal genes, including 851 genes encoding predicted effector proteins, were identified. We used a protoplast transient expression system to assess 42 of the predicted effector proteins for the ability to induce plant cell death. Ectopic expression assays identified five novel effectors that induced host cell death only when they contained the signal peptide for secretion to the extracellular space. Four of them induced cell death in Nicotiana benthamiana. Although the five effectors are highly diverse in their sequences, the physiological basis of cell death induced by each was similar. This study demonstrates that our integrative genomic approach is effective for the identification of in planta-expressed cell death-inducing effectors from M. oryzae that may play an important role facilitating colonization and fungal growth during infection.

Molecular dissection of atrazine-responsive transcriptome and gene networks in rice by high-throughput sequencing

The residue of atrazine (a herbicide) has become hazards in environments due to its intensive use. However, its molecular toxicity to on plants and human beings is not fully understood. In this study, we performed high-throughput sequencing of atrazine-exposed rice (Oryza sativa) to analyze global expression and complexity of genes in the crop. Four libraries were constructed from shoots and roots with or without atrazine exposure. We sequenced 5,751,861, 5,790,013, 5,375,999 and 6,039,618 clean tags that corresponded to 220,806, 111,301, 248,802 and 114,338 distinct tags for Root-Atr (root control, atrazine-free), Shoot-Atr (shoot control, atrazine-free), Root+Atr (root treated with atrazine) and Shoot+Atr (shoot treated with atrazine) libraries, respectively. Mapping the clean tags to gene databases generated 18,833-21,007 annotated genes for each library. Most of annotated genes were differentially expressed among the libraries. The most 40 differentially expressed genes were associated with resistance to environmental stress, degradation of xenobiotics and molecular metabolism. Validation of gene expression by quantitative RT-PCR confirmed the deep-sequencing results. The transcriptome sequences were further subjected to Gene Orthology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and showed modified biological functions and metabolic pathways. Our results not only highlight the transcriptional complexity in rice with atrazine but also represent a major improvement for analyzing transcriptional changes on a large scale in xenobiotics-responsive toxicology.

Genome-wide analysis of the complex transcriptional networks of rice developing seeds

Background

The development of rice (Oryza sativa) seed is closely associated with assimilates storage and plant yield, and is fine controlled by complex regulatory networks. Exhaustive transcriptome analysis of developing rice embryo and endosperm will help to characterize the genes possibly involved in the regulation of seed development and provide clues of yield and quality improvement.

Principal Findings

Our analysis showed that genes involved in metabolism regulation, hormone response and cellular organization processes are predominantly expressed during rice development. Interestingly, 191 transcription factor (TF)-encoding genes are predominantly expressed in seed and 59 TFs are regulated during seed development, some of which are homologs of seed-specific TFs or regulators of Arabidopsis seed development. Gene co-expression network analysis showed these TFs associated with multiple cellular and metabolism pathways, indicating a complex regulation of rice seed development. Further, by employing a cold-resistant cultivar Hanfeng (HF), genome-wide analyses of seed transcriptome at normal and low temperature reveal that rice seed is sensitive to low temperature at early stage and many genes associated with seed development are down-regulated by low temperature, indicating that the delayed development of rice seed by low temperature is mainly caused by the inhibition of the development-related genes. The transcriptional response of seed and seedling to low temperature is different, and the differential expressions of genes in signaling and metabolism pathways may contribute to the chilling tolerance of HF during seed development.

Conclusions

These results provide informative clues and will significantly improve the understanding of rice seed development regulation and the mechanism of cold response in rice seed.