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

Subtilisin-like protease 4 regulates cold tolerance through cell wall modification in rice

Rice is susceptible to cold temperatures, especially during the seedling stage. Despite extensive research into the cold tolerance mechanisms of rice, the number of cloned genes remains limited. Plant subtilisin-like proteases (SUBs or SBTs) are protein-hydrolyzing enzymes which play important roles in various aspects of plant growth as well as the plant response to biotic and abiotic stress. The rice SUB gene family consists of 62 members, but it is unknown whether they are involved in the response to cold stress. In this study, we observed that a loss-of-function SUB4 mutant exhibited enhanced cold tolerance at the seedling stage. The sub4 mutant seedlings exhibited improved survival rates and related physiological parameters, including relative electrolyte conductivity, chlorophyll content, malondialdehyde content, and antioxidant enzyme activity. Transcriptomic analysis revealed that differentially expressed genes responsive to cold stress in the sub4 mutants were primarily associated with metabolism and signal transduction. Notably, the majority of cold-responsive genes were associated with cell wall functions, including those related to cell wall organization, chitin catabolic processes, and oxidoreductases. Our findings suggest that SUB4 negatively regulates the cold response in rice seedlings, possibly by modifying the properties of the cell wall.

Comparative transcriptome analysis and meta-QTLs mapping reveal the regulatory mechanism of cold tolerance in rice at the budding stage

Rice (Oryza sativa L.) is one of the most extensively farmed food crops, but its development and productivity are significantly impacted by cold stress during the budding period. In this study, transcriptome sequencing was conducted on two types of rice: the cold-sensitive indica rice A117 and the substantially cold-tolerant japonica rice B106 under control and cold treatments. Differentially expressed genes between the two materials under cold conditions were analyzed using GO and KEGG enrichment analyses. The results revealed that processes such as the TCA cycle, glycolysis/glycogenesis, oxidative phosphorylation, and glutathione metabolism contribute to B106's cold tolerance. Additionally, an enrichment analysis of cold-induced genes in each material and shared genes identified significant enrichment in pathways such as glutathione metabolism, phenylpropanoid biosynthesis, and photosynthesis-antenna proteins. Initial cold tolerance QTLs at the rice bud stage were collected from published literature, and meta-QTL mapping identified 9 MQTLs. Gene expression profiling led to the identification of 75 potential DEGs within the 9 MQTLs region, from which four candidate genes (Os02g0194100, Os03g0802500, Os05g0129000, and Os07g0462000) were selected using qRT-PCR and gene annotation. These findings provide genetic resources for further research on the molecular mechanisms underlying rice's response to cold stress during the bud stage.

The Molecular Mechanism of Cold-Stress Tolerance: Cold Responsive Genes and Their Mechanisms in Rice (Oryza sativa L.)

Rice (Oryza sativa L.) production is highly susceptible to temperature fluctuations, which can significantly reduce plant growth and development at different developmental stages, resulting in a dramatic loss of grain yield. Over the past century, substantial efforts have been undertaken to investigate the physiological, biochemical, and molecular mechanisms of cold stress tolerance in rice. This review aims to provide a comprehensive overview of the recent developments and trends in this field. We summarized the previous advancements and methodologies used for identifying cold-responsive genes and the molecular mechanisms of cold tolerance in rice. Integration of new technologies has significantly improved studies in this era, facilitating the identification of essential genes, QTLs, and molecular modules in rice. These findings have accelerated the molecular breeding of cold-resistant rice varieties. In addition, functional genomics, including the investigation of natural variations in alleles and artificially developed mutants, is emerging as an exciting new approach to investigating cold tolerance. Looking ahead, it is imperative for scientists to evaluate the collective impacts of these novel genes to develop rice cultivars resilient to global climate change.

Harnessing weedy rice as functional food and source of novel traits for crop improvement

A relative of cultivated rice (Oryza sativa L.), weedy or red rice (Oryza spp.) is currently recognized as the dominant weed, leading to a drastic loss of yield of cultivated rice due to its highly competitive abilities like producing more tillers, panicles, and biomass with better nutrient uptake. Due to its high nutritional value, antioxidant properties (anthocyanin and proanthocyanin), and nutrient absorption ability, weedy rice is gaining immense research attentions to understand its genetic constitution to augment future breeding strategies and to develop nutrition-rich functional foods. Consequently, this review focuses on the unique gene source of weedy rice to enhance the cultivated rice for its crucial features like water use efficiency, abiotic and biotic stress tolerance, early flowering, and the red pericarp of the seed. It explores the debating issues on the origin and evolution of weedy rice, including its high diversity, signalling aspects, quantitative trait loci (QTL) mapping under stress conditions, the intricacy of the mechanism in the expression of the gene flow, and ecological challenges of nutrient removal by weedy rice. This review may create a foundation for future researchers to understand the gene flow between cultivated crops and weedy traits and support an improved approach for the applicability of several models in predicting multiomics variables.

Advances and opportunities in unraveling cold-tolerance mechanisms in the world's primary staple food crops

Temperatures below or above optimal growth conditions are among the major stressors affecting productivity, end-use quality, and distribution of key staple crops including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays L.). Among temperature stresses, cold stress induces cellular changes that cause oxidative stress and slowdown metabolism, limit growth, and ultimately reduce crop productivity. Perception of cold stress by plant cells leads to the activation of cold-responsive transcription factors and downstream genes, which ultimately impart cold tolerance. The response triggered in crops to cold stress includes gene expression/suppression, the accumulation of sugars upon chilling, and signaling molecules, among others. Much of the information on the effects of cold stress on perception, signal transduction, gene expression, and plant metabolism are available in the model plant Arabidopsis but somewhat lacking in major crops. Hence, a complete understanding of the molecular mechanisms by which staple crops respond to cold stress remain largely unknown. Here, we make an effort to elaborate on the molecular mechanisms employed in response to low-temperature stress. We summarize the effects of cold stress on the growth and development of these crops, the mechanism of cold perception, and the role of various sensors and transducers in cold signaling. We discuss the progress in cold tolerance research at the genome, transcriptome, proteome, and metabolome levels and highlight how these findings provide opportunities for designing cold-tolerant crops for the future.

Genetic dissection of cold tolerance at the budding stage of rice in an indica-japonica recombination inbred line population

Rice is highly cold-sensitive, and thus, the promotion of cold resistance in buds is essential. In this study, we conducted a mapping analysis to identify quantitative trait loci (QTLs) associated with cold tolerance in buds. The analysis was performed using a recombinant inbred line (RIL) population consisting of 192 lines derived from the cold-tolerant strain 02428 and the cold-sensitive strain YZX. Seven additive loci on chromosomes 1, 4, 5, and 6 were identified, of which loci 3 and 7 were found in two crop seasons, indicating stability. Three epistatic interactions, one present over two seasons, were found. Loci 3 and 7 pyramided with two main-effect QTLs observed to control the rate of low-temperature germination in our previous study. Two materials with good cold resistance at the germination and bud stages were obtained, namely, G93 and G146. Transcriptome sequencing analysis of the two parent buds after cold treatment found that genes expressed differentially between the two parents were related to photosynthesis, energy metabolism, and reactive oxygen scavenging. Five candidate genes, namely, Os01g0385400, Os01g0388000, Os06g0287700, Os06g0289200, and Os06g0291100, were selected in the two stable intervals based on gene expression profiles and annotations. These genetic loci exhibit strong potential as targets for breeding cold tolerance in buds and require additional investigation. In conclusion, this work provides valuable genetic resources that can be utilized to improve the cold tolerance of rice.

Gene regulation and ionome homeostasis in rice plants in response to arsenite stress: potential connection between transcriptomics and ionomics

Ionomics and transcriptomics were applied to demonstrate response of rice to arsenite [As(III)] stress in the current study. Rice plants were cultured in nutrient solutions treated with 0, 100 and 500 μg/L As(III) coded as CK, As1 and As5, respectively. The rice ionomes exhibited discriminatory response to environmental disturbances. Solid evidence of the effects of As(III) stress on binding, transport or metabolism of P, K, Ca, Zn and Cu was obtained in this work. Differentially expressed genes (DEGs) in the shoots were identified in three datasets: As1 vs CK, As5 vs CK and As5 vs As1. DEGs identified simultaneously in two or three datasets were selected for subsequent interaction and enrichment analyses. Upregulation of genes involved in protein kinase activity, phosphorus metabolic process and phosphorylation were detected in the rice treated with As(III), resulting in the maintenance of P homeostasis in the shoots. Zn and Ca binding genes were up-regulated since excess As inhibited the translocation of Zn and Ca from roots to shoots. Increased expression of responsive genes including HMA, WRKY, NAC and PUB genes conferred As tolerance in the rice plants to cope with external As(III) stress. The results suggested that As(III) stress could disturb the uptake and translocation of macro and essential elements by rice. Plants could regulate the expression of corresponding genes to maintain mineral nutrient homeostasis for essential metabolic processes.

Genome-wide transcriptional profiling provides clues to molecular mechanisms underlying cold tolerance in chickpea

Chickpea is an important food legume cultivated in several countries. A sudden drop in autumn temperature, freezing winter temperature, and late spring cold events result in significant losses in chickpea production. The current study used RNA sequencing of two cold tolerant (Saral) and sensitive (ILC533) Kabuli chickpea genotypes to identify cold tolerance-associated genes/pathways. A total of 200.85 million raw reads were acquired from the leaf samples by Illumina sequencing, and around 86% of the clean reads (199 million) were mapped to the chickpea reference genome. The results indicated that 3710 (1980 up- and 1730 down-regulated) and 3473 (1972 up- and 1501 down-regulated) genes were expressed differentially under cold stress in the tolerant and sensitive genotypes, respectively. According to the GO enrichment analysis of uniquely down-regulated genes under cold stress in ILC533, photosynthetic membrane, photosystem II, chloroplast part, and photosystem processes were enriched, revealing that the photosynthesis is severely sensitive to cold stress in this sensitive genotype. Many remarkable transcription factors (CaDREB1E, CaMYB4, CaNAC47, CaTCP4, and CaWRKY33), signaling/regulatory genes (CaCDPK4, CaPP2C6, CaMKK2, and CaHSFA3), and protective genes (CaCOR47, CaLEA3, and CaGST) were identified among the cold-responsive genes of the tolerant genotype. These findings would help improve cold tolerance across chickpea genotypes by molecular breeding or genetic engineering.

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.

Transcriptomic and metabolomic insights on the molecular mechanisms of flower buds in responses to cold stress in two Camellia oleifera cultivars

INTRODUCTION

The Camellia oleifera (C. oleifera) cultivars 'Huashuo' (HS) and 'Huaxin' (HX) are new high-yielding and economically valuable cultivars that frequently encounter prolonged cold weather during the flowering period, resulting in decreased yields and quality. The flower buds of HS sometimes fail to open or open incompletely under cold stress, whereas the flower buds of HX exhibit delayed opening but the flowers and fruits rarely drop.

METHODS

In this study, flower buds at the same development stage of two C. oleifera cultivars were used as test materials for a combination of physiological, transcriptomic and metabolomic analyses, to unravel the different cold regulatory mechanisms between two cultivars of C. oleifera.

RESULTS AND DISCUSSION

Key differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) involved in sugar metabolism, phenylpropanoid biosynthesis, and hormone signal transduction were significantly higher in HX than in HS, which is consistent with phenotypic observations from a previous study. The results indicate that the flower buds of HX are less affected by long-term cold stress than those of HS, and that cold resistance in C. oleifera cultivars varies among tissues or organs.This study will provide a basis for molecular markers and molecular breeding of C. oleifera.

Integrated transcriptomics and metabolomics analysis reveals key regulatory network that response to cold stress in common Bean (Phaseolus vulgaris L.)

Cold temperatures can be detrimental to crop survival and productivity. Breeding progress can be improved by understanding the molecular basis of low temperature tolerance. We investigated the key routes and critical metabolites related to low temperature resistance in cold-tolerant and -sensitive common bean cultivars 120 and 093, respectively. Many potential genes and metabolites implicated in major metabolic pathways during the chilling stress response were identified through transcriptomics and metabolomics research. Under chilling stress, the expression of many genes involved in lipid, amino acid, and flavonoid metabolism, as well as metabolite accumulation increased in the two bean types. Malondialdehyde (MDA) content was lower in 120 than in 093. Regarding amino acid metabolism, 120 had a higher concentration of acidic amino acids than 093, whereas 093 had a higher concentration of basic amino acids. Methionine accumulation was clearly higher in 120 than in 093. In addition, 120 had a higher concentration of many types of flavonoids than 093. Flavonoids, methionine and malondialdehyde could be used as biomarkers of plant chilling injury. Transcriptome analysis of hormone metabolism revealed considerably greater, expression of abscisic acid (ABA), gibberellin (GA), and jasmonic acid (JA) in 093 than in 120 during chilling stress, indicating that hormone regulation modes in 093 and 120 were different. Thus, chilling stress tolerance is different between 093 and 120 possibly due to transcriptional and metabolic regulation.

Progress and Prospects of the Molecular Basis of Soybean Cold Tolerance

Cold stress is a major factor influencing the geographical distribution of soybean growth and causes immense losses in productivity. Understanding the molecular mechanisms that the soybean has undergone to survive cold temperatures will have immense value in improving soybean cold tolerance. This review focuses on the molecular mechanisms involved in soybean response to cold. We summarized the recent studies on soybean cold-tolerant quantitative trait loci (QTLs), transcription factors, associated cold-regulated (COR) genes, and the regulatory pathways in response to cold stress. Cold-tolerant QTLs were found to be overlapped with the genomic region of maturity loci of E1, E3, E4, pubescence color locus of T, stem growth habit gene locus of Dt1, and leaf shape locus of Ln, indicating that pleiotropic loci may control multiple traits, including cold tolerance. The C-repeat responsive element binding factors (CBFs) are evolutionarily conserved across species. The expression of most GmDREB1s was upregulated by cold stress and overexpression of GmDREB1B;1 in soybean protoplast, and transgenic Arabidopsis plants can increase the expression of genes with the DRE core motif in their promoter regions under cold stress. Other soybean cold-responsive regulators, such as GmMYBJ1, GmNEK1, GmZF1, GmbZIP, GmTCF1a, SCOF-1 and so on, enhance cold tolerance by regulating the expression of COR genes in transgenic Arabidopsis. CBF-dependent and CBF-independent pathways are cross-talking and work together to activate cold stress gene expression. Even though it requires further dissection for precise understanding, the function of soybean cold-responsive transcription factors and associated COR genes studied in Arabidopsis shed light on the molecular mechanism of cold responses in soybeans and other crops. Furthermore, the findings may also provide practical applications for breeding cold-tolerant soybean varieties in high-latitude and high-altitude regions.

Transcriptomic profiling of wheat stem during meiosis in response to freezing stress

Low temperature injury in spring has seriously destabilized the production and grain quality of common wheat. However, the molecular mechanisms underlying spring frost tolerance remain elusive. In this study, we investigated the response of a frost-tolerant wheat variety Zhongmai8444 to freezing stress at the meiotic stage. Transcriptome profiles over a time course were subsequently generated by high-throughput sequencing. Our results revealed that the prolonged freezing temperature led to the significant reductions in plant height and seed setting rate. Cell wall thickening in the vascular tissue was also observed in the stems. RNA-seq analyses demonstrated the identification of 1010 up-regulated and 230 down-regulated genes shared by all time points of freezing treatment. Enrichment analysis revealed that gene activity related to hormone signal transduction and cell wall biosynthesis was significantly modulated under freezing. In addition, among the identified differentially expressed genes, 111 transcription factors belonging to multiple gene families exhibited dynamic expression pattern. This study provided valuable gene resources beneficial for the breeding of wheat varieties with improved spring frost tolerance.

Global whole-genome comparison and analysis to classify subpopulations and identify resistance genes in weedy rice relevant for improving crops

Common weedy rice plants are important genetic resources for modern breeding programs because they are the closest relatives to rice cultivars and their genomes contain elite genes. Determining the utility and copy numbers of WRKY and nucleotide-binding site (NBS) resistance-related genes may help to clarify their variation patterns and lead to crop improvements. In this study, the weedy rice line LM8 was examined at the whole-genome level. To identify the Oryza sativa japonica subpopulation that LM8 belongs to, the single nucleotide polymorphisms (SNPs) of 180 cultivated and 23 weedy rice varieties were used to construct a phylogenetic tree and a principal component analysis and STRUCTURE analysis were performed. The results indicated that LM8 with admixture components from japonica (GJ) and indica (XI) belonged to GJ-admixture (GJ-adm), with more than 60% of its genetic background derived from XI-2 (22.98%), GJ-tropical (22.86%), and GJ-subtropical (17.76%). Less than 9% of its genetic background was introgressed from weedy rice. Our results also suggested LM8 may have originated in a subtropical or tropical geographic region. Moreover, the comparisons with Nipponbare (NIP) and Shuhui498 (R498) revealed many specific structure variations (SVs) in the LM8 genome and fewer SVs between LM8 and NIP than between LM8 and R498. Next, 96 WRKY and 464 NBS genes were identified and mapped on LM8 chromosomes to eliminate redundancies. Three WRKY genes (ORUFILM02g002693, ORUFILM05g002725, and ORUFILM05g001757) in group III and one RNL [including the resistance to powdery mildew 8 (RPW8) domain, NBS, and leucine rich repeats (LRRs)] type NBS gene (ORUFILM12g000772) were detected in LM8. Among the NBS genes, the RPW8 domain was detected only in ORUFILM12g000772. This gene may improve plant resistance to pathogens as previously reported. Its classification and potential utility imply LM8 should be considered as a germplasm resource relevant for rice breeding programs.

Identification of Rice Accessions Having Cold Tolerance at the Seedling Stage and Development of Novel Genotypic Assays for Predicting Cold Tolerance

Rice is susceptible to cold stress at the seedling stage, which can delay growth and decrease yield. We evaluated 187 rice accessions for cold tolerance at the seedling stage and developed genotypic assays for three markers. All japonica (20/20) and 20/140 indica accessions were highly cold tolerant. Two SNP markers specific for COLD1 and LOC_Os10g34840 were practical to use by normal agarose gel. The SNP marker specific for COLD1 was highly specific for predicting cold tolerance. However, the sensitivity of this marker was low as several cold-tolerant indica accessions lacked the cold-tolerant allele. The LOC_Os10g34840 marker was slightly more sensitive than the COLD1 marker for predicting highly cold-tolerant accessions. An insertion/deletion variant in the NAC6 gene was identified as a novel cold tolerance marker. The NAC6 marker predicted more highly cold-tolerant accessions compared with the other two markers. The SNP marker specific for LOC_Os10g34840 and the NAC6 marker were present in several tested subgroups, suggesting their wide effects and distribution. The three markers combined predicted the most highly cold-tolerant accessions, indicating that the marker combination is superior for applications such as marker-assisted breeding. The cold-tolerant accessions and the genotypic marker assays will be useful for future rice breeding.

ERF subfamily transcription factors and their function in plant responses to abiotic stresses

Ethylene Responsive Factor (ERF) subfamily comprise the largest number of proteins in the plant AP2/ERF superfamily, and have been most extensively studied on the biological functions. Members of this subfamily have been proven to regulate plant resistances to various abiotic stresses, such as drought, salinity, chilling and some other adversities. Under these stresses, ERFs are usually activated by mitogen-activated protein kinase induced phosphorylation or escape from ubiquitin-ligase enzymes, and then form complex with nucleic proteins before binding to cis-element in promoter regions of stress responsive genes. In this review, we will discuss the phylogenetic relationships among the ERF subfamily proteins, summarize molecular mechanism how the transcriptional activity of ERFs been regulated and how ERFs of different subgroup regulate the transcription of stress responsive genes, such as high-affinity K⁺ transporter gene PalHKT1;2, reactive oxygen species related genes LcLTP, LcPrx, and LcRP, flavonoids synthesis related genes FtF3H and LhMYBSPLATTER, etc. Though increasing researches demonstrate that ERFs are involved in various abiotic stresses, very few interact proteins and target genes of them have been comprehensively annotated. Hence, future research prospects are described on the mechanisms of how stress signals been transited to ERFs and how ERFs regulate the transcriptional expression of stress responsive genes.

Genome-wide transcriptome profiling revealed biological macromolecules respond to low temperature stress in Brassica napus L

Brassica napus L. (B. napus) is a vital oilseed crop cultivated worldwide; low temperature (LT) is one of the major stress factors that limit its growth, development, distribution, and production. Even though processes have been developed to characterize LT-responsive genes, only limited studies have exploited the molecular response mechanisms in B. napus. Here the transcriptome data of an elite B. napus variety with LT adaptability was acquired and applied to investigate the gene expression profiles of B. napus in response to LT stress. The bioinformatics study revealed a total of 79,061 unigenes, of which 3,703 genes were differentially expressed genes (DEGs), with 2,129 upregulated and 1,574 downregulated. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis pinpointed that the DEGs were enriched in LT-stress-responsive biological functions and metabolic pathways, which included sugar metabolism, antioxidant defense system, plant hormone signal transduction, and photosynthesis. Moreover, a group of LT-stress-responsive transcription factors with divergent expression patterns under LT was summarized. A combined protein interaction suggested that a complex interconnected regulatory network existed in all detected pathways. RNA-seq data was verified using real-time quantitative polymerase chain reaction analysis. Based on these findings, we presented a hypothesis model illustrating valuable information for understanding the LT response mechanisms in B. napus.

Combination of Genomics, Transcriptomics Identifies Candidate Loci Related to Cold Tolerance in Dongxiang Wild Rice

Rice, a cold-sensitive crop, is a staple food for more than 50% of the world’s population. Low temperature severely compromises the growth of rice and challenges China’s food safety. Dongxiang wild rice (DXWR) is the most northerly common wild rice in China and has strong cold tolerance, but the genetic basis of its cold tolerance is still unclear. Here, we report quantitative trait loci (QTLs) analysis for seedling cold tolerance (SCT) using a high-density single nucleotide polymorphism linkage map in the backcross recombinant inbred lines that were derived from a cross of DXWR, and an indica cultivar, GZX49. A total of 10 putative QTLs were identified for SCT under 4 °C cold treatment, each explaining 2.0–6.8% of the phenotypic variation in this population. Furthermore, transcriptome sequencing of DXWR seedlings before and after cold treatment was performed, and 898 and 3413 differentially expressed genes (DEGs) relative to 0 h in cold-tolerant for 4 h and 12 h were identified, respectively. Gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) analysis were performed on these DEGs. Using transcriptome data and genetic linkage analysis, combined with qRT-PCR, sequence comparison, and bioinformatics, LOC_Os08g04840 was putatively identified as a candidate gene for the major effect locus qSCT8. These findings provided insights into the genetic basis of SCT for the improvement of cold stress potential in rice breeding programs.

Identifying conserved genes involved in crop tolerance to cold stress

Low temperature is a limiting factor for crop productivity in tropical and subtropical climates. Cold stress response in plants involves perceiving and relaying the signal through a transcriptional cascade composed of different transduction components, resulting in altered gene activity. We performed a meta-analysis of four previously published datasets of cold-tolerant and cold-sensitive crops to better understand the gene regulatory networks and identify key genes involved in cold stress tolerance conserved across phylogenetically distant species. Re-analysing the raw data with the same bioinformatics pipeline, we identified common cold tolerance-related genes. We found 236 and 242 commonly regulated genes in sensitive and tolerant genotypes, respectively. Gene enrichment analysis showed that protein modifications, hormone metabolism, cell wall, and secondary metabolism are the most conserved pathways involved in cold tolerance. Upregulation of the abiotic stress (heat and drought/salt) related genes [heat shock N -terminal domain-containing protein, 15.7kDa class I-related small heat shock protein-like, DNAJ heat shock N -terminal domain-containing protein, and HYP1 (HYPOTHETICAL PROTEIN 1)] in sensitive genotypes and downregulation of the abiotic stress (heat and drought/salt) related genes (zinc ion binding and pollen Ole e 1 allergen and extensin family protein) in tolerant genotypes was observed across the species. Almost all development-related genes were upregulated in tolerant and downregulated in sensitive genotypes. Moreover, protein-protein network analysis identified highly interacting proteins linked to cold tolerance. Mapping of abiotic stress-related genes on analysed species genomes provided information that could be essential to developing molecular markers for breeding and building up genetic improvement strategies using CRISPR/Cas9 technologies.

Analysis of Thioredoxins and Glutaredoxins in Soybean: Evidence of Translational Regulation under Water Restriction

Soybean (Glycine max (L.) Merr.) establishes symbiosis with rhizobacteria, developing the symbiotic nodule, where the biological nitrogen fixation (BNF) occurs. The redox control is key for guaranteeing the establishment and correct function of the BNF process. Plants have many antioxidative systems involved in ROS homeostasis and signaling, among them a network of thio- and glutaredoxins. Our group is particularly interested in studying the differential response of nodulated soybean plants to water-deficit stress. To shed light on this phenomenon, we set up an RNA-seq experiment (for total and polysome-associated mRNAs) with soybean roots comprising combined treatments including the hydric and the nodulation condition. Moreover, we performed the initial identification and description of the complete repertoire of thioredoxins (Trx) and glutaredoxins (Grx) in soybean. We found that water deficit altered the expression of a greater number of differentially expressed genes (DEGs) than the condition of plant nodulation. Among them, we identified 12 thioredoxin (Trx) and 12 glutaredoxin (Grx) DEGs, which represented a significant fraction of the detected GmTrx and GmGrx in our RNA-seq data. Moreover, we identified an enriched network in which a GmTrx and a GmGrx interacted with each other and associated through several types of interactions with nitrogen metabolism enzymes.

Transcriptome-based gene regulatory network analyses of differential cold tolerance of two tobacco cultivars

BACKGROUND

Cold is one of the main abiotic stresses that severely affect plant growth and development, and crop productivity as well. Transcriptional changes during cold stress have already been intensively studied in various plant species. However, the gene networks involved in the regulation of differential cold tolerance between tobacco varieties with contrasting cold resistance are quite limited.

RESULTS

Here, we conducted multiple time-point transcriptomic analyses using Tai tobacco (TT, cold susceptibility) and Yan tobacco (YT, cold resistance) with contrasting cold responses. We identified similar DEGs in both cultivars after comparing with the corresponding control (without cold treatment), which were mainly involved in response to abiotic stimuli, metabolic processes, kinase activities. Through comparison of the two cultivars at each time point, in contrast to TT, YT had higher expression levels of the genes responsible for environmental stresses. By applying Weighted Gene Co-Expression Network Analysis (WGCNA), we identified two main modules: the pink module was similar while the brown module was distinct between the two cultivars. Moreover, we obtained 100 hub genes, including 11 important transcription factors (TFs) potentially involved in cold stress, 3 key TFs in the brown module and 8 key TFs in the pink module. More importantly, according to the genetic regulatory networks (GRNs) between TFs and other genes or TFs by using GENIE3, we identified 3 TFs (ABI3/VP1, ARR-B and WRKY) mainly functioning in differential cold responses between two cultivars, and 3 key TFs (GRAS, AP2-EREBP and C2H2) primarily involved in cold responses.

CONCLUSION

Collectively, our study provides valuable resources for transcriptome- based gene network studies of cold responses in tobacco. It helps to reveal how key cold responsive TFs or other genes are regulated through network. It also helps to identify the potential key cold responsive genes for the genetic manipulation of tobacco cultivars with enhanced cold tolerance in the future.

OsWRKY115 on qCT7 links to cold tolerance in rice

qCT7, a novel QTL for increasing seedling cold tolerance in rice, was fine-mapped to a 70.9-kb region on chromosome 7, and key OsWRKY115 was identified in transgenic plants. Cold stress caused by underground cold-water irrigation seriously limits rice productivity. We systemically measured the cold-responsive traits of 2,570 F2 individuals derived from two widely cultivated rice cultivars, Kong-Yu-131 and Dong-Nong-422, to identify the major genomic regions associated with cold tolerance. A novel major QTL, qCT7, was mapped on chromosome 7 associated with the cold tolerance and survival, using whole-genome re-sequencing with bulked segregant analysis. Local QTL linkage analysis with F₂ and fine mapping with recombinant plant revealed a 70.9-kb core region on qCT7 encoding 13 protein-coding genes. Only the LOC_Os07g27670 expression level encoding the OsWRKY115 transcription factor on the locus was specifically induced by cold stress in the cold-tolerant cultivar. Moreover, haplotype analysis and the KASP8 marker indicated that OsWRKY115 was significantly associated with cold tolerance. Overexpression and knockout of OsWRKY115 significantly affected cold tolerance in seedlings. Our experiments identified OsWRKY115 as a novel regulatory gene associated with cold response in rice, and the Kong-Yu-131 allele with specific cold-induced expression may be an important molecular variant.

Low Temperature Stress Tolerance: An Insight Into the Omics Approaches for Legume Crops

The change in climatic conditions is the major cause for decline in crop production worldwide. Decreasing crop productivity will further lead to increase in global hunger rate. Climate change results in environmental stress which has negative impact on plant-like deficiencies in growth, crop yield, permanent damage, or death if the plant remains in the stress conditions for prolonged period. Cold stress is one of the main abiotic stresses which have already affected the global crop production. Cold stress adversely affects the plants leading to necrosis, chlorosis, and growth retardation. Various physiological, biochemical, and molecular responses under cold stress have revealed that the cold resistance is more complex than perceived which involves multiple pathways. Like other crops, legumes are also affected by cold stress and therefore, an effective technique to mitigate cold-mediated damage is critical for long-term legume production. Earlier, crop improvement for any stress was challenging for scientific community as conventional breeding approaches like inter-specific or inter-generic hybridization had limited success in crop improvement. The availability of genome sequence, transcriptome, and proteome data provides in-depth sight into different complex mechanisms under cold stress. Identification of QTLs, genes, and proteins responsible for cold stress tolerance will help in improving or developing stress-tolerant legume crop. Cold stress can alter gene expression which further leads to increases in stress protecting metabolites to cope up the plant against the temperature fluctuations. Moreover, genetic engineering can help in development of new cold stress-tolerant varieties of legume crop. This paper provides a general insight into the "omics" approaches for cold stress in legume crops.

Comparing Early Transcriptomic Responses of 18 Soybean (Glycine max) Genotypes to Iron Stress

Iron deficiency chlorosis (IDC) is an abiotic stress that negatively affects soybean (Glycine max [L.] Merr.) production. Much of our knowledge of IDC stress responses is derived from model plant species. Gene expression, quantitative trait loci (QTL) mapping, and genome-wide association studies (GWAS) performed in soybean suggest that stress response differences exist between model and crop species. Our current understanding of the molecular response to IDC in soybeans is largely derived from gene expression studies using near-isogenic lines differing in iron efficiency. To improve iron efficiency in soybeans and other crops, we need to expand gene expression studies to include the diversity present in germplasm collections. Therefore, we collected 216 purified RNA samples (18 genotypes, two tissue types [leaves and roots], two iron treatments [sufficient and deficient], three replicates) and used RNA sequencing to examine the expression differences of 18 diverse soybean genotypes in response to iron deficiency. We found a rapid response to iron deficiency across genotypes, most responding within 60 min of stress. There was little evidence of an overlap of specific differentially expressed genes, and comparisons of gene ontology terms and transcription factor families suggest the utilization of different pathways in the stress response. These initial findings suggest an untapped genetic potential within the soybean germplasm collection that could be used for the continued improvement of iron efficiency in soybean.

Transcriptome profiling of Malus sieversii under freezing stress after being cold-acclimated

Background

Freezing temperatures are an abiotic stress that has a serious impact on plant growth and development in temperate regions and even threatens plant survival. The wild apple tree (Malus sieversii) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in wild apple trees.

Results

Phytohormone and physiology profiles and transcriptome analysis were used to elaborate on the dynamic response mechanism. We determined that JA, IAA, and ABA accumulated in the cold acclimation stage and decreased during freezing stress in response to freezing stress. To elucidate the molecular mechanisms of freezing stress after cold acclimation, we employed single molecular real-time (SMRT) and RNA-seq technologies to study genome-wide expression profiles in wild apple. Using the PacBio and Illumina platform, we obtained 20.79G subreads. These reads were assembled into 61,908 transcripts, and 24,716 differentially expressed transcripts were obtained. Among them, 4410 transcripts were differentially expressed during the whole process of freezing stress, and these were examined for enrichment via GO and KEGG analyses. Pathway analysis indicated that “plant hormone signal transduction”, “starch and sucrose metabolism”, “peroxisome” and “photosynthesis” might play a vital role in wild apple responses to freezing stress. Furthermore, the transcription factors DREB1/CBF, MYC2, WRKY70, WRKY71, MYB4 and MYB88 were strongly induced during the whole stress period.

Conclusions

Our study presents a global survey of the transcriptome profiles of wild apple trees in dynamic response to freezing stress after two days cold acclimation and provides insights into the molecular mechanisms of freezing adaptation of wild apple plants for the first time. The study also provides valuable information for further research on the antifreezing reaction mechanism and genetic improvement of M. sieversii after cold acclimation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07998-0.

Tissue-level transcriptomic responses to local and distal chilling reveal potential chilling survival mechanisms in maize

Chilling is a major stress to plants of subtropical and tropical origins including maize (Zea mays L.). To reveal molecular mechanisms underlying chilling tolerance and survival, we investigated transcriptomic responses to chilling stress in differentiated leaves and roots as well as in crowns with meristem activity in maize. Chilling stress on shoots and roots is found to each contributes to seedling lethality in maize. Comparison of maize lines with different chilling tolerance capacities reveals that chilling survival is highly associated with upregulation of abscisic acid biosynthesis and response as well as transcriptional regulators in leaves and crowns. It is also associated with the downregulation of translation in leaves and heat response in crowns. Chilling treatment on whole or part of the plants reveals that response to distal-chilling is very distinct from, and sometimes opposite to, response to local- or whole-plant chilling in both leaves and roots, suggesting a communication between shoots and roots in environmental response. This study thus provides transcriptomic responses in leaves, roots and crowns under differential chilling stresses in maize and reveals potential chilling tolerance and survival mechanisms which lays ground for improving chilling tolerance in crop plants.

Transcriptome analysis of Sonneratia caseolaris seedlings under chilling stress

Sonneratia caseolaris is a native mangrove species found in China. It is fast growing and highly adaptable for mangrove afforestation, but suffered great damage by chilling event once introduced to high latitude area. To understand the response mechanisms under chilling stress, physiological and transcriptomic analyses were conducted. The relative electrolyte conductivity, malondialdehyde (MDA) content, soluble sugar content and soluble protein content increased significantly under chilling stress. This indicated that S. caseolaris suffered great damage and increased the levels of osmoprotectants in response to the chilling stress. Gene expression comparison analysis of S. caseolaris leaves after 6 h of chilling stress was performed at the transcriptional scale using RNA-Seq. A total of 168,473 unigenes and 3,706 differentially expressed genes (DEGs) were identified. GO and KEGG enrichment analyses showed that the DEGs were mainly involved in carbohydrate metabolism, antioxidant enzyme, plant hormone signal transduction, and transcription factors (TFs). Sixteen genes associated with carbohydrate metabolism, antioxidant enzyme, phytohormones and TFs were selected for qRT-PCR verification, and they indicated that the transcriptome data were reliable. Our work provided a comprehensive review of the chilling response of S. caseolaris at both physiological and transcriptomic levels, which will prove useful for further studies on stress-responses in mangrove plants.

Comparative transcriptomic and physiological analyses of weedy rice and cultivated rice to identify vital differentially expressed genes and pathways regulating the ABA response

Weedy rice is a valuable germplasm resource characterized by its high tolerance to both abiotic and biotic stresses. Abscisic acid (ABA) serves as a regulatory signal in plant cells as part of their adaptive response to stress. However, a global understanding of the response of weedy rice to ABA remains to be elucidated. In the present study, the sensitivity to ABA of weedy rice (WR04-6) was compared with that of temperate japonica Shennong9816 (SN9816) in terms of seed germination and post-germination growth via the application of exogenous ABA and diniconazole, an inhibitor of ABA catabolism. Physiological analysis and a transcriptomic comparison allowed elucidation of the molecular and physiological mechanisms associated with continuous ABA and diniconazole treatment. WR04-6 was found to display higher ABA sensitivity than SN9816, resulting in the rapid promotion of antioxidant enzyme activity. Comparative transcriptomic analyses indicated that the number of differentially expressed genes (DEGs) in WR04-6 seedlings treated with 2 μM ABA or 10 μM diniconazole was greater than that in SN9816 seedlings. Genes involved in stress defense, hormone signal transduction, and glycolytic and citrate cycle pathways were highly expressed in WR04-6 in response to ABA and diniconazole. These findings provide new insight into key processes mediating the ABA response between weedy and cultivated rice.

Comparative transcriptome profile analysis of rice varieties with different tolerance to zinc deficiency

Zinc (Zn) is an indispensable element for rice growth. Zn deficiency results in brown blotches and streaks 2-3 weeks after transplanting, as well as stunting, reduced tillering, and low productivity of rice plants. These processes are controlled by different families of expressed genes. A comparative transcriptome profile analysis was conducted using the roots of two Zn deficiency tolerant varieties (UCP122 and KALIBORO26) and two sensitive varieties (IR26 and IR64) by merging data from untreated control (CK) and Zn deficiency treated samples. Results revealed a total of 4,688 differentially expressed genes (DEGs) between the normal Zn and deficient conditions, with 2,702 and 1,489 unique DEGs upregulated and downregulated, respectively. Functional enrichment analysis identified transcription factors (TFs), such as WRKY, MYB, ERF, and bHLH which are important in the regulation of the Zn deficiency response. Furthermore, chitinases, jasmonic acid, and phenylpropanoid pathways were found to be important in the Zn deficiency response. The metal tolerance protein (MTP) genes also appeared to play an important role in conferring tolerance to Zn deficiency. A heavy metal-associated domain-containing protein 7 was associated with tolerance to Zn deficiency and negatively regulated downstream genes. Collectively, our findings provide valuable expression patterns and candidate genes for the study of molecular mechanisms underlying the response to Zn deficiency and for improvements in breeding for tolerance to Zn deficiency in rice.

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.

Transcriptomic Profiling of Young Cotyledons Response to Chilling Stress in Two Contrasting Cotton (Gossypium hirsutum L.) Genotypes at the Seedling Stage

Young cotyledons of cotton seedlings are most susceptible to chilling stress. To gain insight into the potential mechanism of cold tolerance of young cotton cotyledons, we conducted physiological and comparative transcriptome analysis of two varieties with contrasting phenotypes. The evaluation of chilling injury of young cotyledons among 74 cotton varieties revealed that H559 was the most tolerant and YM21 was the most sensitive. The physiological analysis found that the ROS scavenging ability was lower, and cell membrane damage was more severe in the cotyledons of YM21 than that of H559 under chilling stress. RNA-seq analysis identified a total of 44,998 expressed genes and 19,982 differentially expressed genes (DEGs) in young cotyledons of the two varieties under chilling stress. Weighted gene coexpression network analysis (WGCNA) of all DEGs revealed four significant modules with close correlation with specific samples. The GO-term enrichment analysis found that lots of genes in H559-specific modules were involved in plant resistance to abiotic stress. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that pathways such as plant hormone signal transduction, MAPK signaling, and plant–pathogen interaction were related to chilling stress response. A total of 574 transcription factors and 936 hub genes in these modules were identified. Twenty hub genes were selected for qRT-PCR verification, revealing the reliability and accuracy of transcriptome data. These findings will lay a foundation for future research on the molecular mechanism of cold tolerance in cotyledons of cotton.

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.

Understanding the early cold response mechanism in IR64 indica rice variety through comparative transcriptome analysis

Background

Cellular reprogramming in response to environmental stress involves alteration of gene expression, changes in the protein and metabolite profile for ensuring better stress management in plants. Similar to other plant species originating in tropical and sub-tropical areas, indica rice is highly sensitive to low temperature that adversely affects its growth and grain productivity. Substantial work has been done to understand cold induced changes in gene expression in rice plants. However, adequate information is not available for early gene expression, especially in indica variety. Therefore, a transcriptome profile was generated for cold shock treated seedlings of IR64 variety to identify early responsive genes.

Results

The functional annotation of early DEGs shows enrichment of genes involved in altered membrane rigidity and electrolytic leakage, the onset of calcium signaling, ROS generation and activation of stress responsive transcription factors in IR64. Gene regulatory network suggests that cold shock induced Ca2+ signaling activates DREB/CBF pathway and other groups of transcription factors such as MYB, NAC and ZFP; for activating various cold-responsive genes. The analysis also indicates that cold induced signaling proteins like RLKs, RLCKs, CDPKs and MAPKK and ROS signaling proteins. Further, several late-embryogenesis-abundant (LEA), dehydrins and low temperature-induced-genes were upregulated under early cold shock condition, indicating the onset of water-deficit conditions. Expression profiling in different high yielding cultivars shows high expression of cold-responsive genes in Heera and CB1 indica varieties. These varieties show low levels of cold induced ROS production, electrolytic leakage and high germination rate post-cold stress, compared to IR36 and IR64. Collectively, these results suggest that these varieties may have improved adaptability to cold stress.

Conclusions

The results of this study provide insights about early responsive events in Oryza sativa l.ssp. indica cv IR64 in response to cold stress. Our data shows the onset of cold response is associated with upregulation of stress responsive TFs, hydrophilic proteins and signaling molecules, whereas, the genes coding for cellular biosynthetic enzymes, cell cycle control and growth-related TFs are downregulated. This study reports that the generation of ROS is integral to the early response to trigger the ROS mediated signaling events during later stages.

Comparative Transcriptome-Based Mining and Expression Profiling of Transcription Factors Related to Cold Tolerance in Peanut

Plants tolerate cold stress by regulating gene networks controlling cellular and physiological traits to modify growth and development. Transcription factor (TF)-directed regulation of transcription within these gene networks is key to eliciting appropriate responses. Identifying TFs related to cold tolerance contributes to cold-tolerant crop breeding. In this study, a comparative transcriptome analysis was carried out to investigate global gene expression of entire TFs in two peanut varieties with different cold-tolerant abilities. A total of 87 TF families including 2328 TF genes were identified. Among them, 445 TF genes were significantly differentially expressed in two peanut varieties under cold stress. The TF families represented by the largest numbers of differentially expressed members were bHLH (basic helix—loop—helix protein), C2H2 (Cys2/His2 zinc finger protein), ERF (ethylene-responsive factor), MYB (v-myb avian myeloblastosis viral oncogene homolog), NAC (NAM, ATAF1/2, CUC2) and WRKY TFs. Phylogenetic evolutionary analysis, temporal expression profiling, protein–protein interaction (PPI) network, and functional enrichment of differentially expressed TFs revealed the importance of plant hormone signal transduction and plant-pathogen interaction pathways and their possible mechanism in peanut cold tolerance. This study contributes to a better understanding of the complex mechanism of TFs in response to cold stress in peanut and provides valuable resources for the investigation of evolutionary history and biological functions of peanut TFs genes involved in cold tolerance.

Comprehensive Evaluation and Analysis of the Mechanism of Cold Tolerance Based on the Transcriptome of Weedy Rice Seedlings

In this study, the cold-tolerance capacity of 133 varieties of weedy rice was evaluated based on the comprehensive evaluation index D, with Kongyu 131 used as a cold-tolerant control. A total of 39.8% of the 133 varieties were considered 'strong', indicating that weedy rice populations indeed have relatively strong cold-tolerance capacity as a whole, and the robust cold-tolerant varieties WR29 and WR157 were identified. Regression analysis showed that the metrics including the nitrogen recovery index, superoxide dismutase (SOD) content and malondialdehyde (MDA) content correlated significantly (P < 0.05) with cold tolerance and could be used as indicators of cold tolerance. On the basis of a transcriptome analysis of WR157, a robust cold-tolerant variety identified in this study, a total of 4645 putative DEGs were identified in treated groups compared to the control groups, with 2123 upregulated DEGs and 2522 downregulated DEGs. All upregulated DEGs were enriched on 1388 terms, all downregulated DEGs were enriched on 1566 terms; 911 of the 2123 upregulated DEGs fell into 98 KEGG categories and 1103 of the 2522 downregulated DEGs were in 115 categories. Further analysis showed that GO:0019740 and GO:0006808 are involved in nitrogen utilization; GO:0009269 and GO:0009414 are related to the stress response; and GO:0016491 and GO:0016614 are related to oxidoreductase activity. BACKGROUND: Weedy rice (Oryza) is a related pest species of cultivated rice (Oryza sativa L.) that has strong abiotic stress resistance; however, the comprehensive mechanism governing its cold tolerance is poorly understood. CONCLUSION: Our comprehensive evaluation based on five morphological indices and nine physiological indicators revealed outstanding levels of cold-tolerance capacity among weedy rice varieties from different regions and revealed some terms related to cold tolerance via transcriptome analysis. Our results underscored the reliable evaluation methods for additional cold tolerance studies and revealed several genes related to cold tolerance, which will help researchers breed cultivated rice varieties to increase their cold-tolerance capacity. These traits have the ability to increase seedling survival rate and growth, as well as future yields.

Genome-wide analysis of the bHLH gene family in Chinese jujube (Ziziphus jujuba Mill.) and wild jujube

Background

The bHLH (basic helix-loop-helix) transcription factor is one of the largest families of transcription factors in plants, containing a large number of members with diverse functions. Chinese jujube (Ziziphus jujuba Mill.) is the species with the highest economic value in the family Rhamnaceae. However, the characteristics of the bHLH family in the jujube genome are still unclear. Hence, ZjbHLHs were first searched at a genome-wide level, their expression levels under various conditions were investigated systematically, and their protein-protein interaction networks were predicted.

Results

We identified 92 ZjbHLHs in the jujube genome, and these genes were classified into 16 classes according to bHLH domains. Ten ZjbHLHs with atypical bHLH domains were found. Seventy ZjbHLHs were mapped to but not evenly distributed on 12 pseudo- chromosomes. The domain sequences among ZjbHLHs were highly conserved, and their conserved residues were also identified. The tissue-specific expression of 37 ZjbHLH genes in jujube and wild jujube showed diverse patterns, revealing that these genes likely perform multiple functions. Many ZjbHLH genes were screened and found to be involved in flower and fruit development, especially in earlier developmental stages. A few genes responsive to phytoplasma invasion were also verified. Based on protein-protein interaction prediction and homology comparison, protein-protein interaction networks composed of 92 ZjbHLHs were also established.

Conclusions

This study provides a comprehensive bioinformatics analysis of 92 identified ZjbHLH genes. We explored their expression patterns in various tissues, the flowering process, and fruit ripening and under phytoplasma stress. The protein-protein interaction networks of ZjbHLHs provide valuable clues toward further studies of their biological functions.

Electronic supplementary material

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