Transcriptomic and metabolomic profiling of Camellia sinensis L. cv. 'Suchazao' exposed to temperature stresses reveals modification in protein synthesis and photosynthetic and anthocyanin biosynthetic pathways

Physiological and transcriptomic analyses reveal the regulatory mechanisms of Anoectochilus roxburghii in response to high-temperature stress

BACKGROUND

High temperatures significantly affect the growth, development, and yield of plants. Anoectochilus roxburghii prefers a cool and humid environment, intolerant of high temperatures. It is necessary to enhance the heat tolerance of A. roxburghii and breed heat-tolerant varieties. Therefore, we studied the physiological indexes and transcriptome of A. roxburghii under different times of high-temperature stress treatments.

RESULTS

Under high-temperature stress, proline (Pro), H2O2 content increased, then decreased, then increased again, catalase (CAT) activity increased continuously, peroxidase (POD) activity decreased rapidly, then increased, then decreased again, superoxide dismutase (SOD) activity, malondialdehyde (MDA), and soluble sugars (SS) content all decreased, then increased, and chlorophyll and soluble proteins (SP) content increased, then decreased. Transcriptomic investigation indicated that a total of 2740 DEGs were identified and numerous DEGs were notably enriched for "Plant-pathogen interaction" and "Plant hormone signal transduction". We identified a total of 32 genes in these two pathways that may be the key genes for resistance to high-temperature stress in A. roxburghii.

CONCLUSIONS

To sum up, the results of this study provide a reference for the molecular regulation of A. roxburghii's tolerance to high temperatures, which is useful for further cultivation of high-temperature-tolerant A. roxburghii varieties.

Physiological and Transcriptome Analyses Reveal the Protective Effect of Exogenous Trehalose in Response to Heat Stress in Tea Plant (Camellia sinensis)

It is well known that application of exogenous trehalose can enhance the heat resistance of plants. To investigate the underlying molecular mechanisms by which exogenous trehalose induces heat resistance in C. sinensis, a combination of physiological and transcriptome analyses was conducted. The findings revealed a significant increase in the activity of superoxide dismutase (SOD) and peroxidase (POD) upon treatment with 5.0 mM trehalose at different time points. Moreover, the contents of proline (PRO), endogenous trehalose, and soluble sugar exhibited a significant increase, while malondialdehyde (MDA) content decreased following treatment with 5.0 mM trehalose under 24 h high-temperature stress (38 °C/29 °C, 12 h/12 h). RNA-seq analysis demonstrated that the differentially expressed genes (DEGs) were significantly enriched in the MAPK pathway, plant hormone signal transduction, phenylpropanoid biosynthesis, flavone and flavonol biosynthesis, flavonoid biosynthesis, and the galactose metabolism pathway. The capability to scavenge free radicals was enhanced, and the expression of a heat shock factor gene (HSFB2B) and two heat shock protein genes (HSP18.1 and HSP26.5) were upregulated in the tea plant. Consequently, it was concluded that exogenous trehalose contributes to alleviating heat stress in C. sinensis. Furthermore, it regulates the expression of genes involved in diverse pathways crucial for C. sinensis under heat-stress conditions. These findings provide novel insights into the molecular mechanisms underlying the alleviation of heat stress in C. sinensis with trehalose.

For a Colorful Life: Recent Advances in Anthocyanin Biosynthesis during Leaf Senescence

Leaf senescence is the last stage of leaf development, and it is accompanied by a leaf color change. In some species, anthocyanins are accumulated during leaf senescence, which are vital indicators for both ornamental and commercial value. Therefore, it is essential to understand the molecular mechanism of anthocyanin accumulation during leaf senescence, which would provide new insight into autumn coloration and molecular breeding for more colorful plants. Anthocyanin accumulation is a surprisingly complex process, and significant advances have been made in the past decades. In this review, we focused on leaf coloration during senescence. We emphatically discussed several networks linked to genetic, hormonal, environmental, and nutritional factors in regulating anthocyanin accumulation during leaf senescence. This paper aims to provide a regulatory model for leaf coloration and to put forward some prospects for future development.

The effect of intercropping leguminous green manure on theanine accumulation in the tea plant: A metagenomic analysis

Intercropping is a widely recognised technique that contributes to agricultural sustainability. While intercropping leguminous green manure offers advantages for soil health and tea plants growth, the impact on the accumulation of theanine and soil nitrogen cycle are largely unknown. The levels of theanine, epigallocatechin gallate and soluble sugar in tea leaves increased by 52.87% and 40.98%, 22.80% and 6.17%, 22.22% and 29.04% in intercropping with soybean-Chinese milk vetch rotation and soybean alone, respectively. Additionally, intercropping significantly increased soil amino acidnitrogen content, enhanced extracellular enzyme activities, particularly β-glucosidase and N-acetyl-glucosaminidase, as well as soil multifunctionality. Metagenomics analysis revealed that intercropping positively influenced the relative abundances of several potentially beneficial microorganisms, including Burkholderia, Mycolicibacterium and Paraburkholderia. Intercropping resulted in lower expression levels of nitrification genes, reducing soil mineral nitrogen loss and N2 O emissions. The expression of nrfA/H significantly increased in intercropping with soybean-Chinese milk vetch rotation. Structural equation model analysis demonstrated that the accumulation of theanine in tea leaves was directly influenced by the number of intercropping leguminous green manure species, soil ammonium nitrogen and amino acid nitrogen. In summary, the intercropping strategy, particularly intercropping with soybean-Chinese milk vetch rotation, could be a novel way for theanine accumulation.

Comprehensive analysis of Dendrobium catenatum HSP20 family genes and functional characterization of DcHSP20-12 in response to temperature stress

Heat shock proteins (HSPs) are a class of protective proteins in response to abiotic stress in plants, and HSP20 plays an essential role in response to temperature stress. However, there are few studies on HSP20 in Dendrobium catenatum. In this study, 18 DcHSP20 genes were identified from the D. catenatum genome. Phylogenetic analysis showed that DcHSP20s could be classified into six subgroups, each member of which has similar conserved motifs and gene structures. Gene expression analysis of 18 DcHSP20 genes revealed that they exhibited variable expression patterns in different plant tissues. Meanwhile, all 18 DcHSP20 genes were induced to be up-regulated under high temperature, while six genes (DcHSP20-2/9/10/12/16/17) were significantly up-regulated under low temperature. Moreover, combining gene expression under high and low temperature stress, the DcHSP20-12 gene was cloned for functional analysis. The germination ratios, fresh weights, root lengths of two DcHSP20-12-overexpressing transgenic Arabidopsis thaliana lines were significantly higher, but MDA contents were lower than that of wild-type (WT) plants under heat and cold stresses, displayed enhanced thermotolerance and cold-resistance. These results lay a foundation for the functional characterization of DcHSP20s and provide a candidate gene, DcHSP20-12, for improving the tolerance of D. catenatum to temperature stress in the future.

The light-harvesting chlorophyll a/b-binding proteins of photosystem II family members are responsible for temperature sensitivity and leaf color phenotype in albino tea plant

INTRODUCTION

Light-harvesting chlorophyll a/b-binding (LHCB) protein complexes of photosystem II are integral to the formation of thylakoid structure and the photosynthetic process. They play an important role in photoprotection, a crucial process in leaf development under low-temperature stress. Nonetheless, potential key genes directly related to low-temperature response and albino phenotype have not been precisely identified in tea plant. Moreover, there are no studies simultaneously investigating multiple albino tea cultivars with different temperature sensitivity.

OBJECTIVES

The study aimed to clarify the basic characteristics of CsLHCB gene family members, and identify critical CsLHCB genes potentially influential in leaf color phenotypic variation and low-temperature stress response by contrasting green and albino tea cultivars. Concurrently, exploring the differential expression of the CsLHCB gene family across diverse temperature-sensitive albino tea cultivars.

METHODS

We identified 20 putative CsLHCB genes according to phylogenetic analysis. Evolutionary relationships, gene duplication, chromosomal localization, and structures were analyzed by TBtools; the physiological and biochemical characteristics were analyzed by protein analysis websites; the differences in coding sequences and protein accumulation in green and albino tea cultivars, gene expression with maturity were tested by molecular biology technology; and protein interaction was analyzed in the STRING database.

RESULTS

All genes were categorized into seven groups, mapping onto 7 chromosomes, including three tandem and one segmental duplications. They all own a conserved chlorophyll A/B binding protein domain. The expression of CsLHCB genes was tissue-specific, predominantly in leaves. CsLHCB5 may play a key role in the process of leaf maturation and senescence. In contrast to CsLHCB5, CsLHCB1.1, CsLHCB2, and CsLHCB3.2 were highly conserved in amino acid sequence between green and albino tea cultivars. In albino tea cultivars, unlike in green cultivars, the expression of CsLHCB1.1, CsLHCB1.2, and CsLHCB2 was down-regulated under low-temperature stress. The accumulation of CsLHCB1 and CsLHCB5 proteins was lower in albino tea cultivars. Greater accumulation of CsLHCB2 protein was detected in RX1 and RX2 compared to other albino cultivars.

CONCLUSIONS

CsLHCB1.1, CsLHCB1.2, and CsLHCB2 played a role in the response to low-temperature stress. The amino acid sequence site mutation of CsLHCB5 would distinguish the green and albino tea cultivars. The less accumulation of CsLHCB1 and CsLHCB5 had a potential influence on albino leaves. Albino cultivars more sensitive to temperature exhibited lower CsLHCB gene expression. CsLHCB2 may serve as an indicator of temperature sensitivity differences in albino tea cultivars. This study could provide a reference for further studies of the functions of the CsLHCB family and contribute to research on the mechanism of the albino in tea plant.

Comprehensive applications of metabolomics on tea science and technology: Opportunities, hurdles, and perspectives

With the development of metabolomics analytical techniques, relevant studies have increased in recent decades. The procedures of metabolomics analysis mainly include sample preparation, data acquisition and pre-processing, multivariate statistical analysis, as well as maker compounds' identification. In the present review, we summarized the published articles of tea metabolomics regarding different analytical tools, such as mass spectrometry, nuclear magnetic resonance, ultraviolet-visible spectrometry, and Fourier transform infrared spectrometry. The metabolite variation of fresh tea leaves with different treatments, such as biotic/abiotic stress, horticultural measures, and nutritional supplies was reviewed. Furthermore, the changes of chemical composition of processed tea samples under different processing technologies were also profiled. Since the identification of critical or marker metabolites is a complicated task, we also discussed the procedure of metabolite identification to clarify the importance of omics data analysis. The present review provides a workflow diagram for tea metabolomics research and also the perspectives of related studies in the future.

Physiological and transcriptomic analysis reveals the potential mechanism of Morinda officinalis How in response to freezing stress

BACKGROUND

Morinda officinalis How (MO) is a vine shrub distributed in tropical and subtropical regions, known as one of the "Four Southern Herbal Medicines" in China. The unclear responsive mechanism by which MO adapt to freezing stress limits progress in molecular breeding for MO freezing tolerance.

RESULTS

In this study, morphological, physiological and microstructure changes in MO exposed to -2℃ for 0 h, 3 h, 8 h and 24 h were comprehensively characterized. The results showed that freezing stress caused seedling dehydration, palisade cell and spongy mesophyll destruction. A significant increase in the content of proline, soluble protein and soluble sugars, as well as the activity of superoxide dismutase and peroxidase was observed. Subsequently, we analyzed the transcriptomic changes of MO leaves at different times under freezing treatment by RNA-seq. A total of 24,498 unigenes were annotated and 3252 unigenes were identified as differentially expressed genes (DEGs). Most of these DEGs were annotated in starch and sucrose metabolism, plant hormone signal transduction and MAPK signaling pathways. Family Enrichment analysis showed that the glucosyl/glucuronosyl transferases, oxidoreductase, chlorophyll a/b binding protein and calcium binding protein families were significantly enriched. We also characterized 7 types of transcription factors responding to freezing stress, among which the most abundant family was the MYBs, followed by the AP2/ERFs and NACs. Furthermore, 10 DEGs were selected for qRT-PCR analysis, which validated the reliability and accuracy of RNA-seq data.

CONCLUSIONS

Our results provide an overall view of the dynamic changes in physiology and insight into the molecular regulation mechanisms of MO in response to freezing stress. This study will lay a foundation for freezing tolerance molecular breeding and improving the quality of MO.

Flavonoid metabolites in tea plant (Camellia sinensis) stress response: Insights from bibliometric analysis

In the context of global climate change, tea plants are at risk from elevating environmental stress factors. Coping with this problem relies upon the understanding of tea plant stress response and its underlying mechanisms. Over the past two decades, research in this field has prospered with the contributions of scientists worldwide. Aiming in providing a comprehensive perspective of the research field related to tea plant stress response, we present a bibliometric analysis of the this area. Our results demonstrate the most studied stresses, global contribution, authorship and collaboration, and trending research topics. We highlight the importance of flavonoid metabolites in tea plant stress response, particularly their role in maintaining redox homeostasis, yield, and adjusting tea quality under stress conditions. Further research on the flavonoid response under various stress conditions can promote the development of cultivation measures, thereby improving stress resistance and tea quality.

Anthocyanin metabolism and its differential regulation in purple tea (Camellia sinensis)

Tea plants (Camellia sinensis) typically contain high-flavonoid phytochemicals like catechins. Recently, new tea cultivars with unique purple-colored leaves have gained attention. These purple tea cultivars are enriched with anthocyanin, which provides an interesting perspective for studying the metabolic flux of the flavonoid pathway. An increasing number of studies are focusing on the leaf color formation of purple tea and this review aims to summarize the latest progress made on the composition and accumulation of anthocyanins in tea plants. In addition, the regulation mechanism in its synthesis will be discussed and a hypothetical regulation model for leaf color transformation during growth will be proposed. Some novel insights are presented to facilitate future in-depth studies of purple tea to provide a theoretical basis for targeted breeding programs in leaf color.

The CsHSFA-CsJAZ6 module-mediated high temperature regulates flavonoid metabolism in Camellia sinensis

High temperatures (HTs) seriously affect the yield and quality of tea. Catechins, derived from the flavonoid pathway, are characteristic compounds that contribute to the flavour of tea leaves. In this study, we first showed that the flavonoid content of tea leaves was significantly reduced under HT conditions via metabolic profiles; and then demonstrated that two transcription factors, CsHSFA1b and CsHSFA2 were activated by HT and negatively regulate flavonoid biosynthesis during HT treatment. Jasmonate (JA), a defensive hormone, plays a key role in plant adaption to environmental stress. However, little has been reported on its involvement in HT response in tea. Herein, we demonstrated that CsHSFA1b and CsHSFA2 activate CsJAZ6 expression through directly binding to heat shock elements in its promoter, and thereby repress the JA pathway. Most secondary metabolites are regulated by JA, including catechin in tea. Our study reported that CsJAZ6 directly interacts with CsEGL3 and CsTTG1 and thereby reduces catechin accumulation. From this, we proposed a CsHSFA-CsJAZ6-mediated HT regulation model of catechin biosynthesis. We also determined that negative regulation of the JA pathway by CsHSFAs and its homologues is conserved in Arabidopsis. These findings broaden the applicability of the regulation of JAZ by HSF transcription factors and further suggest the JA pathway as a valuable candidate for HT-resistant breeding and cultivation.

Anthocyanins Profiling Analysis and RNA-Seq Revealed the Dominating Pigments and Coloring Mechanism in Cyclamen Flowers

Pigments in cyclamen (Cyclamen purpurascens) endows flowers with great ornamental and medicinal values. However, little is known about the biosynthetic pathways of pigments, especially anthocyanins, in cyclamen flowers. Herein, anthocyanins profiling and RNA-Seq were used to decipher the molecular events using cyclamen genotypes of red (HXK) or white (BXK) flowers. We found that red cyclamen petals are rich in cyanidin-3-O-rutinoside, cyanidin-3-O-glucoside, delphinidin-3-O-glucoside, malvidin-3-O-glucoside, peonidin-3-O-rutinoside, quercetin-3-O-glucoside, and ruti. In addition, our transcriptomics data revealed 3589 up-regulated genes and 2788 down-regulated genes comparing the BXK to HXK. Our rich dataset also identified eight putative key genes for anthocyanin synthesis, including four chalcone synthase (CHS, g13809_i0, g12097_i0, g18851_i0, g36714_i0), one chalcone isomerase (CHI, g26337_i0), two flavonoid 3-hydroxylase (F3'H, g14710_i0 and g15005_i0), and one anthocyanidin synthase (ANS, g18981_i0). Importantly, we found a 2.5 order of magnitude higher expression of anthocyanin 3-O-glucosyltransferase (g8206_i0), which encodes a key gene in glycosylation of anthocyanins, in HXK compared to BXK. Taken together, our multiomics approach demonstrated massive changes in gene regulatory networks and anthocyanin metabolism in controlling cyclamen flower color.

Integrative Transcriptomic and Phytohormonal Analyses Provide Insights into the Cold Injury Recovery Mechanisms of Tea Leaves

Tea plant is susceptible to low temperature, while the cold injury recovery mechanisms of tea leaves are still unclear. Windbreak has an effective and gradient range of protecting tea plants. Tea plants with increasing cold damage degree have varying recovery status accordingly, which are the ideal objects for investigating the cold injury recovery mechanisms of tea leaves. Here, we investigated the transcriptome and phytohormone profiles of tea leaves with different cold injury degrees in recovery (adjacent to the windbreak), and the levels of chlorophylls, malondialdehyde, major phytohormones as well as the activities of peroxidase (POD) and superoxide dismutase (SOD) were also measured. The results showed the content of total chlorophylls and the activity of POD in mature tea leaves gradually decreased with the distance to windbreak, while SOD showed the opposite. The major phytohormones were highly accumulated in the moderately cold-injured tea leaves. The biosynthesis of abscisic acid (ABA) was enhanced in the moderate cold damaged tea leaves, suggesting that ABA plays an important role in the cold response and resistance of tea plants. The transcriptomic result showed that the samples in different rows were well discriminated, and the pathways of plant-pathogen interaction and flavonoid biosynthesis were enriched based on KEGG analysis. WRKY, GRAS and NAC were the top classes of transcription factors differentially expressed in the different cold-injured tea leaves. Thus, windbreak is effective to protect adjacent tea plants from cold wave, and phytohormones importantly participate in the cold injury recovery of tea leaves.

Genome-Wide Identification and Expression Analysis of Isopentenyl transferase Family Genes during Development and Resistance to Abiotic Stresses in Tea Plant (Camellia sinensis)

The tea plant is an important economic crop and is widely cultivated. Isopentenyl transferase (IPT) is the first and rate-limiting enzyme of cytokinin (CK) signaling, which plays key roles in plant development and abiotic stress. However, the IPT gene family in tea plants has not been systematically investigated until now. The phylogenetic analyses, gene structures, and conserved domains were predicted here. The results showed that a total of 13 CsIPT members were identified from a tea plant genome database and phylogenetically classified into four groups. Furthermore, 10 CsIPT members belonged to plant ADP/ATP-IPT genes, and 3 CsIPTs were tRNA-IPT genes. There is a conserved putative ATP/GTP-binding site (P-loop motif) in all the CsIPT sequences. Based on publicly available transcriptome data as well as through RNA-seq and qRT-PCR analysis, the CsIPT genes which play key roles in the development of different tissues were identified, respectively. Furthermore, CsIPT6.2 may be involved in the response to different light treatments. CsIPT6.4 may play a key role during the dormancy and flush of the lateral buds. CsIPT5.1 may play important regulatory roles during the development of the lateral bud, leaf, and flower. CsIPT5.2 and CsIPT6.2 may both play key roles for increased resistance to cold-stress, whereas CsIPT3.2 may play a key role in improving resistance to high-temperature stress as well as drought-stress and rewatering. This study could provide a reference for further studies of CsIPT family’s functions and could contribute to tea molecular breeding.

Anthocyanins, Carotenoids and Chlorophylls in Edible Plant Leaves Unveiled by Tandem Mass Spectrometry

Natural pigments are a quite relevant group of molecules that are widely distributed in nature, possessing a significant role in our daily lives. Besides their colors, natural pigments are currently recognized as having relevant biological properties associated with health benefits, such as anti-tumor, anti-atherogenicity, anti-aging and anti-inflammatory activities, among others. Some of these compounds are easily associated with specific fruits (such as blueberries with anthocyanins, red pitaya with betalain or tomato with lycopene), vegetables (carrots with carotenoids), plant leaves (chlorophylls in green leaves or carotenoids in yellow and red autumn leaves) and even the muscle tissue of vertebrates (such as myoglobin). Despite being less popular as natural pigment sources, edible plant leaves possess a high variety of chlorophylls, as well as a high variety of carotenoids and anthocyanins. The purpose of this review is to critically analyze the whole workflow employed to identify and quantify the most common natural pigments (anthocyanin, carotenoids and chlorophylls) in edible plant leaves using tandem mass spectrometry. Across the literature there, is a lack of consistency in the methods used to extract and analyze these compounds, and this review aims to surpass this issue. Additionally, mass spectrometry has stood out in the context of metabolomics, currently being a widely employed technique in this field. For the three pigments classes, the following steps will be scrutinized: (i) sample pre-preparation, including the solvents and extraction conditions; (ii) details of the chromatographic separation and mass spectrometry experiments (iii) pigment identification and quantification.

Integrated sRNA-seq and RNA-seq Analyses Reveal a microRNA Regulation Network Involved in Cold Response in Pisum sativum L

(1) Background: Cold stress affects growth and development in plants and is a major environmental factor that decreases productivity. Over the past two decades, the advent of next generation sequencing (NGS) technologies has opened new opportunities to understand the molecular bases of stress resistance by enabling the detection of weakly expressed transcripts and the identification of regulatory RNAs of gene expression, including microRNAs (miRNAs). (2) Methods: In this study, we performed time series sRNA and mRNA sequencing experiments on two pea (Pisum sativum L., Ps) lines, Champagne frost-tolerant and Térèse frost-sensitive, during a low temperature treatment versus a control condition. (3) Results: An integrative analysis led to the identification of 136 miRNAs and a regulation network composed of 39 miRNA/mRNA target pairs with discordant expression patterns. (4) Conclusions: Our findings indicate that the cold response in pea involves 11 miRNA families as well as their target genes related to antioxidative and multi-stress defense mechanisms and cell wall biosynthesis.

Comparative Transcriptome and Proteome Analysis Provides New Insights Into the Mechanism of Protein Synthesis in Kenaf (Hibiscus cannabinus L.) Leaves

Given the rising domestic demand and increasing global prices of corn and soybean, China is looking for alternatives for these imports to produce animal fodder. Kenaf (Hibiscus cannabinus L.) has great potential as a new forage source, due to abundant proteins, phenols and flavonoids in its leaves. However, few studies have evaluated the mechanism of protein synthesis in kenaf leaves. In the current work, compared with kenaf material "L332," the percentage of crude protein content in leaves of material "Q303" increased by 6.13%; combined with transcriptome and proteome data, the kenaf samples were systematically studied to obtain mRNA-protein correlation. Then, the genes/proteins related to protein synthesis in the kenaf leaves were obtained. Moreover, this work detected mRNA expression of 20 differentially expressed genes (DEGs). Meanwhile, 20 differentially expressed proteins (DEPs) related to protein synthesis were performed parallel reaction monitoring. Fructose-1,6-bisphosphatase (FBP), nitrite reductase (NirA), prolyl tRNA synthase (PARS) and glycine dehydrogenase (GLDC) presented increased mRNA and protein levels within kenaf leaves with high protein content. Based on the obtained findings, FBP, NirA, PARS, and GLDC genes may exert a vital function in the protein synthesis of kenaf leaves. The results provide a new idea for further studying the potential genes affecting the quality trait of protein content in kenaf leaves and provide gene resources and a theoretical foundation for further cultivating high protein kenaf varieties.

Exogenous Melatonin Enhances Photosynthetic Capacity and Related Gene Expression in A Dose-Dependent Manner in the Tea Plant (Camellia sinensis (L.) Kuntze)

The enhancement of photosynthesis of tea leaves can increase tea yield. In order to explore the regulation mechanism of exogenous melatonin (MT) on the photosynthetic characteristics of tea plants, tea variety 'Zhongcha 108' was used as the experimental material in this study. The effects of different concentrations (0, 0.2, 0.3, 0.4 mM) of melatonin on the chlorophyll (Chl) content, stomatal opening, photosynthetic and fluorescence parameters, antioxidant enzyme activity, and related gene expression of tea plants were detected and analyzed. The results showed that under 0.2-mM MT treatment, chlorophyll (Chl) content, photosynthetic rate (P_n), stomatal conductance (G_s), intercellular CO₂ concentration (C_i), and transpiration rate (T_r) improved, accompanied by a decrease in stomata density and increase in stomata area. Zero point two millimolar MT increased Chl fluorescence level and superoxide dismutase (SOD) activity, and reduced hydrogen peroxide (H₂O₂) and malondialdehyde (MDA) contents, indicating that MT alleviated PSII inhibition and improved photochemical efficiency. At the same time, 0.2 mM MT induced the expression of genes involved in photosynthesis and chlorophyll metabolism to varying degrees. The study demonstrated that MT can effectively enhance the photosynthetic capacity of tea plants in a dose-dependent manner. These results may promote a comprehensive understanding of the potential regulatory mechanism of exogenous MT on photosynthesis in tea plants.

Tea (Camellia sinensis): A Review of Nutritional Composition, Potential Applications, and Omics Research

Tea (Camelliasinensis) is the world’s most widely consumed non-alcoholic beverage with essential economic and health benefits since it is an excellent source of polyphenols, catechins, amino acids, flavonoids, carotenoids, vitamins, and polysaccharides. The aim of this review is to summarize the main secondary metabolites in tea plants, and the content and distribution of these compounds in six different types of tea and different organs of tea plant were further investigated. The application of these secondary metabolites on food processing, cosmetics industry, and pharmaceutical industry was reviewed in this study. With the rapid advancements in biotechnology and sequencing technology, omics analyses, including genome, transcriptome, and metabolome, were widely used to detect the main secondary metabolites and their molecular regulatory mechanisms in tea plants. Numerous functional genes and regulatory factors have been discovered, studied, and applied to improve tea plants. Research advances, including secondary metabolites, applications, omics research, and functional gene mining, are comprehensively reviewed here. Further exploration and application trends are briefly described. This review provides a reference for basic and applied research on tea plants.

Molecular and Metabolic Changes under Environmental Stresses: The Biosynthesis of Quality Components in Preharvest Tea Shoots

Severe environments impose various abiotic stresses on tea plants. Although much is known about the physiological and biochemical responses of tea (Camellia sinensis L.) shoots under environmental stresses, little is known about how these stresses impact the biosynthesis of quality components. This review summarizes and analyzes the changes in molecular and quality components in tea shoots subjected to major environmental stresses during the past 20 years, including light (shade, blue light, green light, and UV-B), drought, high/low temperature, CO2, and salinity. These studies reveal that carbon and nitrogen metabolism is critical to the downstream biosynthesis of quality components. Based on the molecular responses of tea plants to stresses, a series of artificial methods have been suggested to treat the pre-harvest tea plants that are exposed to inhospitable environments to improve the quality components in shoots. Furthermore, many pleiotropic genes that are up- or down-regulated under both single and concurrent stresses were analyzed as the most effective genes for regulating multi-resistance and quality components. These findings deepen our understanding of how environmental stresses affect the quality components of tea, providing novel insights into strategies for balancing plant resistance, growth, and quality components in field-based cultivation and for breeding plants using pleiotropic genes.

Differential Physiological, Transcriptomic, and Metabolomic Responses of Paspalum wettsteinii Under High-Temperature Stress

Global warming has far-reaching effects on plant growth and development. As a warm-season forage grass, Paspalum wettsteinii is highly adaptable to high temperatures. However, the response mechanism of P. wettsteinii under high-temperature stress is still unclear. Therefore, we investigated the physiological indicators, transcriptome and metabolome of P. wettsteinii under different heat stress treatments. Plant height, the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), and the contents of soluble sugar, proline, chlorophyll a, and chlorophyll b increased and then decreased, while the malondialdehyde (MDA) content decreased and then increased with increasing heat stress. Transcriptomic analysis revealed that genes related to energy and carbohydrate metabolism, heat shock proteins (HSPs), and transcription factors (TFs), secondary metabolite biosynthesis and the antioxidant system significantly changed to varying degrees. Metabolomic analysis showed that only free fatty acids were downregulated, while amino acids and their derivatives, organic acids, flavonoids, and sugars were both up- and downregulated under heat stress. These combined analyses revealed that growth was promoted at 25-40°C, while at 45°C, excess reactive oxygen species (ROS) damage reduced antioxidant and osmoregulatory effects and inactivated genes associated with the light and electron transport chains (ETCs), as well as damaged the PS II system and inhibited photosynthesis. A small number of genes and metabolites were upregulated to maintain the basic growth of P. wettsteinii. The physiological and biochemical changes in response to high-temperature stress were revealed, and the important metabolites and key genes involved in the response to high temperature were identified, providing an important reference for the physiological and molecular regulation of high-temperature stress in plants.

Physiological and Molecular Analysis Reveals the Differences of Photosynthesis between Colored and Green Leaf Poplars

Leaf coloration changes evoke different photosynthetic responses among different poplar cultivars. The aim of this study is to investigate the photosynthetic difference between a red leaf cultivar (ZHP) and a green leaf (L2025) cultivar of Populus deltoides. In this study, ‘ZHP’ exhibited wide ranges and huge potential for absorption and utilization of light energy and CO₂ concentration which were similar to those in ‘L2025’ and even showed a stronger absorption for weak light. However, with the increasing light intensity and CO₂ concentration, the photosynthetic capacity in both ‘L2025’ and ‘ZHP’ was gradually restricted, and the net photosynthetic rate (Pn) in ‘ZHP’ was significantly lower than that in ‘L2025’under high light or high CO₂ conditions, which was mainly attributed to stomatal regulation and different photosynthetic efficiency (including the light energy utilization efficiency and photosynthetic CO₂ assimilation efficiency) in these two poplars. Moreover, the higher anthocyanin content in ‘ZHP’ than that in ‘L2025’ was considered to be closely related to the decreased photosynthetic efficiency in ‘ZHP’. According to the results from the JIP-test, the capture efficiency of the reaction center for light energy in ‘L2025’ was significantly higher than that in ‘ZHP’. Interestingly, the higher levels of light quantum caused relatively higher accumulation of Q_A^- in ‘L2025’, which blocked the electron transport and weakened the photosystem II (PSII) performance as compared with ‘ZHP’; however, the decreased capture of light quantum also could not promote the utilization of light energy, which was the key to the low photosynthetic efficiency in ‘ZHP’. The differential expressions of a series of photosynthesis-related genes further promoted these specific photosynthetic processes between ‘L2025’ and ‘ZHP’.

Full-length transcriptome profiling reveals insight into the cold response of two kiwifruit genotypes (A. arguta) with contrasting freezing tolerances

BACKGROUND

Kiwifruit (Actinidia Lindl.) is considered an important fruit species worldwide. Due to its temperate origin, this species is highly vulnerable to freezing injury while under low-temperature stress. To obtain further knowledge of the mechanism underlying freezing tolerance, we carried out a hybrid transcriptome analysis of two A. arguta (Actinidi arguta) genotypes, KL and RB, whose freezing tolerance is high and low, respectively. Both genotypes were subjected to - 25 °C for 0 h, 1 h, and 4 h.

RESULTS

SMRT (single-molecule real-time) RNA-seq data were assembled using the de novo method, producing 24,306 unigenes with an N50 value of 1834 bp. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs showed that they were involved in the 'starch and sucrose metabolism', the 'mitogen-activated protein kinase (MAPK) signaling pathway', the 'phosphatidylinositol signaling system', the 'inositol phosphate metabolism', and the 'plant hormone signal transduction'. In particular, for 'starch and sucrose metabolism', we identified 3 key genes involved in cellulose degradation, trehalose synthesis, and starch degradation processes. Moreover, the activities of beta-GC (beta-glucosidase), TPS (trehalose-6-phosphate synthase), and BAM (beta-amylase), encoded by the abovementioned 3 key genes, were enhanced by cold stress. Three transcription factors (TFs) belonging to the AP2/ERF, bHLH (basic helix-loop-helix), and MYB families were involved in the low-temperature response. Furthermore, weighted gene coexpression network analysis (WGCNA) indicated that beta-GC, TPS5, and BAM3.1 were the key genes involved in the cold response and were highly coexpressed together with the CBF3, MYC2, and MYB44 genes.

CONCLUSIONS

Cold stress led various changes in kiwifruit, the 'phosphatidylinositol signaling system', 'inositol phosphate metabolism', 'MAPK signaling pathway', 'plant hormone signal transduction', and 'starch and sucrose metabolism' processes were significantly affected by low temperature. Moreover, starch and sucrose metabolism may be the key pathway for tolerant kiwifruit to resist low temperature damages. These results increase our understanding of the complex mechanisms involved in the freezing tolerance of kiwifruit under cold stress and reveal a series of candidate genes for use in breeding new cultivars with enhanced freezing tolerance.

TeaAS: a comprehensive database for alternative splicing in tea plants (Camellia sinensis)

Alternative splicing (AS) increases the diversity of transcripts and proteins through the selection of different splice sites and plays an important role in the growth, development and stress tolerance of plants. With the release of the reference genome of the tea plant (Camellia sinensis) and the development of transcriptome sequencing, researchers have reported the existence of AS in tea plants. However, there is a lack of a platform, centered on different RNA-seq datasets, that provides comprehensive information on AS.To facilitate access to information on AS and reveal the molecular function of AS in tea plants, we established the first comprehensive AS database for tea plants (TeaAS, http://www.teaas.cn/index.php ). In this study, 3.96 Tb reads from 66 different RNA-seq datasets were collected to identify AS events. TeaAS supports four methods of retrieval of AS information based on gene ID, gene name, annotation (non-redundant/Kyoto encyclopedia of genes and genomes/gene ontology annotation or chromosomal location) and RNA-seq data. It integrates data pertaining to genome annotation, type of AS event, transcript sequence, and isoforms expression levels from 66 RNA-seq datasets. The AS events resulting from different environmental conditions and that occurring in varied tissue types, and the expression levels of specific transcripts can be clearly identified through this online database. Moreover, it also provides two useful tools, Basic Local Alignment Search Tool and Generic Genome Browser, for sequence alignment and visualization of gene structure.The features of the TeaAS database make it a comprehensive AS bioinformatics platform for researchers, as well as a reference for studying AS events in woody crops. It could also be helpful for revealing the novel biological functions of AS in gene regulation in tea plants.

Metabolic profiling and gene expression analysis provides insights into flavonoid and anthocyanin metabolism in poplar

Poplar, a woody perennial model, is a common and widespread tree genus. We cultivated two red leaf poplar varieties from bud mutation of Populus sp. Linn. '2025' (also known as Zhonglin 2025, L2025 for shot): Populus deltoides varieties with bright red leaves (LHY) and completely red leaves (QHY). After measuring total contents of flavonoid, anthocyanin, chlorophyll and carotenoid metabolites, a liquid chromatography-electrospray ionization-tandem mass spectrometry system was used for the relative quantification of widely targeted metabolites in leaves of three poplar varieties. A total of 210 flavonoid metabolites (89 flavones, 40 flavonols, 25 flavanones, 18 anthocyanins, 16 isoflavones, 7 dihydroflavonols, 7 chalcones, 5 proanthocyanidins and 3 other flavonoid metabolites) were identified. Compared with L2025, 48 and 8 flavonoids were more and less abundant, respectively, in LHY, whereas 51 and 9 flavonoids were more and less abundant in QHY, respectively. On the basis of a comprehensive analysis of the metabolic network, gene expression levels were analyzed by deep sequencing to screen for potential reference genes for the red leaves. Most phenylpropanoid biosynthesis pathway-involved genes were differentially expressed among the examined varieties. Gene expression analysis also revealed several potential anthocyanin biosynthesis regulators including three MYB genes. The study results provide new insights into poplar flavonoid metabolites and represent the theoretical basis for future studies on leaf coloration in this model tree species.

Differences in Anthocyanin Accumulation Profiles between Teinturier and Non-Teinturier Cultivars during Ripening

Anthocyanins are vital components of plant secondary metabolites, and are also the most important coloring substances in wine. Teinturier cultivars are rich in anthocyanins. However, the differences in anthocyanin accumulation and profiles between teinturier and non-teinturier cultivars have not been reported. In this study, Yan 73 and Dunkelfelder were selected as the experimental materials, and three non-teinturier cultivars were used for comparison. LC-MS and qRT-PCR were used to determine the individual anthocyanin contents and the relative gene expression. The results show that the total anthocyanin content of the teinturier cultivars was considerably higher than that in non-teinturier cultivars, and the levels of individual anthocyanins increased gradually during ripening. Lower ratios of modified anthocyanins were found in the teinturier cultivars, which was not only due to the high expression level of VvUFGT and VvGST4, but also due to the relatively low expression of VvOMT in these cultivars. Cluster analysis of gene expression and anthocyanin accumulation showed that VvUFGT is related to anthocyanin accumulation, and that AM1 is related to the synthesis and transport of methylated anthocyanins. Our results will be useful for further clarifying the pathways of anthocyanin synthesis, modification, and transport in teinturier cultivars.

Chromatin accessibility and translational landscapes of tea plants under chilling stress

Plants have evolved regulatory mechanisms at multiple levels to regulate gene expression in order to improve their cold adaptability. However, limited information is available regarding the stress response at the chromatin and translational levels. Here, we characterize the chromatin accessibility, transcriptional, and translational landscapes of tea plants in vivo under chilling stress for the first time. Chilling stress significantly affected both the transcription and translation levels as well as the translation efficiency of tea plants. A total of 3010 genes that underwent rapid and independent translation under chilling stress were observed, and they were significantly enriched in the photosynthesis-antenna protein and phenylpropanoid biosynthesis pathways. A set of genes that were significantly responsive to cold at the transcription and translation levels, including four (+)-neomenthol dehydrogenases (MNDs) and two (E)-nerolidol synthases (NESs) arranged in tandem on the chromosomes, were also found. We detected potential upstream open reading frames (uORFs) on 3082 genes and found that tea plants may inhibit the overall expression of genes by enhancing the translation of uORFs under chilling stress. In addition, we identified distal transposase hypersensitive sites (THSs) and proximal THSs and constructed a transcriptional regulatory network for tea plants under chilling stress. We also identified 13 high-confidence transcription factors (TFs) that may play a crucial role in cold regulation. These results provide valuable information regarding the potential transcriptional regulatory network in plants and help to clarify how plants exhibit flexible responses to chilling stress.

Up-regulation of phenylpropanoid biosynthesis system in peach species by peach aphids produces anthocyanins that protect the aphids against UVB and UVC radiation

Conspicuous color is a common trait of foliar galls, but their relationship with gall-inducing insects is unknown. Red and green galls were taken from sunny or shady parts of peach species Prunus persica (L.) Batsch. f. rubro-plena Schneid with peach aphid Tuberocephalus momonis (Matsumura) infestation. We found that the loss of photosynthetic pigments was associated with the conspicuous coloration of green gall tissues. The concentrations of anthocyanins significantly increased following ultraviolet (UV) irradiation of green gall tissues, suggesting that accumulation of anthocyanins in red galls is related to ultraviolet B and C (UVB and UVC) radiation. The expression of structural genes related to the biosynthesis of chlorogenic acid and malic acid benzoate was increased in all gall tissues and negatively correlated with the expression profiles of certain genes associated with photosynthetic biosynthesis, indicating that the increased transcript levels of the phenylpropanoid pathway might cause loss of photosynthetic efficiency in the gall tissues. Transcriptome and quantitative reverse transcription PCR analyses revealed that MYB transcription factors that up-regulate the biosynthesis of anthocyanins in red gall tissues might be activated by both UVB and UVC exposure. Comet assays suggest that green and red gall tissues have similar DNA damage following UV irradiation. No obvious effect of the up-regulated compounds on the growth of the peach aphid was observed. Interestingly, peach aphids under leaves painted with anthocyanins had lower mortality following UV irradiation than those in controls. These results suggest that the anthocyanins in red gall tissues have a defensive function for the peach aphid, protecting it against UV radiation.

The Role of Polyphenols in Abiotic Stress Response: The Influence of Molecular Structure

Abiotic stressors such as extreme temperatures, drought, flood, light, salt, and heavy metals alter biological diversity and crop production worldwide. Therefore, it is important to know the mechanisms by which plants cope with stress conditions. Polyphenols, which are the largest group of plant-specialized metabolites, are generally recognized as molecules involved in stress protection in plants. This diverse group of metabolites contains various structures, from simple forms consisting of one aromatic ring to more complex ones consisting of large number of polymerized molecules. Consequently, all these molecules, depending on their structure, may show different roles in plant growth, development, and stress protection. In the present review, we aimed to summarize data on how different polyphenol structures influence their biological activity and their roles in abiotic stress responses. We focused our review on phenolic acids, flavonoids, stilbenoids, and lignans.

Genome-Wide Investigation of N6-Methyladenosine Regulatory Genes and Their Roles in Tea (Camellia sinensis) Leaves During Withering Process

N6-methyladenosine (m6A), one of the internal modifications of RNA molecules, can directly influence RNA abundance and function without altering the nucleotide sequence, and plays a pivotal role in response to diverse environmental stresses. The precise m6A regulatory mechanism comprises three types of components, namely, m6A writers, erasers, and readers. To date, the research focusing on m6A regulatory genes in plant kingdom is still in its infancy. Here, a total of 34 m6A regulatory genes were identified from the chromosome-scale genome of tea plants. The expansion of m6A regulatory genes was driven mainly by whole-genome duplication (WGD) and segmental duplication, and the duplicated gene pairs evolved through purifying selection. Gene structure analysis revealed that the sequence variation contributed to the functional diversification of m6A regulatory genes. Expression pattern analysis showed that most m6A regulatory genes were differentially expressed under environmental stresses and tea-withering stage. These observations indicated that m6A regulatory genes play essential roles in response to environmental stresses and tea-withering stage. We also found that RNA methylation and DNA methylation formed a negative feedback by interacting with each other's methylation regulatory genes. This study provided a foundation for understanding the m6A-mediated regulatory mechanism in tea plants under environmental stresses and tea-withering stage.

Genome-Wide Characterization of the HSP20 Gene Family Identifies Potential Members Involved in Temperature Stress Response in Apple

Apple (Malus domestica Borkh.), an economically important tree fruit worldwide, frequently suffers from temperature stress during growth and development, which strongly affects the yield and quality. Heat shock protein 20 (HSP20) genes play crucial roles in protecting plants against abiotic stresses. However, they have not been systematically investigated in apple. In this study, we identified 41 HSP20 genes in the apple 'Golden Delicious' genome. These genes were unequally distributed on 15 different chromosomes and were classified into 10 subfamilies based on phylogenetic analysis and predicted subcellular localization. Chromosome mapping and synteny analysis indicated that three pairs of apple HSP20 genes were tandemly duplicated. Sequence analysis revealed that all apple HSP20 proteins reflected high structure conservation and most apple HSP20 genes (92.6%) possessed no introns, or only one intron. Numerous apple HSP20 gene promoter sequences contained stress and hormone response cis-elements. Transcriptome analysis revealed that 35 of 41 apple HSP20 genes were nearly unchanged or downregulated under normal temperature and cold stress, whereas these genes exhibited high-expression levels under heat stress. Subsequent qRT-PCR results showed that 12 of 29 selected apple HSP20 genes were extremely up-regulated (more than 1,000-fold) after 4 h of heat stress. However, the heat-upregulated genes were barely expressed or downregulated in response to cold stress, which indicated their potential function in mediating the response of apple to heat stress. Taken together, these findings lay the foundation to functionally characterize HSP20 genes to unravel their exact role in heat defense response in apple.

Genome-Wide Characterization of the HSP20 Gene Family Identifies Potential Members Involved in Temperature Stress Response in Apple

Apple (Malus domestica Borkh.), an economically important tree fruit worldwide, frequently suffers from temperature stress during growth and development, which strongly affects the yield and quality. Heat shock protein 20 (HSP20) genes play crucial roles in protecting plants against abiotic stresses. However, they have not been systematically investigated in apple. In this study, we identified 41 HSP20 genes in the apple ‘Golden Delicious’ genome. These genes were unequally distributed on 15 different chromosomes and were classified into 10 subfamilies based on phylogenetic analysis and predicted subcellular localization. Chromosome mapping and synteny analysis indicated that three pairs of apple HSP20 genes were tandemly duplicated. Sequence analysis revealed that all apple HSP20 proteins reflected high structure conservation and most apple HSP20 genes (92.6%) possessed no introns, or only one intron. Numerous apple HSP20 gene promoter sequences contained stress and hormone response cis-elements. Transcriptome analysis revealed that 35 of 41 apple HSP20 genes were nearly unchanged or downregulated under normal temperature and cold stress, whereas these genes exhibited high-expression levels under heat stress. Subsequent qRT-PCR results showed that 12 of 29 selected apple HSP20 genes were extremely up-regulated (more than 1,000-fold) after 4 h of heat stress. However, the heat-upregulated genes were barely expressed or downregulated in response to cold stress, which indicated their potential function in mediating the response of apple to heat stress. Taken together, these findings lay the foundation to functionally characterize HSP20 genes to unravel their exact role in heat defense response in apple.

Transcriptomic analysis of tea plant (Camellia sinensis) revealed the co-expression network of 4111 paralogous genes and biosynthesis of quality-related key metabolites under multiple stresses

The tea plant is an essential economic plant in many countries. However, its growing season renders them vulnerable to stresses. To understand the transcriptomic influences of these stresses on tea plants, we sequenced and analyzed the transcriptomes under drought, high-temperature, and pest. Paralogs were identified by comparing 14 evolutionarily close genomes. The differentially expressed paralog (DEPs) genes were analyzed regarding single or multiple stresses, and 1075 of the 4111 DEPs were commonly found in all the stresses. The co-expression network of the DEPs and TFs indicated that genes of catechin biosynthesis were associated with most transcription factors specific to each stress. The genes playing a significant role in the late response to drought and pest stress mainly functioned in the early response to high-temperature. This study revealed the relationship between stress and regulation of QRM synthesis and the role of QRMs in response to these (a)biotic stresses.

Carrot Anthocyanins Genetics and Genomics: Status and Perspectives to Improve Its Application for the Food Colorant Industry

Purple or black carrots (Daucus carota ssp. sativus var. atrorubens Alef) are characterized by their dark purple- to black-colored roots, owing their appearance to high anthocyanin concentrations. In recent years, there has been increasing interest in the use of black carrot anthocyanins as natural food dyes. Black carrot roots contain large quantities of mono-acylated anthocyanins, which impart a measure of heat-, light- and pH-stability, enhancing the color-stability of food products over their shelf-life. The genetic pathway controlling anthocyanin biosynthesis appears well conserved among land plants; however, different variants of anthocyanin-related genes between cultivars results in tissue-specific accumulations of purple pigments. Thus, broad genetic variations of anthocyanin profile, and tissue-specific distributions in carrot tissues and organs, can be observed, and the ratio of acylated to non-acylated anthocyanins varies significantly in the purple carrot germplasm. Additionally, anthocyanins synthesis can also be influenced by a wide range of external factors, such as abiotic stressors and/or chemical elicitors, directly affecting the anthocyanin yield and stability potential in food and beverage applications. In this study, we critically review and discuss the current knowledge on anthocyanin diversity, genetics and the molecular mechanisms controlling anthocyanin accumulation in carrots. We also provide a view of the current knowledge gaps and advancement needs as regards developing and applying innovative molecular tools to improve the yield, product performance and stability of carrot anthocyanin for use as a natural food colorant.

Comparative transcriptomic analysis reveals gene expression associated with cold adaptation in the tea plant Camellia sinensis

Background

Low temperature restricts the planting range of all crops, but cold acclimation induces adaption to cold stress in many plants. Camellia sinensis, a perennial evergreen tree that is the source of tea, is mainly grown in warm areas. Camellia sinensis var. sinensis (CSS) has greater cold tolerance than Camellia sinensis var. assamica (CSA). To gain deep insight into the molecular mechanisms underlying cold adaptation, we investigated the physiological responses and transcriptome profiles by RNA-Seq in two tea varieties, cold resistant SCZ (classified as CSS) and cold susceptible YH9 (classified as CSA), during cold acclimation.

Results

Under freezing stress, lower relative electrical conductivity and higher chlorophyll fluorescence (Fv/Fm) values were detected in SCZ than in YH9 when subjected to freezing acclimation. During cold treatment, 6072 and 7749 DEGs were observed for SCZ and YH9, respectively. A total of 978 DEGs were common for both SCZ and YH9 during the entire cold acclimation process. DEGs were enriched in pathways of photosynthesis, hormone signal transduction, and transcriptional regulation of plant-pathogen interactions. Further analyses indicated that decreased expression of Lhca2 and higher expression of SnRK2.8 are correlated with cold tolerance in SCZ.

Conclusions

Compared with CSA, CSS was significantly more resistant to freezing after cold acclimation, and this increased resistance was associated with an earlier expression of cold-induced genes. Because the greater transcriptional differentiation during cold acclimation in SCZ may contribute to its greater cold tolerance, our studies identify specific genes involved in photoinhibition, ABA signal conduction, and plant immunity that should be studied for understanding the processes involved in cold tolerance. Marker-assisted breeding focused on the allelic variation at these loci provides an avenue for the possible generation of CSA cultivars that have CSS-level cold tolerance.

Electronic supplementary material

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