A transcriptomic analysis of Chrysanthemum nankingense provides insights into the basis of low temperature tolerance

Analysis of the Mechanism of Wood Vinegar and Butyrolactone Promoting Rapeseed Growth and Improving Low-Temperature Stress Resistance Based on Transcriptome and Metabolomics

Rapeseed is an important oil crop in the world. Wood vinegar could increase the yield and abiotic resistance of rapeseed. However, little is known about the underlying mechanisms of wood vinegar or its valid chemical components on rapeseed. In the present study, wood vinegar and butyrolactone (γ-Butyrolactone, one of the main components of wood vinegar) were applied to rapeseed at the seedling stage, and the molecular mechanisms of wood vinegar that affect rapeseed were studied by combining transcriptome and metabolomic analyses. The results show that applying wood vinegar and butyrolactone increases the biomass of rapeseed by increasing the leaf area and the number of pods per plant, and enhances the tolerance of rapeseed under low temperature by reducing membrane lipid oxidation and improving the content of chlorophyll, proline, soluble sugar, and antioxidant enzymes. Compared to the control, 681 and 700 differentially expressed genes were in the transcriptional group treated with wood vinegar and butyrolactone, respectively, and 76 and 90 differentially expressed metabolites were in the metabolic group. The combination of transcriptome and metabolomic analyses revealed the key gene-metabolic networks related to various pathways. Our research shows that after wood vinegar and butyrolactone treatment, the amino acid biosynthesis pathway of rapeseed may be involved in mediating the increase in rapeseed biomass, the proline metabolism pathway of wood vinegar treatment may be involved in mediating rapeseed's resistance to low-temperature stress, and the sphingolipid metabolism pathway of butyrolactone treatment may be involved in mediating rapeseed's resistance to low-temperature stress. It is suggested that the use of wood vinegar or butyrolactone are new approaches to increasing rapeseed yield and low-temperature resistance.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family

BACKGROUND

Terpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated.

RESULTS

We identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality.

CONCLUSIONS

We systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Comparative transcriptome analysis provides molecular insights into heterosis of waterlogging tolerance in Chrysanthemum indicum

BACKGROUND

Heterosis breeding is one of the most important breeding methods for chrysanthemum. To date, the genetic mechanisms of heterosis for waterlogging tolerance in chrysanthemum are still unclear. This study aims to analyze the expression profiles and potential heterosis-related genes of two hybrid lines and their parents with extreme differences in waterlogging tolerance under control and waterlogging stress conditions by RNA-seq.

RESULTS

A population of 140 F1 progeny derived from Chrysanthemum indicum (Nanchang) (waterlogging-tolerant) and Chrysanthemum indicum (Nanjing) (waterlogging-sensitive) was used to characterize the extent of genetic variation in terms of seven waterlogging tolerance-related traits across two years. Lines 98 and 95, respectively displaying positive and negative overdominance heterosis for the waterlogging tolerance traits together with their parents under control and waterlogging stress conditions, were used for RNA-seq. In consequence, the maximal number of differentially expressed genes (DEGs) occurred in line 98. Gene ontology (GO) enrichment analysis revealed multiple stress-related biological processes for the common up-regulated genes. Line 98 had a significant increase in non-additive genes under waterlogging stress, with transgressive up-regulation and paternal-expression dominant patterns being the major gene expression profiles. Further, GO analysis identified 55 and 95 transgressive up-regulation genes that overlapped with the up-regulated genes shared by two parents in terms of responses to stress and stimulus, respectively. 6,640 genes in total displaying maternal-expression dominance patterns were observed in line 95. In addition, 16 key candidate genes, including SAP12, DOX1, and ERF017 which might be of significant importance for the formation of waterlogging tolerance heterosis in line 98, were highlighted.

CONCLUSION

The current study provides a comprehensive overview of the root transcriptomes among F1 hybrids and their parents under waterlogging stress. These findings lay the foundation for further studies on molecular mechanisms underlying chrysanthemum heterosis on waterlogging tolerance.

Rubisco and sucrose synthesis and translocation are involved in the regulation of photosynthesis in wheat with different source-sink relationships

Source-sink relationships influence photosynthesis. So far, the limiting factors for photosynthesis of wheat cultivars with different source-sink relationships have not been determined. We aimed to determine the variation patterns of photosynthetic characteristics of wheat cultivars with different source-sink relationships. In this study, two wheat cultivars with different source-sink relationships were selected for photosynthetic physiological analyses. The results showed that YM25 (source-limited cultivar) had higher photosynthetic efficiency compared to YM1 (sink-limited cultivar). This is mainly due to a stronger photochemical efficiency, electron transfer capacity, and Rubisco carboxylation capacity of YM25. YM25 accumulated less soluble carbohydrates in flag leaves than YM1. This is mainly due to the stronger sucrose synthesis and transport capacity of YM25 by presenting higher sucrose-related enzyme activities and gene expression. A PCA analysis showed that Rubisco was the main factor limiting the photosynthetic capacity of YM25. The soluble sugar accumulation in flag leaves and sink limitation decreased the photosynthetic activity of YM1. Increased N application improved source-sink relationships and increased grain yield and source leaf photosynthetic capacity in both two wheat cultivars. Taken together, our findings suggest that Rubisco and sucrose synthesis and translocation are involved in the regulation of photosynthesis of wheat cultivars with different source-sink relationships and that source and sink limitation effects should be considered in photosynthesis.

Research advance on cold tolerance in chrysanthemum

Chrysanthemums are one of the top ten most well-known traditional famous flowers in China and one of the top four cut flowers worldwide, holding a significant position in landscape gardening. The cold temperatures of winter restrict the cultivation, introduction, and application of chrysanthemum, resulting in high costs for year-round production. This severely impacts the ornamental and economic value of chrysanthemum. Therefore, research on cold tolerance is of vital importance for guiding chrysanthemum production and application. With the development of genomics, transcriptomics, metabolomics, and other omics approaches, along with high-throughput molecular marker technologies, research on chrysanthemum cold tolerance has been continuously advancing. This article provides a comprehensive overview of the progress in cold tolerance research from various aspects, including chrysanthemum phenotype, physiological mechanisms, the forward genetics, molecular mechanisms, and breeding. The aim is to offer insights into the mechanisms of cold tolerance in chrysanthemum and provide reference for in-depth research and the development of new cold tolerance chrysanthemum varieties.

Detection of Transcription Factors Related to Axillary Bud Development after Exposure to Cold Conditions in Hexaploid Chrysanthemum morifolium Using Arabidopsis Information

Chrysanthemum is one of the most commercially used ornamental flowering plants in the world. As chrysanthemum is self-incompatible, the propagation of identical varieties is carried out through cuttings rather than through seed. Axillary bud development can be controlled by changing the temperature; for instance, axillary bud development in some varieties is suppressed at high temperatures. In this study, we focused on the simultaneous axillary bud growth from multiple lines of chrysanthemum upon changing conditions from low to normal temperature. Transcriptome analysis was conducted on the Chrysanthemum morifolium cultivar 'Jinba' to identify the important genes for axillary bud development seen when moved from low-temperature treatment to normal cultivation temperature. We performed RNA-Seq analysis on plants after cold conditions in two-day time-course experiments. Under these settings, we constructed a transcriptome of 415,923 C. morifolium and extracted 7357 differentially expressed genes. Our understanding of Arabidopsis axillary meristem development and growth showed that at least 101 genes in our dataset were homologous to transcription factors involved in the biological process. In addition, six genes exhibited statistically significant variations in expression throughout conditions. We hypothesized that these genes were involved in the formation of axillary buds in C. morifolium after cold conditions.

Comprehensive Transcriptome Analysis of Responses during Cold Stress in Wheat (Triticum aestivum L.)

Wheat production is often impacted by pre-winter freezing damage and cold spells in later spring. To study the influences of cold stress on wheat seedlings, unstressed Jing 841 was sampled once at the seedling stage, followed by 4 °C stress treatment for 30 days and once every 10 days. A total of 12,926 differentially expressed genes (DEGs) were identified from the transcriptome. K-means cluster analysis found a group of genes related to the glutamate metabolism pathway, and many genes belonging to the bHLH, MYB, NAC, WRKY, and ERF transcription factor families were highly expressed. Starch and sucrose metabolism, glutathione metabolism, and plant hormone signal transduction pathways were found. Weighted Gene Co-Expression Network Analysis (WGCNA) identified several key genes involved in the development of seedlings under cold stress. The cluster tree diagram showed seven different modules marked with different colors. The blue module had the highest correlation coefficient for the samples treated with cold stress for 30 days, and most genes in this module were rich in glutathione metabolism (ko00480). A total of eight DEGs were validated using quantitative real-time PCR. Overall, this study provides new insights into the physiological metabolic pathways and gene changes in a cold stress transcriptome, and it has a potential significance for improving freezing tolerance in wheat.

Genome-wide identification and expression analysis of SBP-box gene family reveal their involvement in hormone response and abiotic stresses in Chrysanthemum nankingense

SQUAMOSA promoter-binding-protein (SBP)-box family proteins are a class of plant-specific transcription factors, and widely regulate the development of floral and leaf morphology in plant growth and involve in environment and hormone signal response. In this study, we isolated and identified 21 non-redundant SBP-box genes in Chrysanthemum nankingense with bioinformatics analysis. Sequence alignments of 21 CnSBP proteins discovered a highly conserved SBP domain including two zinc finger-like structures and a nuclear localization signal region. According to the amino acid sequence alignments, 67 SBP-box genes from Arabidopsis thaliana, rice, Artemisia annua and C. nankingense were clustered into eight groups, and the motif and gene structure analysis also sustained this classification. The gene evolution analysis indicated the CnSBP genes experienced a duplication event about 10 million years ago (Mya), and the CnSBP and AtSPL genes occurred a divergence at 24 Mya. Transcriptome data provided valuable information for tissue-specific expression profiles of the CnSBPs, which highly expressed in floral tissues and differentially expressed in leaf, root and stem organs. Quantitative Real-time Polymerase Chain Reaction data showed expression patterns of the CnSBPs under exogenous hormone and abiotic stress treatments, separately abscisic acid, salicylic acid, gibberellin A3, methyl jasmonate and ethylene spraying as well as salt and drought stresses, indicating that the candidate CnSBP genes showed differentiated spatiotemporal expression patterns in response to hormone and abiotic stresses. Our study provides a systematic genome-wide analysis of the SBP-box gene family in C. nankingense. In general, it provides a fundamental theoretical basis that SBP-box genes may regulate the resistance of stress physiology in chrysanthemum via exogenous hormone pathways.

Genome-Wide Identification and Expression Analysis of CsCaM/CML Gene Family in Response to Low-Temperature and Salt Stresses in Chrysanthemum seticuspe

Calmodulin (CaM) and calmodulin-like proteins (CML) act as significant Ca²⁺ sensors binding Ca²⁺ with EF-hand motifs and have been reported to be involved in various environmental stresses in plants. In this study, calmodulin CsCaM/CML gene family members were identified based on the genome of Chrysanthemum seticuspe published recently; a phylogenetic tree was constructed; gene structures and chromosomal locations of CsCaM/CML were depicted; cis-acting regulatory elements were predicted; collinearity and duplicate events of CaM/CML were analyzed using MCScanX software; and the expression levels of CsCaM/CML in response to abiotic stress were analyzed, based on the published RNA-seq data. We identified 86 CsCaM/CML (4 CsCaMs and 82 CsCMLs) genes in total. Promoter sequences of CsCaM/CML contained elements related to abiotic stresses (including low-temperature and anaerobic stresses) and plant hormones (including abscisic acid (ABA), MeJA, and salicylic acid). CsCaM/CML genes were distributed on nine chromosomes unevenly. Collinearity analysis indicated that recent segmental duplications significantly enlarged the scale of the CML family in C. seticuspe. Four CsCMLs (CsCML14, CsCML50, CsCML65, and CsCML79) were statistically differentially regulated under low-temperature and salt stress compared with those in the normal condition. These results indicate diverse roles of CsCaM/CML in plant development and in response to environmental stimuli in C. seticuspe.

Expression Analysis of TCP Transcription Factor Family in Autopolyploids of Chrysanthemum nankingense

Autopolyploids often exhibit plant characteristics different from their diploid ancestors and are frequently associated with altered genes expression controlling growth and development. TCP is a unique transcription factor family in plants that is closely related to plant growth and development. Based on transcriptome sequencing of Chrysanthemum nankingense, 23 full-length TCP genes were cloned. The expression of CnTCP9 was most variable in tetraploids, at least threefold greater than diploids. Due to the lack of a C. nankingense transgenic system, we overexpressed CnTCP9 in Arabidopsis thaliana (Col-0) and Chrysanthemum morifolium. Overexpression of CnTCP9 caused enlargement of leaves in A. thaliana and petals in C. morifolium, and the expression of genes downstream of the GA pathway in C. morifolium were increased. Our results suggest that autopolyploidization of C. nankingense led to differential expression of TCP family genes, thereby affecting plant characteristics by the GA pathway. This study improves the understanding of enlarged plant size after autopolyploidization.

Biogeographic and metabolic studies support a glacial radiation hypothesis during Chrysanthemum evolution

Chrysanthemum (Chrysanthemum morifolium Ramat.) is an economically important plant species growing worldwide. However, its origin, especially as revealed by biogeographic and metabolomics research, remains unclear. To understand the geographic distribution of species diversity and metabolomics in three genera (Chrysanthemum, Ajania, and Phaeostigma), geographic information systems and gas chromatography-mass spectrometry were used in 19, 15, and 4 species respectively. China and Japan were two potential panbiogeographic nodes and diverse hotspots of Chrysanthemum, with species richness ratios of 58.97 and 33.33%. We studied different species from two hotspots which in similar geographical environments had closer chemotaxonomic relationships under the same cultivation conditions based on a cluster of 30 secondary metabolites. The average distribution altitude (ADA) differed significantly among Chrysanthemum, Ajania, and Phaeostigma in which it was 1227.49, 2400.12, and 3760.53 m.a.s.l. respectively, and the presence/absence of ray florets (RF) was significantly correlated with ADA (-0.62). Mountain landform was an important contributor to global Chrysanthemum diversity, playing a key role in the divergence and distribution pattern of Chrysanthemum and its allies. The Hengduan Mountains-Qinling Mountains (HDQ) in China was a potential secondary radiation and evolution center of Chrysanthemum and its related genera in the world. During the Quaternary glacial-interglacial cycles, this region became their refuge, and they radiated and evolved from this center.

Characterization of the TCP Gene Family in Chrysanthemum nankingense and the Role of CnTCP4 in Cold Tolerance

Plant-specific TCP transcription factors play a key role in plant development and stress responses. Chrysanthemum nankingense shows higher cold tolerance than its ornamental polyploid counterpart. However, whether the TCP gene family plays a role in conferring cold tolerance upon C. nankingense remains unknown. Here, we identified 23 CnTCP genes in C. nankingense, systematically analyzed their phylogenetic relationships and synteny with TCPs from other species, and evaluated their expression profiles at low temperature. Phylogenetic analysis of the protein sequences suggested that CnTCP proteins fall into two classes and three clades, with a typical bHLH domain. However, differences between C. nankingense and Arabidopsis in predicted protein structure and binding sites suggested a unique function of CnTCPs in C. nankingense. Furthermore, expression profiles showed that expression of most CnTCPs were downregulated under cold conditions, suggesting their importance in plant responses to cold stress. Notably, expression of miR319 and of its predicted target genes, CnTCP2/4/14, led to fast responses to cold. Overexpression of Arabidopsis CnTCP4 led to hypersensitivity to cold, suggesting that CnTCP4 might play a negative role in C. nankingense responses to cold stress. Our results provide a foundation for future functional genomic studies on this gene family in chrysanthemum.

Transcriptomic analysis of Vigna radiata in response to chilling stress and uniconazole application

Background

Chilling injury of mung bean (Vigna radiata (L.)) during the blooming and podding stages is a major agricultural threat in Northeast China. Uniconazole (UNZ) can alleviate water deficit stress in soybean and waterlogging stress in mung bean. However, there has been no report on the effect of UNZ application on the growth and transcriptomic profile of mung bean under chilling stress.

Results

UNZ application before chilling stress at the R1 stage alleviated the decline in mung bean yield. UNZ delayed the decrease in leaf chlorophyll content under chilling stress at the R1 stage and accelerated the increase in leaf chlorophyll content during the recovery period. Eighteen separate RNA-Seq libraries were generated from RNA samples collected from leaves exposed to six different treatment schemes. The numbers of DEGs specific for UNZ treatment between D1 + S vs. D1 and D4 + S vs. D4 were 708 and 810, respectively. GO annotations showed that photosynthesis genes were obviously enriched among the genes affected by chilling stress and UNZ application. KEGG pathway enrichment analysis indicated that 4 pathways (cutin, suberin and wax biosynthesis; photosynthesis; porphyrin and chlorophyll metabolism; and ribosome) were downregulated, while plant–pathogen interaction was upregulated, by chilling stress. UNZ application effectively prevented the further downregulation of the gene expression of members of these 4 KEGG pathways under chilling stress.

Conclusions

UNZ application effectively delayed the decrease in photosynthetic pigment content under chilling stress and accelerated the increase in photosynthetic pigment content during the recovery period, thus effectively limiting the decline in mung bean yield. UNZ application effectively prevented the further downregulation of the gene expression of members of 4 KEGG pathways under chilling stress and increased mung bean tolerance to chilling stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08443-6.

Comparison of chrysanthemum flowers grown under hydroponic and soil-based systems: yield and transcriptome analysis

BACKGROUND

Flowers of Chrysanthemum × morifolium Ramat. are used as tea in traditional Chinese cuisine. However, with increasing population and urbanization, water and land availability have become limiting for chrysanthemum tea production. Hydroponic culture enables effective, rapid nutrient exchange, while requiring no soil and less water than soil cultivation. Hydroponic culture can reduce pesticide residues in food and improve the quantity or size of fruits, flowers, and leaves, and the levels of active compounds important for nutrition and health. To date, studies to improve the yield and active compounds of chrysanthemum have focused on soil culture. Moreover, the molecular effects of hydroponic and soil culture on chrysanthemum tea development remain understudied.

RESULTS

Here, we studied the effects of soil and hydroponic culture on yield and total flavonoid and chlorogenic acid contents in chrysanthemum flowers (C. morifolium 'wuyuanhuang'). Yield and the total flavonoids and chlorogenic acid contents of chrysanthemum flowers were higher in the hydroponic culture system than in the soil system. Transcriptome profiling using RNA-seq revealed 3858 differentially expressed genes (DEGs) between chrysanthemum flowers grown in soil and hydroponic conditions. Gene Ontology (GO) enrichment annotation revealed that these differentially transcribed genes are mainly involved in "cytoplasmic part", "biosynthetic process", "organic substance biosynthetic process", "cell wall organization or biogenesis" and other processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed enrichment in "metabolic pathways", "biosynthesis of secondary metabolites", "ribosome", "carbon metabolism", "plant hormone signal transduction" and other metabolic processes. In functional annotations, pathways related to yield and formation of the main active compounds included phytohormone signaling, secondary metabolism, and cell wall metabolism. Enrichment analysis of transcription factors also showed that under the hydroponic system, bHLH, MYB, NAC, and ERF protein families were involved in metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction.

CONCLUSIONS

Hydroponic culture is a simple and effective way to cultivate chrysanthemum for tea production. A transcriptome analysis of chrysanthemum flowers grown in soil and hydroponic conditions. The large number of DEGs identified confirmed the difference of the regulatory machinery under two culture system.

Transcriptome and Metabolome Analyses of the Flowers and Leaves of Chrysanthemum dichrum

Chrysanthemum dichrum is an important wild species in the family Asteraceae. However, because of a lack of genetic information, there has been relatively little research conducted on the molecular mechanisms in C. dichrum. There is no report describing the transcriptome and metabolome of C. dichrum flowers and leaves at different developmental stages. In this study, high-throughput sequencing and RNA-seq analyses were used to investigate the transcriptome of C. dichrum leaves, flower buds, and blooming flowers. Additionally, these three tissues also underwent a metabolomics analysis. A total of 447,313,764 clean reads were assembled into 77,683 unigenes, with an average length of 839 bp. Of the 44,204 annotated unigenes, 42,189, 28,531, 23,420, and 17,599 were annotated using the Nr, Swiss-Prot, KOG, and KEGG databases, respectively. Furthermore, 31,848 differentially expressed genes (DEGs) were detected between the leaves and flower buds, whereas 23,197 DEGs were detected between the leaves and blooming flowers, and 11,240 DEGs were detected between the flower buds and blooming flowers. Finally, a quantitative real-time Polymerase Chain Reaction (qRT-PCR) assay was conducted to validate the identified DEGs. The metabolome data revealed several abundant metabolites in C. dichrum leaves, flower buds, and blooming flowers, including raffinose, 1-kestose, asparagine, glutamine, and other medicinal compounds. The expression patterns of significant DEGs revealed by the transcriptome analysis as well as the data for the differentially abundant metabolites in three C. dichrum tissues provide important genetic and metabolic information relevant for future investigations of the molecular mechanisms in C. dichrum. Moreover, the results of this study may be useful for the molecular breeding, development, and application of C. dichrum resources.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Physiological and Transcripts Analyses Reveal the Mechanism by Which Melatonin Alleviates Heat Stress in Chrysanthemum Seedlings

Heat stress limits the growth and development of chrysanthemum seedlings. Although melatonin (MT) has been linked to the heat stress response in plants, research on the underlying molecular mechanisms is scarce. In this study, the regulatory networks of MT on heat stress in chrysanthemum seedlings were explored. Physiological measurements suggested that MT not only reduced malondialdehyde accumulation, hydrogen peroxide content, and superoxide anion free radical generation rate, but also significantly promoted osmotic regulation substance synthesis (proline and soluble protein), antioxidant accumulation (GSH and AsA), and the antioxidant enzyme activities (SOD, POD, CAT, and APX) in chrysanthemum leaves under heat stress. Furthermore, MT increased the fresh weight, dry weight, chlorophyll content, photosynthesis rate, and gas exchange indexes. Further, RNA-seq results revealed 33,497 and 36,740 differentially expressed genes in the S/Con and SMT/ConMT comparisons, respectively. The differences in the comparisons revealed that MT regulated heat shock transcription factors (HSFs) and heat shock proteins (HSPs), and the genes involved in Ca²⁺ signal transduction (CNGCs and CAM/CMLs), starch and sucrose metabolism (EDGL, BGLU, SuS, and SPS), hormone (PP2Cs, AUX/IAAs, EBFs, and MYC2), chlorophyll metabolism (HEMA and PORA), flavonoid biosynthesis (CHS, DFR, and FNS), and carotenoid biosynthesis (DXPS, GGDP, and PSY). MT effectively improved chrysanthemum seedling heat-resistance. Our study, thus, provides novel evidence of a gene network regulated by MT under heat stress.

Time-series transcriptomic analysis reveals novel gene modules that control theanine biosynthesis in tea plant (Camellia sinensis)

Theanine (thea) is a unique non-protein amino acid in tea plant (Camellia sinensis) and one of the most important small molecular compounds for tea quality and health effects. The molecular mechanism that maintains thea biosynthesis is not clear but may be reflected in complicated biological networks as other secondary metabolites in plants. We performed an integrative transcriptomic analysis of tea seedlings bud and leave over the time-course of ethylamine (EA) treatment that activated thea pathway. We identified 54 consistent differentially expressed genes (cDEGs, 25 upregulated and 29 downregulated) during thea activation. Gene Ontology (GO) functional enrichment analysis of upregulated genes and downregulated genes showed that they may function as a cascade of biological events during their cooperative contribution to thea biosynthesis. Among the total cDEGs, a diversity of functional genes (e.g., enzymes, transcription factors, transport and binding proteins) were identified, indicating a hierarchy of gene control network underlying thea biosynthesis. A gene network associated with thea biosynthesis was modeled and three interconnected gene functional modules were identified. Among the gene modules, several topologically important genes (e.g., CsBCS-1, CsRP, CsABC2) were experimentally validated using a combined thea content and gene expression analysis. Collectively, we presented here for the first time a comprehensive landscape of the biosynthetic mechanism of thea controlled by a underling gene network, which might provide a theoretical basis for the identification of key genes that contribute to thea biosynthesis.

Transcriptomic analysis of Litopenaeus vannamei hepatopancreas under cold stress in cold-tolerant and cold-sensitive cultivars

Litopenaeus vannamei (L. vannamei) is one of the most widely cultured shrimp species in the world. The species often suffers from cold stress. To understand the molecular mechanism of cold tolerance, we performed transcriptomic analysis on two contrasting cultivars of L. vannamei, namely, cold-tolerant Guihai 2 (GH2) and cold-sensitive Guihai1 (GH1), under a control temperature (28 °C), cold stress (16 °C), and recovery to 28 °C. A total of 84.5 Gb of sequences were generated from 12 L. vannamei hepatopancreas libraries. The de-novo assembly generated a total of 143,029 unigenes with a mean size of 1,052 bp and an N50 of 2,604 bp, of which 34.08% were annotated in the Nr database. We analyzed the differentially expressed genes (DEGs) between nine comparison groups and detected a total of 21,026 DEGs. KEGG pathways, including lysosome, sphingolipid metabolism and nitrogen metabolism, were significantly enriched by DEGs between different temperatures in GH2. Furthermore, eight of the most significantly DEGs under cold stress from the transcriptomic analysis were selected for quantitative real-time PCR (qPCR) validation. Overall, we compared gene expression changes under cold stress in cold-tolerant and cold-sensitive L. vannamei for the first time. The results may further extend our understanding of the cold stress-response mechanism in L. vannamei.

CmBES1 is a regulator of boundary formation in chrysanthemum ray florets

Chrysanthemum (Chrysanthemum morifolium) is an ideal model species for studying petal morphogenesis because of the diversity in the flower form across varieties; however, the molecular mechanisms underlying petal development are poorly understood. Here, we show that the brassinosteroid transcription factor BRI1-EMS-SUPPRESSOR 1 (CmBES1) in chrysanthemum (C. morifolium cv. Jinba) is important for organ boundary formation because it represses organ boundary identity genes. Chrysanthemum plants overexpressing CmBES1 displayed increased fusion of the outermost ray florets due to the loss of differentiation of the two dorsal petals, which developed simultaneously with the ventral petals. RNA-seq analysis of the overexpression lines revealed potential genes and pathways involved in petal development, such as CUP-SHAPED COTYLEDON (CUC2), CYCLOIDEA 4 (CYC4), genes encoding MADS-box transcription factors and homeodomain-leucine zippers (HD-Zips) and auxin pathway-related genes. This study characterizes the role of CmBES1 in ray floret development by its modulation of flower development and boundary identity genes in chrysanthemum.

Responses of Manila Grass (Zoysia matrella) to chilling stress: From transcriptomics to physiology

Manila grass (Zoysia matrella), a warm-season turfgrass, usually wilts and browns by late autumn because of low temperature. To elucidate the molecular mechanisms regarding Manila grass responses to cold stress, we performed transcriptome sequencing of leaves exposed to 4°C for 0 (CK), 2h (2h_CT) and 72h (72h_CT) by Illumina technology. Approximately 250 million paired-end reads were obtained and de novo assembled into 82,605 unigenes. A total of 34,879 unigenes were annotated by comparing their sequence to public protein databases. At the 2h- and 72h-cold time points, 324 and 5,851 differentially expressed genes (DEGs) were identified, respectively. Gene ontology (GO) and metabolism pathway (KEGG) enrichment analyses of DEGs indicated that auxin, gibberellins, ethylene and calcium took part in the cold signal transduction in the early period. And in the late cold period, electron transport activities, photosynthetic machinery and activity, carbohydrate and nitrogen metabolism, redox equilibrium and hormone metabolism were disturbed. Low temperature stress triggered high light, drought and oxidative stress. At the physiological level, cold stress induced a decrease in water content, an increase in levels of total soluble sugar, free proline and MDA, and changes in bioactive gibberellins levels, which supported the changes in gene expression. The results provided a large set of sequence data of Manila grass as well as molecular mechanisms of the grass in response to cold stress. This information will be helpful for future study of molecular breeding and turf management.

TeaCoN: a database of gene co-expression network for tea plant (Camellia sinensis)

Background

Tea plant (Camellia sinensis) is one of the world’s most important beverage crops due to its numerous secondary metabolites conferring tea quality and health effects. However, only a small fraction of tea genes (especially for those metabolite-related genes) have been functionally characterized to date. A cohesive bioinformatics platform is thus urgently needed to aid in the functional determination of the remaining genes.

Description

TeaCoN, a database of gene co-expression network for tea plant, was established to provide genome-wide associations in gene co-expression to survey gene modules (i.e., co-expressed gene sets) for a function of interest. TeaCoN featured a comprehensive collection of 261 high-quality RNA-Seq experiments that covered a wide range of tea tissues as well as various treatments for tea plant. In the current version of TeaCoN, 31,968 (94% coverage of the genome) tea gene models were documented. Users can retrieve detailed co-expression information for gene(s) of interest in four aspects: 1) co-expressed genes with the corresponding Pearson correlation coefficients (PCC-values) and statistical P-values, 2) gene information (gene ID, description, symbol, alias, chromosomal location, GO and KEGG annotation), 3) expression profile heatmap of co-expressed genes across seven main tea tissues (e.g., leaf, bud, stem, root), and 4) network visualization of co-expressed genes. We also implemented a gene co-expression analysis, BLAST search function, GO and KEGG enrichment analysis, and genome browser to facilitate use of the database.

Conclusion

The TeaCoN project can serve as a beneficial platform for candidate gene screening and functional exploration of important agronomical traits in tea plant. TeaCoN is freely available at http://teacon.wchoda.com.

Characterization of cold stress responses in different rapeseed ecotypes based on metabolomics and transcriptomics analyses

The winter oilseed ecotype is more tolerant to low temperature than the spring ecotype. Transcriptome and metabolome analyses of leaf samples of five spring Brassica napus L. (B. napus) ecotype lines and five winter B. napus ecotype lines treated at 4 °C and 28 °C were performed. A total of 25,460 differentially expressed genes (DEGs) of the spring oilseed ecotype and 28,512 DEGs of the winter oilseed ecotype were identified after cold stress; there were 41 differentially expressed metabolites (DEMs) in the spring and 47 in the winter oilseed ecotypes. Moreover, more than 46.2% DEGs were commonly detected in both ecotypes, and the extent of the changes were much more pronounced in the winter than spring ecotype. By contrast, only six DEMs were detected in both the spring and winter oilseed ecotypes. Eighty-one DEMs mainly belonged to primary metabolites, including amino acids, organic acids and sugars. The large number of specific genes and metabolites emphasizes the complex regulatory mechanisms involved in the cold stress response in oilseed rape. Furthermore, these data suggest that lipid, ABA, secondary metabolism, signal transduction and transcription factors may play distinct roles in the spring and winter ecotypes in response to cold stress. Differences in gene expression and metabolite levels after cold stress treatment may have contributed to the cold tolerance of the different oilseed ecotypes.

Transcriptomic response to cold of thermophilous medicinal plant Marsdenia tenacissima

Extracts from Marsdenia tenacissima, involving tenacissoside H, I and G, have been used as remedies of cancer, inflammation and asthma. Low temperature serves as one of the main factors constrain the planting expansion and quality of M. tenacissima, but its functional mechanism has been known scarcely for the lack of genomic information and transcriptional profile. Here we investigated the transcriptomic responses of M. tenacissima under cold stress to gain insight into the molecular mechanism of low temperature sensitivity. Total RNAs were collected from samples obtained at 4-time points (after 0, 3, 6 and 48 h cold treatments with 4 °C, respectively), then used for library construction and sequenced on the Illumina Hiseq™ 4000 platform. Passing quality assessments, 500794 transcripts, and 206137 unigenes were de novo assembly out in Trinity v2.4.0, holding contig N50 of 2566 bp and unigene mean length of 754 bp. 44.20% of assembled unigenes were annotated to the well-known public protein database on a basis of sequence similarity. Using statistical comparison of the fragments per kilo base of transcript per million reads mapped (FPKM) values between conditions, 6082 group-specific differentially expressed genes (DEGs) were identified and considered as cold-responsive genes, which contained copious transcription factors and active secondary metabolism. Among them, 43 unigenes were constantly up-regulated expression along with cold time, which mainly implicated in the biosynthesis of secondary metabolites, carbon metabolism, RNA and DNA metabolism. Conversely, 21 unigenes involved in photosynthesis, cell wall, protein degradation, and transporters were downregulated continually with cold timescale. Experimentally, MtEF1α was chosen as the best housekeeping gene. Functional enrichments found that damaging of cold stress on M. tenacissima may be ascribed to inability of photosynthesis, ribsome processing, flavonoid biosynthesis and terpenoids degradation. Correlation analysis between cold induced transcription factors and tenacissoside biosynthesis-related genes indicated that 3β-HSD significant positively correlated with bHLH51, and 4-MSO with NF-YB, GRAS3, Trihelix, FAR1, MYB60, MYBS1, bZIP43. Further promoter clone found MYB-binding site in the promoter of 4-MSO. In view of the reported cold tolerance of MYB60, it is recommended as a potential candidate suitable for future molecular design of exaptation cultivation with high bioactive constituents.

Cold-inducible expression of an Arabidopsis thaliana AP2 transcription factor gene, AtCRAP2, promotes flowering under unsuitable low-temperatures in chrysanthemum

Flowering time is regulated by biotic and abiotic stresses and affected by the ambient temperature. For chrysanthemum, a low ambient growth temperature can cause a flowering delay, which limits the annual commercial production. Therefore, it is important to improve the low-temperature flowering capability of chrysanthemum through genetic modifications. Here, we isolated a natural variation of a CRT/DRE-binding factor (CBF/DREB) 3 gene, CRAP2, from the Arabidopsis thaliana accession Condara (190AV) that encodes a stop codon at position 151 of the CBF3 protein. Unlike AtCBF3, the overexpression AtCRAP2 in Arabidopsis did not cause detectable growth retardation nor delayed flowering and it conferred cold tolerance. The cold-inducible expression of AtCRAP2 in chrysanthemum promoted flowering under short-day conditions with a low 15 °C nighttime temperature. RNA-sequencing of rd29A:AtCRAP2 and qRT-PCR assays of flowering time-related genes and AtCRAP2 expressed at an ambient temperature revealed that AtCRAP2 positively affected SOC1 and FTL3, thereby promoting flowering under low temperature stress and short-day conditions. These results indicate that DREB genes can be used in the genetic engineering of crop plants without accompanying negative effects by modifying the encoded proteins' C termini.

Gene expression profile indicates involvement of uniconazole in Coix lachryma-jobi L. seedlings at low temperature

Uniconazole (UNZ) can alleviate a variety of abiotic stresses such as low temperature. With application of UNZ on Coix lachryma-jobi L. (coix) under low-temperature stress, growth and physiological parameters were investigated in seedlings. Meanwhile, transcriptome profile in coix seedlings was characterized as well. The results showed an increase of 11.90%, 13.59%, and 10.98% in stem diameter, the aboveground and belowground biomass in 5 mg/L uniconazole application group (U3), compared with control check low-temperature group (CKL). Some anti-oxidase activities also show significant difference between CKL and U3 (p < .05). Transcriptome results showed that 3,901 and 1,040 genes had different expression level at control check (CK) and CKL, CKL and U3. A considerable number of different expressing genes (DEGs) related to the plant hormone signal transduction, photosynthesis, reactive oxygen species (ROS)-related genes, and secondary metabolism in response to uniconazole application were identified in this study. The transcriptomic gene expression profiles present a valuable genomic tool to improve studying the molecular mechanisms underlying low-temperature tolerance in coix. At the same time, it would provide a certain basis for the application of UNZ in the production of coix resistance under low temperature.

Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single-nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south-western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.

The CmTCP20 gene regulates petal elongation growth in Chrysanthemum morifolium

Chrysanthemum morifolium is one of the most popular ornamental species worldwide, with high ornamental and economic value. Petal size is an important factor that influences the ornamental value. CmTCP20 is a member of TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORs (TCPs) gene family, which is closely associated with the growth and development of plants. Our previous study found that the expression of CmTCP20 was obviously down-regulated during chrysanthemum petal elongation, but its function in petal elongation has not yet been revealed. We show here that the overexpression CmTCP20 in Arabidopsis and chrysanthemum leads to similar phenotypes, including larger flower buds (or inflorescences) and longer petals. Interestingly, ectopic expression in Schizosaccharomyces pombe yeast cells showed that CmTCP20 could repress cell division and promote cell elongation. Moreover, the yeast two-hybrid, BiFC and pull-down experimental results indicated that CmTCP20 may regulate petal size via interacting with CmJAZ1-like and inducing down-regulation of CmBPE2 gene expression. This study preliminarily clarifies the function of CmTCP20 on chrysanthemum petal elongation, providing the basic theory for improving the ornamental characteristic of chrysanthemum.

Comparative transcriptome analysis to elucidate the enhanced thermotolerance of tea plants (Camellia sinensis) treated with exogenous calcium

The molecular mechanisms regulating calcium-mediated thermotolerance in Camellia sinensis were revealed by RNA-Sequencing. Heat stress is one of the most remarkable abiotic factors limiting the growth and productivity of Camellia sinensis plants. Calcium helps regulate plant responses to various adverse environmental conditions, including heat stress. In this study, the effects of exogenous calcium on the physiological characteristics of heat-stressed C. sinensis were investigated. A calcium pretreatment increased the proline, soluble sugar, Ca²⁺, and chlorophyll contents, but decreased the malondialdehyde content and relative electrical conductivity in C. sinensis leaves under heat stress. Further analysis of the ultra-structure of chloroplasts indicated that heat stress induced accumulation of starch granules and destruction of the stroma lamella in C. sinensis. However, calcium pretreatment counteracted the adverse effects of heat stress on the structure of the photosynthetic apparatus. These results imply that the calcium pretreatment increased C. sinensis thermotolerance. Moreover, RNA-sequencing was applied to characterize the calcium-mediated transcript-level responses to heat stress. A total of 923 differentially expressed genes (DEGs) including 299 up-regulated and 624 down-regulated genes were identified. Functional annotations indicated that these DEGs were primarily related to signal transduction, transcriptional regulation, and post-translational modification. In addition, a C. sinensis gene [CsCML45 (GenBank: KY652927)] encoding a calmodulin-like protein was isolated. The heterologous expression of CsCML45 enhanced the thermotolerance of transgenic Arabidopsis thaliana plants. These results may be useful for characterizing the calcium-mediated molecular mechanism responsible for C. sinensis thermotolerance.

RNA Sequencing Characterizes Transcriptomes Differences in Cold Response Between Northern and Southern Alternanthera philoxeroides and Highlight Adaptations Associated With Northward Expansion

Alternanthera philoxeroides recently expanded its range northwards in China. It is unknown if the range expansion has a genetic and/or epigenetic basis, or merely an environmental basis due to a warming climate. To test these possibilities, we used an RNAseq approach with a common greenhouse design to examine gene expression in individuals from the northern edge and central portion of alligator weed range from China to determine if there were differences in their responses to cold temperatures. We hypothesized that if the recent range expansion was primarily environmental, we would observe few differences or only differences unrelated to low-temperature adaptations. We assembled over 75,000 genes of which over 65,000 had long open reading frames with similarity to sequences from arabidopsis. Differences in expression between northern and southern populations that were both exposed to low temperatures showed similar expression among genes in the C-REPEAT/DRE BINDING FACTOR (CBF) regulon. However, gene set and sub-network enrichment analysis indicated differences in the response of photosynthetic processes and oxidative stress responses were different between the two populations and we relate these differences to cold adaptation. The transcriptome differences in response to cold between the individuals from the two populations is consistent with adaptations potentiating or resulting from selection after expansion into colder environments and may indicate that genetic changes have accompanied the recent northward expansion of A. philoxeroides in China. However, we cannot rule out the possibility of epigenetic changes may have a role in this expansion.

Transcriptome Analysis Reveals the Mechanism of Fluoride Treatment Affecting Biochemical Components in Camellia sinensis

Tea (Camellia sinensis (L.) O. Kuntze), one of the main crops in China, is high in various bioactive compounds including flavonoids, catechins, caffeine, theanine, and other amino acids. C. sinensis is also known as an accumulator of fluoride (F), and the bioactive compounds are affected by F, however, the mechanism remains unclear. Here, the effects of F treatment on the accumulation of F and major bioactive compounds and gene expression were investigated, revealing the molecular mechanisms affecting the accumulation of bioactive compounds by F treatment. The results showed that F accumulation in tea leaves gradually increased under exogenous F treatments. Similarly, the flavonoid content also increased in the F treatment. In contrast, the polyphenol content, free amino acids, and the total catechins decreased significantly. Special amino acids, such as sulfur-containing amino acids and proline, had the opposite trend of free amino acids. Caffeine was obviously induced by exogenous F, while the theanine content peaked after two day-treatment. These results suggest that the F accumulation and content of bioactive compounds were dramatically affected by F treatment. Furthermore, differentially expressed genes (DEGs) related to the metabolism of main bioactive compounds and amino acids, especially the pivotal regulatory genes of catechins, caffeine, and theanine biosynthesis pathways, were identified and analyzed using high-throughput Illumina RNA-Seq technology and qRT-PCR. The expression of pivotal regulatory genes is consistent with the changes of the main bioactive compounds in C. sinensis leaves, indicating a complicated molecular mechanism for the above findings. Overall, these data provide a reference for exploring the possible molecular mechanism of the accumulation of major bioactive components such as flavonoid, catechins, caffeine, theanine and other amino acids in tea leaves in response to fluoride treatment.

Transcriptomic responses to low temperature stress in the Nile tilapia, Oreochromis niloticus

The Nile tilapia, Oreochromis niloticus, is a species of high economic value and extensively cultured. The limited stress tolerance of this species to a low temperature usually leads to mass mortality and great loss. Nevertheless, there is limited information on the molecular mechanisms underlying the susceptibility to low temperature in the tilapia. In this study, tilapia was treated at 28 °C to a lethal temperature of 8 °C by a gradual decrement. Transcriptomic response of the immune organ, kidney, in tilapia was characterized using RNA-seq. In total, 2191 genes were annotated for significant expression, which were mainly associated with metabolism and immunity. Pathway analysis showed that immune-related pathways of phagosome and cell adhesion molecules (CAMs) pathway were significantly down-regulated under low temperature. Moreover, ferroptosis, a significantly changed pathway involved in tissue damage and acute renal failure, is reported here for the first time. The levels of serum parameters associated with kidney damage such as urea and uric acid (UA) increased significantly under low temperature. The immunofluorescence staining of the kidney showed that cell apoptosis occurred at low temperature. The results of the present study indicate that exposure to low temperature can cause kidney disfunction and down-regulate the immune-related pathway in the kidney of tilapia. This study provides new insight into the mechanism of kidney damage in fish under low temperature.

Analysis of the transcriptome data in Litopenaeus vannamei reveals the immune basis and predicts the hub regulation-genes in response to high-pH stress

Soil salinization erodes the farmlands and poses a serious threat to human life, reuse of the saline-alkali lands as cultivated resources becomes increasingly prominent. Pacific white shrimp (Litopenaeus vannamei) is an important farmed aquatic species for the development and utilization of the saline-alkali areas. However, little is known about the adaptation mechanism of this species in terms of high-pH stress. In the present study, a transcriptome analysis on the gill tissues of L. vannamei in response to high-pH stress (pH 9.3 ± 0.1) was conducted. After analyzing, the cyclic nucleotide gated channel-Ca2+ (CNGC-Ca2+) and patched 1 (Ptc1) were detected as the majority annotated components in the cAMP signaling pathway (KO04024), indicating that the CNGC-Ca2+ and Ptc1 might be the candidate components for transducing and maintaining the high-pH stress signals, respectively. The immunoglobulin superfamily (IgSF), heat shock protein (HSP), glutathione s-transferase (GST), prophenoloxidase/phenoloxidase (proPO/PO), superoxide dismutase (SOD), anti-lipopolysaccharide factor (ALF) and lipoprotein were discovered as the major transcribed immune factors in response to high-pH stress. To further detect hub regulation-genes, protein-protein interaction (PPI) networks were constructed; the genes/proteins "Polymerase (RNA) II (DNA directed) polypeptide A" (POLR2A), "Histone acetyltransferase p300" (EP300) and "Heat shock 70kDa protein 8" (HSPA8) were suggested as the top three hub regulation-genes in response to acute high-pH stress; the genes/proteins "Heat shock 70kDa protein 4" (HSPA4), "FBJ murine osteosarcoma viral oncogene homolog" (FOS) and "Nucleoporin 54kDa" (NUP54) were proposed as the top three hub regulation-genes involved in adapting endurance high-pH stress; the protein-interactions of "EP300-HSPA8" and "HSPA4-NUP54" were detected as the most important biological interactions in response to the high-pH stress; and the HSP70 family genes might play essential roles in the adaptation of the high-pH stress environment in L. vannamei. These findings provide the first insight into the molecular and immune basis of L. vannamei in terms of high-pH environments, and the construction of a PPI network might improve our understanding in revealing the hub regulation-genes in response to abiotic stress in shrimp species and might be beneficial for further studies.

The Chrysanthemum nankingense Genome Provides Insights into the Evolution and Diversification of Chrysanthemum Flowers and Medicinal Traits

The Asteraceae (Compositae), a large plant family of approximately 24 000-35 000 species, accounts for ∼10% of all angiosperm species and contributes a lot to plant diversity. The most representative members of the Asteraceae are the economically important chrysanthemums (Chrysanthemum L.) that diversified through reticulate evolution. Biodiversity is typically created by multiple evolutionary mechanisms such as whole-genome duplication (WGD) or polyploidization and locally repetitive genome expansion. However, the lack of genomic data from chrysanthemum species has prevented an in-depth analysis of the evolutionary mechanisms involved in their diversification. Here, we used Oxford Nanopore long-read technology to sequence the diploid Chrysanthemum nankingense genome, which represents one of the progenitor genomes of domesticated chrysanthemums. Our analysis revealed that the evolution of the C. nankingense genome was driven by bursts of repetitive element expansion and WGD events including a recent WGD that distinguishes chrysanthemum from sunflower, which diverged from chrysanthemum approximately 38.8 million years ago. Variations of ornamental and medicinal traits in chrysanthemums are linked to the expansion of candidate gene families by duplication events including paralogous gene duplication. Collectively, our study of the assembled reference genome offers new knowledge and resources to dissect the history and pattern of evolution and diversification of chrysanthemum plants, and also to accelerate their breeding and improvement.

Cellular and molecular characteristics of pollen abortion in chrysanthemum cv. Kingfisher

Microspore degeneration at the tetrad stage is associated with tapetum degeneration retardation. Some genes and proteins related to cell senescence and death are the key factors for pollen abortion. Chrysanthemum (Chrysanthemum morifolium) is a major floriculture crop in the world, but pollen contamination is an urgent problem to be solved in chrysanthemum production. C. morifolium 'Kingfisher' is a chrysanthemum cultivar that does not contain any pollen in mature anthers, thus it is a very important material for developing chrysanthemum without pollen contamination. However, the mechanism of its pollen abortion remains unclear. In this study, the cellular and molecular mechanisms of 'Kingfisher' pollen abortion were investigated using transmission electron microscopy, RNA sequencing, isobaric tags for relative and absolute quantitation, and bioinformatics. It was found that the meiosis of microspore mother cells was normal before the tetrad stage, the microspores began to degenerate at the tetrad stage, and no microspores were observed in the anthers after the tetrad stage. In addition, transcriptomic and proteomic analyses showed that some genes and proteins related to cell senescence and death were identified to be implicated in chrysanthemum pollen abortion. These results indicated that the tetrad stage was the main period of pollen abortion, and the genes and proteins related to cell senescence and death contributed to pollen abortion. These add to our understanding of chrysanthemum pollen abortion and will be helpful for development of flowers without pollen contamination in the future.

sRNAome and transcriptome analysis provide insight into chilling response of cowpea pods

Cowpea is an important horticultural crop in tropical and subtropical areas of Asia, Africa, and Latin America, as well as parts of southern Europe and Central and South America. Chilling injury is a common physiological hazard of cowpea in cold chain logistics which reduce the cowpea pods nutritional quality and product value. However, the molecular mechanism involved in chilling injury remains unclear in cowpea pods. RNA-Seq and sRNA-Seq technologies were employed to decipher the miRNAs and mRNAs expression profiles and their regulatory networks in cowpea pods involved in chilling stress. Differentially expressed miRNAs and mRNA profiles were obtained based on cluster analysis, miRNAs and target genes were found to show coherent relationships in the regulatory networks of chilling injury. Furthermore, we found that numerous miRNAs and nat-siRNAs' targets were predicted to be key enzymes involved in the redox reactions such as POD, CAT, AO and LOX, energy metabolism such as ATPase, FAD and NAD related enzymes and different transcription factors such as WRKY, bHLH, MYB, ERF and NAC which play important roles in chilling injury.

Comparative Transcriptome Analysis of Waterlogging-Sensitive and Waterlogging-Tolerant Chrysanthemum morifolium Cultivars under Waterlogging Stress and Reoxygenation Conditions

Waterlogging stress is among the most severe abiotic stressors in the breeding and the production of Chrysanthemum morifolium. However, the mechanism underlying the response to waterlogging and post-waterlogging reoxygenation in C. morifolium remains unknown. In this study, we compared the differences between the transcriptomes of two chrysanthemum cultivars, i.e., the waterlogging-tolerant cultivar "Nannongxuefeng" and the waterlogging-sensitive cultivar "Qinglu", by performing RNA-seq to elucidate the possible mechanism of waterlogging and reoxygenation in C. morifolium. "Nannongxuefeng" had a higher ethylene production under the waterlogging and reoxygenation conditions. Furthermore, the expression of transcription factors and genes that are involved in the hormone response, N-end rule pathway and ROS signaling significantly differed between the two cultivars. "Nannongxuefeng" and "Qinglu" significantly differed in their response to waterlogging and reoxygenation, providing a deeper understanding of the mechanism underlying the response to waterlogging and guidance for the breeding of C. morifolium.

Overexpression of Phosphate Transporter Gene CmPht1;2 Facilitated Pi Uptake and Alternated the Metabolic Profiles of Chrysanthemum Under Phosphate Deficiency

Low availability of phosphorus (P) in the soil is the principal limiting factor for the growth of cut chrysanthemum. Plant phosphate transporters (PTs) facilitate acquisition of inorganic phosphate (Pi) and its homeostasis within the plant. In the present study, CmPht1;2 of the Pht1 family was cloned from chrysanthemum. CmPht1;2 is composed of 12 transmembrane domains and localized to the plasma membrane. Expression of CmPht1;2 in roots was induced by Pi starvation. Chrysanthemum plants with overexpression of CmPht1;2 (Oe) showed higher Pi uptake, as compared to the wild type (WT), both under Pi-starvation and Pi-sufficient conditions, and also showed a higher root biomass compared to WT in the Pi-starvation conditions. Seven days after the P-deficiency treatment, 85 distinct analytes were identified in the roots and 27 in the shoots between the Oe1 plant and WT, in which sophorose, sorbitol (sugars), hydroxybutyric acid (organic acids), and ornithine (amino acid) of CmPht1;2 overexpressing chrysanthemum are specific responses to P-starvation.

Genomic and transcriptomic alterations following intergeneric hybridization and polyploidization in the Chrysanthemum nankingense×Tanacetum vulgare hybrid and allopolyploid (Asteraceae)

Allopolyploid formation involves two major events: interspecific hybridization and polyploidization. A number of species in the Asteraceae family are polyploids because of frequent hybridization. The effects of hybridization on genomics and transcriptomics in Chrysanthemum nankingense×Tanacetum vulgare hybrids have been reported. In this study, we obtained allopolyploids by applying a colchicine treatment to a synthesized C. nankingense×T. vulgare hybrid. Sequence-related amplified polymorphism (SRAP), methylation-sensitive amplification polymorphism (MSAP), and high-throughput RNA sequencing (RNA-Seq) technologies were used to investigate the genomic, epigenetic, and transcriptomic alterations in both the hybrid and allopolyploids. The genomic alterations in the hybrid and allopolyploids mainly involved the loss of parental fragments and the gain of novel fragments. The DNA methylation level of the hybrid was reduced by hybridization but was restored somewhat after polyploidization. There were more significant differences in gene expression between the hybrid/allopolyploid and the paternal parent than between the hybrid/allopolyploid and the maternal parent. Most differentially expressed genes (DEGs) showed down-regulation in the hybrid/allopolyploid relative to the parents. Among the non-additive genes, transgressive patterns appeared to be dominant, especially repression patterns. Maternal expression dominance was observed specifically for down-regulated genes. Many methylase and methyltransferase genes showed differential expression between the hybrid and parents and between the allopolyploid and parents. Our data indicate that hybridization may be a major factor affecting genomic and transcriptomic changes in newly formed allopolyploids. The formation of allopolyploids may not simply be the sum of hybridization and polyploidization changes but also may be influenced by the interaction between these processes.

Overexpression of GhDof1 improved salt and cold tolerance and seed oil content in Gossypium hirsutum

A homologous GhDof1, which belongs to a large family of plant-specific transcription factor DOF, was isolated from Upland cotton (Gossypium hirsutum L.). GhDof1 protein was located in the nucleus of onion epidermal cells, the core domain of transcriptional activity existed in the C-terminal, and the activity elements of GhDof1 promoter existed in the regions of -645∼ -469bp and -286∼ -132bp of transcriptional start codon. GhDof1 constitutively expressed in leaves, roots and stems, accumulated highest in leaves. The salinity and cold treatments induced GhDof1 transcript accumulation. The GhDof1-overexpressed cotton showed significantly higher salt and cold tolerance over the wild-type plants. Under salt stress, the root growth of overexpressed GhDof1 lines was promoted. The expression levels of stress-responsive genes, GhP5CS, GhSOD and GhMYB, were differently up-regulated in transgenic lines. Oil contents increased in some transgenic plants, and protein contents reduced compared with transformed receptor. These results suggested that GhDof1 was a functional transcription factor for improving the abiotic tolerance and seed oil content in Upland cotton.

De novo transcriptome sequencing of Isaria cateniannulata and comparative analysis of gene expression in response to heat and cold stresses

Isaria cateniannulata is a very important and virulent entomopathogenic fungus that infects many insect pest species. Although I. cateniannulata is commonly exposed to extreme environmental temperature conditions, little is known about its molecular response mechanism to temperature stress. Here, we sequenced and de novo assembled the transcriptome of I. cateniannulata in response to high and low temperature stresses using Illumina RNA-Seq technology. Our assembly encompassed 17,514 unigenes (mean length = 1,197 bp), in which 11,445 unigenes (65.34%) showed significant similarities to known sequences in NCBI non-redundant protein sequences (Nr) database. Using digital gene expression analysis, 4,483 differentially expressed genes (DEGs) were identified after heat treatment, including 2,905 up-regulated genes and 1,578 down-regulated genes. Under cold stress, 1,927 DEGs were identified, including 1,245 up-regulated genes and 682 down-regulated genes. The expression patterns of 18 randomly selected candidate DEGs resulting from quantitative real-time PCR (qRT-PCR) were consistent with their transcriptome analysis results. Although DEGs were involved in many pathways, we focused on the genes that were involved in endocytosis: In heat stress, the pathway of clathrin-dependent endocytosis (CDE) was active; however at low temperature stresses, the pathway of clathrin-independent endocytosis (CIE) was active. Besides, four categories of DEGs acting as temperature sensors were observed, including cell-wall-major-components-metabolism-related (CWMCMR) genes, heat shock protein (Hsp) genes, intracellular-compatible-solutes-metabolism-related (ICSMR) genes and glutathione S-transferase (GST). These results enhance our understanding of the molecular mechanisms of I. cateniannulata in response to temperature stresses and provide a valuable resource for the future investigations.

'Bartlett' pear fruit (Pyrus communis L.) ripening regulation by low temperatures involves genes associated with jasmonic acid, cold response, and transcription factors

Low temperature (LT) treatments enhance ethylene production and ripening rate in the European pear (Pyrus communis L.). However, the underlying molecular mechanisms are not well understood. This study aims to identify genes responsible for ripening enhancement by LT. To this end, the transcriptome of 'Bartlett' pears treated with LT (0°C or 10°C for up to 14 d), which results in faster ripening, and control pears without conditioning treatment was analyzed. LT conditioned pears reached eating firmness (18N) in 6 d while control pears took about 12 d when left to ripen at 20°C. We identified 8,536 differentially expressed (DE) genes between the 0°C-treated and control fruit, and 7,938 DE genes between the 10°C-treated and control fruit. In an attempt to differentiate temperature-induced vs. ethylene-responsive pathways, we also monitored gene expression in fruit sequentially treated with 1-MCP then exposed to low temperature. This analysis revealed that genes associated with jasmonic acid biosynthesis and signaling, as well as the transcription factors TCP9a, TCP9b, CBF1, CBF4, AGL24, MYB1R1, and HsfB2b could be involved in the LT-mediated enhancement of ripening independently or upstream of ethylene.

Gibberellic Acid Signaling Is Required to Induce Flowering of Chrysanthemums Grown under Both Short and Long Days

Flower bud formation and flowering in chrysanthemums occur under short day conditions (SD), but the molecular basis for the switch to reproductive growth is less well understood than in model plants. Here, a spontaneous mutant able to flower under long days is described. In an attempt to reveal the pathway(s) involved in the formation of flower buds under contrasting daylengths, transcriptome sequencing was carried out in plants grown both under SD and long day conditions (LD). A number of differentially transcribed genes involved in the various known flowering pathways were identified. Both circadian clock genes and Chrysanthemum FLOWERING LOCUS T Like3 (CmFTL3) were up-regulated under SD, thereby inducing floral bud formation and flowering. The gibberellin (GA) signaling pathway-related genes Gibberellin 20-oxidase (GA20ox) and Gibberellin receptor (GID1) were up-regulated in the mutant under LD, while the catabolic genes Gibberellin 2-oxidase (GA2ox) and GA-INSENSITIVE (GAI) were down-regulated, thereby inducing the transcription of SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) and LEAFY (LFY). The GA content of the leaf was higher in the mutant than in the wild type (WT) under LD and SD, and the mutant has more branching than WT plants under LD or SD. When treated with GA, the mutant flowered earlier under both SD and LD relative to WT, but there was no detectable phenotype difference between the two lines. The indication was that the photoperiod pathway majorly regulates flower bud formation and flowering time in chrysanthemums under SD. The GA signaling pathway only plays a subsidiary role for flowering. However, the GA signaling pathway predominated for flowering under LD.

Transcriptomic analyses on muscle tissues of Litopenaeus vannamei provide the first profile insight into the response to low temperature stress

The Pacific white shrimp (Litopenaeus vannamei) is an important cultured crustacean species worldwide. However, little is known about the molecular mechanism of this species involved in the response to cold stress. In this study, four separate RNA-Seq libraries of L. vannamei were generated from 13°C stress and control temperature. Total 29,662 of Unigenes and overall of 19,619 annotated genes were obtained. Three comparisons were carried out among the four libraries, in which 72 of the top 20% of differentially-expressed genes were obtained, 15 GO and 5 KEGG temperature-sensitive pathways were fished out. Catalytic activity (GO: 0003824) and Metabolic pathways (ko01100) were the most annotated GO and KEGG pathways in response to cold stress, respectively. In addition, Calcium, MAPK cascade, Transcription factor and Serine/threonine-protein kinase signal pathway were picked out and clustered. Serine/threonine-protein kinase signal pathway might play more important roles in cold adaptation, while other three signal pathway were not widely transcribed. Our results had summarized the differentially-expressed genes and suggested the major important signaling pathways and related genes. These findings provide the first profile insight into the molecular basis of L. vannamei response to cold stress.

Transcriptomic and hormone analyses reveal mechanisms underlying petal elongation in Chrysanthemum morifolium 'Jinba'

Auxin regulates chrysanthemum petal elongation by promoting cell elongation. Transcriptomic analysis shows that auxin signal transduction may connect with other transcription factors by TCPs to regulate chrysanthemum petal elongation. As an ornamental species, Chrysanthemum morifolium has high ornamental and economic value. Petal size is the primary factor that influences the ornamental value of chrysanthemum, but the mechanism underlying the development of C. morifolium petals remains unclear. In our study, we tracked the growth of petals and found that the basal region of 'Jinba' petals showed a higher elongation rate, exhibiting rapid cell elongation during petal growth. During petal elongation growth, auxin was demonstrated to promote cell elongation and an increase in cell numbers in the petal basal region. To further study the molecular mechanisms underlying petal growth, the RNA-seq (high-throughput cDNA sequencing) technique was employed. Four cDNA libraries were assembled from petals in the budding, bud breaking, early blooming and full blooming stages of 'Jinba' flower development. Analysis of differentially expressed genes (DEGs) showed that auxin was the most important regulator in controlling petal growth. The TEOSINTEBRANCHED 1, CYCLOIDEA and PCF transcription factor genes (TCPs), basic helix-loop-helix-encoding gene (bHLH), glutaredoxin-C (GRXC) and other zinc finger protein genes exhibited obvious up-regulation and might have significant effects on the growth of 'Jinba' petals. Given the interaction between these genes in Arabidopsis thaliana, we speculated that auxin signal transduction might exhibit a close relationship with transcription factors through TCPs. In summary, we present the first comprehensive transcriptomic and hormone analyses of C. morifolium petals. The results offer direction in identifying the mechanism underlying the development of chrysanthemum petals in the elongated phase and have great significance in improving the ornamental characteristics of C. morifolium via molecular breeding.

Transcriptomic analysis of differentially expressed genes in the floral transition of the summer flowering chrysanthemum

Background

Chrysanthemum is a leading cut flower species. Most conventional cultivars flower during the fall, but the Chrysanthemum morifolium ‘Yuuka’ flowers during the summer, thereby filling a gap in the market. To date, investigations of flowering time determination have largely focused on fall-flowering types. Little is known about molecular basis of flowering time in the summer-flowering chrysanthemum. Here, the genome-wide transcriptome of ‘Yuuka’ was acquired using RNA-Seq technology, with a view to shedding light on the molecular basis of the shift to reproductive growth as induced by variation in the photoperiod.

Results

Two sequencing libraries were prepared from the apical meristem and leaves of plants exposed to short days, three from plants exposed to long days and one from plants sampled before any photoperiod treatment was imposed. From the ~316 million clean reads obtained, 115,300 Unigenes were assembled. In total 70,860 annotated sequences were identified by reference to various databases. A number of transcription factors and genes involved in flowering pathways were found to be differentially transcribed. Under short days, genes acting in the photoperiod and gibberellin pathways might accelerate flowering, while under long days, the trehalose-6-phosphate and sugar signaling pathways might be promoted, while the phytochrome B pathway might block flowering. The differential transcription of eight of the differentially transcribed genes was successfully validated using quantitative real time PCR.

Conclusions

A transcriptome analysis of the summer-flowering cultivar ‘Yuuka’ has been described, along with a global analysis of floral transition under various daylengths. The large number of differentially transcribed genes identified confirmed the complexity of the regulatory machinery underlying floral transition.

Electronic supplementary material

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

Comparative Analysis of the Chrysanthemum Leaf Transcript Profiling in Response to Salt Stress

Salt stress has some remarkable influence on chrysanthemum growth and productivity. To understand the molecular mechanisms associated with salt stress and identify genes of potential importance in cultivated chrysanthemum, we carried out transcriptome sequencing of chrysanthemum. Two cDNA libraries were generated from the control and salt-treated samples (Sample_0510_control and Sample_0510_treat) of leaves. By using the Illumina Solexa RNA sequencing technology, 94 million high quality sequencing reads and 161,522 unigenes were generated and then we annotated unigenes through comparing these sequences to diverse protein databases. A total of 126,646 differentially expressed transcripts (DETs) were identified in leaf. Plant hormones, amino acid metabolism, photosynthesis and secondary metabolism were all changed under salt stress after the complete list of GO term and KEGG enrichment analysis. The hormone biosynthesis changing and oxidative hurt decreasing appeared to be significantly related to salt tolerance of chrysanthemum. Important protein kinases and major transcription factor families involved in abiotic stress were differentially expressed, such as MAPKs, CDPKs, MYB, WRKY, AP2 and HD-zip. In general, these results can help us to confirm the molecular regulation mechanism and also provide us a comprehensive resource of chrysanthemum under salt stress.

Whole-transcriptome analysis of differentially expressed genes in the ray florets and disc florets of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most important global cut flower and pot plants, and has been cultivated worldwide. However, limited genomic resources are available and the molecular mechanisms involved in the two morphologically distinct floret developmental cycles in chrysanthemum remain unclear.

RESULTS

The transcriptomes of chrysanthemum ray florets, disc florets and leaves were sequenced using Illumina paired-end sequencing technology. In total, 16.9 G reads were assembled into 93,138 unigenes with an average length of 738 bp, of which 44,364 unigenes showed similarity to known proteins in the Swissprot or NCBI non-redundant protein databases. Additionally, 26,320, 22,304 and 13,949 unigenes were assigned to 54 gene ontology (GO) categories, 25 EuKaryotic Orthologous Groups (KOG) categories, and 280 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. A total of 1863 differentially expressed genes (DEGs) (1210 up-regulated and 653 down-regulated) were identified between ray florets and disc florets, including genes encoding transcription factors and protein kinases. GO and KEGG pathway enrichment analyses were performed on the DEGs to identify differences in the biological processes and pathways between ray florets and disc florets. The important regulatory genes controlling flower development and flower organ determination, as well as important functional genes in the anthocyanin biosynthetic pathway, were identified, of which two leucoanthocyanidin dioxygenase-encoding genes showed specific expression in ray florets. Lastly, reverse transcription quantitative PCR was conducted to validate the DEGs identified in our study.

CONCLUSIONS

Comparative transcriptome analysis revealed significant differences in patterns of gene expression and signaling pathways between ray florets and disc florets in Chrysanthemum morifolium. This study provided the first step to understanding the molecular mechanism of the differential development of ray florets and disc florets in chrysanthemum, and also provided valuable genomic resources for candidate genes applicable for the breeding of novel varieties in chrysanthemum.