De novo transcriptome sequencing and gene expression profiling of Elymus nutans under cold stress

Tissue-specific transcriptomic analysis reveals the molecular mechanisms responsive to cold stress in Poa crymophila, and development of EST-SSR markers linked to cold tolerance candidate genes

BACKGROUND

Poa crymophila is a perennial, cold-tolerant, native grass species, widely distributed in the Qinghai-Tibet Plateau. However, the tissue-specific regulatory mechanisms and key regulatory genes underlying its cold tolerance remain poorly characterized. Therefore, in this study, based on the screening and evaluation of cold tolerance of four Poa species, the cold tolerance mechanism of P. crymophila's roots, stems, and leaves and its cold tolerance candidate genes were investigated through physiological and transcriptomic analyses.

RESULTS

Results of the present study suggested that the cold tolerance of the four Poa species was in the following order: P. crymophila > P. botryoides > P. pratensis var. anceps > P. pratensis. Cold stress significantly changed the physiological characteristics of roots, stems, and leaves of P. crymophila in this study. In addition, the transcriptome results showed that 4434, 8793, and 14,942 differentially expressed genes (DEGs) were identified in roots, stems, and leaves, respectively; however, 464 DEGs were commonly identified in these three tissues. KEGG enrichment analysis showed that these DEGs were mainly enriched in the phenylpropanoid biosynthesis pathway (roots), photosynthesis pathway (stems and leaves), circadian rhythm-plant pathway (stems and leaves), starch and sucrose metabolism pathway (roots, stems, and leaves), and galactose metabolism pathway (roots, stems, and leaves). A total of 392 candidate genes involved in Ca2+ signaling, ROS scavenging system, hormones, circadian clock, photosynthesis, and transcription factors (TFs) were identified in P. crymophila. Weighted gene co-expression network analysis (WGCNA) identified nine hub genes that may be involved in P. crymophila cold response. A total of 200 candidate gene-based EST-SSRs were developed and characterized. Twenty-nine polymorphic EST-SSRs primers were finally used to study genetic diversity of 40 individuals from four Poa species with different cold tolerance characteristics. UPGMA cluster and STRUCTURE analysis showed that the 40 Poa individuals were clustered into three major groups, individual plant with similar cold tolerance tended to group together. Notably, markers P37 (PcGA2ox3) and P148 (PcERF013) could distinguish P. crymophila from P. pratensis var. anceps, P. pratensis, and P. botryoides.

CONCLUSIONS

This study provides new insights into the molecular mechanisms underlying the cold tolerance of P. crymophila, and also lays a foundation for molecular marker-assisted selection for cold tolerance improvement in Poa species.

Effects of appropriate low-temperature treatment on the yield and quality of pigmented potato (Solanum tuberosum L.) tubers

Temperature is one of the important environmental factors affecting plant growth, yield and quality. Moreover, appropriately low temperature is also beneficial for tuber coloration. The red potato variety Jianchuanhong, whose tuber color is susceptible to temperature, and the purple potato variety Huaxinyangyu, whose tuber color is stable, were used as experimental materials and subjected to 20 °C (control check), 15 °C and 10 °C treatments during the whole growth period. The effects of temperature treatment on the phenotype, the expression levels of structural genes related to anthocyanins and the correlations of each indicator were analyzed. The results showed that treatment at 10 °C significantly inhibited the potato plant height, and the chlorophyll content and photosynthetic parameters in the leaves were reduced, and the enzyme activities of SOD and POD were significantly increased, all indicating that the leaves were damaged. Treatment at 10 °C also affected the tuberization of Huaxinyangyu and reduced the tuberization and coloring of Jianchuanhong, while treatment at 15 °C significantly increased the stem diameter, root-to-shoot ratio, yield and content of secondary metabolites, especially anthocyanins. Similarly, the expression of structural genes were enhanced in two pigmented potatoes under low-temperature treatment conditions. In short, proper low temperature can not only increase yield but also enhance secondary metabolites production. Previous studies have not focused on the effects of appropriate low-temperature treatment during the whole growth period of potato on the changes in metabolites during tuber growth and development, these results can provide a theoretical basis and technical guidance for the selection of pigmented potatoes with better nutritional quality planting environment and the formulation of cultivation measures.

Metabolomics combined with physiology and transcriptomics reveal how Nicotiana tabacum leaves respond to cold stress

Low temperature-induced cold stress is a major threat to plant growth, development and distribution. Unraveling the responses of temperature-sensitive crops to cold stress and the mechanisms of cold acclimation are critical for food demand. In this study, combined physiological, transcriptomic, and metabolomic analyses were conducted on Nicotiana tabacum suffering short-term 4 °C cold stress. Our results showed that cold stress destroyed cellular membrane stability, decreased the chlorophyll (Chl) and carotenoid contents, and closed stomata, resulting in lipid peroxidation and photosynthesis restriction. Chl fluorescence measurements revealed that primary photochemistry, photoelectrochemical quenching and photosynthetic electron transport in Nicotiana tabacum leaves were seriously suppressed upon exposer to cold stress. Enzymatic and nonenzymatic antioxidants, including superoxide dismutase, catalase, peroxidase, reduced glutathione, proline, and soluble sugar, were all profoundly increased to trigger the cold acclimation defense against oxidative damage. A total of 178 metabolites and 16,204 genes were differentially expressed in cold-stressed Nicotiana tabacum leaves. MEturquoise and MEblue modules identified by WGCNA were highly correlated with physiological indices, and the corresponding hub genes were significantly enriched in pathways related to photosynthesis - antenna proteins and flavonoid biosynthesis. Untargeted metabolomic analysis identified specific metabolites, including sucrose, phenylalanine, glutamine, glutamate, and proline, that enhance plant cold acclimation. Combined transcriptomics and metabolomic analysis highlight the vital roles of carbohydrate and amino acid metabolism in enhancing the cold tolerance of Nicotiana tabacum. Our comprehensive investigation provides novel insights for efforts to alleviate low temperature-induced oxidative damage to Nicotiana tabacum plants and proposes a breeding target for cold stress-tolerant cultivars.

The Overexpression of Peanut (Arachis hypogaea L.) AhALDH2B6 in Soybean Enhances Cold Resistance

Soybeans are the main source of oils and protein for humans and animals; however, cold stress jeopardizes their growth and limits the soybean planting area. Aldehyde dehydrogenases (ALDH) are conserved enzymes that catalyze aldehyde oxidation for detoxification in response to stress. Additionally, transgenic breeding is an efficient method for producing stress-resistant germplasms. In this study, the peanut ALDH gene AhALDH2B6 was heterologously expressed in soybean, and its function was tested. We performed RNA-seq using transgenic and wild-type soybeans with and without cold treatment to investigate the potential mechanism. Transgenic soybeans developed stronger cold tolerance, with longer roots and taller stems than P3 soybeans. Biochemically, the transgenic soybeans exhibited a decrease in malondialdehyde activity and an increase in peroxidase and catalase content, both of which are indicative of stress alleviation. They also possessed higher levels of ALDH enzyme activity. Two phenylpropanoid-related pathways were specifically enriched in up-regulated differentially expressed genes (DEGs), including the phenylpropanoid metabolic process and phenylpropanoid biosynthetic process. Our findings suggest that AhALDH2B6 specifically up-regulates genes involved in oxidoreductase-related functions such as peroxidase, oxidoreductase, monooxygenase, and antioxidant activity, which is partially consistent with our biochemical data. These findings established the function of AhALDH2B6, especially its role in cold stress processes, and provided a foundation for molecular plant breeding, especially plant-stress-resistance breeding.

Integrated Transcriptomics and Metabolomics Analysis of the Fructan Metabolism Response to Low-Temperature Stress in Garlic

As the main reserve carbohydrate in garlic, fructan contributes to garlic's yield and quality formation. Numerous studies have shown that plant fructan metabolism induces a stress response to adverse environments. However, the transcriptional regulation mechanism of garlic fructan in low-temperature environments is still unknown. In this study, the fructan metabolism of garlic seedlings under low-temperature stress was revealed by transcriptome and metabolome approaches. With the extension of stress time, the number of differentially expressed genes and metabolites increased. Using weighted gene co-expression network analysis (WGCNA), three key enzyme genes related to fructan metabolism were screened (a total of 12 transcripts): sucrose: sucrose 1-fructosyltransferase (1-SST) gene; fructan: fructan 6G fructosyltransferase (6G-FFT) gene; and fructan 1-exohydrolase (1-FEH) gene. Finally, two hub genes were obtained, namely Cluster-4573.161559 (6G-FFT) and Cluster-4573.153574 (1-FEH). The correlation network and metabolic heat map analysis between fructan genes and carbohydrate metabolites indicate that the expression of key enzyme genes in fructan metabolism plays a positive promoting role in the fructan response to low temperatures in garlic. The number of genes associated with the key enzyme of fructan metabolism in trehalose 6-phosphate was the highest, and the accumulation of trehalose 6-phosphate content may mainly depend on the key enzyme genes of fructan metabolism rather than the enzyme genes in its own synthesis pathway. This study not only obtained the key genes of fructan metabolism in garlic seedlings responding to low temperatures but also preliminarily analyzed its regulatory mechanism, providing an important theoretical basis for further elucidating the cold resistance mechanism of garlic fructan metabolism.

Metabolic profiles in the xylem sap of Brassica juncea exposed to cadmium

Metabolic profiles in xylem sap are considered a fundamental mechanism for Cadmium (Cd) detoxification in plants. However, the metabolic mechanism of Brassica juncea xylem sap in response to Cd is still unclear. Here, we investigated the effects on the metabolomics of B. juncea xylem sap treated with Cd at different times by utilizing a nontargeted liquid chromatography-mass spectrometry (LC-MS)-based metabolomics method for further elucidating the response mechanism of Cd exposure. The findings indicated that 48 h and 7 days Cd exposure caused significant differences in metabolic profiles of the B. juncea xylem sap. Those differential metabolites are primarily involved in amino acids, organic acids, lipids, and carbohydrates, and most of them were downregulated, which played essential roles in response to Cd stress. Furthermore, B. juncea xylem sap resisted 48-h Cd exposure via regulation of glycerophospholipid metabolism, carbon metabolism, aminoacyl-tRNA biosynthesis, glyoxylate and dicarboxylate metabolism, linoleic acid metabolism, C5-branched dibasic acid metabolism, alpha-linolenic acid metabolism, cyanoamino acid metabolism, ABC transporters, biosynthesis of amino acids, and pyrimidine metabolism; whereas alpha-linolenic acid metabolism, glycerophospholipid metabolism, photosynthesis, and oxidative phosphorylation were regulated for resisting 7-day Cd exposure.

Improving abiotic stress tolerance of forage grasses - prospects of using genome editing

Due to an increase in the consumption of food, feed, and fuel and to meet global food security needs for the rapidly growing human population, there is a necessity to obtain high-yielding crops that can adapt to future climate changes. Currently, the main feed source used for ruminant livestock production is forage grasses. In temperate climate zones, perennial grasses grown for feed are widely distributed and tend to suffer under unfavorable environmental conditions. Genome editing has been shown to be an effective tool for the development of abiotic stress-resistant plants. The highly versatile CRISPR-Cas system enables increasingly complex modifications in genomes while maintaining precision and low off-target frequency mutations. In this review, we provide an overview of forage grass species that have been subjected to genome editing. We offer a perspective view on the generation of plants resilient to abiotic stresses. Due to the broad factors contributing to these stresses the review focuses on drought, salt, heat, and cold stresses. The application of new genomic techniques (e.g., CRISPR-Cas) allows addressing several challenges caused by climate change and abiotic stresses for developing forage grass cultivars with improved adaptation to the future climatic conditions. Genome editing will contribute towards developing safe and sustainable food systems.

Molecular Mechanism of Cold Tolerance of Centipedegrass Based on the Transcriptome

Low temperature is an important limiting factor in the environment that affects the distribution, growth and development of warm-season grasses. Transcriptome sequencing has been widely used to mine candidate genes under low-temperature stress and other abiotic stresses. However, the molecular mechanism of centipedegrass in response to low-temperature stress was rarely reported. To understand the molecular mechanism of centipedegrass in response to low-temperature stress, we measured physiological indicators and sequenced the transcriptome of centipedegrass under different stress durations. Under cold stress, the SS content and APX activity of centipedegrass increased while the SOD activity decreased; the CAT activity, POD activity and flavonoid content first increased and then decreased; and the GSH-Px activity first decreased and then increased. Using full-length transcriptome and second-generation sequencing, we obtained 38.76 G subreads. These reads were integrated into 177,178 isoforms, and 885 differentially expressed transcripts were obtained. The expression of AUX_IAA and WRKY transcription factors and HSF transcription-influencing factors increased during cold stress. Through KEGG enrichment analysis, we determined that arginine and proline metabolism, plant circadian rhythm, plant hormone signal transduction and the flavonoid biosynthesis pathways played important roles in the cold stress resistance of centipedegrass. In addition, by using weighted gene coexpression network analysis (WGCNA), we determined that the turquoise module was significantly correlated with SS content and APX activity, while the blue module was significantly negatively correlated with POD and CAT activity. This paper is the first to report the response of centipedegrass to cold stress at the transcriptome level. Our results help to clarify the molecular mechanisms underlying the cold tolerance of warm-season grasses.

Outcomes of Low-Temperature Stress on Biological Alterations within Pothos (Epipremnum aureum) Leaves

Pothos (Epipremnum aureum) is a commonly used indoor ornamental foliage, particularly in the middle and lower regions of the Yangtze River in China. It typically grows in the tropical area, and it is yet unclear whether prolonged winter temperatures cause plant damage and impact its development. In this study, the E. aureum chilling injury response was explored by maintaining it at 1 °C. Based on the acquired results, low-temperature stress (LTS) induced wilting and yellowing of leaves and diminished chloroplast pigment concentrations, particularly the chlorophyll b content. LTS also induced overproduction of reactive oxygen species (ROS) within E. aureum and enhanced the relative electrical conductivity and superoxide dismutase activity. In addition, with prolonged LTS, the anatomical structure of E. aureum was severely damaged, resulting in a marked reduction in the photochemical activity of the photosystem Ⅱ reaction center and suppressed photosynthesis. Moreover, results of the transcriptomic analysis revealed that LTS induced the expression of genes involved in the α-linolenic acid metabolic pathway, plant hormone network, host plant–pathogen association, and MAPK axis, suggesting that LTS would activate its resistant response to cold stress. These results unraveled the physiological and transcriptomical response of E. aureum to chilling injury, which would lay a theoretical foundation for the cultivation of low-temperature-tolerant varieties of E. aureum.

A De Novo Transcriptome Analysis Identifies Cold-Responsive Genes in the Seeds of Taxillus chinensis (DC.) Danser

Taxillus chinensis (DC.) Danser, a parasitic plant of the Loranthaceae family, grows by attacking other plants. It has a long history of being used in Chinese medicine to treat multiple chronic diseases. We previously observed that T. chinensis seeds are sensitive to cold. In this study, we performed transcriptome sequencing for T. chinensis seeds treated by cold (0°C) for 0 h, 12 h, 24 h, and 36 h. TRINITY assembled 257,870 transcripts from 223,512 genes. The GC content and N50 were calculated as 42.29% and 1,368, respectively. Then, we identified 42,183 CDSs and 35,268 likely proteins in the assembled transcriptome, which contained 1,622 signal peptides and 6,795 transmembrane domains. Next, we identified 17,217 genes (FPKM > 5) and 2,333 differentially expressed genes (DEGs) in T. chinensis seeds under cold stress. The MAPK pathway, as an early cold response, was significantly enriched by the DEGs in the T. chinensis seeds after 24 h of cold treatment. Known cold-responsive genes encoding abscisic acid-associated, aquaporin, C-repeat binding factor (CBF), cold-regulated protein, heat shock protein, protein kinase, ribosomal protein, transcription factor (TF), zinc finger protein, and ubiquitin were deregulated in the T. chinensis seeds under cold stress. Notably, the upregulation of CBF gene might be the consequences of the downregulation of MYB and GATA TFs. Additionally, we identified that genes encoding CDC20, YLS9, EXORDIUM, and AUX1 and wound-responsive family protein might be related to novel mechanisms of T. chinensis seeds exposed to cold. This study is first to report the differential transcriptional induction in T. chinensis seeds under cold stress. It will improve our understanding of parasitic plants in response to cold and provide a valuable resource for future studies.

Meta-analysis of transcriptomic responses to cold stress in plants

Transcriptomic analyses are needful tools to gain insight into the molecular mechanisms underlying plant responses to abiotic stresses. The aim of this study was to identify key genes differentially regulated in response to chilling stress in various plant species with different levels of tolerance to low temperatures. A meta-analysis was performed using the RNA-Seq data of published studies whose experimental conditions were comparable. The results confirmed the importance of ethylene in the hormonal cross-talk modulating the defensive responses against chilling stress, especially in sensitive species. The transcriptomic activity of five Ethylene Response Factors genes and a REDOX Responsive Transcription Factor 1 involved in hormone-related pathways belonging to ethylene metabolism and signal transduction were induced. Transcription activity of two genes encoding for heat shock factors was enhanced, together with various genes associated with developmental processes. Several transcription factor families showed to be commonly induced between different plant species. Protein-protein interaction networks highlighted the role of the photosystems I and II, as well as genes encoding for HSF and WRKY transcription factors. A model of gene regulatory network underlying plant responses to chilling stress was developed, allowing the delivery of new candidate genes for genetic improvement of crops towards low temperatures tolerance.

Transcriptome sequencing and gene expression profiling of Pinus sibirica under different cold stresses

Cold stress is a major abiotic factor that affects plant growth and geographical distribution. Pinus sibirica is extremely frigostable tree species. To understand the molecular mechanisms of cold tolerance by P. sibirica, physiological responses were analyzed and transcriptome profiling was conducted to the plants treated by cold stress. The physiological data showed that membrane permeability relative conductivity (REC), reactive oxygen species (ROS), malonaldehyde (MDA) content, peroxidase (POD) and catalase (CAT) activity, soluble sugar, soluble protein and proline contents were increased significantly (p < 0.05) in response to cold stress. Transcriptome analysis identified a total of 871, 1397 and 872 differentially expressed genes (DEGs) after cold treatment for 6 h, 24 h and 48 h at -20°C, respectively. The signaling pathway mediated by Ca²⁺ as a signaling molecule and abscisic acid pathways were the main cold signal transduction pathways in P. sibirica. The APETALA2/Ethylene-Responsive Factor (AP2/ERF) and MYB transcription factor families also play an important role in the transcriptional regulation of P. sibirica. In addition, many genes related to photosynthesis were differentially expressed under cold stress. We also validated the reliability of transcriptome data with quantitative real-time PCR. This study lays the foundation for understanding the molecular mechanisms related to cold responses in P. sibirica.

Physiological and transcriptome analysis of Magnolia denudata leaf buds during long-term cold acclimation

BACKGROUND

Magonlia denudata is an important perennial tree species of the Magnoliaceae family, known for its ornamental value, resistance to smoke pollution and wind, role in air purification, and robust cold tolerance. In this study, a high-throughput transcriptome analysis of leaf buds was performed, and gene expression following artificial acclimation 22 °C, 4 °C and 0 °C, was compared by RNA sequencing.

RESULTS

Over 426 million clean reads were produced from three libraries (22 °C, 4 °C and 0 °C). A total of 74,503 non-redundant unigenes were generated, with an average length of 1173.7 bp (N50 = 1548). Based on transcriptional results, 357 and 235 unigenes were identified as being upregulated and downregulated under cold stress conditions, respectively. Differentially expressed genes were annotated using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway analyses. The transcriptomic analysis focused on carbon metabolism and plant hormone signal transduction associated with cold acclimation. Transcription factors such as those in the basic helix-loop-helix and AP2/ERF families were found to play an important role in M. denudata cold acclimation.

CONCLUSION

M. denudata exhibits responses to non-freezing cold temperature (4 °C) to increase its cold tolerance. Cold resistance was further strengthened with cold acclimation under freezing conditions (0 °C). Cold tolerance genes, and cold signaling transcriptional pathways, and potential functional key components for the regulation of the cold response were identified in M. denudata. These results provide a basis for further studies, and the verification of key genes involved in cold acclimation responses in M. denudata lays a foundation for developing breeding programs for Magnoliaceae species.

The Effect of Low Temperature Stress on the Leaves and MicroRNA Expression of Potato Seedlings

In recent years, with the wanton destruction of the ecological environment by humans and the frequent occurrence of extreme bad weather, many places that should have been warm and blooming in spring have instead experienced the phenomenon of the “April blizzard,” which has seriously affected China's crops, especially spring potato production in most areas. Potato cultivars, especially potato seedlings, are sensitive to frost, and low temperature frost has become one of the most important abiotic stresses affecting potato production. Potato cold tolerance is regulated by a complex gene network. Although some low temperature resistant microRNAs have been identified, little is known about the role of miRNAs in response to low temperature stress in potato. Therefore, the objective of this study is to clarify the influence of low temperature stress on the miRNA expression of potato by comparing the expression differences of miRNA in potato which was treated with different low temperatures. For the study, 307 known miRNAs belonging to 73 small RNA families and 211 novel miRNAs were obtained. When the temperature decreased, the number of both known and novel miRNA decreased, and the minimum temperature was −2°C. Most of the miRNAs respond to low temperature, drought, and disease stress; some conserved miRNAs were first found to respond to low temperature stress in potato, such as stu-miR530, stu-miR156d, and stu-miR167b. The Gene Ontology, Kyoto Encyclopedia of Genes, and Genomes pathway enrichment analysis of 442 different expression miRNAs target genes indicated that there existed diversified low temperature responsive pathways, but Abscisic Acid was found likely to play a central coordinating role in response to low temperature stress in many metabolism pathways. Quantitative real-time PCR assays indicated that the related targets were negatively regulated by the tested different expression miRNAs during low temperature stress. The results indicated that miRNAs may play an important coordination role in response to low temperature stress in many metabolic pathways by regulating abscisic acid and gibberellin, which provided insight into the roles of miRNAs during low temperature stress and would be helpful for alleviating low temperature stress and promoting low temperature resistant breeding in potatoes.

Leaf transcriptomic response mediated by cold stress in two maize inbred lines with contrasting tolerance levels

Maize (Zea mays L.) is a thermophilic plant and a minor drop in temperature can prolong the maturity period. Plants respond to cold stress through structural and functional modification in cell membranes as well as changes in the photosynthesis and energy metabolism. In order to understand the molecular mechanisms underlying cold tolerance and adaptation, we employed leaf transcriptome sequencing together with leaf microstructure and relative electrical conductivity measurements in two maize inbred lines, having different cold stress tolerance potentials. The leaf physiological and transcriptomic responses of maize seedlings were studied after growing both inbred lines at 5 °C for 0, 12 and 24 h. Differentially expressed genes were enriched in photosynthesis antenna proteins, MAPK signaling pathway, plant hormone signal transduction, circadian rhythm, secondary metabolites related pathways, ribosome, and proteasome. The seedlings of both genotypes employed common stress responsive pathways to respond to cold stress. However, the cold tolerant line B144 protected its photosystem II from photooxidation by upregulating D1 proteins. The sensitive line Q319 was unable to close its stomata. Collectively, B144 exhibited a cold tolerance owing to its ability to mediate changes in stomata opening as well as protecting photosystem. These results increase our understanding on the cold stress tolerance in maize seedlings and propose multiple key regulators of stress responses such as modifications in photosystem II, stomata guard cell opening and closing, changes in secondary metabolite biosynthesis, and circadian rhythm. This study also presents the signal transduction related changes in MAPK and phytohormone signaling pathways in response to cold stress during seedling stage of maize.

Integration of Transcriptional and Post-transcriptional Analysis Revealed the Early Response Mechanism of Sugarcane to Cold Stress

Cold stress causes major losses to sugarcane production, yet the precise molecular mechanisms that cause losses due to cold stress are not well-understood. To survey miRNAs and genes involved in cold tolerance, RNA-seq, miRNA-seq, and integration analyses were performed on Saccharum spontaneum. Results showed that a total of 118,015 genes and 6,034 of these differentially expressed genes (DEGs) were screened. Protein–protein interaction (PPI) analyses revealed that ABA signaling via protein phosphatase 2Cs was the most important signal transduction pathway and late embryogenesis abundant protein was the hub protein associated with adaptation to cold stress. Furthermore, a total of 856 miRNAs were identified in this study and 109 of them were differentially expressed in sugarcane responding to cold stress. Most importantly, the miRNA–gene regulatory networks suggested the complex post-transcriptional regulation in sugarcane under cold stress, including 10 miRNAs−42 genes, 16 miRNAs−70 genes, and three miRNAs−18 genes in CT vs. LT0.5, CT vs. LT1, and CT0.5 vs. LT1, respectively. Specifically, key regulators from 16 genes encoding laccase were targeted by novel-Chr4C_47059 and Novel-Chr4A_40498, while five LRR-RLK genes were targeted by Novel-Chr6B_65233 and Novel-Chr5D_60023, 19 PPR repeat proteins by Novel-Chr5C_57213 and Novel-Chr5D_58065. Our findings suggested that these miRNAs and cell wall-related genes played vital regulatory roles in the responses of sugarcane to cold stress. Overall, the results of this study provide insights into the transcriptional and post-transcriptional regulatory network underlying the responses of sugarcane to cold stress.

Transcriptome Sequencing Analysis of Birch (Betula platyphylla Sukaczev) under Low-Temperature Stress

Low temperature is one of the common abiotic stresses that adversely affect the growth and development of plants. In this study, we used RNA-Seq to identify low-temperature-responsive genes in birch and further analyzed the underlying molecular mechanism. Birch seedlings were treated by the low temperature (6 °C) for 0, 1, 1.5, 2, 2.5, and 3 h, respectively. A total of 3491 genes were differentially expressed after low-temperature stress. Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) functional enrichment analysis were performed for the differentially expressed genes (DEGs). GO analysis indicated that 3491 DEGs were distributed into 1002 categories, and these DEGs were enriched in “cell process”, “metabolic process”, and “stimulus response”, under the “biological process” category; in “organelles” and “cell components”, under the “cell component” category; and in “catalytic activity” and “adhesion”, under the “molecular function” category. The KEGG enrichment indicated that 119 DEGs were involved in Ca2+ and plant hormone signal transduction; 205 DEGs were involved in secondary metabolic processes, such as lipid metabolism and phenylpropanoid biosynthesis pathway; and 20 DEGs were involved in photosynthesis. In addition, a total of 362 transcription factors (TFs) were differentially expressed under low-temperature stress, including AP2/ERF, C2H2, MYB-HB-like, WRKY, bHLH, WD40-like, and GRAS families. Gene Bpev01.c0480.g0081 (calmodulin-like CML38), Bpev01.c1074.g0005 (calmodulin-like CML25), Bpev01.c1074.g0001 (Calcium-binding EF-hand family protein), Bpev01.c2029.g0005 (calmodulin-like protein), Bpev01.c0154.g0008 (POD), Bpev01.c0015.g0143 (N-acetyl-1-glutamate synthase), and Bpev01.c0148.g0010 (branched chain amino acid transferase) were up-regulated at a high level, under low-temperature stress.

Physiological and transcriptome analysis of Poa pratensis var. anceps cv. Qinghai in response to cold stress

BACKGROUND

Low temperature limits the growth and development and geographical distribution of plants. Poa pratensis is a cool-season turfgrass mainly grown in urban areas. However, low winter temperature or cold events in spring and autumn may cause P.pratensis mortality, affecting the appearance of lawns. P.pratensis var. anceps cv. Qinghai (PQ) is widely distributed in the Qinghai-Tibet Plateau above 3000 m. PQ has greater cold tolerance than the commercially cultivated P.pratensis varieties. However, existing studies on the response mechanism of PQ to low temperatures have mainly focused on physiological and biochemical perspectives, while changes in the PQ transcriptome during the response to cold stress have not been reported.

RESULTS

To investigate the molecular mechanism of the PQ cold response and identify genes to improve the low-temperature tolerance of P.pratensis, we analyzed and compared the transcriptomes of PQ and the cold-sensitive P.pratensis cv. 'Baron' (PB) under cold stress using RNA sequencing. We identified 5996 and 3285 differentially expressed genes (DEGs) between the treatment vs control comparison of PQ and PB, respectively, with 5612 DEGs specific to PQ. Based on the DEGs, important Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, such as "starch and sucrose metabolism", "protein processing in endoplasmic reticulum", "phenylalanine metabolism" and "glycolysis/gluconeogenesis" were significantly enriched in PQ, and "starch and sucrose metabolism", "phenylpropanoid biosynthesis", "galactose metabolism" and "glutathione metabolism" were significantly enriched in PB. In addition, the "glycolysis" and "citrate cycle (TCA cycle)" pathways were identified as involved in cold tolerance of P.pratensis.

CONCLUSIONS

As we know, this is the first study to explore the transcriptome of P.pratensis var. anceps cv. Qinghai. Our study not noly provides important insights into the molecular mechanisms of P.pratensis var. anceps cv. Qinghai responds to cold stress, but also systematically reveals the changes of key genes and products of glycolysis and TCA cycle in response to cold stress, which is conductive to the breeding of cold-tolerance P.pratensis genotype.

Transcriptome profiles identify the common responsive genes to drought stress in two Elymus species

Elymus, the largest genus of the Triticeae Dumort, is a forage grass in the Qinghai-Tibetan Plateau, where the climate has gradually become increasingly dry in recent years. To understand the mechanisms of the response to drought stress in Elymus species, we first investigated physiological and biochemical responses to polyethylene glycol (PEG-6000) simulated drought stress in two Elymus species, Elymus nutans and Elymus sibiricus, and found that E. nutans was more tolerant to drought stress than E. sibiricus. De novo transcriptome analysis of these two Elymus species treated with or without 10 % PEG-6000 revealed that a total of 1695 unigenes were commonly regulated by drought treatment in these two Elymus species, with 1614 unigenes up-regulated and 81 unigenes down-regulated. The coexpressed differentially expressed genes (DEGs) were enriched in regulation of transcription and gene expression in the GO database. KEGG pathway analysis indicated plant hormone signaling transduction were mostly enriched in co-expressed DEGs. Furthermore, genes annotated in the plant hormone signaling transduction were screened from co-expressed DEGs, and found that abscisic acid plays the major role in the drought stress tolerance of Elymus. Meanwhile, transcription factors screened from co-expressed DEGs were mainly classified into the ERF subfamily and WRKY, DREB, and HSF family members. Our results provide further reference for studying the response mechanism and culturing highly tolerant grasses of the Elymus species under drought stress.

Correlation analysis of cold-related gene expression with physiological and biochemical indicators under cold stress in oil palm

Oil palm (Elaeis guineensis Jacq.) is a representative tropical oil crop that is sensitive to low temperature. Oil palm can experience cold damage when exposed to low temperatures for a long period. During these unfavorable conditions, a series of gene induction/repression and physico-chemical changes occur in oil palm. To better understand the link between these events, we investigated the expression levels of various genes (including COR410, COR413, CBF1, CBF2, CBF3, ICE1-1, ICE1-2, ICE1-4, SIZ1-1, SIZ1-2, ZAT10, ZAT12) and the accumulation of osmolytes (proline, malondialdehyde and sucrose). Likewise, the activity of superoxide dismutase (SOD) in oil palm under cold stress (4°C, 8°C and 12°C) was examined. The results showed a clear link among the expression of CBFs (especially CBF1 and CBF3) and the all genes examined under cold stress (12°C). The expression of CBF1 and CBF2 also exhibited a positive link with the accumulation of sucrose and proline under cold stress in oil palm. At 4°C, the proline content exhibited a very significant correlation with electrolyte leakage in oil palm. The results of this study provide necessary information regarding the mechanism of the response and adaption of oil palm to cold stress. Additionally, they offer clues for the selection or development of cold-tolerant cultivars from the available germplasms of oil palm.

Comparative transcriptome analysis of Rosa chinensis 'Slater's crimson China' provides insights into the crucial factors and signaling pathways in heat stress response

Heat stress limits the growth of roses and adversely affects the yield and the quality of the rose cut-flowers. To investigate the heat stress response (HSR) mechanisms of rose, we compared the transcriptome profiling generated from Rosa chinensis 'Slater's crimson China' exposed to heat stress for five different time duration (0, 0.5, 2, 6, 12 h). Overall, 6175 differentially expressed genes (DGEs) were identified and exhibited different temporal expression patterns. Up-regulated genes related to chaperone-mediated protein folding, signal transduction and ROS scavenging were rapidly induced after 0.5-2 h of heat treatment, which provides evidence for the early adjustments of heat stress response in R. chinensis. While the down-regulated genes related to light reaction, sucrose biosynthesis, starch biosynthesis and cell wall biosynthesis were identified after as short as 6 h of heat stress, which indicated the ongoing negative effects on the physiology of R. chinensis. Using weighted gene co-expression network analysis, we found that different heat stress stages could be delineated by several modules. Based on integrating the transcription factors with upstream enriched DNA motifs of co-expressed genes in these modules, the gene regulation networks were predicted and several regulators of HSR were identified. Of particular importance was the discovery of the module associated with rapid sensing and signal transduction, in which numerous co-expressed genes related to chaperones, Ca²⁺ signaling pathways and transcription factors were identified. The results of this study provided an important resource for further dissecting the role of candidate genes governing the transcriptional regulatory network of HSR in Rose.

Transcriptome profiling reveals candidate flavonol-related genes of Tetrastigma hemsleyanum under cold stress

Background

Tetrastigma hemsleyanum Diels et Gilg is a valuable medicinal herb, whose main bioactive constituents are flavonoids. Chilling sensitivity is the dominant environmental factor limiting growth and development of the plants. But the mechanisms of cold sensitivity in this plant are still unclear. Also, not enough information on genes involved in flavonoid biosynthesis in T. hemsleyanum is available to understand the mechanisms of its physiological and pharmaceutical effects.

Results

The electrolyte leakage, POD activity, soluble protein, and MDA content showed a linear sustained increase under cold stress. The critical period of cold damage in T. hemsleyanum was from 12 h to 48 h. Expression profiles revealed 18,104 differentially expressed genes (DEGs) among these critical time points. Most of the cold regulated DEGs were early-response genes. A total of 114 unigenes were assigned to the flavonoid biosynthetic pathway. Fourteen genes most likely to encode flavonoid biosynthetic enzymes were identified. Flavonols of T. hemsleyanum might play a crucial role in combating cold stress. Genes encoding PAL, 4CL, CHS, ANR, FLS, and LAR were significantly up-regulated by cold stress, which could result in a significant increase in crucial flavonols (catechin, epicatechin, rutin, and quercetin) in T. hemsleyanum.

Conclusions

Overall, our results show that the expression of genes related to flavonol biosynthesis as well as flavonol content increased in T. hemsleyanum under cold stress. These findings provide valuable information regarding the transcriptome changes in response to cold stress and give a clue for identifying candidate genes as promising targets that could be used for improving cold tolerance via molecular breeding. The study also provides candidate genes involved in flavonoid biosynthesis and may be useful for clarifying the biosynthetic pathway of flavonoids in T. hemsleyanum.

Electronic supplementary material

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

De novo transcriptome sequencing and gene expression profiling of Magnolia wufengensis in response to cold stress

BACKGROUND

Magnolia wufengensis is a new species of Magnolia L. and has considerable ornamental and economic value due to its unique characteristics. However, because of its characteristic of poor low temperature resistance, M. wufengensis is hardly popularization and application in the north of China. Furthermore, the mechanisms of gene regulation and signaling pathways involved in the cold-stress response remained unclear in this species. In order to solve the above-mentioned problems, we performed de novo transcriptome assembly and compared the gene expression under the natural (25 °C) and cold (4 °C) conditions for M. wufengensis seedlings.

RESULTS

More than 46 million high-quality clean reads were produced from six samples (RNA was extracted from the leaves) and were used for performing de novo transcriptome assembly. A total of 59,764 non-redundant unigenes with an average length of 899 bp (N50 = 1,110) were generated. Among these unigenes, 31,038 unigenes exhibited significant sequence similarity to known genes, as determined by BLASTx searches (E-value ≤1.0E-05) against the Nr, SwissProt, String, GO, KEGG, and Cluster of COG databases. Based on a comparative transcriptome analysis, 3,910 unigenes were significantly differentially expressed (false discovery rate [FDR] < 0.05 and |log2FC (CT/CK)| ≥ 1) in the cold-treated samples, and 2,616 and 1,294 unigenes were up- and down-regulated by cold stress, respectively. Analysis of the expression patterns of 16 differentially expressed genes (DEGs) by quantitative real-time RT-PCR (qRT-PCR) confirmed the accuracy of the RNA-Seq results. Gene Ontology and KEGG pathway functional enrichment analyses allowed us to better understand these differentially expressed unigenes. The most significant transcriptomic changes observed under cold stress were related to plant hormone and signal transduction pathways, primary and secondary metabolism, and photosynthesis. In addition, 113 transcription factors, including members of the AP2-EREBP, bHLH, WRKY, MYB, NAC, HSF, and bZIP families, were identified as cold responsive.

CONCLUSION

We generated a genome-wide transcript profile of M. wufengensis and a de novo-assembled transcriptome that can be used to analyze genes involved in biological processes. In this study, we provide the first report of transcriptome sequencing of cold-stressed M. wufengensis. Our findings provide important clues not only for understanding the molecular mechanisms of cold stress in plants but also for introducing cold hardiness into M. wufengensis.

Comparative transcriptome analysis reveals an early gene expression profile that contributes to cold resistance in Hevea brasiliensis (the Para rubber tree)

The rubber tree (Hevea brasiliensis Muell. Arg) is a tropical, perennial, woody plant that is susceptible to cold stress. In China, cold stress has been found to severely damage rubber plants in plantations in past decades. Although several Hevea clones that are resistant to cold have been developed, their cold hardiness mechanism has yet to be elucidated. For the study reported herein, we subjected the cold-resistant clone CATAS93-114 and the cold-sensitive clone Reken501 to chilling stress, and characterized their transcriptomes at 0, 2, 8 and 24 h after the start of chilling. We found that 7870 genes were differentially expressed in the transcriptomes of the two clones. In CATAS93-114, a greater number of genes were found to be up- or downregulated between 2 h and 8 h than in Reken501, which indicated a more rapid and intensive response by CATAS93-114 than by Reken501. The differentially expressed genes were grouped into seven major clusters, according to their Gene Ontology terms. The expression profiles for genes involved in abscisic acid metabolism and signaling, in an abscisic acid-independent pathway, and in early signal perception were found to have distinct expression patterns for the transcriptomes of the two clones. The differential expression of 22 genes that appeared to have central roles in response to chilling was confirmed by quantitative real-time PCR.

Shifts in functional trait-species abundance relationships over secondary subalpine meadow succession in the Qinghai-Tibetan Plateau

Although trait-based processes of community assembly during secondary succession invokes multiple factors that ultimately determine the presence or absence of a species, little is known regarding the impacts of functional traits on species abundance in successional plant communities. Here in species-rich subalpine secondary successional meadows of the Qinghai-Tibetan Plateau, we measured photosynthesis rate and leaf proline content that are related to plant growth and abiotic stress resistance, respectively, and seed germination rate that is closely correlated with plant germination strategy to test their influence on species abundance during succession. We used a linear mixed effects model framework to examine the shifts in trait-abundance relationships and the correlations among these three traits in successional communities. We observed significant shifts in trait-abundance relationships during succession, e.g., abundant species in early-successional meadows exhibited relatively high photosynthesis rates and leaf proline content, but showed low seed germination rates, whereas the converse were true in late successional communities. However, the correlations among the three traits were insignificant in most meadow communities. Our results show that functional traits associated with plant growth, stress resistance, and reproduction impose strong influence on species abundance during secondary subalpine meadow succession in the Qinghai-Tibetan Plateau.

De novo transcriptome sequencing of Flammulina velutipes uncover candidate genes associated with cold-induced fruiting

To understand molecular mechanism of cold-induced fruiting in Flammulina velutipes, which is one of most popular edible fungi in east Asia, de novo assembly of the F. velutipes transcriptome was carried out. There were 26,888,494 and 26,275,146 clean reads obtained from mycelium and primordia of F. velutipes, respectively. A total of 20,157 unigenes were de novo assembled and 15,058 of them were annotated. Moreover, 7935 unigenes were differentially expressed between mycelium and primordia, 4025 of them were up-regulated and 3910 were down-regulated. GO and KEGG pathway analysis of the differentially expressed unigenes indicated that functional groups associated with two-component signaling pathway, calcium signaling, mitogen-actived protein kinase pathway, molecular chaperones, cell wall and membrane system, play an important role in cold-induced fruiting of F. velutipes. In this work 643 EST-SSRs were identified in 20,157 unigenes and 1560 EST-SSRs primers pairs were designed. Moreover, 5548 and 5955 SNPs were detected in mycelium and primordia, respectively. Consequently, results of this work can serve as a valuable resource for functional genomics study of cold-induced fruiting in F. velutipes.

Plant Glycine-Rich Proteins in Stress Response: An Emerging, Still Prospective Story

Seed plants are sessile organisms that have developed a plethora of strategies for sensing, avoiding, and responding to stress. Several proteins, including the glycine-rich protein (GRP) superfamily, are involved in cellular stress responses and signaling. GRPs are characterized by high glycine content and the presence of conserved segments including glycine-containing structural motifs composed of repetitive amino acid residues. The general structure of this superfamily facilitates division of GRPs into five main subclasses. Although the participation of GRPs in plant stress response has been indicated in numerous model and non-model plant species, relatively little is known about the key physiological processes and molecular mechanisms in which those proteins are engaged. Class I, II, and IV members are known to be involved in hormone signaling, stress acclimation, and floral development, and are crucial for regulation of plant cells growth. GRPs of class IV [RNA-binding proteins (RBPs)] are involved in alternative splicing or regulation of transcription and stomatal movement, seed, pollen, and stamen development; their accumulation is regulated by the circadian clock. Owing to the fact that the overexpression of GRPs can confer tolerance to stress (e.g., some are involved in cold acclimation and may improve growth at low temperatures), these proteins could play a promising role in agriculture through plant genetic engineering. Consequently, isolation, cloning, characterization, and functional validation of novel GRPs expressed in response to the diverse stress conditions are expected to be growing areas of research in the coming years. According to our knowledge, this is the first comprehensive review on participation of plant GRPs in the response to diverse stress stimuli.

Comparison of Gene Expression Profiles in Nonmodel Eukaryotic Organisms with RNA-Seq

With recent advances of next-generation sequencing technology, RNA-Sequencing (RNA-Seq) has emerged as a powerful approach for the transcriptomic profiling. RNA-Seq has been used in almost every field of biological studies, and has greatly extended our view of transcriptomic complexity in different species. In particular, for nonmodel organisms which are usually without high-quality reference genomes, the de novo transcriptome assembly from RNA-Seq data provides a solution for their comparative transcriptomic study. In this chapter, we focus on the comparative transcriptomic analysis of nonmodel organisms. Two analysis strategies (without or with reference genome) are described step-by-step, with the differentially expressed genes explored.

Cold-responsive miRNAs and their target genes in the wild eggplant species Solanum aculeatissimum

Background

Low temperature is an important abiotic stress in plant growth and development, especially for thermophilic plants. Eggplants are thermophilic vegetables, although the molecular mechanism of their response to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although the role of many plant miRNAs in facilitating chilling tolerance has been verified, little is known about the mechanisms of eggplant chilling tolerance.

Results

Here, we used high-throughput sequencing to extract the miRNA and target genes expression profiles of Solanum aculeatissimum (S. aculeatissimum) under low temperature stress at different time periods(0 h, 2 h, 6 h, 12 h, 24 h). Differentially regulated miRNAs and their target genes were analyzed by comparing the small RNA (sRNA) and miRBase 20.0 databases using BLAST or BOWTIE, respectively. Fifty-six down-regulated miRNAs and 28 up-regulated miRNAs corresponding to 220 up-regulated mRNAs and 94 down-regulated mRNAs, respectively, were identified in S. aculeatissimum. Nine significant differentially expressed miRNAs and twelve mRNAs were identified by quantitative Real-time PCR and association analysis, and analyzed for their GO function enrichment and KEGG pathway association.

Conclusions

In summary, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, which provides a theoretical basis for the further study of low temperature stress-related miRNAs and the regulation of cold-tolerance mechanisms of eggplant at the miRNA level.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4341-y) contains supplementary material, which is available to authorized users.

Genome-Wide Identification and Analysis of Genes, Conserved between japonica and indica Rice Cultivars, that Respond to Low-Temperature Stress at the Vegetative Growth Stage

Cold stress is very detrimental to crop production. However, only a few genes in rice have been identified with known functions related to cold tolerance. To meet this agronomic challenge more effectively, researchers must take global approaches to select useful candidate genes and find the major regulatory factors. We used five Gene expression omnibus series data series of Affymetrix array data, produced with cold stress-treated samples from the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/), and identified 502 cold-inducible genes common to both japonica and indica rice cultivars. From them, we confirmed that the expression of two randomly chosen genes was increased by cold stress in planta. In addition, overexpression of OsWRKY71 enhanced cold tolerance in 'Dongjin,' the tested japonica cultivar. Comparisons between japonica and indica rice, based on calculations of plant survival rates and chlorophyll fluorescence, confirmed that the japonica rice was more cold-tolerant. Gene Ontology enrichment analysis indicate that the 'L-phenylalanine catabolic process,' within the Biological Process category, was the most highly overrepresented under cold-stress conditions, implying its significance in that response in rice. MapMan analysis classified 'Major Metabolic' processes and 'Regulatory Gene Modules' as two other major determinants of the cold-stress response and suggested several key cis-regulatory elements. Based on these results, we proposed a model that includes a pathway for cold stress-responsive signaling. Results from our functional analysis of the main signal transduction and transcription regulation factors identified in that pathway will provide insight into novel regulatory metabolism(s), as well as a foundation by which we can develop crop plants with enhanced cold tolerance.