Transcriptome profiling of rice seedlings under cold stress

Identification of candidate genes controlling cold tolerance at the early seedling stage from Dongxiang wild rice by QTL mapping, BSA-Seq and RNA-Seq

BACKGROUND

The cold tolerance of rice is closely related to its production and geographic distribution. The identification of cold tolerance-related genes is of important significance for developing cold-tolerant rice. Dongxiang wild rice (Oryza rufipogon Griff.) (DXWR) is well-adapted to the cold climate of northernmost-latitude habitats ever found in the world, and is one of the most valuable rice germplasms for cold tolerance improvement.

RESULTS

Transcriptome analysis revealed genes differentially expressed between Xieqingzao B (XB; a cold sensitive variety) and 19H19 (derived from an interspecific cross between DXWR and XB) in the room temperature (RT), low temperature (LT), and recovery treatments. The results demonstrated that chloroplast genes might be involved in the regulation of cold tolerance in rice. A high-resolution SNP genetic map was constructed using 120 BC5F2 lines derived from a cross between 19H19 and XB based on the genotyping-by-sequencing (GBS) technique. Two quantitative trait loci (QTLs) for cold tolerance at the early seedling stage (CTS), qCTS12 and qCTS8, were detected. Moreover, a total of 112 candidate genes associated with cold tolerance were identified based on bulked segregant analysis sequencing (BSA-seq). These candidate genes were divided into eight functional categories, and the expression trend of candidate genes related to 'oxidation-reduction process' and 'response to stress' differed between XB and 19H19 in the RT, LT and recovery treatments. Among these candidate genes, the expression level of LOC_Os12g18729 in 19H19 (related to 'response to stress') decreased in the LT treatment but restored and enhanced during the recovery treatment whereas the expression level of LOC_Os12g18729 in XB declined during recovery treatment. Additionally, XB contained a 42-bp deletion in the third exon of LOC_Os12g18729, and the genotype of BC5F2 individuals with a survival percentage (SP) lower than 15% was consistent with that of XB. Weighted gene coexpression network analysis (WGCNA) and modular regulatory network learning with per gene information (MERLIN) algorithm revealed a gene interaction/coexpression network regulating cold tolerance in rice. In the network, differentially expressed genes (DEGs) related to 'oxidation-reduction process', 'response to stress' and 'protein phosphorylation' interacted with LOC_Os12g18729. Moreover, the knockout mutant of LOC_Os12g18729 decreased cold tolerance in early rice seedling stage signifcantly compared with that of wild type.

CONCLUSIONS

In general, study of the genetic basis of cold tolerance of rice is important for the development of cold-tolerant rice varieties. In the present study, QTL mapping, BSA-seq and RNA-seq were integrated to identify two CTS QTLs qCTS8 and qCTS12. Furthermore, qRT-PCR, genotype sequencing and knockout analysis indicated that LOC_Os12g18729 could be the candidate gene of qCTS12. These results are expected to further exploration of the genetic mechanism of CTS in rice and improve cold tolerance of cultivated rice by introducing the cold tolerant genes from DXWR through marker-assisted selection.

Whole-transcriptome profiling and identification of cold tolerance-related ceRNA networks in japonica rice varieties

Introduction

Low-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined.

Methods

We conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress.

Results

We identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice.

Discussion and conclusion

We identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.

Comparative transcriptomic analysis of germinating rice seedlings to individual and combined anaerobic and cold stress

BACKGROUND

Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination.

RESULTS

The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein).

CONCLUSION

A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.

Transcriptomic profiling of the cold stress and recovery responsiveness of two contrasting Guizhou HE rice genotypes

BACKGROUND

At the seed germination stage, rice is sensitive to cold stress, which adversely affects its growth and development. Guizhou HE rice comprises several different landraces, most of which are cold tolerant.

OBJECTIVE

To identify differentially expressed genes and molecular mechanism underlying the cold tolerance of Guizhou HE.

METHODS

Two Guizhou HE genotypes, AC44 (cold-sensitive) and AC96 (cold-tolerant), which exhibit opposite phenotypes in response to cold treatment at the seed germination stage were used. Comprehensive gene expressions of AC44 and AC96 under 4 °C cold treatment and subsequent recovery conditions were comparatively analyzed by RNA sequencing.

RESULTS

Overall, 11,082 and 7749 differentially expressed genes were detected in AC44 and AC96, respectively. Comparative transcriptome analysis demonstrated that, compared with AC44, AC96 presented fewer upregulated and downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses demonstrated that AC96 presented more upregulated GO terms, especially terms associated with biological processes. However, AC44 presented more terms related to cellular components, mainly chloroplasts. Moreover, DEGs related to the auxin signaling pathway (including ARF and IAA family members) and transcription factors (including members of the F-box, bZIP, basic helix-loop-helix [bHLH], and MYB-like transcription factor families) were found to be expressed specifically in AC96; thus, these DEGs may be responsible for the cold tolerance of AC96.

CONCLUSIONS

These findings present information about the cold tolerance mechanism of Guizhou HE rice at the germination stage, providing valuable resources and candidate genes for breeding cold-tolerant rice genotypes.

A meta-analysis of microarray data revealed hub genes and transcription factors involved in drought stress response in rice (Oryza sativa L.)

Although there are various studies attempted to clarify the genetic mechanism of plant response to drought stress that reduces crop yield, a meta-analysis can integrate the results of them to provide a better picture of the issue. Therefore, in this study, several microarray datasets of rice were meta-analysed under drought stress and normal condition using the R packages. Accordingly, differentially expressed genes (meta-DEGs) were identified. The results showed 643 and 677 upregulated and downregulated genes, respectively. The significant common Gene Ontology (GO) terms between the up- and downregulated genes were responses to abiotic stimulus , water deprivation , oxygen-containing compound and abscisic acid . The transcription factors (TF) survey showed that bHLH under drought stress activates up genes 42% more than down genes while bzip Homeodomain activates down genes 54% more than up genes. The hub downregulated genes obtained from this study were mainly related to photosynthesis and the hub upregulated genes were mainly related to stress tolerance which include heat shock proteins (HSPs), late embryogenesis abundant (LEAs), calmodulin-like protein (CML), phosphatase 2C (PP2Cs) and IAA genes. Moreover, this meta-analysis data were compared with other experimental data and the results confirmed the up and down expression of them. Our findings can provide novel insights into the molecular mechanism of rice (Oryza sativa L.) response to drought stress.

Full-length transcriptome analysis of maize root tips reveals the molecular mechanism of cold stress during the seedling stage

BACKGROUND

As maize originated in tropical or subtropical zones, most maize germplasm is extremely sensitive to low temperatures during the seedling stage. Clarifying the molecular mechanism of cold acclimation would facilitate the breeding of cold tolerant maize varieties, which is one of the major sustainability factors for crop production. To meet this goal, we investigated two maize inbred lines with contrasting levels of cold tolerance at the seedling stage (IL85, a cold tolerant line; B73, a cold sensitive line), and performed full-length transcriptome sequencing on the root tips of seedlings before and after 24 h of cold treatment.

RESULTS

We identified 152,263 transcripts, including 20,993 novel transcripts, and determined per-transcript expression levels. A total of 1,475 transcripts were specifically up-regulated in the cold tolerant line IL85 under cold stress. GO enrichment analysis revealed that 25 transcripts were involved in reactive oxygen species (ROS) metabolic processes and 15 transcripts were related to the response to heat. Eight genes showed specific differential alternative splicing (DAS) in IL85 under cold stress, and were mainly involved in amine metabolism. A total of 1,111 lncRNAs were further identified, 62 of which were up-regulated in IL85 or B73 under cold stress, and their corresponding target genes were enriched in protein phosphorylation.

CONCLUSIONS

These results provide new insights into the molecular mechanism of cold acclimation during the seedling stage in maize, and will facilitate the development of cultivars with improved cold stress tolerance.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, and NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the nested association mapping pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long-term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

Comparative Transcriptome Analysis Revealing the Potential Mechanism of Low-Temperature Stress in Machilus microcarpa

Machilus microcarpa is a rare national tree species in China and possesses important ornamental and ecological value. M. microcarpa can be planted in low-temperature areas, depending on whether its seedlings can withstand the harm. To face this problem, the annual seedlings of M. microcarpa were subjected to five temperature treatments, and eight physiological indicators were measured. Furthermore, comparative transcriptome analysis was performed between M. microcarpa leaves treated at 25°C and -2.8°C. A total of 9,385 differentially expressed genes (DEGs) were involved in low-temperature stress in M. microcarpa. An upregulated (cobA) and five downregulated (HEM, CHLM, CRD, CLH, and PORA) genes associated with the porphyrin and chlorophyll metabolism pathway may reduce chlorophyll synthesis under low-temperature stress. Upregulation of six DEGs (two GAPDHs, PFK, PGAM, PDC, and PK) involved in the glycolysis/gluconeogenesis pathway provided energy for M. microcarpa under adverse cold conditions. Thirteen upregulated and seven downregulated genes related to antioxidant enzymes were also observed under low-temperature stress. Candidate transcription factors (TFs) played key roles in signal transduction under low-temperature stress in M. microcarpa, and quantitative real-time PCR (qRT-PCR) analysis validated the RNA-seq data. The results provide valuable information for further studies on the cold response mechanisms for low-temperature stress in M. microcarpa.

Comparative Proteomics Combined with Morphophysiological Analysis Revealed Chilling Response Patterns in Two Contrasting Maize Genotypes

Maize yield is significantly influenced by low temperature, particularly chilling stress at the maize seedling stage. Various physiological approaches have been established to resist chilling stress; however, the detailed proteins change patterns underlying the maize chilling stress response at the seedling stage remain unknown, preventing the development of breeding-based methods to resist chilling stress in maize. Thus, we performed comprehensive physiological, comparative proteomics and specific phytohormone abscisic acid (ABA) assay on different maize inbred lines (tolerant-line KR701 and sensitive-line hei8834) at different seedling stages (the first leaf stage and third leaf stage) under chilling stress. The results revealed several signalling proteins and pathways in response to chilling stress at the maize seedling stage. Meanwhile, we found ABA pathway was important for chilling resistance of tolerant-line KR701 at the first leaf stage. Related chilling-responsive proteins were further catalogued and analysed, providing a resource for further investigation and maize breeding.

Transcriptome Analysis of Near-Isogenic Lines Provides Novel Insights into Genes Associated with Seed Low-Temperature Germination Ability in Maize (Zea mays L.)

Maize originated from tropical regions and is extremely sensitive to low temperature during germination. Our previous work identified a major QTL, qp1ER1-1, for low temperature germination ability (LTGA) of maize. Here, we introgressed qp1ER1-1 from the tolerant line L220 into the sensitive line PH4CV to generate two near isogenic lines NIL^220-3 and NIL^220-25. When germinated under cold temperature for 25 days (Cold-25), the NILs showed similar seedling root length and shoot length to L220, but significantly higher than PH4CV. However, when germinated under cold temperature for 15 days (Cold-15) or under normal temperature (25 °C) for 3 days (CK-3), all lines showed similar seedling performance, indicating that introgression of qp1ER1-1 from L220 into PH4CV could improve LTGA of NIL^220-3 and NIL^220-25. The whole seedlings, including root and shoot, of Cold-15 and CK-3 were harvested for transcriptome analysis, when both stayed at a similar developmental stage. Dry seed embryo was sequenced as a non-germination control (CK-0). Compared with PH4CV, the tolerant line (L220, NIL^220-3, and NIL^220-25) specifically expressed genes (different expressed genes, DEGs) were identified for CK-0, Cold-15, and CK-3. Then, DEGs identified from Cold-15, but not from CK-0 or CK-3, were defined as tolerant line specifically expressed LTGA genes. Finally, 1786, 174, and 133 DEGs were identified as L220, NIL^220-3, and NIL^220-25 specifically expressed LTGA genes, respectively. Of them, 27 were common LTGA genes that could be identified from all three tolerant lines, with two (Zm00001d031209 and Zm00001d031292) locating in the confidence interval of qp1ER1-1. In addition, GO analysis revealed that L220 specifically expressed LTGA genes were majorly enriched in the cell division process and plasma membrane related categories. Taken together, these results provided new insight into the molecular mechanism of maize seed LTGA and facilitated the cloning of the qp1ER1-1 gene.

Selection of transcripts related to low-temperature tolerance using RNA sequencing from F2 plants between japonica and indica rice (Oryza sativa L.) cultivars

At low temperatures (18°C), seedlings of an indica rice (Oryza sativa L.) cultivar Kasalath showed symptoms of chlorosis, although the leaves of a japonica cultivar Arroz da Terra remained green. In this study, transcripts related to the chlorophyll content of rice seedlings grown at 18°C were investigated using RNA-sequencing (RNA-Seq) data for F2 crosses between cultivars Arroz da Terra and Kasalath, as well as their parental cultivars. Differential expression analysis revealed that gene ontology terms related to 'photosynthesis' were significantly enriched in lowly expressed genes at 18°C than at 25°C in Kasalath. However, the gene ontology terms related to 'response to stress' were significantly enriched in highly expressed genes at 18°C than at 25°C in Kasalath. When the F2 plants were grown at 18°C, their chlorophyll contents varied. Transcripts with expression levels related to chlorophyll content were statistically selected using RNA-Seq data from 21 F2 plants. In regression models, frequently selected genes included four photosynthetic and two stress-responsive genes. The expression values of four photosynthetic and two stress-responsive genes in high-frequency selected genes were significantly correlated with chlorophyll content not only in plants analysed using RNA-Seq but also in 95 F2 plants.

De novo transcriptome assembly and comparative transcriptomic analysis provide molecular insights into low temperature stress response of Canarium album

A de novo transcriptome analysis was performed in C. album, a temperature sensitive fruit tree in China, after treatment with varied temperatures. A total number of 168,385 transcripts were assembled, comprising of 109,439 unigenes, of which 70,530 were successfully annotated. Compared with control check group (CK), which was treated under 25 °C, the chilling stress (4 °C) treated group (CT), showed about 2810 up-regulated and 2567 down-regulated genes. Whereas, group treated under freezing (- 3 °C) stress (FT) showed an up-regulation and a down-regulation of 1748 and 1459 genes, respectively. GO classification analysis revealed that DEGs related to metabolic processes, single-organism metabolic process, and catalytic activity are significantly enriched in both CT and FT conditions. KEGG pathway enrichment analysis for both CT and FT treatments showed an enrichment of genes encoding or related to glycine/serine and threonine metabolism, alpha-linolenic acid metabolism, carotenoid biosynthesis, photosynthesis-antenna proteins, and circadian rhythm. However, genes related to photosynthesis, carbon fixation in photosynthetic organisms, glutathione metabolism, pyruvate metabolism, nicotinate and nicotinamide metabolism were specifically enriched in CT condition. Nevertheless, FT treatment induced genes related to plant-pathogen interaction, linoleic acid metabolism, plant hormone signal transduction and pentose phosphate pathway. Many of the genes involved in plant hormone signal transduction showed significantly different expression in both FT and CT conditions. However, the change was more evident in FT. Here we present the first of the reports for a de novo transcriptomic analysis in C. album, suggesting that the plant shows differential responses in chilling and freezing temperatures, where the hormone signaling and transduction contribute greatly to FT responses. Our study thus paves way for future research regarding functions of these potentially identified genes.

When rice gets the chills: comparative transcriptome profiling at germination shows WRKY transcription factor responses

Rice is vital for food security. Due to its tropical origin, rice suffers from cold temperatures that affect its entire life cycle. Key genes have been identified involved in cold tolerance. WRKYs are generally downstream of the MAPK cascade and can act together with VQ proteins to regulate stress-responsive genes. Chilling treatment was applied at germination to two rice genotypes (tolerant and sensitive). Shoots at S3 stage were collected for RNA-seq to identify OsWRKY, OsMAPKs and OsVQs expression. Relationships among MAPKs, WRKYs and VQs were predicted through correlation analysis. OsWRKY transcriptional regulation was predicted by in silico analysis of cis-regulatory elements. A total of 39 OsWRKYs were differentially expressed. OsWRKY21, OsWRK24 and OsWRKY69 are potential positive regulators, while OsWRKY10, OsWRK47, OsWRKY62, OsWRKY72 and OsWRKY77 are potential negative regulators, of chilling tolerance. 12 OsMAPKs were differentially expressed. OsMAPKs were downregulated and negatively correlated with the upregulated OsWRKYs in the tolerant genotype. 19 OsVQs were differentially expressed, three and six OsVQs were positively correlated with OsWRKYs in the tolerant and sensitive genotypes, respectively. Seven differentially expressed OsWRKYs have cold-responsive elements in their promoters and five upregulated OsWRKYs in the tolerant genotype contained the W-box motif. Chilling causes changes in OsWRKY, OsMAPK and OsVQ gene expression at germination. OsWRKYs may not act downstream of the MAPK cascade to coordinate chilling tolerance, but OsWRKYs may act with VQs to regulate chilling tolerance. Candidate OsWRKYs are correlated and have a W-box in the promoter, suggesting an auto-regulation mechanism.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

Transcriptomic profiling of germinating seeds under cold stress and characterization of the cold-tolerant gene LTG5 in rice

BACKGROUND

Low temperature is a limiting factor of rice productivity and geographical distribution. Wild rice (Oryza rufipogon Griff.) is an important germplasm resource for rice improvement. It has superior tolerance to many abiotic stresses, including cold stress, but little is known about the mechanism underlying its resistance to cold.

RESULTS

This study elucidated the molecular genetic mechanisms of wild rice in tolerating low temperature. Comprehensive transcriptome profiles of two rice genotypes (cold-sensitive ce 253 and cold-tolerant Y12-4) at the germinating stage under cold stress were comparatively analyzed. A total of 42.44-68.71 million readings were obtained, resulting in the alignment of 29,128 and 30,131 genes in genotypes 253 and Y12-4, respectively. Many common and differentially expressed genes (DEGs) were analyzed in the cold-sensitive and cold-tolerant genotypes. Results showed more upregulated DEGs in the cold-tolerant genotype than in the cold-sensitive genotype at four stages under cold stress. Gene ontology enrichment analyses based on cellular process, metabolic process, response stimulus, membrane part, and catalytic activity indicated more upregulated genes than downregulated ones in the cold-tolerant genotype than in the cold-sensitive genotype. Quantitative real-time polymerase chain reaction was performed on seven randomly selected DEGs to confirm the RNA Sequencing (RNA-seq) data. These genes showed similar expression patterns corresponding with the RNA-Seq method. Weighted gene co-expression network analysis (WGCNA) revealed Y12-4 showed more positive genes than 253 under cold stress. We also explored the cold tolerance gene LTG5 (Low Temperature Growth 5) encoding a UDP-glucosyltransferase. The overexpression of the LTG5 gene conferred cold tolerance to indica rice.

CONCLUSION

Gene resources related to cold stress from wild rice can be valuable for improving the cold tolerance of crops.

Integrated transcriptomic and metabolomic analyses of yellow horn (Xanthoceras sorbifolia) in response to cold stress

Xanthoceras sorbifolia, a medicinal and oil-rich woody plant, has great potential for biodiesel production. However, little study explores the link between gene expression level and metabolite accumulation of X. sorbifolia in response to cold stress. Herein, we performed both transcriptomic and metabolomic analyses of X. sorbifolia seedlings to investigate the regulatory mechanism of resistance to low temperature (4 °C) based on physiological profile analyses. Cold stress resulted in a significant increase in the malondialdehyde content, electrolyte leakage and activity of antioxidant enzymes. A total of 1,527 common differentially expressed genes (DEGs) were identified, of which 895 were upregulated and 632 were downregulated. Annotation of DEGs revealed that amino acid metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, galactose metabolism, fructose and mannose metabolism, and the citrate cycle (TCA) were strongly affected by cold stress. In addition, DEGs within the plant mitogen-activated protein kinase (MAPK) signaling pathway and TF families of ERF, WRKY, NAC, MYB, and bHLH were transcriptionally activated. Through metabolomic analysis, we found 51 significantly changed metabolites, particularly with the analysis of primary metabolites, such as sugars, amino acids, and organic acids. Moreover, there is an overlap between transcript and metabolite profiles. Association analysis between key genes and altered metabolites indicated that amino acid metabolism and sugar metabolism were enhanced. A large number of specific cold-responsive genes and metabolites highlight a comprehensive regulatory mechanism, which will contribute to a deeper understanding of the highly complex regulatory program under cold stress in X. sorbifolia.

Genome-wide lncRNAs identification and association analysis for cold-responsive genes at the booting stage in rice (Oryza sativa L.)

Long non-coding RNAs (lncRNAs) are essential regulators of a broad range of biological processes in plants. The spectacular progress made in next-generation sequencing technologies has enabled a genome-wide identification of lncRNAs in multiple plant species. In this study, a genome-wide lncRNA sequencing technology was used to identify cold-responsive lncRNAs at the booting stage in rice by comparing a tolerant variety, Kongyu131 (KY131) and a sensitive variety, Dongnong422 (DN422). A total of 1485 lncRNAs were identified, and 566 of these lncRNAs were defined as differential lncRNAs by comparing four samples. GO and KEGG enrichment analyses were performed, focusing on the cis- and trans- target genes of the differential lncRNAs. To identify cold-responsive genes, a meta-analysis was used to integrate 35 cold-tolerant QTLs at the booting stage. In summary, 12 candidate genes and their target lncRNAs were identified by qRT-PCR. LncTar was used to identify the interaction between lncRNAs and the candidate genes. In addition, 130 rice cultivars with rich genetic diversity were collected to verify the association of candidate genes with cold-resistance. The results revealed that five SNPs in LOC_Os07g42940, three SNP and one InDel in LOC_Os02g03410 were associated with cold-resistance at a significant level using association analysis. This study provides new gene resources and insights into cold-resistance research for rice.

Integrated RNA-Seq Analysis and Meta-QTLs Mapping Provide Insights into Cold Stress Response in Rice Seedling Roots

Rice (Oryza sativa L.) is a widely cultivated food crop around the world, especially in Asia. However, rice seedlings often suffer from cold stress, which affects their growth and yield. Here, RNA-seq analysis and Meta-QTLs mapping were performed to understand the molecular mechanisms underlying cold tolerance in the roots of 14-day-old seedlings of rice (RPY geng, cold-tolerant genotype). A total of 4779 of the differentially expressed genes (DEGs) were identified, including 2457 up-regulated and 2322 down-regulated DEGs. The GO, COG, KEEG, and Mapman enrichment results of DEGs revealed that DEGs are mainly involved in carbohydrate transport and metabolism, signal transduction mechanisms (plant hormone signal transduction), biosynthesis, transport and catabolism of secondary metabolites (phenylpropanoid biosynthesis), defense mechanisms, and large enzyme families mechanisms. Notably, the AP2/ERF-ERF, NAC, WRKY, MYB, C2H2, and bHLH transcription factors participated in rice’s cold–stress response and tolerance. On the other hand, we mapped the identified DEGs to 44 published cold–stress-related genes and 41 cold-tolerant Meta-QTLs regions. Of them, 12 DEGs were the published cold–stress-related genes and 418 DEGs fell into the cold-tolerant Meta-QTLs regions. In this study, the identified DEGs and the putative molecular regulatory network can provide insights for understanding the mechanism of cold stress tolerance in rice. In addition, DEGs in KEGG term-enriched terms or cold-tolerant Meta-QTLs will help to secure key candidate genes for further functional studies on the molecular mechanism of cold stress response in rice.

A Meta-Analysis of Comparative Transcriptomic Data Reveals a Set of Key Genes Involved in the Tolerance to Abiotic Stresses in Rice

Several environmental factors, such as drought, salinity, and extreme temperatures, negatively affect plant growth and development, which leads to yield losses. The tolerance or sensitivity to abiotic stressors are the expression of a complex machinery involving molecular, biochemical, and physiological mechanisms. Here, a meta-analysis on previously published RNA-Seq data was performed to identify the genes conferring tolerance to chilling, osmotic, and salt stresses, by comparing the transcriptomic changes between tolerant and susceptible rice genotypes. Several genes encoding transcription factors (TFs) were identified, suggesting that abiotic stress tolerance involves upstream regulatory pathways. A gene co-expression network defined the metabolic and signalling pathways with a prominent role in the differentiation between tolerance and susceptibility: (i) the regulation of endogenous abscisic acid (ABA) levels, through the modulation of genes that are related to its biosynthesis/catabolism, (ii) the signalling pathways mediated by ABA and jasmonic acid, (iii) the activity of the "Drought and Salt Tolerance" TF, involved in the negative regulation of stomatal closure, and (iv) the regulation of flavonoid biosynthesis by specific MYB TFs. The identified genes represent putative key players for conferring tolerance to a broad range of abiotic stresses in rice; a fine-tuning of their expression seems to be crucial for rice plants to cope with environmental cues.

Transcriptomics profiling in response to cold stress in cultivated rice and weedy rice

Weedy rice is an important germplasm resource for rice improvement because it has useful genes for many abiotic stresses including cold tolerance. We identified the cold tolerance and cold sensitivity of two weedy rice lines (WR 03-35 and WR 03-26) and two cultivated rice lines (Kongyu 131 and 9311). During the seedling stage of these lines, we used RNA-seq to measure changes in weedy rice and cultivated rice whole-genome transcriptome before and after cold treatment. We identified 14,213 and 14,730 differentially expressed genes (DEGs) in cold-tolerant genotypes (WR 03-35, Kongyu 131), and 9219 and 720 DEGs were observed in two cold-sensitive genotypes (WR 03-26, 9311). Many common and special DEGs were analyzed in cold-tolerant and cold-sensitive genotypes, respectively. Some typical genes related to cold stress such as the basic helix-loop-helix (bHLH) gene and leucine-rich repeat (LRR) domain gene etc. The number of these DEGs in cold-tolerant genotypes is more than those found in cold-sensitive genotypes. The gene ontology (GO) enrichment analyses showed significantly enriched terms for biological processes, cellular components and molecular functions. In addition, some genes related to several plant hormones such as abscisic acid (ABA), gibberellic acid (GA), auxin and ethylene were identified. To confirm the RNA-seq data, semi-quantitative RT-PCR and qRT-PCR were performed on 12 randomly selected DEGs. The expression patterns of RNA-seq on these genes corresponded with the semi-quantitative RT-PCR and qRT-PCR method. This study suggests the gene resources related to cold stress from weedy rice could be valuable for understanding the mechanisms involved in cold stress and rice breeding for improving cold tolerance.

Transcriptome profiling of short-term response to chilling stress in tolerant and sensitive Oryza sativa ssp. Japonica seedlings

Low temperature is a major factor limiting rice growth and yield, and seedling is one of the developmental stages at which sensitivity to chilling stress is higher. Tolerance to chilling is a complex quantitative trait, so one of the most effective approaches to identify genes and pathways involved is to compare the stress-induced expression changes between tolerant and sensitive genotypes. Phenotypic responses to chilling of 13 Japonica cultivars were evaluated, and Thaibonnet and Volano were selected as sensitive and tolerant genotypes, respectively. To thoroughly profile the short-term response of the two cultivars to chilling, RNA-Seq was performed on Thaibonnet and Volano seedlings after 0 (not stressed), 2, and 10 h at 10 °C. Differential expression analysis revealed that the ICE-DREB1/CBF pathway plays a primary role in chilling tolerance, mainly due to some important transcription factors involved (some of which had never been reported before). Moreover, the expression trends of some genes that were radically different between Thaibonnet and Volano (i.e., calcium-dependent protein kinases OsCDPK21 and OsCDPK23, cytochrome P450 monooxygenase CYP76M8, etc.) suggest their involvement in low temperature tolerance too. Density of differentially expressed genes along rice genome was determined and linked to the position of known QTLs: remarkable co-locations were reported, delivering an overview of genomic regions determinant for low temperature response at seedling stage. Our study contributes to a better understanding of the molecular mechanisms underlying rice response to chilling and provides a solid background for development of low temperature-tolerant germplasm.