An RNA sequencing transcriptome analysis of the high-temperature stressed tall fescue reveals novel insights into plant thermotolerance

Adequate water supply enhances seedling growth and metabolism in Festuca kryloviana: insights from physiological and transcriptomic analys

BACKGROUND

Festuca kryloviana is a significant native grass species in the Qinghai Lake region, and its low emergence rate is a primary factor limiting the successful establishment of cultivated grasslands. The region's arid and low-rainfall climate characteristics result in reduced soil moisture content at the surface. Despite the recognized impact of water availability on plant growth, the specific role of moisture in seedling development remains not fully elucidated. This study aims to investigate the germination rate and seedling growth velocity of F. kryloviana seeds under varying moisture conditions, and to integrate physiological and transcriptomic analyses of seedlings under these conditions to reveal the mechanisms by which water influences seedling development.

RESULTS

The emergence rate of F. kryloviana seedlings exhibited an initial increase followed by a decrease with increasing moisture content. The highest emergence rate, reaching 75%, was observed under 20% soil moisture conditions. By the eighth day of the experiment, the lengths of the plumules and radicles under the optimal emergence rate (full water, FW) were 21.82% and 10.87% longer, respectively, than those under closely matching the soil moisture content during the background survey (stress water, SW). The differential development of seedlings under varying moisture regimes is attributed to sugar metabolism within the seeds and the accumulation of abscisic acid (ABA). At FW conditions, enhanced sugar metabolism, which generates more energy for seedling development, is facilitated by higher activities of α-amylase, sucrose synthase, and trehalose-6-phosphate synthase compared to SW conditions. This is reflected at the transcriptomic level with upregulated expression of the α-amylase (AMY2) gene and trehalose-6-phosphate synthase (TPS6), while genes associated with ABA signaling and transduction are downregulated. Additionally, under FW conditions, the expression of genes related to the chloroplast thylakoid photosystems, such as photosystem II (PSII) and photosystem I (PSI), is upregulated, enhancing the seedlings' light-capturing ability and photosynthetic efficiency, thereby improving their autotrophic capacity. Furthermore, FW treatment enhances the expression of the non-enzymatic antioxidant system, promoting metabolism within the seeds. In contrast, SW treatment increases the activity of the enzymatic antioxidant system, including peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT), to cope with water stress.

CONCLUSIONS

Our experiment systematically evaluated the impact of moisture conditions on the growth and development of F. kryloviana seedlings. Physiological and transcriptomic data collectively indicate that adequate water (20%) supply enhances seedling growth and development by reducing ABA levels and increasing α-amylase activity within seeds, thereby boosting sugar metabolism and promoting the growth of seedling, which in turn leads to an improved emergence rate. Considering water management in future cultivation practices may be a crucial strategy for enhancing the successful establishment of F. kryloviana in grassland ecosystems.

Heat stress during male meiosis impairs cytoskeletal organization, spindle assembly and tapetum degeneration in wheat

The significance of heat stress in agriculture is ever-increasing with the progress of global climate changes. Due to a negative effect on the yield of staple crops, including wheat, the impairment of plant reproductive development triggered by high ambient temperature became a restraint in food production. Although the heat sensitivity of male meiosis and the following gamete development in wheat has long been recognized, a detailed structural characterization combined with a comprehensive gene expression analysis has not been done about this phenomenon. We demonstrate here that heat stress severely alters the cytoskeletal configuration, triggers the failure of meiotic division in wheat. Moreover, it changes the expression of genes related to gamete development in male meiocytes and the tapetum layer in a genotype-dependent manner. 'Ellvis', a heat-tolerant winter wheat cultivar, showed high spikelet fertility rate and only scarce structural aberrations upon exposure to high temperature. In addition, heat shock genes and genes involved in scavenging reactive oxygen species were significantly upregulated in 'Ellvis', and the expression of meiosis-specific and major developmental genes showed high stability in this cultivar. In the heat-sensitive 'Mv 17-09', however, genes participating in cytoskeletal fiber nucleation, the spindle assembly checkpoint genes, and tapetum-specific developmental regulators were downregulated. These alterations may be related to the decreased cytoskeleton content, frequent micronuclei formation, and the erroneous persistence of the tapetum layer observed in the sensitive genotype. Our results suggest that understanding the heat-sensitive regulation of these gene functions would be an essential contribution to the development of new, heat-tolerant cultivars.

Deficiency in NDH-cyclic electron transport retards heat acclimation of photosynthesis in tobacco over day and night shift

In order to cope with the impact of global warming and frequent extreme weather, thermal acclimation ability is particularly important for plant development and growth, but the mechanism behind is still not fully understood. To investigate the role of NADH dehydrogenase-like complex (NDH) mediated cyclic electron flow (CEF) contributing to heat acclimation, wild type (WT) tobacco (Nicotiana tabacum) and its NDH-B or NDH-C, J, K subunits deficient mutants (ΔB or ΔCJK) were grown at 25/20°C before being shifted to a moderate heat stress environment (35/30°C). The photosynthetic performance of WT and ndh mutants could all eventually acclimate to the increased temperature, but the acclimation process of ndh mutants took longer. Transcriptome profiles revealed that ΔB mutant exhibited distinct photosynthetic-response patterns and stress-response genes compared to WT. Metabolite analysis suggested over-accumulated reducing power and production of more reactive oxygen species in ΔB mutant, which were likely associated with the non-parallel recovery of CO2 assimilation and light reactions shown in ΔB mutant during heat acclimation. Notably, in the warm night periods that could happen in the field, NDH pathway may link to the re-balance of excess reducing power accumulated during daytime. Thus, understanding the diurnal cycle contribution of NDH-mediated CEF for thermal acclimation is expected to facilitate efforts toward enhanced crop fitness and survival under future climates.

Comparative transcriptomic analysis of the gene expression and underlying molecular mechanism of submergence stress response in orchardgrass roots

INTRODUCTION

Submergence stress creates a hypoxic environment. Roots are the first plant organ to face these low-oxygen conditions, which causes damage and affects the plant growth and yield. Orchardgrass (Dactylis glomerata L.) is one of the most important cold-season forage grasses globally. However, their submergence stress-induced gene expression and the underlying molecular mechanisms of orchardgrass roots are still unknown.

METHODS

Using the submergence-tolerant 'Dianbei' and submergence-sensitive 'Anba', the transcriptomic analysis of orchardgrass roots at different time points of submergence stress (0 h, 8 h, and 24 h) was performed.

RESULTS

We obtained 118.82Gb clean data by RNA-Seq. As compared with the control, a total of 6663 and 9857 differentially expressed genes (DEGs) were detected in Dianbei, while 7894 and 11215 DEGs were detected in Anba at 8 h and 24 h post-submergence-stress, respectively. Gene Ontology (GO) enrichment analysis obtained 986 terms, while Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis obtained 123 pathways. Among them, the DEGs in plant hormones, mitogen-activated protein kinase (MAPK) and Ca²⁺ signal transduction were significantly differentially expressed in Dianbei, but not in Anba.

DISCUSSION

This study was the first to molecularly elucidate the submergence stress tolerance in the roots of two orchardgrass cultivars. These findings not only enhanced our understanding of the orchardgrass submergence tolerance, but also provided a theoretical basis 36 for the cultivation of submergence-tolerant forage varieties.

Transcriptome analysis of perennial ryegrass reveals the regulatory role of Aspergillus aculeatus under salt stress

Perennial ryegrass (Lolium perenne) is an important turf grass and forage grass with moderately tolerant to salinity stress. Aspergillus aculeatus has been documented to involved in salt stress response of perennial ryegrass, while the A. aculeatus-mediated molecular mechanisms are unclear. Therefore, to investigate the molecular mechanisms underlying A. aculeatus-mediated salt tolerance, the comprehensive transcriptome analysis of the perennial ryegrass roots was performed. Twelve cDNA libraries from roots were constructed after 12 h of plant-fungus cocultivation under 300 mM salt stress concentrations. A total of 21,915 differentially expressed genes (DEGs) were identified through pairwise comparisons. Enrichment analysis revealed that potentially important A. aculeatus-induced salt responsive genes belonging to specific categories, such as hormonal metabolism (auxin and salicylic acid metabolism related genes), secondary metabolism (flavonoid's metabolism related genes) and transcription factors (MYB, HSF and AP2/EREBP family). In addition, weighted gene co-expression network analysis (WGCNA) showed that blue and black modules were significantly positively correlated with the peroxidase activity and proline content, then the hub genes within these two modules were further identified. Taken together, we found the categories of A. aculeatus-induced salt responsive genes, revealing underlying fungus-induced molecular mechanisms of salt stress response in perennial ryegrass roots. Besides, fungus-induced salt-tolerant hub genes represent a foundation for further exploring the molecular mechanisms.

Morpho-histology, endogenous hormone dynamics, and transcriptome profiling in Dacrydium pectinatum during female cone development

Dacrydium pectinatum de Laubenfels is a perennial dioeciously gymnosperm species dominant in tropical montane rain forests. Due to deforestation, natural disasters, long infancy, and poor natural regeneration ability, the population of this species has been significantly reduced and listed as an endangered protected plant. To better understand the female cone development in D. pectinatum, we examined the morphological and anatomical changes, analyzed the endogenous hormone dynamics, and profiled gene expression. The female reproductive structures were first observed in January. The morpho-histological observations suggest that the development of the D. pectinatum megaspore can be largely divided into six stages: early flower bud differentiation, bract primordium differentiation, ovule primordium differentiation, dormancy, ovule maturity, and seed maturity. The levels of gibberellins (GA), auxin (IAA), abscisic acid (ABA), and cytokinin (CTK) fluctuate during the process of female cone development. The female cones of D. pectinatum need to maintain a low level of GA3-IAA-ABA steady state to promote seed germination. The first transcriptome database for female D. pectinatum was generated, revealing 310,621 unigenes. Differential expression analyses revealed several floral (MADS2, AGL62, and LFY) and hormone biosynthesis and signal transduction (CKX, KO, KAO, ABA4, ACO, etc.) genes that could be critical for female cone development. Our study provides new insights into the cone development in D. pectinatum and the foundation for female cone induction with hormones.

Identification of Heat-Tolerant Genes in Non-Reference Sequences in Rice by Integrating Pan-Genome, Transcriptomics, and QTLs

The availability of large-scale genomic data resources makes it very convenient to mine and analyze genes that are related to important agricultural traits in rice. Pan-genomes have been constructed to provide insight into the genome diversity and functionality of different plants, which can be used in genome-assisted crop improvement. Thus, a pan-genome comprising all genetic elements is crucial for comprehensive variation study among the heat-resistant and -susceptible rice varieties. In this study, a rice pan-genome was firstly constructed by using 45 heat-tolerant and 15 heat-sensitive rice varieties. A total of 38,998 pan-genome genes were identified, including 37,859 genes in the reference and 1141 in the non-reference contigs. Genomic variation analysis demonstrated that a total of 76,435 SNPs were detected and identified as the heat-tolerance-related SNPs, which were specifically present in the highly heat-resistant rice cultivars and located in the genic regions or within 2 kbp upstream and downstream of the genes. Meanwhile, 3214 upregulated and 2212 downregulated genes with heat stress tolerance-related SNPs were detected in one or multiple RNA-seq datasets of rice under heat stress, among which 24 were located in the non-reference contigs of the rice pan-genome. We then mapped the DEGs with heat stress tolerance-related SNPs to the heat stress-resistant QTL regions. A total of 1677 DEGs, including 990 upregulated and 687 downregulated genes, were mapped to the 46 heat stress-resistant QTL regions, in which 2 upregulated genes with heat stress tolerance-related SNPs were identified in the non-reference sequences. This pan-genome resource is an important step towards the effective and efficient genetic improvement of heat stress resistance in rice to help meet the rapidly growing needs for improved rice productivity under different environmental stresses. These findings provide further insight into the functional validation of a number of non-reference genes and, especially, the two genes identified in the heat stress-resistant QTLs in rice.

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.

Integrated physiological and transcriptomic analyses of two warm- and cool-season turfgrass species in response to heat stress

Warm- and cool-season turfgrasses were originated from different locations with contrasting heat tolerance. The molecular mechanisms of heat tolerance have not been extensively studied in turfgrass species. In this study, transcriptomic analysis showed that bermudagrass was more tolerant to heat stress as evidenced by lower contents of H₂O₂, proline and glutathione than those in tall fescue after heat treatment. RNA sequencing analysis revealed that 32.7% and 17.7% more genes were changed in tall fescue than in bermudagrass after 2 and 12h heat treatment, respectively. GO terms of redox were enriched in bermudagrass whereas metabolite transportation ones were over-represented in tall fescue after 2h treatment. Ubiquitin dependent degradation pathways were commonly regulated in both grass species. CdF-box and FaF-box transgenic Arabidopsis exhibited improved tolerance to heat stress. Regulatory elements analysis revealed that four ABA-responsive elements present in CdF-box promoter, indicating CdF-box could be potentially regulated by ABRE binding factors (ABFs). All these findings provide evidences for understanding heat stress response in warm- and cool-season grass species.

FaHSP17.8-CII orchestrates lead tolerance and accumulation in shoots via enhancing antioxidant enzymatic response and PSII activity in tall fescue

Tall fescue (Festuca arundinacea Schreb.) shows huge potential for lead (Pb) phytoremediation, while little is known on the molecular mechanisms involved in Pb tolerance and accumulation. Here, genetic engineering strategy was firstly used to investigate Pb tolerance and accumulation in tall fescue. The transgenic tall fescue overexpressing a class II (CII) sHSP gene FaHSP17.8-CII was generated. After exposure to 1000 mg/L Pb(NO₃)₂, two FaHSP17.8-CII overexpressing lines, OE#3 and OE#7, showed higher tolerance to Pb as illustrated by the reduced levels of electrolyte leakage (EL) and malondialdehyde (MDA) as compared to the wild-type (WT) plants under Pb stress. Moreover, the FaHSP17.8-CII overexpression lines, OE#3 and OE#7, exhibited 36.3% and 46.6% higher shoot Pb accumulation relative to the WT grasses. When the grasses were exposed to Pb stress, the two OE lines had higher CAT, POD and SOD activities as compared to WT. Additionally, overexpression of FaHSP17.8-CII improved the synthesis of chlorophyll and transcript abundance of FapsbC, FapsbD and FapsbE, and alleviated the photoinhibition of PSII in tall fescue under Pb stress. This study provides an initial genetic engineering strategy to improve Pb phytoremediation efficiency in tall fescue by FaHSP17.8-CII overexpression.

Overexpression of FaHSP17.8-CII improves cadmium accumulation and tolerance in tall fescue shoots by promoting chloroplast stability and photosynthetic electron transfer of PSII

Genetic improvement could play a significant role in enhancing the Cd accumulation, translocation and tolerance in plants. In this study, for the first time, we constructed transgenic tall fescue overexpressing a class II (CII) sHSP gene FaHSP17.8-CII, which enhanced Cd tolerance and the root-to-shoot Cd translocation. After exposed to 400 μM CdCl₂, two FaHSP17.8-CII overexpressing lines (OE#3 and OE#7) exhibited 30% and 40% more shoot fresh weight, respectively, relative to the wild-type (WT). Both transgenic lines showed higher tolerance to Cd, as evidenced by lower levels of electrolyte leakage and malondialdehyde compared to the WT plants under Cd stress. FaHSP17.8-CII overexpression increased shoot Cd contents 49-59% over the WT plants. The Cd translocation factor of root-to-shoot in OE grasses was 69-85% greater than WT under Cd stress. Furthermore, overexpression of FaHSP17.8-CII reduced Cd-induced damages of chloroplast ultra-structure and chlorophyll synthesis, and then improved photosystem II (PSII) function under Cd stress, which resulted in less reactive oxygen species (ROS) accumulation in OE grasses than that in WT exposed to Cd stress. The study suggests a novel FaHSP17.8-CII-PSII-ROS module to understand the mechanisms of Cd detoxification and tolerance, which provides a new strategy to improve phytoremediation efficiency in Cd-stressed grasses.

Simultaneous gene editing of three homoeoalleles in self-incompatible allohexaploid grasses

Clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) system has been widely used for precise gene editing in plants. However, simultaneous gene editing of multiple homoeoalleles remains challenging, especially in self-incompatible polyploid plants. Here, we simultaneously introduced targeted mutations in all three homoeoalleles of two genes in the self-incompatible allohexaploid tall fescue, using both CRISPR/Cas9 and LbCas12a (LbCpf1) systems. Loss-of-function mutants of FaPDS exhibited albino leaves, while knockout of FaHSP17.9 resulted in impaired heat resistance in T0 generation of tall fescue. Moreover, these mutations were inheritable. Our findings demonstrate the feasibility of generating loss-of-function mutants in T0 generation polyploid perennial grasses using CRISPR/Cas systems.

Dissecting the molecular responses of lentil to individual and combined drought and heat stresses by comparative transcriptomic analysis

Lentil cultivation could be challenged by combined heat and drought stress in semi-arid regions. We used RNA-seq approach to profile transcriptome changes of Lens culinaris exposed to individual and combined heat and drought stresses. It was determined that most of the differentially expressed genes observed in response to combined stress, could not be identified by analysis of transcriptome exposed to corresponding individual stresses. Interestingly, this study results revealed that the expression of ribosome generation and protein biosynthesis and starch degradation pathways related genes were uniquely up-regulated under the combined stress. Although multiple genes related to antioxidant activity were up-regulated in response to all stresses, variation in types and expression levels of these genes under the combined stress were higher than that of individual stresses. Using this comparative approach, for the first time, we reported up-regulation of several TF, CDPK, CYP, and antioxidant genes in response to combined stress in plants.

High temperature damage to fatty acids and carbohydrate metabolism in tall fescue by coupling deep transcriptome and metabolome analysis

High temperature damage impairs the growth of tall fescue by inhibiting secondary metabolites. Little is known about the regulation pattern of the fatty acids and carbohydrate metabolism at the whole-transcriptome level in tall fescue under high temperature stress. Here, two tall fescue genotypes, heat tolerant PI578718 and heat sensitive PI234881 were subjected to high temperature stress for 36 h. PI 578718 showed higher SPAD chloroplast value, lower EL and leaf injury than PI 234881 during the first 36 h high-temperature stress. Furthermore, by transcriptomic analysis, 121 genes were found to be induced during the second energy production phase in tall fescue exposed to high-temperature conditions, indicating that there may be one energy-sensing system in cool-season turfgrass to adapt high-temperature conditions. PI 578718 showed higher differentially expressed unigenes involved in fatty acids and carbohydrate metabolism compared with PI 234881 for 36 h heat stress. Interestingly, a metabolomic analysis using GC-MS uncovered that the sugars and sugar alcohol accounted for more than 65.06% of the total 41 metabolites content and high-temperature elevated the rate to 82.89-91.16% in PI 578718. High-temperature damage decreased the rate of fatty acid in the total 41 metabolites content and PI 578718 showed lower content than in PI 234881, which might be attributed to the down-regulated genes in fatty acid biosynthesis pathway in tall fescue. The integration of deep transcriptome and metabolome analyses provides systems-wide datasets to facilitate the identification of crucial regulation factors in cool-season turfgrass in response to high-temperature damage.

Genome-wide small RNA profiling reveals tiller development in tall fescue (Festuca arundinacea Schreb)

Background

Tall fescue (Festuca arundinacea Schreb.) is a major cool-season forage and turfgrass species. The low tiller density and size dramatically limits its turf performance and forage yield. MicroRNAs (miRNA)-genes modules play critical roles in tiller development in plants. In this study, a genome-wide small RNA profiling was carried out in two tall fescue genotypes contrasting for tillering production (‘Ch-3’, high tiller production rate and ‘Ch-5’, low tiller production rate) and two types of tissue samples at different tillering development stage (Pre-tillering, grass before tillering; Tillering, grass after tillering). ‘Ch-3’, ‘Ch-5’, Pre-tillering, and Tillering samples were analyzed using high-throughput RNA sequencing.

Results

A total of 222 million high-quality clean reads were generated and 208 miRNAs were discovered, including 148 known miRNAs belonging to 70 families and 60 novel ones. Furthermore, 18 miRNAs were involved in tall fescue tiller development process. Among them, 14 miRNAs displayed increased abundance in both Ch-3 and Tillering plants compared with that in Ch-5 and Pre-tillering plants and were positive with tillering, while another four miRNAs were negative with tiller development. Out of the three miRNAs osa-miR156a, zma-miR528a-3p and osa-miR444b.2, the rest of 15 miRNAs were newfound and associated with tiller development in plants. Based on our previous full-length transcriptome analysis in tall fescue, 28,927 potential target genes were discovered for all identified miRNAs. Most of the 212 target genes of the 18 miRNAs were dominantly enriched into “ubiquitin-mediated proteolysis”, “phagosome”, “fatty acid biosynthesis”, “oxidative phosphorylation”, and “biosynthesis of unsaturated fatty acids” KEGG pathways. In addition, bdi-miR167e-3p targets two kinase proteins EIF2AK4 and IRAK4, and osa-miR397a targets auxin response factor 5, which may be the significant miRNA-genes controllers in tillering development.

Conclusions

This is the first genome-wide miRNA profiles analysis to identify regulators involved in tiller development in cool-season turfgrass. Tillering related 18 miRNAs and their 212 target genes provide novel information for the understanding of the molecular mechanisms of miRNA-genes mediated tiller development in cool-season turfgrass.

SMRT and Illumina RNA sequencing reveal novel insights into the heat stress response and crosstalk with leaf senescence in tall fescue

BACKGROUND

As a cool-season grass species, tall fescue (Festuca arundinacea) is challenged by increasing temperatures. Heat acclimation or activation of leaf senescence, are two main strategies when tall fescue is exposed to heat stress (HS). However, lacking a genome sequence, the complexity of hexaploidy nature, and the short read of second-generation sequencing hinder a comprehensive understanding of the mechanism. This study aims to characterize the molecular mechanism of heat adaptation and heat-induced senescence at transcriptional and post-transcriptional levels.

RESULTS

Transcriptome of heat-treated (1 h and 72 h) and senescent leaves of tall fescue were generated by combining single-molecular real-time and Illumina sequencing. In total, 4076; 6917, and 11,918 differentially expressed genes (DEGs) were induced by short- and long-term heat stress (HS), and senescence, respectively. Venn and bioinformatics analyses of DEGs showed that short-term HS strongly activated heat shock proteins (Hsps) and heat shock factors (Hsfs), as well as specifically activated FK506-binding proteins (FKBPs), calcium signaling genes, glutathione S-transferase genes, photosynthesis-related genes, and phytohormone signaling genes. By contrast, long-term HS shared most of DEGs with senescence, including the up-regulated chlorophyll catabolic genes, phytohormone synthesis/degradation genes, stress-related genes, and NACs, and the down-regulated photosynthesis-related genes, FKBPs, and catalases. Subsequently, transient overexpression in tobacco showed that FaHsfA2a (up-regulated specifically by short-term HS) reduced cell membrane damages caused by HS, but FaNAC029 and FaNAM-B1 (up-regulated by long-term HS and senescence) increased the damages. Besides, alternative splicing was widely observed in HS and senescence responsive genes, including Hsps, Hsfs, and phytohormone signaling/synthesis genes.

CONCLUSIONS

The short-term HS can stimulate gene responses and improve thermotolerance, but long-term HS is a damage and may accelerate leaf senescence. These results contribute to our understanding of the molecular mechanism underlying heat adaptation and heat-induced senescence.

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.

Third-Generation Sequencing Reveals LncRNA-Regulated HSP Genes in the Populus x canadensis Moench Heat Stress Response

Long non-coding RNAs (lncRNAs) regulate plant responses to abiotic stresses. However, the short reads produced by second-generation sequencing technology make it difficult to accurately explore full-length transcripts, limiting the study of lncRNAs. In this study, we used third-generation long-read sequencing technology with the PacBio Sequel and Illumina platform to explore the role of lncRNAs in the heat stress response of Populus x canadensis Moench trees. We using 382,034,416 short reads to correct 4,297,179 long reads by resulted in 66,657 full-length transcripts, representing 33,840 genes. Then, 753 putative lncRNAs were identified, including 658 sense lncRNAs (87.38%), 41 long intervening/intergenic non-coding RNAs (lincRNAs) (5.44%), 12 antisense lncRNAs (1.59%), and 42 sense intronic lncRNAs (5.58%). Using the criteria | log₂FC| ≥ 1 and q-value < 0.05, 3,493 genes and 78 lncRNAs were differentially expressed under the heat treatment. Furthermore, 923 genes were detected as targets of 43 differently expressed lncRNAs by cis regulation. Functional annotation demonstrated that these target genes were related to unfolded protein binding, response to stress, protein folding, and response to stimulus. Lastly, we identified a lncRNA–gene interaction network consisting of four lncRNAs and six genes [Heat Shock Protein 82 (HSP82), HSP83, Disease Resistance Protein 27 (DRL27), DnaJ family protein (DNJH), and two other predicted protein-coding genes], which showed that lncRNAs could regulate HSP family genes in response to heat stress in Populus. Therefore, our third-generation sequencing has improved the description of the P. canadensis transcriptome. The potential lncRNAs and HSP family genes identified here present a genetic resource to improve our understanding of the heat-adaptation mechanisms of trees.

Molecular and physiological responses during thermal acclimation of leaf photosynthesis and respiration in rice

To further our understanding of how sustained changes in temperature affect the carbon economy of rice (Oryza sativa), hydroponically grown plants of the IR64 cultivar were developed at 30°C/25°C (day/night) before being shifted to 25/20°C or 40/35°C. Leaf messenger RNA and protein abundance, sugar and starch concentrations, and gas-exchange and elongation rates were measured on preexisting leaves (PE) already developed at 30/25°C or leaves newly developed (ND) subsequent to temperature transfer. Following a shift in growth temperature, there was a transient adjustment in metabolic gene transcript abundance of PE leaves before homoeostasis was reached within 24 hr, aligning with R_dark (leaf dark respiratory CO₂ release) and A_n (net CO₂ assimilation) changes. With longer exposure, the central respiratory protein cytochrome c oxidase (COX) declined in abundance at 40/35°C. In contrast to R_dark , A_n was maintained across the three growth temperatures in ND leaves. Soluble sugars did not differ significantly with growth temperature, and growth was fastest with extended exposure at 40/35°C. The results highlight that acclimation of photosynthesis and respiration is asynchronous in rice, with heat-acclimated plants exhibiting a striking ability to maintain net carbon gain and growth when exposed to heat-wave temperatures, even while reducing investment in energy-conserving respiratory pathways.

Transcriptome characterization and differentially expressed genes under flooding and drought stress in the biomass grasses Phalaris arundinacea and Dactylis glomerata

BACKGROUND AND AIMS

Perennial grasses are a global resource as forage, and for alternative uses in bioenergy and as raw materials for the processing industry. Marginal lands can be valuable for perennial biomass grass production, if perennial biomass grasses can cope with adverse abiotic environmental stresses such as drought and waterlogging.

METHODS

In this study, two perennial grass species, reed canary grass (Phalaris arundinacea) and cocksfoot (Dactylis glomerata) were subjected to drought and waterlogging stress to study their responses for insights to improving environmental stress tolerance. Physiological responses were recorded, reference transcriptomes established and differential gene expression investigated between control and stress conditions. We applied a robust non-parametric method, RoDEO, based on rank ordering of transcripts to investigate differential gene expression. Furthermore, we extended and validated vRoDEO for comparing samples with varying sequencing depths.

KEY RESULTS

This allowed us to identify expressed genes under drought and waterlogging whilst using only a limited number of RNA sequencing experiments. Validating the methodology, several differentially expressed candidate genes involved in the stage 3 step-wise scheme in detoxification and degradation of xenobiotics were recovered, while several novel stress-related genes classified as of unknown function were discovered.

CONCLUSIONS

Reed canary grass is a species coping particularly well with flooding conditions, but this study adds novel information on how its transcriptome reacts under drought stress. We built extensive transcriptomes for the two investigated C3 species cocksfoot and reed canary grass under both extremes of water stress to provide a clear comparison amongst the two species to broaden our horizon for comparative studies, but further confirmation of the data would be ideal to obtain a more detailed picture.

Alfalfa response to heat stress is modulated by microRNA156

Heat stress and extreme temperatures negatively affect plant development by disrupting regular cellular and biochemical functions, ultimately leading to reduced crop production. Alfalfa (Medicago sativa) is an important forage crop grown worldwide as forage for livestock feed. Limiting the effects of abiotic stress by developing alfalfa cultivars that are stress tolerant would help mitigate losses to crop production. Members of the microRNA156 (miR156) family regulate the Squamosa Promoter-Binding Protein-Like (SPL) genes that in turn impact plant growth and development by regulating downstream genes in response to various abiotic stresses. In this study, alfalfa with miR156 overexpression and SPL13 RNAi knockdown show increased tolerance to heat stress (40°C). Transgenic plants show high water potential and increased non-enzymatic antioxidant content under heat stress. Moreover, anthocyanin content and chlorophyll abundance were increased under stress. Expression of some important transcription factors and downstream genes involved in abiotic stress response were altered in miR156-overexpressing genotypes under heat. Taken together, our results demonstrate that the miR156/SPL13 network contributes to improving heat stress tolerance in alfalfa.

Effects of adult temperature on gene expression in a butterfly: identifying pathways associated with thermal acclimation

BACKGROUND

Phenotypic plasticity is a pervasive property of all organisms and considered to be of key importance for dealing with environmental variation. Plastic responses to temperature, which is one of the most important ecological factors, have received much attention over recent decades. A recurrent pattern of temperature-induced adaptive plasticity includes increased heat tolerance after exposure to warmer temperatures and increased cold tolerance after exposure to cooler temperatures. However, the mechanisms underlying these plastic responses are hitherto not well understood. Therefore, we here investigate effects of adult acclimation on gene expression in the tropical butterfly Bicyclus anynana, using an RNAseq approach.

RESULTS

We show that several antioxidant markers (e.g. peroxidase, cytochrome P450) were up-regulated at a higher temperature compared with a lower adult temperature, which might play an important role in the acclamatory responses subsequently providing increased heat tolerance. Furthermore, several metabolic pathways were up-regulated at the higher temperature, likely reflecting increased metabolic rates. In contrast, we found no evidence for a decisive role of the heat shock response.

CONCLUSIONS

Although the important role of antioxidant defence mechanisms in alleviating detrimental effects of oxidative stress is firmly established, we speculate that its potentially important role in mediating heat tolerance and survival under stress has been underestimated thus far and thus deserves more attention.

Transcriptome analysis providing novel insights for Cd-resistant tall fescue responses to Cd stress

Cadmium (Cd) is a severely toxic heavy metal and environmental pollutant. Tall fescue is a cold season turf grass which has high resistance to Cd as well as the ability to enrich it. To investigate the molecular mechanism underlying the adaptability of tall fescue to Cd stress, RNA-Seq was used to examine Cd stress responses of tall fescue at a transcriptional level. A total of 12 cDNA libraries were constructed from the total RNA of roots or leaves of tall fescue with or without Cd treatments. A total of 2594 (1768 up- and 826 down-regulated) differentially expressed genes (DEGs) were detected in the roots of Cd-stressed tall fescue compared with control roots (R_cd vs R_ck), while only 52 (29 up- and 23 down-regulated) DEGs were found in the leaves of Cd-stressed plants versus the controls (L_cd vs L_ck). The genes encoding glutathione S-transferase (GST), transporter proteins including the ABC transporter, ZRT/IRT-like protein, potassium transporter/channel, nitrate transporter, putative iron-phytosiderophore transporter, copper-transporting ATPase or transporter and multidrug and toxic compound extrusion (MATE) proteins, and numerous transcription factors were found to be significantly induced in Cd-treated roots. In addition, pathogenesis/disease-related gene mRNAs were accumulated in Cd-treated roots of tall fescue. Furthermore, the significantly enriched KEGG pathways in roots were related to 'Glutathione metabolism', 'Ribosome', 'alpha-Linolenic acid metabolism', 'Diterpenoid biosynthesis', 'Sulfur metabolism', 'Phenylpropanoid biosynthesis', 'Protein processing in endoplasmic reticulum', 'Protein export' and 'Nitrogen metabolism'. The study provides novel insights for further understanding the molecular mechanisms of tall fescue responses to Cd stress.

Comprehensive Transcriptome Profiling and Identification of Potential Genes Responsible for Salt Tolerance in Tall Fescue Leaves under Salinity Stress

Soil salinity is a serious threat to plant growth and crop productivity. Tall fescue utilization in saline areas is limited by its inferior salt tolerance. Thus, a transcriptome study is a prerequisite for future research aimed at providing deeper insights into the molecular mechanisms of tall fescue salt tolerance as well as molecular breeding. Recent advances in sequencing technology offer a platform to achieve this. Here, Illumina RNA sequencing of tall fescue leaves generated a total of 144,339 raw reads. After de novo assembly, unigenes with a total length of 129,749,938 base pairs were obtained. For functional annotations, the unigenes were aligned to various databases. Further structural analyses revealed 79,352 coding DNA sequences and 13,003 microsatellites distributed across 11,277 unigenes as well as single nucleotide polymorphisms. In total, 1862 unigenes were predicted to encode for 2120 transcription factors among which most were key salt-responsive. We determined differential gene expression and distribution per sample and most genes related to salt tolerance and photosynthesis were upregulated in 48 h vs. 24 h salt treatment. Protein interaction analysis revealed a high interaction of chaperonins and Rubisco proteins in 48 h vs. 24 h salt treatment. The gene expressions were finally validated using quantitative polymerase chain reaction (qPCR), which was coherent with sequencing results.

Transcriptomic study to understand thermal adaptation in a high temperature-tolerant strain of Pyropia haitanensis

Pyropia haitanensis, a high-yield commercial seaweed in China, is currently undergoing increasing levels of high-temperature stress due to gradual global warming. The mechanisms of plant responses to high temperature stress vary with not only plant type but also the degree and duration of high temperature. To understand the mechanism underlying thermal tolerance in P. haitanensis, gene expression and regulation in response to short- and long-term temperature stresses (SHS and LHS) was investigated by performing genome-wide high-throughput transcriptomic sequencing for a high temperature tolerant strain (HTT). A total of 14,164 differential expression genes were identified to be high temperature-responsive in at least one time point by high-temperature treatment, representing 41.10% of the total number of unigenes. The present data indicated a decrease in the photosynthetic and energy metabolic rates in HTT to reduce unnecessary energy consumption, which in turn facilitated in the rapid establishment of acclimatory homeostasis in its transcriptome during SHS. On the other hand, an increase in energy consumption and antioxidant substance activity was observed with LHS, which apparently facilitates in the development of resistance against severe oxidative stress. Meanwhile, ubiquitin-mediated proteolysis, brassinosteroids, and heat shock proteins also play a vital role in HTT. The effects of SHS and LHS on the mechanism of HTT to resist heat stress were relatively different. The findings may facilitate further studies on gene discovery and the molecular mechanisms underlying high-temperature tolerance in P. haitanensis, as well as allow improvement of breeding schemes for high temperature-tolerant macroalgae that can resist global warming.

Transcriptomic profiling of tall fescue in response to heat stress and improved thermotolerance by melatonin and 24-epibrassinolide

Background

Tall fescue is a widely used cool season turfgrass and relatively sensitive to high temperature. Chemical compounds like melatonin (MT) and 24-epibrassinolide (EBL) have been reported to improve plant heat stress tolerance effectively.

Results

In this study, we reported that MT and EBL pretreated tall fescue seedlings showed decreased reactive oxygen species (ROS), electrolyte leakage (EL) and malondialdehide (MDA), but increased chlorophyll (Chl), total protein and antioxidant enzyme activities under heat stress condition, resulting in improved plant growth. Transcriptomic profiling analysis showed that 4311 and 8395 unigenes were significantly changed after 2 h and 12 h of heat treatments, respectively. Among them, genes involved in heat stress responses, DNA, RNA and protein degradation, redox, energy metabolisms, and hormone metabolism pathways were highly induced after heat stress. Genes including FaHSFA3, FaAWPM and FaCYTC2 were significantly upregulated by both MT and EBL treatments, indicating that these genes might function as the putative target genes of MT and EBL.

Conclusions

These findings indicated that heat stress caused extensively transcriptomic reprogramming of tall fescue and exogenous application of MT and EBL effectively improved thermotolerance in tall fescue.

Electronic supplementary material

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

Transcriptomic analysis reveals unique molecular factors for lipid hydrolysis, secondary cell-walls and oxidative protection associated with thermotolerance in perennial grass

Background

Heat stress is the primary abiotic stress limiting growth of cool-season grass species. The objective of this study was to determine molecular factors and metabolic pathways associated with superior heat tolerance in thermal bentgrass (Agrostis scabra) by comparative analysis of transcriptomic profiles with its co-generic heat-sensitive species creeping bentgrass (A. stolonifera).

Results

Transcriptomic profiling by RNA-seq in both heat-sensitive A. stolonifera (cv. ‘Penncross’) and heat-tolerant A. scabra exposed to heat stress found 1393 (675 up- and 718 down-regulated) and 1508 (777 up- and 731 down-regulated) differentially-expressed genes, respectively. The superior heat tolerance in A. scabra was associated with more up-regulation of genes in oxidative protection, proline biosynthesis, lipid hydrolysis, hemicellulose and lignin biosynthesis, compared to heat-sensitive A. stolonifera. Several transcriptional factors (TFs), such as high mobility group B protein 7 (HMGB7), dehydration-responsive element-binding factor 1a (DREB1a), multiprotein-bridging factor 1c (MBF1c), CCCH-domain containing protein 47 (CCCH47), were also found to be up-regulated in A. scabra under heat stress.

Conclusions

The unique TFs and genes identified in thermal A. scabra could be potential candidate genes for genetic modification of cultivated grass species for improving heat tolerance, and the associated pathways could contribute to the transcriptional regulation for superior heat tolerance in bentgrass species.

Electronic supplementary material

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

Exogenous Calcium Enhances the Photosystem II Photochemistry Response in Salt Stressed Tall Fescue

Calcium enhances turfgrass response to salt stress. However, little is known about PSII photochemical changes when exogenous calcium was applied in salinity-stressed turfgrass. Here, we probe into the rearrangements of PSII electron transport and endogenous ion accumulation in tall fescue (Festuca arundinacea Schreber) treated with exogenous calcium under salt stress. Three-month-old seedlings of genotype "TF133" were subjected to the control (CK), salinity (S), salinity + calcium nitrate (SC), and salinity + ethylene glycol tetraacetic acid (SE). Calcium nitrate and ethylene glycol tetraacetic acid was used as exogenous calcium donor and calcium chelating agent respectively. At the end of a 5-day duration treatment, samples in SC regime had better photochemistry performance on several parameters than salinity only. Such as the Area (equal to the plastoquinone pool size), N (number of [Formula: see text] redox turnovers until Fm is reached), ψE0, or δRo (Efficiencdy/probability with which a PSII trapped electron is transferred from QA to QB or PSI acceptors), ABS/RC (Absorbed photon flux per RC). All the above suggested that calcium enhanced the electron transfer of PSII (especially beyond [Formula: see text]) and prevented reaction centers from inactivation in salt-stressed tall fescue. Furthermore, both grass shoot and root tissues generally accumulated more C, N, Ca2+, and K+ in the SC regime than S regime. Interrelated analysis indicated that ψE0, δRo, ABS/RC, C, and N content in shoots was highly correlated to each other and significantly positively related to Ca2+ and K+ content in roots. Besides, high salt increased ATP6E and CAMK2 transcription level in shoot at 1 and 5 day, respectively while exogenous calcium relieved it. In root, CAMK2 level was reduced by Salinity at 5 day and exogenous calcium recovered it. These observations involved in electron transport capacity and ion accumulation assist in understanding better the protective role of exogenous calcium in tall fescue under salt stress.

Ascorbic Acid Alleviates Damage from Heat Stress in the Photosystem II of Tall Fescue in Both the Photochemical and Thermal Phases

L-Ascorbate (Asc) plays important roles in plant development, hormone signaling, the cell cycle and cellular redox system, etc. The higher content of Asc in plant chloroplasts indicates its important role in the photosystem. The objective of this study was to study the roles of Asc in tall fescue leaves against heat stress. After a heat stress treatment, we observed a lower value of the maximum quantum yield for primary photochemistry (φPo), which reflects the inhibited activity of the photochemical phase of photosystem II (PSII). Moreover, we observed a higher value of efficiency of electron transfer from QB to photosystem I acceptors (δR0), which reflects elevated activity of the thermal phase of the photosystem of the tall fescue. The addition of Asc facilitate the behavior of the photochemical phase of the PSII by lowering the ROS content as well as that of the alternative electron donor to provide electron to the tyrosine residue of the D1 protein. Additionally, exogenous Asc reduces the activity of the thermal phase of the photosystem, which could contribute to the limitation of energy input into the photosystem in tall fescue against heat stress. Synthesis of the Asc increased under heat stress treatment. However, under heat stress this regulation does not occur at the transcription level and requires further study.

Genome-wide identification of heat stress-responsive small RNAs in tall fescue (Festuca arundinacea) by high-throughput sequencing

MicroRNAs (miRNAs) play vital roles in the adaptive response of plants to various abiotic and biotic stresses. Tall fescue (Festuca arundinacea Schreb.) is a major cool-season forage and turf grass species which is severely influenced by heat stress. To unravel possible heat stress-responsive miRNAs, high-throughput sequencing was employed for heat-tolerant PI578718 and heat-sensitive PI234881 genotypes growing in presence and absence of heat stress (40°C for 36h). By searching against the miRBase database, among 1421 reference monocotyledon miRNAs, more than 850 were identified in all samples. Among these miRNAs, 1.46% and 2.29% were differentially expressed in PI234881 and PI578718 under heat stress, respectively, and most of them were down-regulated. In addition, a total of 170 novel miRNAs belonging to 145 miRNA families were identified. Furthermore, putative targets of differentially expressed miRNAs were predicted. The regulation of selected miRNAs by heat stress was revalidated through quantitative reverse transcription PCR (qRT-PCR) analysis. Most of these miRNAs shared similar expression patterns; however, some showed distinct expression patterns under heat stress, with their putative targets displaying different transcription levels. This is the first genome-wide miRNA identification in tall fescue. miRNAs specific to PI578718, or those that exhibited differential expression profiles between the two genotypes under high temperature, were probably associated with the variation in thermotolerance of tall fescue. The differentially expressed miRNAs between these two tall fescue genotypes and their putative targeted genes will provide essential information for further study on miRNAs mediating heat response and facilitate to improve turf grass breeding.

Transcriptome profilings of two tall fescue (Festuca arundinacea) cultivars in response to lead (Pb) stress

BACKGROUND

Lead (Pb) is one of the most toxic heavy metal environmental pollutants. Tall fescue is an important cold season turf grass which can tolerate and accumulate substantial amount of Pb. To estimate genes related to Pb response and the molecular mechanism associated with Pb tolerance and accumulation, we analyzed the transcriptome of tall fescue in response to Pb treatment.

RESULTS

RNA-sequencing was performed in two tall fescue cultivars, Pb tolerant Silverado and Pb sensitive AST7001. A total of 810,146 assembled unique transcripts representing 25,415 unigenes were obtained from the tall fescue leaves. Among the panel, 3,696 differentially expressed genes (DEGs) were detected between the Pb treated (1000 mg/L) and untreated samples. Gene ontology (GO) and pathway enrichment analysis demonstrated that the DEGs were mainly implicated in energy metabolism, metabolism of terpenoids and polyketides, and carbohydrate metabolism related pathways. The expression patterns of 16 randomly selected genes were in consistent with that from the Solexa analysis using quantitative reverse-transcription PCR. In addition, compared to the common transcriptional response to Pb stress in both cultivars, the regulation of numerous genes including those involved in zeatin biosynthesis, limonene and pinene degradation, phagosome was exclusive to one cultivar.

CONCLUSIONS

The tall fescue assembled transcriptome provided substantial molecular resources for further genomics analysis of turfgrass in response to heavy metal stress. The significant expression difference of specific unigenes may account for Pb tolerance or accumulation in two different tall fescue cultivars. This study provided new insights for the investigation of the molecular basis of Pb tolerance and accumulation in tall fescue as well as other related turf grass species.

Transcriptome response of Lolium arundinaceum to its fungal endophyte Epichloë coenophiala

Tall fescue (Lolium arundinaceum) is one of the primary forage and turf grasses in temperate regions of the world. A number of favourable characteristics of tall fescue are enhanced by its seed-transmissible fungal symbiont (endophyte) Epichloë coenophiala. Our approach was to assemble the tall fescue transcriptome, then identify differentially expressed genes (DEGs) for endophyte-symbiotic (E+) vs endophyte-free (E-) clones in leaf blades, pseudostems, crowns and roots. RNA-seq reads were used to construct a tall fescue reference transcriptome and compare gene expression profiles. Over all tissues examined, 478 DEGs were identified between the E+ and E- clones for at least one tissue (more than two-fold; P < 0.0001, 238 E+ > E- and 240 E- > E+), although no genes were differentially expressed in all four tissues. Gene ontology (GO) terms, GO:0010200 (response to chitin), GO:0002679 (respiratory burst during defence response) and GO:0035556 (intracellular signal transduction) were significantly overrepresented among 25 E- > E+ DEGs in leaf blade, and a number of other DEGs were associated with defence and abiotic response. In particular, endophyte effects on various WRKY transcription factors may have implications for symbiotic stability, endophyte distribution in the plant, or defence against pathogens.

Transcriptional Profiling and Identification of Heat-Responsive Genes in Perennial Ryegrass by RNA-Sequencing

Perennial ryegrass (Lolium perenne) is one of the most widely used forage and turf grasses in the world due to its desirable agronomic qualities. However, as a cool-season perennial grass species, high temperature is a major factor limiting its performance in warmer and transition regions. In this study, a de novo transcriptome was generated using a cDNA library constructed from perennial ryegrass leaves subjected to short-term heat stress treatment. Then the expression profiling and identification of perennial ryegrass heat response genes by digital gene expression analyses was performed. The goal of this work was to produce expression profiles of high temperature stress responsive genes in perennial ryegrass leaves and further identify the potentially important candidate genes with altered levels of transcript, such as those genes involved in transcriptional regulation, antioxidant responses, plant hormones and signal transduction, and cellular metabolism. The de novo assembly of perennial ryegrass transcriptome in this study obtained more total and annotated unigenes compared to previously published ones. Many DEGs identified were genes that are known to respond to heat stress in plants, including HSFs, HSPs, and antioxidant related genes. In the meanwhile, we also identified four gene candidates mainly involved in C₄ carbon fixation, and one TOR gene. Their exact roles in plant heat stress response need to dissect further. This study would be important by providing the gene resources for improving heat stress tolerance in both perennial ryegrass and other cool-season perennial grass plants.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.

Global Transcriptome Profiles of 'Meyer' Zoysiagrass in Response to Cold Stress

A long green period is essential for a turfgrass species with high ornamental value and a wide area of use. Zoysiagrasses (Zoysia spp. Willd.) are perennial turfgrass species popular in tropical, subtropical and temperate zones, possessing many properties necessary to be economically useful turfgrass. They do not have a long green period because of cold sensitivity. A main focus in zoysiagrass research is to develop cold tolerant cultivars. Understanding the cold response in zoysiagrass is a fundamental area of research. In the present study, ‘Meyer’ zoysiagrass (Zoysia japonica), a widely cultivated variety in the genus, is used. We employed RNA-Seq to investigate genome-wide gene expression profiles in leaves under cold stress (4°C). Using the Illumina sequencing platform, we obtained approximately 206 million high-quality paired-end reads from three libraries (0 h, 2 h, and 72 h cold treatment at 4°C). After de novo assembly and quantitative assessment, 46,412 unigenes were generated with an average length of 998 bp and an N50 of 1,522 bp. A total of 25,644 (55.2%) unigenes were annotated by alignment with public protein databases including NR, SwissProt, KEGG and KOG. Differentially expressed genes (DEGs) were investigated using the RPKM method. A total of 756 DEGs were identified between 0h and 2h-cold treatment, with 522 up-regulated and 234 down-regulated; and 5327 DEGs were identified between 0h and 72h-cold treatment, with 2453 up-regulated and 2874 down-regulated. The expression profile of 15 DEGs selected randomly was confirmed with qRT-PCR. The results suggest that cold stress can induce desiccation and oxidative stress, inhibit photosynthesis and substance transport. In response to the stress, genes involved in proline synthesis, in starch hydrolysis, in methionine and ascorbic acid metabolism, in SOD activity, and in DREBs response pathway were up-regulated. GA metabolism, ABA and JA stimulus response were affected under cold exposure. This is the first transcriptome sequencing of Z. japonica, providing a large set of sequence data as well as gene expression profiles under cold stress. It will improve our current understanding of the cold response of zoysiagrass and be beneficial in breeding research.

Heat Shock Factor Genes of Tall Fescue and Perennial Ryegrass in Response to Temperature Stress by RNA-Seq Analysis

Heat shock factors (Hsfs) are important regulators of stress-response in plants. However, our understanding of Hsf genes and their responses to temperature stresses in two Pooideae cool-season grasses, Festuca arundinacea, and Lolium perenne, is limited. Here we conducted comparative transcriptome analyses of plant leaves exposed to heat or cold stress for 10 h. Approximately, 30% and 25% of the genes expressed in the two species showed significant changes under heat and cold stress, respectively, including subsets of Hsfs and their target genes. We uncovered 74 Hsfs in F. arundinacea and 52 Hsfs in L. perenne, and categorized these genes into three subfamilies, HsfA, HsfB, and HsfC based on protein sequence homology to known Hsf members in model organisms. The Hsfs showed a strong response to heat and/or cold stress. The expression of HsfAs was elevated under heat stress, especially in class HsfA2, which exhibited the most dramatic responses. HsfBs were upregulated by the both temperature conditions, and HsfCs mainly showed an increase in expression under cold stress. The target genes of Hsfs, such as heat shock protein (HSP), ascorbate peroxidase (APX), inositol-3-phosphate synthase (IPS), and galactinol synthase (GOLS1), showed strong and unique responses to different stressors. We comprehensively detected Hsfs and their target genes in F. arundinacea and L. perenne, providing a foundation for future gene function studies and genetic engineering to improve stress tolerance in grasses and other crops.