RNA-seq reveals differentially expressed genes of rice (Oryza sativa) spikelet in response to temperature interacting with nitrogen at meiosis stage

Overexpression of lncRNA22524 from Dongxiang Wild Rice Reduces Drought and Salt Stress Tolerance in Cultivated Rice

Drought and salt stresses are major challenges to rice production, and a deep understanding of the mechanisms for tolerance could help deal with the challenges. Long non-coding RNAs (lncRNAs) play crucial roles in gene regulation. Previously, lncRNA22524 has been identified as a drought stress-responsive lncRNA from Dongxiang wild rice (DXWR). Nevertheless, its reactions to abiotic stresses in genetics and physiology remained unclear. In this study, we employed a rapid amplification of cDNA ends (RACE) to obtain the full-length cDNA of lncRNA22524 from DXWR, analyzed its cellular localization, built an overexpression vector to generate transgenic lines of cultivated rice and evaluated its impact in genetics and physiology. After treated with drought and salt stress, the overexpressed lines exhibited much more injuries and lower rates of survival, more reactive oxygen species (ROS) and malondialdehyde (MDA), lower antioxidant enzymes and lower proline (Pro) and soluble sugar (SS) than their wild-type (WT). Furthermore, transcriptome analysis of overexpressed lines with weaker tolerance than WT revealed 1,233 differentially expressed genes (DEGs), where most DEGs were involved in phenylpropanoid biosynthesis, photosynthesis and glutathione metabolism. These findings demonstrated that lncRNA22524 negatively regulated rice responses to drought and salt stress, which clear way of working from transcription to metabolic products should be worth of further study.

Increased panicle nitrogen application improves rice yield by alleviating high-temperature damage during panicle initiation to anther development

The grain yield is closely associated with spikelet fertility in rice (Oryza sativa L.) under high temperatures, and nitrogen (N) plays a crucial role in yield formation. To investigate the effect of panicle N application on yield formation under high temperatures at the panicle initiation stage, two rice varieties [Liangyoupeijiu (LYPJ, heat susceptible) and Shanyou63 (SY63, heat tolerant)] were grown and exposed to high daytime temperature (HT) and control temperature (Control) during the panicle initiation stage. Low (LPN) and high (HPN) panicle N applications were conducted. HT markedly decreased the yields by 87% at LPN and 48% at HPN in LYPJ and 31% at LPN and 36% at HPN in SY63. The decrease in grain yield under HT was primarily attributed to the decline in spikelet fertility, HPN increased spikelet fertility. HT resulted in the abnormal development of anthers, which included disordered, enlarged, and broken anther wall layers, degraded and irregularly shaped microspores, delayed tapetum degradation, less vacuolated microspores per locule, abnormal and aborted pollen grains; however, HPN improved the development of anthers under HT, particularly in LYPJ. A high rate of evapotranspiration resulted in an approximately 1°C decrease in panicle temperatures at HPN compared with that at LPN in both varieties under HT. Overall, these results demonstrate that the increased panicle N application favors normal anther development in LYPJ by decreasing the panicle temperature, which results in high pollen viability and spikelet fertility, and consequently less yield loss under HT.

Transcriptome analysis of mulberry (Morus alba L.) leaves to identify differentially expressed genes associated with post-harvest shelf-life elongation

Present study deals with molecular expression patterns responsible for post-harvest shelf-life extension of mulberry leaves. Quantitative profiling showed retention of primary metabolite and accumulation of stress markers in NS7 and CO7 respectively. The leaf mRNA profiles was sequenced using the Illumina platform to identify DEGs. A total of 3413 DEGs were identified between the treatments. Annotation with Arabidopsis database has identified 1022 DEGs unigenes. STRING generated protein-protein interaction, identified 1013 DEGs nodes with p < 1.0e-16. KEGG classifier has identified genes and their participating biological processes. MCODE and BiNGO detected sub-networking and ontological enrichment, respectively at p ≤ 0.05. Genes associated with chloroplast architecture, photosynthesis, detoxifying ROS and RCS, and innate-immune response were significantly up-regulated, responsible for extending shelf-life in NS7. Loss of storage sucrose, enhanced activity of senescence-related hormones, accumulation of xenobiotics, and development of osmotic stress inside tissue system was the probable reason for tissue deterioration in CO7. qPCR validation of DEGs was in good agreement with RNA sequencing results, indicating the reliability of the sequencing platform. Present outcome provides a molecular insight regarding involvement of genes in self-life extension, which might help the sericulture industry to overcome their pre-existing problems related to landless farmers and larval feeding during monsoon.

Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars

Identified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.

Sugar metabolism during pre- and post-fertilization events in plants under high temperature stress

High temperature challenges global crop production by limiting the growth and development of the reproductive structures and seed. It impairs the developmental stages of male and female gametogenesis, pollination, fertilization, endosperm formation and embryo development. Among these, the male reproductive processes are highly prone to abnormalities under high temperature at various stages of development. The disruption of source-sink balance is the main constraint for satisfactory growth of the reproductive structures which is disturbed at the level of sucrose import and utilization within the tissue. Seed development after fertilization is affected by modulation in the activity of enzymes involved in starch metabolism. In addition, the alteration in the seed-filling rate and its duration affects the seed weight and quality. The present review critically discusses the role of sugar metabolism in influencing the various stages of gamete and seed development under high temperature stress. It also highlights the interaction of the sugars with hormones that mediate the transport of sugars to sink tissues. The role of transcription factors for the regulation of sugar availability under high temperature has also been discussed. Further, the omics-based systematic investigation has been suggested to understand the synergistic or antagonistic interactions between sugars, hormones and reactive oxygen species at various points of sucrose flow from source to sink under high temperature stress.

Comparative Analysis Based on Transcriptomics and Metabolomics Data Reveal Differences between Emmer and Durum Wheat in Response to Nitrogen Starvation

Mounting evidence indicates the key role of nitrogen (N) on diverse processes in plant, including development and defense. Using a combined transcriptomics and metabolomics approach, we studied the response of seedlings to N starvation of two different tetraploid wheat genotypes from the two main domesticated subspecies: emmer and durum wheat. We found that durum wheat exhibits broader and stronger response in comparison to emmer as seen from the expression pattern of both genes and metabolites and gene enrichment analysis. They showed major differences in the responses to N starvation for transcription factor families, emmer showed differential reduction in the levels of primary metabolites while durum wheat exhibited increased levels of most of them to N starvation. The correlation-based networks, including the differentially expressed genes and metabolites, revealed tighter regulation of metabolism in durum wheat in comparison to emmer. We also found that glutamate and γ-aminobutyric acid (GABA) had highest values of centrality in the metabolic correlation network, suggesting their critical role in the genotype-specific response to N starvation of emmer and durum wheat, respectively. Moreover, this finding indicates that there might be contrasting strategies associated to GABA and glutamate signaling modulating shoot vs. root growth in the two different wheat subspecies.

RNA Sequencing Reveals Rice Genes Involved in Male Reproductive Development under Temperature Alteration

Rice (Oryza sativa L.) is one of the most important food crops, providing food for nearly half of the world population. Rice grain yields are affected by temperature changes. Temperature stresses, both low and high, affect male reproductive development, resulting in yield reduction. Thermosensitive genic male sterility (TGMS) rice is sterile at high temperature and fertile at low temperature conditions, facilitating hybrid production, and is a good model to study effects of temperatures on male development. Semithin sections of the anthers of a TGMS rice line under low (fertile) and high (sterile) temperature conditions showed differences starting from the dyad stage, suggesting that genes involved in male development play a role during postmeiotic microspore development. Using RNA sequencing (RNA-Seq), transcriptional profiling of TGMS rice panicles at the dyad stage revealed 232 genes showing differential expression (DEGs) in a sterile, compared to a fertile, condition. Using qRT-PCR to study expression of 20 selected DEGs using panicles of TGMS and wild type rice plants grown under low and high temperature conditions, revealed that six out of the 20 selected genes may be unique to TGMS, while the other 14 genes showed common responses to temperatures in both TGMS and wild-type rice plants. The results presented here would be useful for further investigation into molecular mechanisms controlling TGMS and rice responses to temperature alteration.

Morpho-physiological integrators, transcriptome and coexpression network analyses signify the novel molecular signatures associated with axillary bud in chrysanthemum

BACKGROUND

Axillary bud is an important agronomic and economic trait in cut chrysanthemum. Bud outgrowth is an intricate process controlled by complex molecular regulatory networks, physio-chemical integrators and environmental stimuli. Temperature is one of the key regulators of bud's fate. However, little is known about the temperature-mediated control of axillary bud at molecular levels in chrysanthemum. A comprehensive study was designed to study the bud outgrowth at normal and elevated temperature in cut chrysanthemum. Leaf morphology, histology, physiological parameters were studied to correlate the leaf activity with bud morphology, sucrose and hormonal regulation and the molecular controllers.

RESULTS

Temperature caused differential bud outgrowth along bud positions. Photosynthetic leaf area, physiological indicators and sucrose utilization were changed considerable due to high temperature. Comparative transcriptome analysis identified a significant proportion of bud position-specific genes.Weighted Gene Co-expression Network Analysis (WGCNA) showed that axillary bud control can be delineated by modules of coexpressed genes; especially, MEtan3, MEgreen2 and MEantiquewhite presented group of genes specific to bud length. A comparative analysis between different bud positions in two temperatures revealed the morpho-physiological traits associated with specific modules. Moreover, the transcriptional regulatory networks were configured to identify key determinants of bud outgrowth. Cell division, organogenesis, accumulation of storage compounds and metabolic changes were prominent during the bud emergence.

CONCLUSIONS

RNA-seq data coupled with morpho-physiological integrators from three bud positions at two temperature regimes brings a robust source to understand bud outgrowth status influenced by high temperature in cut chrysanthemum. Our results provide helpful information for elucidating the regulatory mechanism of temperature on axillary bud growth in chrysanthemum.

Transcriptome Reveals the Rice Response to Elevated Free Air CO2 Concentration and TiO2 Nanoparticles

Increasing CO₂ levels are speculated to change the effects of engineered nanomaterials in soil and on plant growth. How plants will respond to a combination of elevated CO₂ and nanomaterials stress has rarely been investigated, and the underlying mechanism remains largely unknown. Here, we conducted a field experiment to investigate the rice (Oryza sativa L. cv. IIyou) response to TiO₂ nanoparticles (nano-TiO₂, 0 and 200 mg kg^-1) using a free-air CO₂ enrichment system with different CO₂ levels (ambient ∼370 μmol mol^-1 and elevated ∼570 μmol mol^-1). The results showed that elevated CO₂ or nano-TiO₂ alone did not significantly affect rice chlorophyll content and antioxidant enzyme activities. However, in the presence of nano-TiO₂, elevated CO₂ significantly enhanced the rice height, shoot biomass, and panicle biomass (by 9.4%, 12.8%, and 15.8%, respectively). Furthermore, the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that genes involved in photosynthesis were up-regulated while most genes associated with secondary metabolite biosynthesis were down-regulated in combination-treated rice. This indicated that elevated CO₂ and nano-TiO₂ might stimulate rice growth by adjusting resource allocation between photosynthesis and metabolism. This study provides novel insights into rice responses to increasing contamination under climate change.

iTRAQ-based quantitative proteomic and physiological analysis of the response to N deficiency and the compensation effect in rice

Background

The crop growth compensation effect is a naturally biological phenomenon, and nitrogen (N) is essential for crop growth and development, especially for yield formation. Little is known about the molecular mechanism of N deficiency and N compensation in rice. Thus, the N-sensitive stage of rice was selected to study N deficiency at the tillering stage and N compensation at the young panicle differentiation stage. In this study, a proteome analysis was performed to analyze leaf differentially expressed proteins (DEPs), and to investigate the leaf physiological characteristics and yield under N deficiency and after N compensation.

Results

The yield per plant presented an equivalent compensatory effect. The net photosynthetic rate, optimal/maximal quantum yield of photosystem II (Fv/Fm), soil and plant analyzer development (SPAD) value, and glutamic pyruvic transaminase (GPT) activity of T1 (N deficiency at the tillering stage, and N compensation at the young panicle differentiation stage) were lower than those of CK (N at different stages of growth by constant distribution) under N deficiency. However, after N compensation, the net photosynthetic rate, Fv/Fm, SPAD value and GPT activity were increased. Using an iTRAQ-based quantitative approach, a total of 1665 credible proteins were identified in the three 4-plex iTRAQ experiments. Bioinformatics analysis indicated that DEPs were enriched in photosynthesis, photosynthesis-antenna proteins, carbon metabolism and carbon fixation in the photosynthetic organism pathways. Moreover, the photosynthesis-responsive proteins of chlorophyll a-b binding protein, ribulose bisphosphate carboxylase small chain and phosphoglycerate kinase were significantly downregulated under N deficiency. After N compensation, chlorophyll a-b binding protein, NADH dehydrogenase subunit 5, NADH dehydrogenase subunit 7, and peroxidase proteins were significantly upregulated in rice leaves.

Conclusion

Through physiological and quantitative proteomic analysis, we concluded that a variety of metabolic pathway changes was induced by N deficiency and N compensation. GO and KEGG enrichment analysis revealed that DEPs were significantly associated with photosynthesis pathway-, energy metabolism pathway- and stress resistance-related proteins. The DEPs play an important role in the regulation of N deficiency and the compensation effect in rice.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6031-4) contains supplementary material, which is available to authorized users.

Transcriptomic and Metabolomic Analysis of the Heat-Stress Response of Populus tomentosa Carr

Plants have evolved mechanisms of stress tolerance responses to heat stress. However, little is known about metabolic responses to heat stress in trees. In this study, we exposed Populus tomentosa Carr. to control (25 °C) and heat stress (45 °C) treatments and analyzed the metabolic and transcriptomic effects. Heat stress increased the cellular concentration of H2O2 and the activities of antioxidant enzymes. The levels of proline, raffinose, and melibiose were increased by heat stress, whereas those of pyruvate, fumarate, and myo-inositol were decreased. The expression levels of most genes (PSB27, PSB28, LHCA5, PETB, and PETC) related to the light-harvesting complexes and photosynthetic electron transport system were downregulated by heat stress. Association analysis between key genes and altered metabolites indicated that glycolysis was enhanced, whereas the tricarboxylic acid (TCA) cycle was suppressed. The inositol phosphate; galactose; valine, leucine, and isoleucine; and arginine and proline metabolic pathways were significantly affected by heat stress. In addition, several transcription factors, including HSFA2, HSFA3, HSFA9, HSF4, MYB27, MYB4R1, and bZIP60 were upregulated, whereas WRKY13 and WRKY50 were downregulated by heat stress. Interestingly, under heat stress, the expression of DREB1, DREB2, DREB2E, and DREB5 was dramatically upregulated at 12 h. Our results suggest that proline, raffinose, melibiose, and several genes (e.g., PSB27, LHCA5, and PETB) and transcription factors (e.g., HSFAs and DREBs) are involved in the response to heat stress in P. tomentosa.

Comparative transcriptome analysis of two selenium-accumulating genotypes of Aegilops tauschii Coss. in response to selenium

Background

Selenium (Se), an essential micronutrient in both animals and humans, has various biological functions, and its deficiency can lead to various diseases. The most common method for increasing Se uptake is the consumption of Se-rich plants, which transform inorganic Se into organic forms. Wheat is eaten daily by many people. The Se content of Aegilops tauschii (Ae. tauschii), one of the ancestors of hexaploid common wheat, is generally higher than that of wheat. In this study, two genotypes of Ae. tauschii with contrasting Se-accumulating abilities were subjected to different Se treatments followed by high-throughput transcriptome sequencing.

Results

Sequencing of 12 transcriptome libraries of Ae. tauschii grown under different Se treatments produced about a total of 47.72 GB of clean reads. After filtering out rRNA sequences, approximately 19.3 million high-quality clean reads were mapped to the reference genome (ta IWGSC_MIPSv2.1 genome DA). The total number of reference genome gene is 32,920 and about 26,407 known genes were detected in four groups. Functional annotation of these mapped genes revealed a large number of genes and some pivotal pathways that may participate in Se metabolism. The expressions of several genes potentially involved in Se metabolism were confirmed by quantitative real-time PCR.

Conclusions

Our study, the first to examine Se metabolism in Ae. tauschii, has provided a theoretical foundation for future elucidation of the mechanism of Se metabolism in this species.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0700-1) contains supplementary material, which is available to authorized users.

A comparison of heat-stress transcriptome changes between wild-type Arabidopsis pollen and a heat-sensitive mutant harboring a knockout of cyclic nucleotide-gated cation channel 16 (cngc16)

Background

In flowering plants, the male gametophyte (pollen) is one of the most vulnerable cells to temperature stress. In Arabidopsis thaliana, a pollen-specific CyclicNucleotide-Gated cationChannel 16 (cngc16), is required for plant reproduction under temperature-stress conditions. Plants harboring a cncg16 knockout are nearly sterile under conditions of hot days and cold nights. To understand the underlying cause, RNA-Seq was used to compare the pollen transcriptomes of wild type (WT) and cngc16 under normal and heat stress (HS) conditions.

Results

Here we show that a heat-stress response (HSR) in WT pollen resulted in 2102 statistically significant transcriptome changes (≥ 2-fold changes with adjusted p-value ≤0.01), representing approximately 15% of 14,226 quantified transcripts. Of these changes, 89 corresponded to transcription factors, with 27 showing a preferential expression in pollen over seedling tissues. In contrast to WT, cngc16 pollen showed 1.9-fold more HS-dependent changes (3936 total, with 2776 differences between WT and cngc16). In a quantitative direct comparison between WT and cngc16 transcriptomes, the number of statistically significant differences increased from 21 pre-existing differences under normal conditions to 192 differences under HS. Of the 20 HS-dependent changes in WT that were most different in cngc16, half corresponded to genes encoding proteins predicted to impact cell wall features or membrane dynamics.

Conclusions

Results here define an extensive HS-dependent reprogramming of approximately 15% of the WT pollen transcriptome, and identify at least 27 transcription factor changes that could provide unique contributions to a pollen HSR. The number of statistically significant transcriptome differences between WT and cngc16 increased by more than 9-fold under HS, with most of the largest magnitude changes having the potential to specifically impact cell walls or membrane dynamics, and thereby potentiate cngc16 pollen to be hypersensitive to HS. However, HS-hypersensitivity could also be caused by the extensive number of differences throughout the transcriptome having a cumulative effect on multiple cellular pathways required for tip growth and fertilization. Regardless, results here support a model in which a functional HS-dependent reprogramming of the pollen transcriptome requires a specific calcium-permeable Cyclic Nucleotide-Gated cation Channel, CNGC16.

Electronic supplementary material

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

Insight into Genes Regulating Postharvest Aflatoxin Contamination of Tetraploid Peanut from Transcriptional Profiling

Postharvest aflatoxin contamination is a challenging issue that affects peanut quality. Aflatoxin is produced by fungi belonging to the Aspergilli group, and is known as an acutely toxic, carcinogenic, and immune-suppressing class of mycotoxins. Evidence for several host genetic factors that may impact aflatoxin contamination has been reported, e.g., genes for lipoxygenase (PnLOX1 and PnLOX2/PnLOX3 that showed either positive or negative regulation with Aspergillus infection), reactive oxygen species, and WRKY (highly associated with or differentially expressed upon infection of maize with Aspergillus flavus); however, their roles remain unclear. Therefore, we conducted an RNA-sequencing experiment to differentiate gene response to the infection by A. flavus between resistant (ICG 1471) and susceptible (Florida-07) cultivated peanut genotypes. The gene expression profiling analysis was designed to reveal differentially expressed genes in response to the infection (infected vs. mock-treated seeds). In addition, the differential expression of the fungal genes was profiled. The study revealed the complexity of the interaction between the fungus and peanut seeds as the expression of a large number of genes was altered, including some in the process of plant defense to aflatoxin accumulation. Analysis of the experimental data with "keggseq," a novel designed tool for Kyoto Encyclopedia of Genes and Genomes enrichment analysis, showed the importance of α-linolenic acid metabolism, protein processing in the endoplasmic reticulum, spliceosome, and carbon fixation and metabolism pathways in conditioning resistance to aflatoxin accumulation. In addition, coexpression network analysis was carried out to reveal the correlation of gene expression among peanut and fungal genes. The results showed the importance of WRKY, toll/Interleukin1 receptor-nucleotide binding site leucine-rich repeat (TIR-NBS-LRR), ethylene, and heat shock proteins in the resistance mechanism.

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

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

Transcriptomic response of durum wheat to nitrogen starvation

Nitrogen (N) is a key macronutrient representing a limiting factor for plant growth and development and affects productivity in wheat. In this study, durum wheat response to N chronic starvation during grain filling was investigated through a transcriptomic approach in roots, leaves/stems, flag leaf and spikes of cv. Svevo. Nitrogen stress negatively influenced plant height, tillering, flag leaf area, spike and seed traits, and total N content. RNA-seq data revealed 4,626 differentially expressed genes (DEGs). Most transcriptomic changes were observed in roots, with 3,270 DEGs, while 963 were found in leaves/stems, 470 in flag leaf, and 355 in spike tissues. A total of 799 gene ontology (GO) terms were identified, 180 and 619 among the upregulated and downregulated genes, respectively. Among the most addressed GO categories, N compound metabolism, carbon metabolism, and photosynthesis were mostly represented. Interesting DEGs, such as N transporters, genes involved in N assimilation, along with transcription factors, protein kinases and other genes related to stress were highlighted. These results provide valuable information about the transcriptomic response to chronic N stress in durum wheat, which could be useful for future improvement of N use efficiency.

Transcriptome profiling of rice seedlings under cold stress

Rice (Oryza sativa L.) is one of the most important species for food production worldwide, besides being an excellent genetic model among the grasses. Cold is one of the major abiotic factors reducing rice yield, primarily affecting germination and reproduction phases. Currently, the RNAseq technique allows the identification of differential expressed genes in response to a given treatment, such as cold stress. In the present work, a transcriptome (RNAseq) analysis was performed in the V3 phase for contrasting genotypes Oro (tolerant) and Tio Taka (sensitive), in response to cold (13°C). A total of 241 and 244M readings were obtained, resulting in the alignment of 25.703 and 26.963 genes in genotypes Oro and Tio Taka respectively. The analyses revealed 259 and 5579 differential expressed genes in response to cold in the genotypes Oro and Tio Taka respectively. Ontology classes with larger changes were metabolic process ~27%, cellular process ~21%, binding ~30% and catalytic activity ~22%. In the genotype Oro, 141 unique genes were identified, 118 were common between Oro and Tio Taka and 5461 were unique to Tio Taka. Genes involved in metabolic routes of signal transduction, phytohormones, antioxidant system and biotic stress were identified. These results provide an understanding that breeding for a quantitative trait, such as cold tolerance at germination, several gene loci must be simultaneously selected. In general, few genes were identified, but it was not possible to associate only one gene function as responsible for the cultivar tolerance; since different genes from different metabolic routes were identified. The genes described in the present work will be useful for future investigations and for the detailed validation in marker assisted selection projects for cold tolerance in the germination of rice.

Integrated analysis of rice transcriptomic and metabolomic responses to elevated night temperatures identifies sensitivity- and tolerance-related profiles

Transcript and metabolite profiling were performed on leaves from six rice cultivars under high night temperature (HNT) condition. Six genes were identified as central for HNT response encoding proteins involved in transcription regulation, signal transduction, protein-protein interactions, jasmonate response and the biosynthesis of secondary metabolites. Sensitive cultivars showed specific changes in transcript abundance including abiotic stress responses, changes of cell wall-related genes, of ABA signaling and secondary metabolism. Additionally, metabolite profiles revealed a highly activated TCA cycle under HNT and concomitantly increased levels in pathways branching off that could be corroborated by enzyme activity measurements. Integrated data analysis using clustering based on one-dimensional self-organizing maps identified two profiles highly correlated with HNT sensitivity. The sensitivity profile included genes of the functional bins abiotic stress, hormone metabolism, cell wall, signaling, redox state, transcription factors, secondary metabolites and defence genes. In the tolerance profile, similar bins were affected with slight differences in hormone metabolism and transcription factor responses. Metabolites of the two profiles revealed involvement of GABA signaling, thus providing a link to the TCA cycle status in sensitive cultivars and of myo-inositol as precursor for inositol phosphates linking jasmonate signaling to the HNT response specifically in tolerant cultivars.