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

Transcriptomic and metabolomic profiling provide insight into the role of sugars and hormones in leaf senescence of Pinellia ternata

KEY MESSAGE

The interaction network and pathway map uncover the potential crosstalk between sugar and hormone metabolisms as a possible reason for leaf senescence in P. ternata. Pinellia ternata, an environmentally sensitive medicinal plant, undergoes leaf senescence twice a year, affecting its development and yield. Understanding the potential mechanism that delays leaf senescence could theoretically decrease yield losses. In this study, a typical senescent population model was constructed, and an integrated analysis of transcriptomic and metabolomic profiles of P. ternata was conducted using two early leaf senescence populations and two stay-green populations. The result showed that two key gene modules were associated with leaf senescence which were mainly enriched in sugar and hormone signaling pathways, respectively. A network constructed by unigenes and metabolisms related to the obtained two pathways revealed that several compounds such as D-arabitol and 2MeScZR have a higher significance ranking. In addition, a total of 130 hub genes in this network were categorized into 3 classes based on connectivity. Among them, 34 hub genes were further analyzed through a pathway map, the potential crosstalk between sugar and hormone metabolisms might be an underlying reason of leaf senescence in P. ternata. These findings address the knowledge gap regarding leaf senescence in P. ternata, providing candidate germplasms for molecular breeding and laying theoretical basis for the realization of finely regulated cultivation in future.

Transcriptomic Analysis of Three Differentially Senescing Maize (Zea mays L.) Inbred Lines upon Heat Stress

Maize, one of the world's major food crops, is facing the challenge of rising temperature. Leaf senescence is the most significant phenotypic change of maize under heat stress at the seedling stage, but the underlying molecular mechanism is still unknown. Here, we screened for three inbred lines (PH4CV, B73, and SH19B) that showed differentially senescing phenotypes under heat stress. Among them, PH4CV showed no obviously senescing phenotype under heat stress, while SH19B demonstrated a severely senescing phenotype, with B73 being between the two extremes. Subsequently, transcriptome sequencing showed that differentially expressed genes (DEGs) were generally enriched in response to heat stress, reactive oxygen species (ROS), and photosynthesis in the three inbred lines under heat treatment. Notably, ATP synthesis and oxidative phosphorylation pathway genes were only significantly enriched in SH19B. Then, the expression differences of oxidative phosphorylation pathways, antioxidant enzymes, and senescence-related genes in response to heat stress were analyzed in the three inbred lines. In addition, we demonstrated that silencing ZmbHLH51 by virus-induced gene silencing (VIGS) inhibits the heat-stress-induced senescence of maize leaves. This study helps to further elucidate the molecular mechanisms of heat-stress-induced leaf senescence at the seedling stage of maize.

Comprehensive transcriptome analysis reveals heat-responsive genes in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis) using RNA sequencing

Flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee, 2n=20, AA) is a vegetable species in southern parts of China that faces high temperatures in the summer and winter seasons. While heat stress adversely impacts plant productivity and survival, the underlying molecular and biochemical causes are poorly understood. This study investigated the gene expression profiles of heat-sensitive (HS) '3T-6' and heat-tolerant (HT) 'Youlu-501' varieties of flowering Chinese cabbage in response to heat stress using RNA sequencing. Among the 37,958 genes expressed in leaves, 20,680 were differentially expressed genes (DEGs) at 1, 6, and 12 h, with 1,078 simultaneously expressed at all time points in both varieties. Hierarchical clustering analysis identified three clusters comprising 1,958, 556, and 591 down-regulated, up-regulated, and up- and/or down-regulated DEGs (3205 DEGs; 8.44%), which were significantly enriched in MAPK signaling, plant-pathogen interactions, plant hormone signal transduction, and brassinosteroid biosynthesis pathways and involved in stimulus, stress, growth, reproductive, and defense responses. Transcription factors, including MYB (12), NAC (13), WRKY (11), ERF (31), HSF (17), bHLH (16), and regulatory proteins such as PAL, CYP450, and photosystem II, played an essential role as effectors of homeostasis, kinases/phosphatases, and photosynthesis. Among 3205 DEGs, many previously reported genes underlying heat stress were also identified, e.g., BraWRKY25, BraHSP70, BraHSPB27, BraCYP71A23, BraPYL9, and BraA05g032350.3C. The genome-wide comparison of HS and HT provides a solid foundation for understanding the molecular mechanisms of heat tolerance in flowering Chinese cabbage.

Fatty Acid Biosynthesis Pathways Are Downregulated during Stigma Development and Are Critical during Self-Incompatible Responses in Ornamental Kale

In Brassicaceae, the papillary cells of the stigma are the primary site of the self-incompatibility (SI) responses. SI preserves the genetic diversity by selectively rejecting irrelevant or incompatible pollen, thus promoting cross fertilization and species fitness. Mechanisms that regulate SI responses in Brassica have been studied mainly on the mature stigma that often undermines how stigma papillary cells attain the state of SI during development. To understand this, we integrated PacBio SMRT-seq with Illumina RNA-seq to construct a de novo full-length transcriptomic database for different stages of stigma development in ornamental kale. A total of 48,800 non-redundant transcripts, 31,269 novel transcripts, 24,015 genes, 13,390 alternative splicing, 22,389 simple sequence repeats, 21,816 complete ORF sequences, and 4591 lncRNAs were identified and analyzed using PacBio SMRT-seq. The Illumina RNA-seq revealed 15,712 differentially expressed genes (DEGs) and 8619 transcription factors. The KEGG enrichment analysis of 4038 DEGs in the “incompatibility” group revealed that the flavonoid and fatty acid biosynthesis pathways were significantly enriched. The cluster and qRT-PCR analysis indicated that 11 and 14 candidate genes for the flavonoid and fatty acid biosynthesis pathways have the lowest expression levels at stigma maturation, respectively. To understand the physiological relevance of the downregulation of fatty acid biosynthesis pathways, we performed inhibitor feeding assays on the mature stigma. The compatible pollination response was drastically reduced when mature stigmas were pre-treated with a fatty acid synthase inhibitor. This finding suggested that fatty acid accumulation in the stigmas may be essential for compatible pollination and its downregulation during maturity must have evolved as a support module to discourage the mounting of self-incompatible pollen.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Comprehensive Biochemical, Physiological, and Transcriptomic Analyses Provide Insights Into Floral Bud Dormancy in Rhododendron delavayi Franch

Due to a scarcity of relevant data, the ornamental woody flower Rhododendron delavayi Franch. is examined in the current study for its low temperature-induced floral bud dormancy (late October-end December) aspect. This study used transcriptome data profiling and co-expression network analyses to identify the interplay between endogenous hormones and bud dormancy phases such as pre-dormancy, para-dormancy, endo-dormancy, eco-dormancy, and dormancy release. The biochemical and physiological assays revealed the significance of the abundance of phytohormones (abscisic acid, auxin, zeatin, and gibberellins), carbohydrate metabolism, oxidative species, and proteins (soluble proteins, proline, and malondialdehyde) in the regulatory mechanism of floral bud dormancy. The transcriptome sequencing generated 65,531 transcripts, out of which 504, 514, 307, and 240 expressed transcripts were mapped uniquely to pre-, para-, endo-, and eco-phases of dormancy, showing their roles in the stimulation of dormancy. The transcripts related to LEA29, PGM, SAUR family, RPL9e, ATRX, FLOWERING LOCUS T, SERK1, ABFs, ASR2, and GID1 were identified as potential structural genes involved in floral bud dormancy. The transcription factors, including Zinc fingers, CAD, MADS-box family, MYB, and MYC2, revealed their potential regulatory roles concerning floral bud dormancy. The gene co-expression analysis highlighted essential hub genes involved in cold stress adaptations encoding proteins, viz, SERPIN, HMA, PMEI, LEA_2, TRX, PSBT, and AMAT. We exposed the connection among low temperature-induced dormancy in floral buds, differentially expressed genes, and hub genes via strict screening steps to escalate the confidence in selected genes as being truly putative in the pathways regulating bud dormancy mechanism. The identified candidate genes may prove worthy of further in-depth studies on molecular mechanisms involved in floral bud dormancy of Rhododendron species.

Single-Molecule Real-Time and Illumina Sequencing to Analyze Transcriptional Regulation of Flavonoid Synthesis in Blueberry

Blueberries (Vaccinium corymbosum) contain large amounts of flavonoids, which play important roles in the plant's ability to resist stress and can also have beneficial effects on human health when the fruits are eaten. However, the molecular mechanisms that regulate flavonoid synthesis in blueberries are still unclear. In this study, we combined two different transcriptome sequencing platforms, single-molecule real-time (SMRT) and Illumina sequencing, to elucidate the flavonoid synthetic pathways in blueberries. We analyzed transcript quantity, length, and the number of annotated genes. We mined genes associated with flavonoid synthesis (such as anthocyanins, flavonols, and proanthocyanidins) and employed fluorescence quantitative PCR to analyze the expression of these genes and their correlation with flavonoid synthesis. We discovered one R2R3 MYB transcription factor from the sequencing library, VcMYB1, that can positively regulate anthocyanin synthesis in blueberries. VcMYB1 is mainly expressed in colored (mature) fruits. Experiments showed that overexpression and transient expression of VcMYB1 promoted anthocyanin synthesis in Arabidopsis, tobacco (Nicotiana benthamiana) plants and green blueberry fruits. Yeast one-hybrid (Y1H) assay, electrophoretic mobility shift assay, and transient expression experiments showed that VcMYB1 binds to the MYB binding site on the promoter of the structural gene for anthocyanin synthesis, VcMYB1 to positively regulate the transcription of VcDFR, thereby promoting anthocyanin synthesis. We also performed an in-depth investigation of transcriptional regulation of anthocyanin synthesis. This study provides background information and data for studying the synthetic pathways of flavonoids and other secondary metabolites in blueberries.

Chitosan (CTS) Alleviates Heat-Induced Leaf Senescence in Creeping Bentgrass by Regulating Chlorophyll Metabolism, Antioxidant Defense, and the Heat Shock Pathway

Chitosan (CTS) is a deacetylated derivative of chitin that is involved in adaptive response to abiotic stresses. However, the regulatory role of CTS in heat tolerance is still not fully understood in plants, especially in grass species. The aim of this study was to investigate whether the CTS could reduce heat-induced senescence and damage to creeping bentgrass associated with alterations in antioxidant defense, chlorophyll (Chl) metabolism, and the heat shock pathway. Plants were pretreated exogenously with or without CTS (0.1 g L-1) before being exposed to normal (23/18 °C) or high-temperature (38/33 °C) conditions for 15 days. Heat stress induced detrimental effects, including declines in leaf relative water content and photochemical efficiency, but significantly increased reactive oxygen species (ROS) accumulation, membrane lipid peroxidation, and Chl loss in leaves. The exogenous application of CTS significantly alleviated heat-induced damage in creeping bentgrass leaves by ameliorating water balance, ROS scavenging, the maintenance of Chl metabolism, and photosynthesis. Compared to untreated plants under heat stress, CTS-treated creeping bentgrass exhibited a significantly higher transcription level of genes involved in Chl biosynthesis (AsPBGD and AsCHLH), as well as a lower expression level of Chl degradation-related gene (AsPPH) and senescence-associated genes (AsSAG12, AsSAG39, Asl20, and Ash36), thus reducing leaf senescence and enhancing photosynthetic performance under heat stress. In addition, the foliar application of CTS significantly improved antioxidant enzyme activities (SOD, CAT, POD, and APX), thereby effectively reducing heat-induced oxidative damage. Furthermore, heat tolerance regulated by the CTS in creeping bentgrass was also associated with the heat shock pathway, since AsHSFA-6a and AsHSP82 were significantly up-regulated by the CTS during heat stress. The potential mechanisms of CTS-regulated thermotolerance associated with other metabolic pathways still need to be further studied in grass species.

Gene co-expression network analysis of the heat-responsive core transcriptome identifies hub genes in Brassica rapa

Gene co-expression network analysis of the heat-responsive core transcriptome in two contrasting Brassica rapa accessions reveals the main metabolic pathways, key modules and hub genes, are involved in long-term heat stress. Brassica rapa is a widely cultivated and economically important vegetable in Asia. High temperature is a common stress that severely impacts leaf head formation in B. rapa, resulting in reduced quality and production. The purpose of this study was thus to identify candidate heat tolerance genes by comparative transcriptome analysis of two contrasting B. rapa accessions in response to long-term heat stress. Two B. rapa accessions, '268' and '334', which showed significant differences in heat tolerance, were used for RNA sequencing analysis. We identified a total of 11,055 and 8921 differentially expressed genes (DEGs) in '268' and '334', respectively. Functional enrichment analyses of all of the identified DEGs, together with the genes identified from weighted gene co-expression network analyses (WGCNA), revealed that the autophagy pathway, glutathione metabolism, and ribosome biogenesis in eukaryotes were significantly up-regulated, whereas photosynthesis was down-regulated, in the heat resistance of B. rapa '268'. Furthermore, when B. rapa '334' was subjected to long-term high-temperature stress, heat stress caused significant changes in the expression of certain functional genes linked to protein processing in the endoplasmic reticulum and plant hormone signal transduction pathways. Autophagy-related genes might have been induced by persistent heat stress and remained high during recovery. Several hub genes like HSP17.6, HSP17.6B, HSP70-8, CLPB1, PAP1, PYR1, ADC2, and GSTF11 were discussed in this study, which may be potential candidates for further analyses of the response to long-term heat stress. These results should help elucidate the molecular mechanisms of heat stress adaptation in B. rapa.

Mechanisms of specific systemic response in wheat plants under different locally acting heat stimuli

Mechanisms of the specific systemic response of plant to different adverse factors are poorly understood. We studied the mechanisms acting in wheat (Triticum aestivum L.) under the action of local burn and gradual heating. Both stimuli induce a variation potential (VP) propagation and a biphasic (fast and long-term phases) photosynthetic response in non-stimulated zones of plant with stimulus-specific parameters of the latter: the fast phase or long-term phase predominance in responses induced by burn or heating, respectively. The burn-induced VP and photosynthetic response attenuate with distance, while the heating-induced VP and photosynthetic response were of more stable amplitude in distant part of the stimulated plant. VP propagation in both cases induced apoplast alkalization with dynamics well corresponding to such of VP and of the fast phase of photosynthetic response. Gradual heating induced a significant rise in jasmonate production along with a decrease in stomatal conductance with characteristic times well corresponding to the long-term phase of the photosynthetic response. We suppose that the VP-induced pH shift is responsible for in the induction of the fast phase, while jasmonate production for the long-term phase of the photosynthetic response. The revealed differences in the systemic response to stressors studied, apparently, reflect two distinct plant adaptation strategies to fast and slow-growing stimuli. The immediate response in the tissue nearest to the damage zone is the most important under a fast-growing stimulus. The fundamentally different situation is under a slowly-growing stimulus which provokes long-term changes in the plant that ensure the preparation of the whole organism for impending environmental changes.