Comparative transcriptomics analysis reveals difference of key gene expression between banana and plantain in response to cold stress

Discovery of gene regulation mechanisms associated with uniconazole-induced cold tolerance in banana using integrated transcriptome and metabolome analysis

BACKGROUND

The gibberellic acid (GA) inhibitor, uniconazole, is a plant growth regulator commonly used in banana cultivation to promote dwarfing but also enhances the cold resistance in plants. However, the mechanism of this induced cold resistance remains unclear.

RESULTS

We confirmed that uniconazole induced cold tolerance in bananas and that the activities of Superoxide dismutase and Peroxidase were increased in the uniconazole-treated bananas under cold stress when compared with the control groups. The transcriptome and metabolome of bananas treated with or without uniconazole were analyzed at different time points under cold stress. Compared to the control group, differentially expressed genes (DEGs) between adjacent time points in each uniconazole-treated group were enriched in plant-pathogen interactions, MAPK signaling pathway, and plant hormone signal transduction, which were closely related to stimulus-functional responses. Furthermore, the differentially abundant metabolites (DAMs) between adjacent time points were enriched in flavone and flavonol biosynthesis and linoleic acid metabolism pathways in the uniconazole-treated group than those in the control group. Temporal analysis of DEGs and DAMs in uniconazole-treated and control groups during cold stress showed that the different expression patterns in the two groups were enriched in the linoleic acid metabolism pathway. In addition to strengthening the antioxidant system and complex hormonal changes caused by GA inhibition, an enhanced linoleic acid metabolism can protect cell membrane stability, which may also be an important part of the cold resistance mechanism of uniconazole treatment in banana plants.

CONCLUSIONS

This study provides information for understanding the mechanisms underlying inducible cold resistance in banana, which will benefit the production of this economically important crop.

Genome-wide identification of XTH gene family in Musa acuminata and response analyses of MaXTHs and xyloglucan to low temperature

Banana (Musa spp.) production is seriously threatened by low temperature (LT) in tropical and subtropical regions. Xyloglucan endotransglycosylase/hydrolases (XTHs) are considered chief enzymes in cell wall remodelling and play a central role in stress responses. However, whether MaXTHs are involved in the low temperature stress tolerance in banana is not clear. Here, the identification and characterization of MaXTHs were carried out, followed by prediction of their cis-acting elements and protein-protein interactions. In addition, candidate MaXTHs involved in banana tolerance to LT were screened through a comparison of their responses to LT between tolerant and sensitive cultivars using RNA-Seq analysis. Moreover, immunofluorescence (IF) labelling was employed to compare changes in the temporal and spatial distribution of different types of xyloglucan components between these two cultivars upon stress. In total, 53 MaXTHs have been identified, and all were predicted to be located in the cell wall, 14 of them also in the cytoplasm. Only 11 MaXTHs have been found to interact with other proteins. Among 16 MaXTHs with LT responsiveness elements, MaXTH26/29/32/35/50 (Group I/II members) and MaXTH7/8 (Group IIIB members) might be involved in banana tolerance to LT stress. IF results suggested that the content of xyloglucan components recognized by CCRC-M87/103/104/106 antibodies might be negatively related to banana chilling tolerance. In conclusion, we have identified the MaXTH gene family and assessed cell wall re-modelling under LT stress. These results will be beneficial for banana breeding against stresses and enrich the cell wall-mediated resistance mechanism in plants to stresses.

Comparative transcriptomics of the chilling stress response in two Asian mangrove species, Bruguiera gymnorhiza and Rhizophora apiculata

Low temperatures largely determine the geographic limits of plant species by reducing survival and growth. Inter-specific differences in the geographic distribution of mangrove species have been associated with cold tolerance, with exclusively tropical species being highly cold-sensitive and subtropical species being relatively cold-tolerant. To identify species-specific adaptations to low temperatures, we compared the chilling stress response of two widespread Indo-West Pacific mangrove species from Rhizophoraceae with differing latitudinal range limits-Bruguiera gymnorhiza (L.) Lam. ex Savigny (subtropical range limit) and Rhizophora apiculata Blume (tropical range limit). For both species, we measured the maximum photochemical efficiency of photosystem II (Fv/Fm) as a proxy for the physiological condition of the plants and examined gene expression profiles during chilling at 15 and 5 °C. At 15 °C, B. gymnorhiza maintained a significantly higher Fv/Fm than R. apiculata. However, at 5 °C, both species displayed equivalent Fv/Fm values. Thus, species-specific differences in chilling tolerance were only found at 15 °C, and both species were sensitive to chilling at 5 °C. At 15 °C, B. gymnorhiza downregulated genes related to the light reactions of photosynthesis and upregulated a gene involved in cyclic electron flow regulation, whereas R. apiculata downregulated more RuBisCo-related genes. At 5 °C, both species repressed genes related to CO2 assimilation. The downregulation of genes related to light absorption and upregulation of genes related to cyclic electron flow regulation are photoprotective mechanisms that likely contributed to the greater photosystem II photochemical efficiency of B. gymnorhiza at 15 °C. The results of this study provide evidence that the distributional range limits and potentially the expansion rates of plant species are associated with differences in the regulation of photosynthesis and photoprotective mechanisms under low temperatures.

Integration of metabolomics and transcriptomics analyses reveals the mechanism of nano-selenium treated to activate phenylpropanoid metabolism and enhance the antioxidant activity of peach

Foliar spraying to improve the quality of fruits is a general approach nowadays. In this study, 10 ppm nano-selenium (nano-Se) diluted with distilled water was sprayed on peach leaves every 10 days for a total of 7 sprays during the fruit set period. And then their fruit quality was compared with that of control group. It was found that the firmness, soluble solid concentration, total phenol, and proanthocyanidin content of the peaches were raised after the nano-Se treatment. Moreover, the ascorbic acid glutathione loop (ASA-GSH loop) was fully activated in the nano-Se treated group, and the associated antioxidant capacity and enzyme activity were significantly increased. Metabolomics revealed that nano-Se could upregulate some metabolites, such as phenylalanine, naringenin, and pinocembrin, to fully activate the metabolism of phenylpropanoids. Further, based on transcriptomics, nano-Se treatment was found to affect fruit quality by regulating genes related to phenylpropanoid metabolism, such as arogenate/prephenate dehydratase (ADT), genes related to abscisic acid metabolism such as (+)-abscisic acid 8'-hydroxylase (CYP707A), and some transcription factors such as MYB. Based on the comprehensive analysis of physicochemical indicators, metabolomics, and transcriptomics, it was found that nano-Se improved fruit quality by activating phenylpropanoid metabolism and enhancing antioxidant capacity. This work provides insights into the mechanism of the effect of nano-Se fertilizer on peach fruit quality. PRACTICAL APPLICATION: The firmness and soluble solid concentration of peaches are higher after nano-Se treatment, which is more in line with people's demand for hard soluble peaches like "Yingzui." The antioxidant capacity, antioxidant substance content, and antioxidant enzyme activity of nano-Se-treated peaches are higher, with potential storage resistance and health effects on human body. The mechanism of nano-Se affecting peach quality was analyzed by metabolomics and transcriptomics, which is a reference and guide for the research and application of nano-Se.

Transcriptome assessment in 'Red Globe' grapevine zygotic embryos during the cooling and warming phase of the cryopreservation procedure

Cryopreservation has the potential for long-term germplasm storage. The metabolic pathways and gene regulation involved in cryopreservation procedures are still not well documented. Hence, the genetic expression profile was evaluated using RNA-Seq in zygotic embryos of grapevines subjected to cryopreservation by vitrification. Sequencing was performed on the Illumina NextSeq 500. The average alignment of reads was 96% against the reference genome. The expression profiles showed 229 genes differentially expressed (186 repressed and 46 induced). The main biological processes showing upregulated enrichment were related to nucleosome assembly, while downregulated processes were related to organ growth. The most highly repressed processes were associated with the organization of the cell wall and membrane components. The unnamed protein product and 17.3 kDa class II heat shock protein (HSP) were highly induced, while ATPase subunit 1 and expansin-A1 were repressed. The response to the cooling and warming process during cryopreservation probably indicates that the changes occurring in transcription may be related to epigenetics. In addition, the cell exhibits an increase in the reserve of nutrients while seeking to survive modestly using available energy and pausing the plant's development. Additionally, energy containment occurred to cope with the stress caused by the treatment where deactivation of components of the cell membrane was observed, possibly due to changes in fluidity caused by alterations in temperature.

Genome-wide identification and expression reveal the involvement of the FCS-like zinc finger (FLZ) gene family in Gossypium hirsutum at low temperature

FCS-like zinc finger (FLZ) is a plant-specific gene family that plays an important regulatory role in plant growth and development and its response to stress. However, studies on the characteristics and functions of cotton FLZ family genes are still lacking. This study systematically identified members of the cotton FLZ gene family based on cotton genome data. The cotton FLZ family genes were systematically analyzed by bioinformatics, and their expression patterns in different tissues and under low-temperature stress were analyzed by transcriptome and qRT-PCR. The G. hirsutum genome contains 56 FLZ genes distributed on 20 chromosomes, and most of them are located in the nucleus. According to the number and evolution analysis of FLZ family genes, FLZ family genes can be divided into five subgroups in cotton. The G. hirsutum FLZ gene has a wide range of tissue expression types, among which the expression is generally higher in roots, stems, leaves, receptacles and calyx. Through promoter analysis, it was found that it contained the most cis-acting elements related to methyl jasmonate (MeJA) and abscisic acid (ABA). Combined with the promoter and qRT-PCR results, it was speculated that GhFLZ11, GhFLZ25, GhFLZ44 and GhFLZ55 were involved in the response of cotton to low-temperature stress. Taken together, our findings suggest an important role for the FLZ gene family in the cotton response to cold stress. This study provides an important theoretical basis for further research on the function of the FLZ gene family and the molecular mechanism of the cotton response to low temperature.

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.

Recent advances in proteomics and metabolomics in plants

Over the past decade, systems biology and plant-omics have increasingly become the main stream in plant biology research. New developments in mass spectrometry and bioinformatics tools, and methodological schema to integrate multi-omics data have leveraged recent advances in proteomics and metabolomics. These progresses are driving a rapid evolution in the field of plant research, greatly facilitating our understanding of the mechanistic aspects of plant metabolisms and the interactions of plants with their external environment. Here, we review the recent progresses in MS-based proteomics and metabolomics tools and workflows with a special focus on their applications to plant biology research using several case studies related to mechanistic understanding of stress response, gene/protein function characterization, metabolic and signaling pathways exploration, and natural product discovery. We also present a projection concerning future perspectives in MS-based proteomics and metabolomics development including their applications to and challenges for system biology. This review is intended to provide readers with an overview of how advanced MS technology, and integrated application of proteomics and metabolomics can be used to advance plant system biology research.

A chromosome-level reference genome of Ensete glaucum gives insight into diversity and chromosomal and repetitive sequence evolution in the Musaceae

BACKGROUND

Ensete glaucum (2n = 2x = 18) is a giant herbaceous monocotyledonous plant in the small Musaceae family along with banana (Musa). A high-quality reference genome sequence assembly of E. glaucum is a resource for functional and evolutionary studies of Ensete, Musaceae, and the Zingiberales.

FINDINGS

Using Oxford Nanopore Technologies, chromosome conformation capture (Hi-C), Illumina and RNA survey sequence, supported by molecular cytogenetics, we report a high-quality 481.5 Mb genome assembly with 9 pseudo-chromosomes and 36,836 genes. A total of 55% of the genome is composed of repetitive sequences with predominantly LTR-retroelements (37%) and DNA transposons (7%). The single 5S ribosomal DNA locus had an exceptionally long monomer length of 1,056 bp, more than twice that of the monomers at multiple loci in Musa. A tandemly repeated satellite (1.1% of the genome, with no similar sequence in Musa) was present around all centromeres, together with a few copies of a long interspersed nuclear element (LINE) retroelement. The assembly enabled us to characterize in detail the chromosomal rearrangements occurring between E. glaucum and the x = 11 species of Musa. One E. glaucum chromosome has the same gene content as Musa acuminata, while others show multiple, complex, but clearly defined evolutionary rearrangements in the change between x= 9 and 11.

CONCLUSIONS

The advance towards a Musaceae pangenome including E. glaucum, tolerant of extreme environments, makes a complete set of gene alleles, copy number variation, and a reference for structural variation available for crop breeding and understanding environmental responses. The chromosome-scale genome assembly shows the nature of chromosomal fusion and translocation events during speciation, and features of rapid repetitive DNA change in terms of copy number, sequence, and genomic location, critical to understanding its role in diversity and evolution.

Genome-wide identification and expression analysis of the cucumber PYL gene family

Abscisic acid (ABA) is a very important hormone in plants. It regulates growth and development of plants and plays an important role in biotic and abiotic stresses. The Pyrabactin resistance 1-like (PYR/PYL) proteins play a central role in ABA signal transduction pathways. The working system of PYL genes in cucumber, an important economical vegetable (Cucumis sativus L.), has not been fully studied yet. Through bioinformatics, a total of 14 individual PYL genes were identified in Chinese long '9930' cucumber. Fourteen PYL genes were distributed on six chromosomes of cucumber, and their encoded proteins predicted to be distributed in cytoplasm and nucleus. Based on the phylogenetic analysis, the PYL genes of cucumber, Arabidopsis, rice, apple, Brachypodium distachyon and soybeancould be classified into three groups. Genetic structures and conserved domains analysis revealed that CsPYL genes in the same group have similar exons and conserved domains. By predicting cis-elements in the promoters, we found that all CsPYL members contained hormone and stress-related elements. Additionally, the expression patterns of CsPYL genes were specific in tissues. Finally, we further examined the expression of 14 CsPYL genes under ABA, PEG, salt stress. The qRT-PCR results showed that most PYL gene expression levels were up-regulated. Furthermore, with different treatments about 3h, the relative expression of PYL8 was up-regulated and more than 20 times higher than 0h. It indicated that this gene may play an important role in abiotic stress.

Genome-wide identification and expression analysis of the MYC transcription factor family and its response to sulfur stress in cabbage (Brassica oleracea L.)

MYC transcriptional factors are members of the bHLH (basic helix-loop-helix) superfamily, and play important roles in plant growth, biological and abiotic stress. Recent studies have revealed that some MYCs are involved in the synthesis of sulfur-containing secondary metabolites. Cabbage, as a typical sulfur-loving crop and rich in sulfur-containing secondary metabolites, the regulatory relationship between sulfur stress and MYC gene family, related reports are relatively rare. In this study, we conducted the first genome-wide analysis of the MYC transcription factor family of cabbage and identified 17 BoMYC genes. Homology of the 17 BoMYC genes, 12 Arabidopsis, 12 Chinese cabbage, 8 wheat and 21 maize MYC were analyzed using the phylogenetic analysis. Meanwhile, chromosome locations, physical and chemical characteristics, gene structures, conserved motif, cis-element, specific expression in different tissues were studied. Finally, we analyzed the expression of the BoMYC gene under sulfur stress and its GO annotation and KEGG enrichment analysis, determined the expression of the BoMYC gene under hormone treatment and the growth index, photosynthetic capacity and hormone content in the leaves. This study is of great significance for functional identification and revealed the effect of S on BoMYC transcription factors.

Comprehensive Proteome and Lysine Acetylome Analysis Reveals the Widespread Involvement of Acetylation in Cold Resistance of Pepper (Capsicum annuum L.)

Pepper is a typical warmth-loving vegetable that lacks a cold acclimation mechanism and is sensitive to cold stress. Lysine acetylation plays an important role in diverse cellular processes, but limited knowledge is available regarding acetylation modifications in the resistance of pepper plants to cold stress. In this study, the proteome and acetylome of two pepper varieties with different levels of cold resistance were investigated by subjecting them to cold treatments of varying durations followed by recovery periods. In total, 6,213 proteins and 4,574 lysine acetylation sites were identified, and this resulted in the discovery of 3,008 differentially expressed proteins and 768 differentially expressed acetylated proteins. A total of 1,988 proteins were identified in both the proteome and acetylome, and the functional differences in these co-identified proteins were elucidated through GO enrichment. KEGG analysis showed that 397 identified acetylated proteins were involved in 93 different metabolic pathways. The dynamic changes in the acetylated proteins in photosynthesis and the "carbon fixation in the photosynthetic organisms" pathway in pepper under low-temperature stress were further analyzed. It was found that acetylation of the PsbO and PsbR proteins in photosystem II and the PsaN protein in photosystem I could regulate the response of pepper leaves to cold stress. The acetylation levels of key carbon assimilation enzymes, such as ribulose bisphosphate carboxylase, fructose-1,6-bisphosphatase, sedoheptulose-1,7-bisphosphatase, glyceraldehyde 3-phosphate dehydrogenase, phosphoribulokinase, and triosephosphate isomerase decreased, leading to decreases in carbon assimilation capacity and photosynthetic efficiency, reducing the cold tolerance of pepper leaves. This study is the first to identify the acetylome in pepper, and it greatly expands the catalog of lysine acetylation substrates and sites in Solanaceae crops, providing new insights for posttranslational modification studies.

Transcriptomic responses and physiological changes to cold stress among natural populations provide insights into local adaptation of weeping forsythia

Genetic mechanisms of species local adaptation are an emerging topic of great interest in evolutionary biology and molecular ecology. In this study, we compared the changes of physiological and phenotypic indexes and gene expression of four weeping forsythia populations under cold stress through a common garden experiment. Physiological and phenotypic results showed that there were differences in cold tolerance among populations. cold tolerance of high the latitude population (HBWZ) was the strongest, followed by the middle latitude population (SXWL), while the low latitude populations (SXHM) and (SXLJ) expressed the weakest cold tolerance. We identified significant differences in gene expression of cold tolerance related pathways and ontologies, including genes of oxylipin and isoquinoline alkaloid biosynthetic process, galactose, tyrosine and unsaturated fatty acids metabolism, among these populations under the same experimental temperature treatments. Even under the same degree of stress, there were notable differences in gene expression among natural populations. In this study, we present a working model of weeping forsythia populations which evolved in the context of different intensities of cold stress. Our study provides new insights for comprehending the genetic mechanisms of local adaptation for non-model species.

The Endophytic Fungus Piriformospora indica Reprograms Banana to Cold Resistance

Banana (Musa spp.), one of the most important fruits worldwide, is generally cold sensitive. In this study, by using the cold-sensitive banana variety Tianbaojiao (Musa acuminate) as the study material, we investigated the effects of Piriformospora indica on banana cold resistance. Seedlings with and without fungus colonization were subjected to 4 °C cold treatment. The changes in plant phenotypes, some physiological and biochemical parameters, chlorophyll fluorescence parameters, and the expression of eight cold-responsive genes in banana leaves before and after cold treatment were measured. Results demonstrated that P. indica colonization reduced the contents of malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) but increased the activities of superoxide dismutase (SOD) and catalase (CAT) and the contents of soluble sugar (SS) and proline. Noteworthily, the CAT activity and SS content in the leaves of P. indica-colonized banana were significant (p < 0.05). After 24 h cold treatment, the decline in maximum photochemistry efficiency of photosystem II (F_v/F_m), photochemical quenching coefficient (qP), efficient quantum yield [Y(II)], and photosynthetic electron transport rate (ETR) in the leaves of P. indica-colonized banana was found to be lower than in the non-inoculated controls (p < 0.05). Moreover, although the difference was not significant, P. indica colonization increased the photochemical conversion efficiency and electron transport rate and alleviated the damage to the photosynthetic reaction center of banana leaves under cold treatment to some extent. Additionally, the expression of the most cold-responsive genes in banana leaves was significantly induced by P. indica during cold stress (p < 0.05). It was concluded that P. indica confers banana with enhanced cold resistance by stimulating antioxidant capacity, SS accumulation, and the expression of cold-responsive genes in leaves. The results obtained from this study are helpful for understanding the P. indica-induced cold resistance in banana.

MaMAPK3-MaICE1-MaPOD P7 pathway, a positive regulator of cold tolerance in banana

BACKGROUND

Banana is a tropical fruit with a high economic impact worldwide. Cold stress greatly affects the development and production of banana.

RESULTS

In the present study, we investigated the functions of MaMAPK3 and MaICE1 involved in cold tolerance of banana. The effect of RNAi of MaMAPK3 on Dajiao (Musa spp. 'Dajiao'; ABB Group) cold tolerance was evaluated. The leaves of the MaMAPK3 RNAi transgenic plants showed wilting and severe necrotic symptoms, while the wide-type (WT) plants remained normal after cold exposure. RNAi of MaMAPK3 significantly changed the expressions of the cold-responsive genes, and the oxidoreductase activity was significantly changed in WT plants, while no changes in transgenic plants were observed. MaICE1 interacted with MaMAPK3, and the expression level of MaICE1 was significantly decreased in MaMAPK3 RNAi transgenic plants. Over-expression of MaICE1 in Cavendish banana (Musa spp. AAA group) indicated that the cold resistance of transgenic plants was superior to that of the WT plants. The POD P7 gene was significantly up-regulated in MaICE1-overexpressing transgenic plants compared with WT plants, and the POD P7 was proved to interact with MaICE1.

CONCLUSIONS

Taken together, our work provided new and solid evidence that MaMAPK3-MaICE1-MaPOD P7 pathway positively improved the cold tolerance in monocotyledon banana, shedding light on molecular breeding for the cold-tolerant banana or other agricultural species.

Genome-Wide Identification of Banana Csl Gene Family and Their Different Responses to Low Temperature between Chilling-Sensitive and Tolerant Cultivars

The cell wall plays an important role in responses to various stresses. The cellulose synthase-like gene (Csl) family has been reported to be involved in the biosynthesis of the hemicellulose backbone. However, little information is available on their involvement in plant tolerance to low-temperature (LT) stress. In this study, a total of 42 Csls were identified in Musa acuminata and clustered into six subfamilies (CslA, CslC, CslD, CslE, CslG, and CslH) according to phylogenetic relationships. The genomic features of MaCsl genes were characterized to identify gene structures, conserved motifs and the distribution among chromosomes. A phylogenetic tree was constructed to show the diversity in these genes. Different changes in hemicellulose content between chilling-tolerant and chilling-sensitive banana cultivars under LT were observed, suggesting that certain types of hemicellulose are involved in LT stress tolerance in banana. Thus, the expression patterns of MaCsl genes in both cultivars after LT treatment were investigated by RNA sequencing (RNA-Seq) technique followed by quantitative real-time PCR (qPCR) validation. The results indicated that MaCslA4/12, MaCslD4 and MaCslE2 are promising candidates determining the chilling tolerance of banana. Our results provide the first genome-wide characterization of the MaCsls in banana, and open the door for further functional studies.

Integrated Analysis of Metabolome and Transcriptome Reveals Insights for Cold Tolerance in Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is an important oilseed crop in the world. Its productivity is significantly influenced by numerous abiotic stresses, including cold stress (CS). Consequently, enhancement in CS tolerance is becoming an important area for agricultural investigation and crop improvement. Therefore, the current study aimed to identify the stress-responsive genes, metabolites, and metabolic pathways based on a combined transcriptome and metabolome analysis to understand the CS responses and tolerance mechanisms in the cold-tolerant (C18) and cold-sensitive (C6) rapeseed varieties. Based on the metabolome analysis, 31 differentially accumulated metabolites (DAMs) were identified between different comparisons of both varieties at the same time points. From the transcriptome analysis, 2,845, 3,358, and 2,819 differentially expressed genes (DEGs) were detected from the comparison of C6-0 vs. C18-0, C6-1 vs. C18-1, and C6-7 vs. C18-7. By combining the transcriptome and metabolome data sets, we found that numerous DAMs were strongly correlated with several differentially expressed genes (DEGs). A functional enrichment analysis of the DAMs and the correlated DEGs specified that most DEGs and DAMs were mainly enriched in diverse carbohydrates and amino acid metabolisms. Among them, starch and sucrose metabolism and phenylalanine metabolism were significantly enriched and played a vital role in the CS adaption of rapeseed. Six candidate genes were selected from the two pathways for controlling the adaption to low temperature. In a further validation, the T-DNA insertion mutants of their Arabidopsis homologous, including 4cl3, cel5, fruct4, ugp1, axs1, and bam2/9, were characterized and six lines differed significantly in levels of freezing tolerance. The outcome of the current study provided new prospects for the understanding of the molecular basis of CS responses and tolerance mechanisms in rapeseed and present a set of candidate genes for use in improving CS adaptability in the same plant.

Acceleration of Carbon Fixation in Chilling-Sensitive Banana under Mild and Moderate Chilling Stresses

Banana is one of the most important food and fruit crops in the world and its growth is ceasing at 10–17 °C. However, the mechanisms determining the tolerance of banana to mild (>15 °C) and moderate chilling (10–15 °C) are elusive. Furthermore, the biochemical controls over the photosynthesis in tropical plant species at low temperatures above 10 °C is not well understood. The purpose of this research was to reveal the response of chilling-sensitive banana to mild (16 °C) and moderate chilling stress (10 °C) at the molecular (transcripts, proteins) and physiological levels. The results showed different transcriptome responses between mild and moderate chilling stresses, especially in pathways of plant hormone signal transduction, ABC transporters, ubiquinone, and other terpenoid-quinone biosynthesis. Interestingly, functions related to carbon fixation were assigned preferentially to upregulated genes/proteins, while photosynthesis and photosynthesis-antenna proteins were downregulated at 10 °C, as revealed by both digital gene expression and proteomic analysis. These results were confirmed by qPCR and immunofluorescence labeling methods. Conclusion: Banana responded to the mild chilling stress dramatically at the molecular level. To compensate for the decreased photosynthesis efficiency caused by mild and moderate chilling stresses, banana accelerated its carbon fixation, mainly through upregulation of phosphoenolpyruvate carboxylases.

The R2R3-MYB gene family in banana (Musa acuminata): Genome-wide identification, classification and expression patterns

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, defense responses and metabolite accumulation. To date MYB family genes have not yet been comprehensively identified in the major staple fruit crop banana. In this study, we present a comprehensive, genome-wide analysis of the MYB genes from Musa acuminata DH-Pahang (A genome). A total of 285 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Organ- and development-specific expression patterns were determined from RNA-Seq data. For 280 M. acuminata MYB genes for which expression was found in at least one of the analysed samples, a variety of expression patterns were detected. The M. acuminata R2R3-MYB genes were functionally categorised, leading to the identification of seven clades containing only M. acuminata R2R3-MYBs. The encoded proteins may have specialised functions that were acquired or expanded in Musa during genome evolution. This functional classification and expression analysis of the MYB gene family in banana establishes a solid foundation for future comprehensive functional analysis of MaMYBs and can be utilized in banana improvement programmes.

Investigation of an Antioxidative System for Salinity Tolerance in Oenanthe javanica

Abiotic stress, such as drought and salinity, severely affect the growth and yield of many plants. Oenanthe javanica (commonly known as water dropwort) is an important vegetable that is grown in the saline-alkali soils of East Asia, where salinity is the limiting environmental factor. To study the defense mechanism of salt stress responses in water dropwort, we studied two water dropwort cultivars, V11E0022 and V11E0135, based on phenotypic and physiological indexes. We found that V11E0022 were tolerant to salt stress, as a result of good antioxidant defense system in the form of osmolyte (proline), antioxidants (polyphenols and flavonoids), and antioxidant enzymes (APX and CAT), which provided novel insights for salt-tolerant mechanisms. Then, a comparative transcriptomic analysis was conducted, and Gene Ontology (GO) analysis revealed that differentially expressed genes (DEGs) involved in the carbohydrate metabolic process could reduce oxidative stress and enhance energy production that can help in adaptation against salt stress. Similarly, lipid metabolic processes can also enhance tolerance against salt stress by reducing the transpiration rate, H₂O₂, and oxidative stress. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs involved in hormone signals transduction pathway promoted the activities of antioxidant enzymes and reduced oxidative stress; likewise, arginine and proline metabolism, and flavonoid pathways also stimulated the biosynthesis of proline and flavonoids, respectively, in response to salt stress. Moreover, transcription factors (TFs) were also identified, which play an important role in salt stress tolerance of water dropwort. The finding of this study will be helpful for crop improvement under salt stress.

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.

Genome-wide analyses of banana fasciclin-like AGP genes and their differential expression under low-temperature stress in chilling sensitive and tolerant cultivars

Thirty MaFLAs vary in their molecular features. MaFLA14/18/27/29 are likely to be involved in banana chilling tolerance by facilitating the cold signaling pathway and enhancing the cell wall biosynthesis. Although several studies have identified the molecular functions of individual fasciclin-like arabinogalactan protein (FLA) genes in plant growth and development, little information is available on their involvement in plant tolerance to low-temperature (LT) stress, and the related underlying mechanism is far from clear. In this study, the different expression of FLAs of banana (Musa acuminata) (MaFLAs) in the chilling-sensitive (CS) and chilling-tolerant (CT) banana cultivars under natural LT was investigated. Based on the latest banana genome database, a genome-wide identification of this gene family was done and the molecular features were analyzed. Thirty MaFLAs were distributed in 10 out of 11 chromosomes and these clustered into four major phylogenetic groups based on shared gene structure. Twenty-four MaFLAs contained N-terminal signal, 19 possessed predicted glycosylphosphatidylinositol (GPI), while 16 had both. Most MaFLAs were downregulated by LT stress. However, MaFLA14/18/29 were upregulated by LT in both cultivars with higher expression level recorded in the CT cultivar. Interestingly, MaFLA27 was significantly upregulated in the CT cultivar, but the opposite occurred for the CS cultivar. MaFLA27 possessed only N-terminal signal, MaFLA18 contained only GPI anchor, MaFLA29 possessed both, while MaFLA14 had neither. Thus, it was suggested that the accumulation of these FLAs in banana under LT could improve banana chilling tolerance through facilitating cold signal pathway and thereafter enhancing biosynthesis of plant cell wall components. The results provide background information of MaFLAs, suggest their involvement in plant chilling tolerance and their potential as candidate genes to be targeted when breeding CT banana.

MaMYB4, an R2R3-MYB Repressor Transcription Factor, Negatively Regulates the Biosynthesis of Anthocyanin in Banana

Anthocyanins spatiotemporally accumulate in certain tissues of particular species in the banana plant, and MYB transcription factors (TFs) serve as their primary regulators. However, the precise regulatory mechanism in banana remains to be determined. Here, we report the identification and characterization of MaMYB4, an R2R3-MYB repressor TF, characterized by the presence of EAR (ethylene-responsive element binding factor-associated amphiphilic repression) and TLLLFR motifs. MaMYB4 expression was induced by the accumulation of anthocyanins. Transgenic banana plants overexpressing MaMYB4 displayed a significant reduction in anthocyanin compared to wild type. Consistent with the above results, metabolome results showed that there was a decrease in all three identified cyanidins and one delphinidin, the main anthocyanins that determine the color of banana leaves, whereas both transcriptome and reverse transcription-quantitative polymerase chain reaction analysis showed that many key anthocyanin synthesis structural genes and TF regulators were downregulated in MaMYB4 overexpressors. Furthermore, dual-luciferase assays showed that MaMYB4 was able to bind to the CHS, ANS, DFR, and bHLH promoters, leading to inhibition of their expression. Yeast two-hybrid analysis verified that MaMYB4 did not interact with bHLH, which ruled out the possibility that MaMYB4 could be incorporated into the MYB-bHLH-WD40 complex. Our results indicated that MaMYB4 acts as a repressor of anthocyanin biosynthesis in banana, likely due to a two-level repression mechanism that consists of reduced expression of anthocyanin synthesis structural genes and the parallel downregulation of bHLH to interfere with the proper assembly of the MYB-bHLH-WD40 activation complex. To the best of our knowledge, this is the first MYB TF that regulates anthocyanin synthesis that was identified by genetic methods in bananas, which will be helpful for manipulating anthocyanin coloration in banana programs in the future.

Stable reference genes for RT-qPCR analysis of gene expression in the Musa acuminata-Pseudocercospora musae interaction

Leaf pathogens are limiting factors in banana (Musa spp.) production, with Pseudocercospora spp. responsible for the important Sigatoka disease complex. In order to investigate cellular processes and genes involved in host defence responses, quantitative real-time PCR (RT-qPCR) is an analytical technique for gene expression quantification. Reliable RT-qPCR data, however, requires that reference genes for normalization of mRNA levels in samples are validated under the conditions employed for expression analysis of target genes. We evaluated the stability of potential reference genes ACT1, α-TUB, UBQ1, UBQ2, GAPDH, EF1α, APT and RAN. Total RNA was extracted from leaf tissues of Musa acuminata genotypes Calcutta 4 (resistant) and Cavendish Grande Naine (susceptible), both subjected to P. musae infection. Expression stability was determined with NormFinder, BestKeeper, geNorm and RefFinder algorithms. UBQ2 and RAN were the most stable across all M. acuminata samples, whereas when considering inoculated and non-inoculated leaf samples, APT and UBQ2 were appropriate for normalization in Calcutta 4, with RAN and α-TUB most stable in Cavendish Grande Naine. This first study of reference genes for relative quantification of target gene expression in the M. acuminata-P. musae interaction will enable reliable analysis of gene expression in this pathosystem, benefiting elucidation of disease resistance mechanisms.

Survival Factor Gene FgSvf1 Is Required for Normal Growth and Stress Resistance in Fusarium graminearum

Survival factor 1 (Svf1) is a protein involved in cell survival pathways. In Saccharomyces cerevisiae, Svf1 is required for the diauxic growth shift and survival under stress conditions. In this study, we characterized the role of FgSvf1, the Svf1 homolog in the homothallic ascomycete fungus Fusarium graminearum. In the FgSvf1 deletion mutant, conidial germination was delayed, vegetative growth was reduced, and pathogenicity was completely abolished. Although the FgSvf1 deletion mutant produced perithecia, the normal maturation of ascospore was dismissed in deletion mutant. The FgSvf1 deletion mutant also showed reduced resistance to osmotic, fungicide, and cold stress and reduced sensitivity to oxidative stress when compared to the wild-type strain. In addition, we showed that FgSvf1 affects glycolysis, which results in the abnormal vegetative growth in the FgSvf1 deletion mutant. Further, intracellular reactive oxygen species (ROS) accumulated in the FgSvf1 deletion mutant, and this accumulated ROS might be related to the reduced sensitivity to oxidative stress and the reduced resistance to cold stress and fungicide stress. Overall, understanding the role of FgSvf1 in F. graminearum provides a new target to control F. graminearum infections in fields.

Deep expression analysis reveals distinct cold-response strategies in rubber tree (Hevea brasiliensis)

Background

Natural rubber, an indispensable commodity used in approximately 40,000 products, is fundamental to the tire industry. The rubber tree species Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell-Arg., which is native the Amazon rainforest, is the major producer of latex worldwide. Rubber tree breeding is time consuming, expensive and requires large field areas. Thus, genetic studies could optimize field evaluations, thereby reducing the time and area required for these experiments. In this work, transcriptome sequencing was used to identify a full set of transcripts and to evaluate the gene expression involved in the different cold-response strategies of the RRIM600 (cold-resistant) and GT1 (cold-tolerant) genotypes.

Results

We built a comprehensive transcriptome using multiple database sources, which resulted in 104,738 transcripts clustered in 49,304 genes. The RNA-seq data from the leaf tissues sampled at four different times for each genotype were used to perform a gene-level expression analysis. Differentially expressed genes (DEGs) were identified through pairwise comparisons between the two genotypes for each time series of cold treatments.

DEG annotation revealed that RRIM600 and GT1 exhibit different chilling tolerance strategies. To cope with cold stress, the RRIM600 clone upregulates genes promoting stomata closure, photosynthesis inhibition and a more efficient reactive oxygen species (ROS) scavenging system. The transcriptome was also searched for putative molecular markers (single nucleotide polymorphisms (SNPs) and microsatellites) in each genotype. and a total of 27,111 microsatellites and 202,949 (GT1) and 156,395 (RRIM600) SNPs were identified in GT1 and RRIM600. Furthermore, a search for alternative splicing (AS) events identified a total of 20,279 events.

Conclusions

The elucidation of genes involved in different chilling tolerance strategies associated with molecular markers and information regarding AS events provides a powerful tool for further genetic and genomic analyses of rubber tree breeding.

Electronic supplementary material

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

Identification, Expression, and Interaction Network Analyses of the CDPK Gene Family Reveal Their Involvement in the Development, Ripening, and Abiotic Stress Response in Banana

Calcium-dependent protein kinases (CDPKs) play vital roles in the regulation of plant growth, development, and tolerance to various abiotic stresses. However, little information is available for this gene family in banana. In this study, 44 CDPKs were identified in banana and were classified into four groups based on phylogenetic, gene structure, and conserved motif analyses. The majority of MaCDPKs generally exhibited similar expression patterns in the different tissues. Transcriptome analyses revealed that many CDPKs showed strong transcript accumulation at the early stages of fruit development and postharvest ripening in both varieties. Interaction network and co-expression analysis further identified some CDPKs-mediated network that was potentially active at the early stages of fruit development. Comparative expression analysis suggested that the high levels of CDPK expression in FJ might be related to its fast ripening characteristic. CDPK expression following the abiotic stress treatments indicated a significant transcriptional response to osmotic, cold, and salt treatment, as well as differential expression profiles, between BX and FJ. The findings of this study elucidate the transcriptional control of CDPKs in development, ripening, and the abiotic stress response in banana. Some tissue-specific, development/ripening-dependent, and abiotic stress-responsive candidate MaCDPK genes were identified for further genetic improvement of banana.

Evolution of Cold Acclimation and Its Role in Niche Transition in the Temperate Grass Subfamily Pooideae

The grass subfamily Pooideae dominates the grass floras in cold temperate regions and has evolved complex physiological adaptations to cope with extreme environmental conditions like frost, winter, and seasonality. One such adaptation is cold acclimation, wherein plants increase their frost tolerance in response to gradually falling temperatures and shorter days in the autumn. However, understanding how complex traits like cold acclimation evolve remains a major challenge in evolutionary biology. Here, we investigated the evolution of cold acclimation in Pooideae and found that a phylogenetically diverse set of Pooideae species displayed cold acclimation capacity. However, comparing differential gene expression after cold treatment in transcriptomes of five phylogenetically diverse species revealed widespread species-specific responses of genes with conserved sequences. Furthermore, we studied the correlation between gene family size and number of cold-responsive genes as well as between selection pressure on coding sequences of genes and their cold responsiveness. We saw evidence of protein-coding and regulatory sequence evolution as well as the origin of novel genes and functions contributing toward evolution of a cold response in Pooideae. Our results reflect that selection pressure resulting from global cooling must have acted on already diverged lineages. Nevertheless, conservation of cold-induced gene expression of certain genes indicates that the Pooideae ancestor may have possessed some molecular machinery to mitigate cold stress. Evolution of adaptations to seasonally cold climates is regarded as particularly difficult. How Pooideae evolved to transition from tropical to temperate biomes sheds light on how complex traits evolve in the light of climate changes.

Biological Control Agents Against Fusarium Wilt of Banana

In the last century, the banana crop and industry experienced dramatic losses due to an epidemic of Fusarium wilt of banana (FWB), caused by Fusarium oxysporum f.sp. cubense (Foc) race 1. An even more dramatic menace is now feared due to the spread of Foc tropical race 4. Plant genetic resistance is generally considered as the most plausible strategy for controlling effectively such a devastating disease, as occurred for the first round of FWB epidemic. Nevertheless, with at least 182 articles published since 1970, biological control represents a large body of knowledge on FWB. Remarkably, many studies deal with biological control agents (BCAs) that reached the field-testing stage and even refer to high effectiveness. Some selected BCAs have been repeatedly assayed in independent trials, suggesting their promising value. Overall under field conditions, FWB has been controlled up to 79% by using Pseudomonas spp. strains, and up to 70% by several endophytes and Trichoderma spp. strains. Lower biocontrol efficacy (42-55%) has been obtained with arbuscular mycorrhizal fungi, Bacillus spp., and non-pathogenic Fusarium strains. Studies on Streptomyces spp. have been mostly limited to in vitro conditions so far, with very few pot-experiments, and none conducted in the field. The BCAs have been applied with diverse procedures (e.g., spore suspension, organic amendments, bioformulations, etc.) and at different stages of plant development (i.e., in vitro, nursery, at transplanting, post-transplanting), but there has been no evidence for a protocol better than another. Nonetheless, new bioformulation technologies (e.g., nanotechnology, formulation of microbial consortia and/or their metabolites, etc.) and tailor-made consortia of microbial strains should be encouraged. In conclusion, the literature offers many examples of promising BCAs, suggesting that biocontrol can greatly contribute to limit the damage caused by FWB. More efforts should be done to further validate the currently available outcomes, to deepen the knowledge on the most valuable BCAs, and to improve their efficacy by setting up effective formulations, application protocols, and integrated strategies.

Transcriptome Analysis Reveals Key Cold-Stress-Responsive Genes in Winter Rapeseed (Brassica rapa L.)

Low ambient air temperature limits the growth and selection of crops in cold regions, and cold tolerance is a survival strategy for overwintering plants in cold winters. Studies of differences in transcriptional levels of winter rapeseed (Brassica rapa L.) under cold stress can improve our understanding of transcript-mediated cold stress responses. In this study, two winter rapeseed varieties, Longyou-7 (cold-tolerant) and Lenox (cold-sensitive), were used to reveal morphological, physiological, and transcriptome levels after 24 h of cold stress, and 24 h at room temperature, to identify the mechanism of tolerance to cold stress. Compared to Lenox, Longyou-7 has a shorter growth period and greater belowground mass, and exhibits stronger physiological activity after cold stress. Subsequently, more complete genomic annotation was obtained by sequencing. A total of 10,251 and 10,972 differentially expressed genes (DEG) were identified in Longyou-7 and Lenox, respectively. Six terms closely related to cold stress were found by the Gene Ontology (GO) function annotation. Some of these terms had greater upregulated expression in Longyou-7, and the expression of these genes was verified by qRT-PCR. Most of these DEGs are involved in phenylpropanoid biosynthesis, plant hormone signal transduction, ribosome biogenesis, MAPK signaling pathway, basal transcription factors, and photosynthesis. Analysis of the genes involved in the peroxisome pathway revealed that Longyou-7 and Lenox may have different metabolic patterns. Some transcription factors may play an important role in winter rapeseed tolerance to cold stress, and Longyou-7 is slightly slower than Lenox. Our results provide a transcriptome database and candidate genes for further study of winter rapeseed cold stress.

Low-temperature tolerance in land plants: Are transcript and membrane responses conserved?

Plants' tolerance of low temperatures is an economically and ecologically important limitation on geographic distributions and growing seasons. Tolerance for low temperatures varies significantly across different plant species, and different mechanisms likely act in different species. In order to survive low-temperature stress, plant membranes must maintain their fluidity in increasingly cold and oxidative cellular environments. The responses of different species to low-temperature stress include changes to the types and desaturation levels of membrane lipids, though the precise lipids affected tend to vary by species. Regulation of membrane dynamics and other low-temperature tolerance factors are controlled by both transcriptional and post-transcriptional mechanisms. Here, we review low-temperature induced changes in both membrane lipid composition and gene transcription across multiple related plant species with differing degrees of low-temperature tolerance. We attempt to define a core set of changes for transcripts and lipids across species and treatment variations. Some responses appear to be consistent across all species for which data are available, while many others appear likely to be species or family-specific. Potential rationales are presented, including variance in testing, reporting and the importance of considering the level of stress perceived by the plant.

Genome-wide identification and characterization of mRNAs and lncRNAs involved in cold stress in the wild banana (Musa itinerans)

Cold stress seriously affects banana growth, yield and fruit quality. Long noncoding RNAs (lncRNAs) have been demonstrated as key regulators of biotic and abiotic stress in plants, but the identification and prediction of cold responsive mRNAs and lncRNAs in wild banana remains unexplored. In present study, a cold resistant wild banana line from China was used to profile the cold-responsive mRNAs and lncRNAs by RNA-seq under cold stress conditions, i.e. 13°C (critical growth temperature), 4°C (chilling temperature), 0°C (freezing temperature) and normal growing condition, i.e. 28°C (control group). A total of 12,462 lncRNAs were identified in cold-stressed wild banana. In mRNA, much more alternative splicing events occurred in wild banana under the cold stress conditions compared with that in the normal growing condition. The GO analysis of differential expression genes (DEGs) showed the biochemical processes and membrane related genes responded positively to the cold stress. The KEGG pathway enrichment analysis of the DEGs showed that the pathways of photosynthesis, photosynthesis–antenna proteins, circadian rhythm–plant, glutathione metabolism, starch and sucrose metabolism, cutin/suberine/biosynthesis were altered or affected by the cold stress conditions. Our analyses of the generated transcriptome and lncRNAs provide new insights into regulating expression of genes and lncRNAs that respond to cold stress in the wild banana.

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

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

Deep sequencing-based characterization of transcriptome of Pyrus ussuriensis in response to cold stress

Pyrus ussuriensis is extremely cold hardy when fully acclimated, but knowledge relevant to the molecular mechanisms underlying this economically valuable trait is still limited so far. In this study, global transcriptome profiles of Pyrus ussuriensis under cold conditions (4 °C) over a time course were generated by high-throughput sequencing. In total, >57,121,199 high quality clean reads were obtained with approximately 11.0 M raw data for each library. Among them, the values of 66.84%-72.03% of clean reads in the digital transcript abundance measurement could be well mapped to the pear genome database, resulting in the identification of 8544 differentially expressed genes (DEGs) having 43 Gene Ontology (GO) terms and 17 clusters of orthologous groups (COG) involved in 385 Kyoto Encyclopedia of Genes and Genomes (KEGG) defined pathways. These comprised 3124 (1033 up-regulated, 2091 down-regulated), 1243 (729 up-regulated, 514 down-regulated), and 750 (458 up-regulated, 292 down-regulated) genes from the cold-treated samples at 5, 12 and 24 h, respectively. The accuracy of the RNA-Seq derived transcript expression data was validated by analyzing the expression patterns of 16 DGEs by quantitative real-time PCR. Plant-pathogen interaction, plant hormone signal transduction, Photosynthesis, signal transduction, innate immune response and response to biotic stimulus were the most significantly enriched GO categories among in the DEGs. A total of 335 transcription factors were shown to be cold responsive. In addition, a number of genes involved in the catabolism and signaling of hormones were significantly affected by the cold stress. The RNA-Seq and digital expression profiling provides valuable insights into the understanding the molecular events related to cold responses in Pyrus ussuriensis and dataset may help guide future identification and functional analysis of potential genes that are important for enhancing cold hardiness.

Transcriptomic changes in Cucurbita pepo fruit after cold storage: differential response between two cultivars contrasting in chilling sensitivity

Background

Zucchini fruit is susceptible to chilling injury (CI), but the response to low storage temperature is cultivar dependent. Previous reports about the response of zucchini fruit to chilling storage have been focused on the physiology and biochemistry of this process, with little information about the molecular mechanisms underlying it. In this work, we present a comprehensive analysis of transcriptomic changes that take place after cold storage in zucchini fruit of two commercial cultivars with contrasting response to chilling stress.

Results

RNA-Seq analysis was conducted in exocarp of fruit at harvest and after 14 days of storage at 4 and 20 °C. Differential expressed genes (DEGs) were obtained comparing fruit stored at 4 °C with their control at 20 °C, and then specific and common up and down-regulated DEGs of each cultivar were identified. Functional analysis of these DEGs identified similarities between the response of zucchini fruit to low temperature and other stresses, with an important number of GO terms related to biotic and abiotic stresses overrepresented in both cultivars. This study also revealed several molecular mechanisms that could be related to chilling tolerance, since they were up-regulated in cv. Natura (CI tolerant) or down-regulated in cv. Sinatra (CI sensitive). These mechanisms were mainly those related to carbohydrate and energy metabolism, transcription, signal transduction, and protein transport and degradation. Among DEGs belonging to these pathways, we selected candidate genes that could regulate or promote chilling tolerance in zucchini fruit including the transcription factors MYB76-like, ZAT10-like, DELLA protein GAIP, and AP2/ERF domain-containing protein.

Conclusions

This study provides a broader understanding of the important mechanisms and processes related to coping with low temperature stress in zucchini fruit and allowed the identification of some candidate genes that may be involved in the acquisition of chilling tolerance in this crop. These genes will be the basis of future studies aimed to identify markers involved in cold tolerance and aid in zucchini breeding programs.

Electronic supplementary material

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

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.

Early Cold-Induced Peroxidases and Aquaporins Are Associated With High Cold Tolerance in Dajiao (Musa spp. 'Dajiao')

Banana is an important tropical fruit with high economic value. One of the main cultivars ('Cavendish') is susceptible to low temperatures, while another closely related specie ('Dajiao') has considerably higher cold tolerance. We previously reported that some membrane proteins appear to be involved in the cold tolerance of Dajiao bananas via an antioxidation mechanism. To investigate the early cold stress response of Dajiao, here we applied comparative membrane proteomics analysis for both cold-sensitive Cavendish and cold-tolerant Dajiao bananas subjected to cold stress at 10°C for 0, 3, and 6 h. A total of 2,333 and 1,834 proteins were identified in Cavendish and Dajiao, respectively. Subsequent bioinformatics analyses showed that 692 Cavendish proteins and 524 Dajiao proteins were predicted to be membrane proteins, of which 82 and 137 differentially abundant membrane proteins (DAMPs) were found in Cavendish and Dajiao, respectively. Interestingly, the number of DAMPs with increased abundance following 3 h of cold treatment in Dajiao (80) was seven times more than that in Cavendish (11). Gene ontology molecular function analysis of DAMPs for Cavendish and Dajiao indicated that they belong to eight categories including hydrolase activity, binding, transporter activity, antioxidant activity, etc., but the number in Dajiao is twice that in Cavendish. Strikingly, we found peroxidases (PODs) and aquaporins among the protein groups whose abundance was significantly increased after 3 h of cold treatment in Dajiao. Some of the PODs and aquaporins were verified by reverse-transcription PCR, multiple reaction monitoring, and green fluorescent protein-based subcellular localization analysis, demonstrating that the global membrane proteomics data are reliable. By combining the physiological and biochemical data, we found that membrane-bound Peroxidase 52 and Peroxidase P7, and aquaporins (MaPIP1;1, MaPIP1;2, MaPIP2;4, MaPIP2;6, MaTIP1;3) are mainly involved in decreased lipid peroxidation and maintaining leaf cell water potential, which appear to be the key cellular adaptations contributing to the cold tolerance of Dajiao. This membrane proteomics study provides new insights into cold stress tolerance mechanisms of banana, toward potential applications for ultimate genetic improvement of cold tolerance in banana.

Integrated transcriptomics and metabolomics analysis to characterize cold stress responses in Nicotiana tabacum

BACKGROUND

CB-1 and K326 are closely related tobacco cultivars; however, their cold tolerance capacities are different. K326 is much more cold tolerant than CB-1.

RESULTS

We studied the transcriptomes and metabolomes of CB-1 and K326 leaf samples treated with cold stress. Totally, we have identified 14,590 differentially expressed genes (DEGs) in CB-1 and 14,605 DEGs in K326; there was also 200 differentially expressed metabolites in CB-1 and 194 in K326. Moreover, there were many overlapping genes (around 50%) that were cold-responsive in both plant cultivars, although there were also many differences in the cold responsive genes between the two cultivars. Importantly, for most of the overlapping cold responsive genes, the extent of the changes in expression were typically much more pronounced in K326 than in CB-1, which may help explain the superior cold tolerance of K326. Similar results were found in the metabolome analysis, particularly with the analysis of primary metabolites, including amino acids, organic acids, and sugars. The large number of specific responsive genes and metabolites highlight the complex regulatory mechanisms associated with cold stress in tobacco. In addition, our work implies that the energy metabolism and hormones may function distinctly between CB-1 and K326.

CONCLUSIONS

Differences in gene expression and metabolite levels following cold stress treatment seem likely to have contributed to the observed difference in the cold tolerance phenotype of these two tobacco cultivars.

Comparative physiological and transcriptomic analyses provide integrated insight into osmotic, cold, and salt stress tolerance mechanisms in banana

The growth, development, and production of banana plants are constrained by multiple abiotic stressors. However, it remains elusive for the tolerance mechanisms of banana responding to multiple abiotic stresses. In this study, we found that Fen Jiao (FJ) was more tolerant to osmotic, cold, and salt stresses than BaXi Jiao (BX) by phenotypic and physiological analyses. Comparative transcriptomic analyses highlighted stress tolerance genes that either specifically regulated in FJ or changed more than twofold in FJ relative to BX after treatments. In total, 933, 1644, and 133 stress tolerance genes were identified after osmotic, cold, and salt treatments, respectively. Further integrated analyses found that 30 tolerance genes, including transcription factor, heat shock protein, and E3 ubiquitin protein ligase, could be commonly regulated by osmotic, cold, and salt stresses. Finally, ABA and ROS signaling networks were found to be more active in FJ than in BX under osmotic, cold, and salt treatments, which may contribute to the strong stress tolerances of FJ. Together, this study provides new insights into the tolerance mechanism of banana responding to multiple stresses, thus leading to potential applications in the genetic improvement of multiple abiotic stress tolerances in banana.

Comparative Phosphoproteomics Reveals an Important Role of MKK2 in Banana (Musa spp.) Cold Signal Network

Low temperature is one of the key environmental stresses, which greatly affects global banana production. However, little is known about the global phosphoproteomes in Musa spp. and their regulatory roles in response to cold stress. In this study, we conducted a comparative phosphoproteomic profiling of cold-sensitive Cavendish Banana and relatively cold tolerant Dajiao under cold stress. Phosphopeptide abundances of five phosphoproteins involved in MKK2 interaction network, including MKK2, HY5, CaSR, STN7 and kinesin-like protein, show a remarkable difference between Cavendish Banana and Dajiao in response to cold stress. Western blotting of MKK2 protein and its T31 phosphorylated peptide verified the phosphoproteomic results of increased T31 phosphopeptide abundance with decreased MKK2 abundance in Daojiao for a time course of cold stress. Meanwhile increased expression of MKK2 with no detectable T31 phosphorylation was found in Cavendish Banana. These results suggest that the MKK2 pathway in Dajiao, along with other cold-specific phosphoproteins, appears to be associated with the molecular mechanisms of high tolerance to cold stress in Dajiao. The results also provide new evidence that the signaling pathway of cellular MKK2 phosphorylation plays an important role in abiotic stress tolerance that likely serves as a universal plant cold tolerance mechanism.

Breeding approaches and genomics technologies to increase crop yield under low-temperature stress

Improved knowledge about plant cold stress tolerance offered by modern omics technologies will greatly inform future crop improvement strategies that aim to breed cultivars yielding substantially high under low-temperature conditions. Alarmingly rising temperature extremities present a substantial impediment to the projected target of 70% more food production by 2050. Low-temperature (LT) stress severely constrains crop production worldwide, thereby demanding an urgent yet sustainable solution. Considerable research progress has been achieved on this front. Here, we review the crucial cellular and metabolic alterations in plants that follow LT stress along with the signal transduction and the regulatory network describing the plant cold tolerance. The significance of plant genetic resources to expand the genetic base of breeding programmes with regard to cold tolerance is highlighted. Also, the genetic architecture of cold tolerance trait as elucidated by conventional QTL mapping and genome-wide association mapping is described. Further, global expression profiling techniques including RNA-Seq along with diverse omics platforms are briefly discussed to better understand the underlying mechanism and prioritize the candidate gene (s) for downstream applications. These latest additions to breeders' toolbox hold immense potential to support plant breeding schemes that seek development of LT-tolerant cultivars. High-yielding cultivars endowed with greater cold tolerance are urgently required to sustain the crop yield under conditions severely challenged by low-temperature.

The ERF transcription factor family in cassava: genome-wide characterization and expression analyses against drought stress

Cassava (Manihot esculenta) shows strong tolerance to drought stress; however, the mechanisms underlying this tolerance are poorly understood. Ethylene response factor (ERF) family genes play a crucial role in plants responding to abiotic stress. Currently, less information is known regarding the ERF family in cassava. Herein, 147 ERF genes were characterized from cassava based on the complete genome data, which was further supported by phylogenetic relationship, gene structure, and conserved motif analyses. Transcriptome analysis suggested that most of the MeERF genes have similar expression profiles between W14 and Arg7 during organ development. Comparative expression profiles revealed that the function of MeERFs in drought tolerance may be differentiated in roots and leaves of different genotypes. W14 maintained strong tolerance by activating more MeERF genes in roots compared to Arg7 and SC124, whereas Arg7 and SC124 maintained drought tolerance by inducing more MeERF genes in leaves relative to W14. Expression analyses of the selected MeERF genes showed that most of them are significantly upregulated by osmotic and salt stresses, whereas slightly induced by cold stress. Taken together, this study identified candidate MeERF genes for genetic improvement of abiotic stress tolerance and provided new insights into ERF-mediated cassava tolerance to drought stress.

Proteome profiling reveals insights into cold-tolerant growth in sea buckthorn

Background

Low temperature is one of the crucial environmental factors limiting the productivity and distribution of plants. Sea buckthorn (Hippophae rhamnoides L.), a well recognized multipurpose plant species, live successfully in in cold desert regions. But their molecular mechanisms underlying cold tolerance are not well understood.

Methods

Physiological and biochemical responses to low-temperature stress were studied in seedlings of sea buckthorn. Differentially expressed protein spots were analyzed using multiplexing fluorescent two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) coupled with matrix-assisted laser desorption/ionization (MALDI) time-of-flight/time-of-flight (TOF/TOF) mass spectrometry (MS), the concentration of phytohormone was measured using enzyme-linked immunosorbent assay, and a spectrophotometric assay was used to measure enzymatic reactions.

Results

With the increase of cold stress intensity, the photosynthesis rate, transpiration rate, stomatal conductance in leaves and contents of abscisic acid (ABA) and indole acetic acid (IAA) in roots decreased significantly; however, water-use efficiency, ABA and zeatin riboside in leaves increased significantly, while cell membrane permeability, malondialdehyde and IAA in leaves increased at 7 d and then decreased at 14 d. DIGE and MS/MS analysis identified 32 of 39 differentially expressed protein spots under low-temperature stress, and their functions were mainly involved in metabolism, photosynthesis, signal transduction, antioxidative systems and post-translational modification.

Conclusion

The changed protein abundance and corresponding physiological–biochemical response shed light on the molecular mechanisms related to cold tolerance in cold-tolerant plants and provide key candidate proteins for genetic improvement of plants.

Electronic supplementary material

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

Transcriptome response to temperature stress in the wolf spider Pardosa pseudoannulata (Araneae: Lycosidae)

The wolf spider Pardosa pseudoannulata is a dominant predator in paddy ecosystem and an important biological control agent of rice pests. Temperature represents a primary factor influencing its biology and behavior, although the underlying molecular mechanisms remain unknown. To understand the response of P. pseudoannulata to temperature stress, we performed comparative transcriptome analyses of spider adults exposed to 10°C and 40°C for 12 h. We obtained 67,725 assembled unigenes, 21,765 of which were annotated in P. pseudoannulata transcriptome libraries, and identified 905 and 834 genes significantly up- or down-regulated by temperature stress. Functional categorization revealed the differential regulation of transcription, signal transduction, and metabolism processes. Calcium signaling pathway and metabolic pathway involving respiratory chain components played important roles in adapting to low temperature, whereas at high temperature, oxidative phosphorylation and amino acid metabolism were critical. Differentially expressed ribosomal protein genes contributed to temperature stress adaptation, and heat shock genes were significantly up-regulated. This study represents the first report of transcriptome identification related to the Araneae species in response to temperature stress. These results will greatly facilitate our understanding of the physiological and biochemical mechanisms of spiders in response to temperature stress.

Molecular Characterization of MaCCS, a Novel Copper Chaperone Gene Involved in Abiotic and Hormonal Stress Responses in Musa acuminata cv. Tianbaojiao

Copper/zinc superoxide dismutases (Cu/ZnSODs) play important roles in improving banana resistance to adverse conditions, but their activities depend on the copper chaperone for superoxide dismutase (CCS) delivering copper to them. However, little is known about CCS in monocots and under stress conditions. Here, a novel CCS gene (MaCCS) was obtained from a banana using reverse transcription PCR and rapid-amplification of cDNA ends (RACE) PCR. Sequence analyses showed that MaCCS has typical CCS domains and a conserved gene structure like other plant CCSs. Alternative transcription start sites (ATSSs) and alternative polyadenylation contribute to the mRNA diversity of MaCCS. ATSSs in MaCCS resulted in one open reading frame containing two in-frame start codons to form two protein versions, which is supported by the MaCCS subcellular localization of in both cytosol and chloroplasts. Furthermore, MaCCS promoter was found to contain many cis-elements associated with abiotic and hormonal responses. Quantitative real-time PCR analysis showed that MaCCS was expressed in all tested tissues (leaves, pseudostems and roots). In addition, MaCCS expression was significantly induced by light, heat, drought, abscisic acid and indole-3-acetic acid, but inhibited by relatively high concentrations of CuSO₄ and under cold treatment, which suggests that MaCCS is involved in abiotic and hormonal responses.

Differential root transcriptomics in a polyploid non-model crop: the importance of respiration during osmotic stress

To explore the transcriptomic global response to osmotic stress in roots, 18 mRNA-seq libraries were generated from three triploid banana genotypes grown under mild osmotic stress (5% PEG) and control conditions. Illumina sequencing produced 568 million high quality reads, of which 70-84% were mapped to the banana diploid reference genome. Using different uni- and multivariate statistics, 92 genes were commonly identified as differentially expressed in the three genotypes. Using our in house workflow to analyze GO enriched and underlying biochemical pathways, we present the general processes affected by mild osmotic stress in the root and focus subsequently on the most significantly overrepresented classes associated with: respiration, glycolysis and fermentation. We hypothesize that in fast growing and oxygen demanding tissues, mild osmotic stress leads to a lower energy level, which induces a metabolic shift towards (i) a higher oxidative respiration, (ii) alternative respiration and (iii) fermentation. To confirm the mRNA-seq results, a subset of twenty up-regulated transcripts were further analysed by RT-qPCR in an independent experiment at three different time points. The identification and annotation of this set of genes provides a valuable resource to understand the importance of energy sensing during mild osmotic stress.

Genome-Wide Identification, Phylogeny, and Expression Analyses of the 14-3-3 Family Reveal Their Involvement in the Development, Ripening, and Abiotic Stress Response in Banana

Plant 14-3-3 proteins act as critical components of various cellular signaling processes and play an important role in regulating multiple physiological processes. However, less information is known about the 14-3-3 gene family in banana. In this study, 25 14-3-3 genes were identified from the banana genome. Based on the evolutionary analysis, banana 14-3-3 proteins were clustered into ε and non-ε groups. Conserved motif analysis showed that all identified banana 14-3-3 genes had the typical 14-3-3 motif. The gene structure of banana 14-3-3 genes showed distinct class-specific divergence between the ε group and the non-ε group. Most banana 14-3-3 genes showed strong transcript accumulation changes during fruit development and postharvest ripening in two banana varieties, indicating that they might be involved in regulating fruit development and ripening. Moreover, some 14-3-3 genes also showed great changes after osmotic, cold, and salt treatments in two banana varieties, suggested their potential role in regulating banana response to abiotic stress. Taken together, this systemic analysis reveals the involvement of banana 14-3-3 genes in fruit development, postharvest ripening, and response to abiotic stress and provides useful information for understanding the functions of 14-3-3 genes in banana.

Comparative Transcriptome Analysis Reveals Heat-Responsive Genes in Chinese Cabbage (Brassica rapa ssp. chinensis)

Chinese cabbage (Brassica rapa ssp. chinensis) is an economically and agriculturally significant vegetable crop and is extensively cultivated throughout the world. Heat stress disturbs cellular homeostasis and causes visible growth inhibition of shoots and roots, severe retardation in growth and development, and even death. However, there are few studies on the transcriptome profiling of heat stress in non-heading Chinese cabbage. In this study, we investigated the transcript profiles of non-heading Chinese cabbage from heat-sensitive and heat-tolerant varieties "GHA" and "XK," respectively, in response to high temperature using RNA sequencing (RNA seq). Approximately 625 genes were differentially expressed between the two varieties. The responsive genes can be divided into three phases along with the time of heat treatment: response to stimulus, programmed cell death and ribosome biogenesis. Differentially expressed genes (DEGs) were identified in the two varieties, including transcription factors (TFs), kinases/phosphatases, genes related to photosynthesis and effectors of homeostasis. Many TFs were involved in the heat stress response of Chinese cabbage, including NAC069 TF which was up-regulated at all the heat treatment stages. And their expression levels were also validated by quantitative real-time-PCR (qRT-PCR). These candidate genes will provide genetic resources for further improving the heat-tolerant characteristics in non-heading Chinese cabbage.

Transcriptome profiling and validation of gene based single nucleotide polymorphisms (SNPs) in sorghum genotypes with contrasting responses to cold stress

BACKGROUND

Sorghum is a versatile cereal crop, with excellent heat and drought tolerance. However, it is susceptible to early-season cold stress (12-15 °C) which limits stand-establishment and seedling growth. To gain further insights on the molecular mechanism of cold tolerance in sorghum we performed transcriptome profiling between known cold sensitive and tolerant sorghum lines using RNA sequencing technology under control and cold stress treatments.

RESULTS

Here we report on the identification of differentially expressed genes (DEGs) between contrasting sorghum genotypes, HongkeZi (cold tolerant) and BTx623 (cold sensitive) under cool and control temperatures using RNAseq approach to elucidate the molecular basis of sorghum response to cold stress. Furthermore, we validated bi-allelic variants in the form of single nucleotide polymorphism (SNPs) between the cold susceptible and tolerant lines of sorghum. An analysis of transcriptome profile showed that in response to cold, a total of 1910 DEGs were detected under cold and control temperatures in both genotypes. We identified a subset of genes under cold stress for downstream analysis, including transcription factors that exhibit differential abundance between the sensitive and tolerant genotypes. We identified transcription factors including Dehydration-responsive element-binding factors, C-repeat binding factors, and Ethylene responsive transcription factors as significantly upregulated during cold stress in cold tolerant HKZ. Additionally, specific genes such as plant cytochromes, glutathione s-transferases, and heat shock proteins were found differentially regulated under cold stress between cold tolerant and susceptible genotype of sorghum. A total of 41,603 SNP were identified between the cold sensitive and tolerant genotypes with minimum read of four. Approximately 89 % of the 114 SNP sites selected for evaluation were validated using endpoint genotyping technology.

CONCLUSION

A new strategy which involved an integrated analysis of differential gene expression and identification of bi-allelic single nucleotide polymorphism (SNP) was conducted to determine and analyze differentially expressed genes and variation involved in cold stress response of sorghum. The results gathered provide an insight into the complex mechanisms associated with cold response in sorghum, which involve an array of transcription factors and genes which were previously related to abiotic stress response. This study also offers resource for gene based SNP that can be applied towards targeted genomic studies of cold tolerance in sorghum and other cereal crops.