Transcriptomic network analyses of leaf dehydration responses identify highly connected ABA and ethylene signaling hubs in three grapevine species differing in drought tolerance

Gene Co-expression Network Analysis of the Comparative Transcriptome Identifies Hub Genes Associated With Resistance to Aspergillus flavus L. in Cultivated Peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.), a cosmopolitan oil crop, is susceptible to a variety of pathogens, especially Aspergillus flavus L., which not only vastly reduce the quality of peanut products but also seriously threaten food safety for the contamination of aflatoxin. However, the key genes related to resistance to Aspergillus flavus L. in peanuts remain unclear. This study identifies hub genes positively associated with resistance to A. flavus in two genotypes by comparative transcriptome and weighted gene co-expression network analysis (WGCNA) method. Compared with susceptible genotype (Zhonghua 12, S), the rapid response to A. flavus and quick preparation for the translation of resistance-related genes in the resistant genotype (J-11, R) may be the drivers of its high resistance. WGCNA analysis revealed that 18 genes encoding pathogenesis-related proteins (PR10), 1-aminocyclopropane-1-carboxylate oxidase (ACO1), MAPK kinase, serine/threonine kinase (STK), pattern recognition receptors (PRRs), cytochrome P450, SNARE protein SYP121, pectinesterase, phosphatidylinositol transfer protein, and pentatricopeptide repeat (PPR) protein play major and active roles in peanut resistance to A. flavus. Collectively, this study provides new insight into resistance to A. flavus by employing WGCNA, and the identification of hub resistance-responsive genes may contribute to the development of resistant cultivars by molecular-assisted breeding.

Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families

The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.

Grapevine aquaporins: Diversity, cellular functions, and ecophysiological perspectives

High-scored premium wines are typically produced under moderate drought stress, suggesting that the water status of grapevine is crucial for wine quality. Aquaporins greatly influence the plant water status by facilitating water diffusion across the plasma membrane in a tightly regulated manner. They adjust the hydraulic conductance of the plasma membrane rapidly and reversibly, which is essential in specific physiological events, including adaptation to soil water scarcity. The comprehension of the sophisticated plant-water relations at the molecular level are thus important to optimize agricultural practices or to assist plant breeding programs. This review explores the recent progresses in understanding the water transport in grapevine at the cellular level through aquaporins and its regulation. Important aspects, including aquaporin structure, diversity, cellular localization, transport properties, and regulation at the cellular and whole plant level are addressed. An ecophysiological perspective about the roles of grapevine aquaporins in plant response to drought stress is also provided.

Comparative transcriptome analyses between cultivated and wild grapes reveal conservation of expressed genes but extensive rewiring of co-expression networks

The transcriptomes of wild and cultivated grapes consists of similar expressed genes but distinct wiring of co-expressed modules associated with environmental conditions. Grapevine is an important fruit crop worldwide, with high economic value and widespread distribution. Commercial production is based on Vitis vinifera, and, to a lesser extent, on hybrids with American grapes, such as V. labrusca. Wild grape relatives are important sources of resistance against biotic and abiotic factors; however, their global gene expression patterns remain poorly characterized. We associated genome-wide transcript profiling to phenotypic analyses to investigate the responses of cultivated and wild vines to vineyard conditions. The expressed genes in the Vitis reference transcriptome are largely shared by wild grapes, V. labrusca hybrids and vinifera cultivars. In contrast, significant differential regulation between wild and vinifera genotypes represents 80% of gene expression variation, regardless of the environment. In wild grapes, genes associated to regulatory processes are downregulated, whereas those involved in metabolic pathways are upregulated, in comparison to vinifera. Photosynthesis-related ontologies are overrepresented in the induced genes, in agreement with higher contents of chlorophyll in wild grapes. Co-regulated gene network analyses provide evidence of more complex transcriptome organization in vinifera. In wild grapes, genes involved in signaling pathways of stress-related hormones are overrepresented in modules associated with the environment. Consensus network analyses revealed high preservation within co-regulated gene modules between cultivated and wild grapes, but divergent relationships among the expression clusters. In conclusion, the distinct phenotypes of wild and cultivated grapes are underlain by differences in gene expression, but also by distinct higher-order organization of the transcriptome and contrasting association of co-expressed gene clusters with the environment.

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.

Comparative Phenotypic and Transcriptomic Analysis Reveals Key Responses of Upland Cotton to Salinity Stress During Postgermination

To understand the molecular mechanisms of salinity tolerance during seed germination and post-germination stages, this study characterized phenotypic and transcriptome responses of two cotton cultivars during salinity stress. The two cultivars were salt-tolerant (ST) LMY37 and salt-sensitive (SS) ZM12, with the former exhibiting higher germination rate, growth, and primary-root fresh weight under salinity stress. Transcriptomic comparison revealed that up-regulation of differentially expressed genes (DEGs) was the main characteristic of transcriptional regulation in ST, while SS DEGs were mainly down-regulated. GO and KEGG analyses uncovered both common and specific responses in ST and SS. Common processes, such as reactive oxygen species (ROS) metabolism and cell wall biosynthesis, may be general responses to salinity in cotton. In contrast, DEGs involved in MAPK-signaling pathway activated by ROS, carotenoid biosynthesis pathway and cysteine and methionine metabolism pathway [producing the precursors of stress hormone abscisic acid (ABA) and ethylene (ET), respectively] as well as stress tolerance related transcription factor genes, showed significant expression differences between ST and SS. These differences might be the molecular basis leading to contrasting salinity tolerance. Silencing of GhERF12, an ethylene response factor gene, caused higher salinity sensitivity and increased ROS accumulation after salinity stress. In addition, peroxidase (POD) and superoxide dismutase (SOD) activity obviously declined after silencing GhERF12. These results suggest that GhERF12 is involved in salinity tolerance during early development. This study provides a novel and comprehensive perspective to understand key mechanisms of salinity tolerance and explores candidate genes that may be useful in developing stress-tolerant crops through biotechnology.

Imperative role of sugar signaling and transport during drought stress responses in plants

Cellular sugar status is essentially maintained during normal growth conditions but is impacted negatively during various environmental perturbations. Drought presents one such unfavorable environmental cue that hampers the photosynthetic fixation of carbon into sugars and affects their transport by lowering the cellular osmotic potential. The transport of cellular sugar is facilitated by a specific set of proteins known as sugar transporters. These transporter proteins are the key determinant of influx/ efflux of various sugars and their metabolite intermediates that support the plant growth and developmental process. Abiotic stress and especially drought stress-mediated injury results in reprogramming of sugar distribution across the cellular and subcellular compartments. Here, we have reviewed the imperative role of sugar accumulation, signaling, and transport under typical and atypical stressful environments. We have discussed the physiological effects of drought on sugar accumulation and transport through different transporter proteins involved in monosaccharide and disaccharide sugar transport. Further, we have illustrated sugar-mediated signaling and regulation of sugar transporter proteins along with the overall crosstalk of this signaling with the phytohormone module of abiotic stress response under osmotic stress. Overall, the present review highlights the critical role of sugar transport, distribution and signaling in plants under drought stress conditions.

Molecular Tools for Adapting Viticulture to Climate Change

Adaptation of viticulture to climate change includes exploration of new geographical areas, new training systems, new management practices, or new varieties, both for rootstocks and scions. Molecular tools can be defined as molecular approaches used to study DNAs, RNAs, and proteins in all living organisms. We present here the current knowledge about molecular tools and their potential usefulness in three aspects of grapevine adaptation to the ongoing climate change. (i) Molecular tools for understanding grapevine response to environmental stresses. A fine description of the regulation of gene expression is a powerful tool to understand the physiological mechanisms set up by the grapevine to respond to abiotic stress such as high temperatures or drought. The current knowledge on gene expression is continuously evolving with increasing evidence of the role of alternative splicing, small RNAs, long non-coding RNAs, DNA methylation, or chromatin activity. (ii) Genetics and genomics of grapevine stress tolerance. The description of the grapevine genome is more and more precise. The genetic variations among genotypes are now revealed with new technologies with the sequencing of very long DNA molecules. High throughput technologies for DNA sequencing also allow now the genetic characterization at the same time of hundreds of genotypes for thousands of points in the genome, which provides unprecedented datasets for genotype-phenotype associations studies. We review the current knowledge on the genetic determinism of traits for the adaptation to climate change. We focus on quantitative trait loci and molecular markers available for developmental stages, tolerance to water stress/water use efficiency, sugar content, acidity, and secondary metabolism of the berries. (iii) Controlling the genome and its expression to allow breeding of better-adapted genotypes. High-density DNA genotyping can be used to select genotypes with specific interesting alleles but genomic selection is also a powerful method able to take into account the genetic information along the whole genome to predict a phenotype. Modern technologies are also able to generate mutations that are possibly interesting for generating new phenotypes but the most promising one is the direct editing of the genome at a precise location.

Dynamic changes of transcriptome of fifth-instar spodoptera litura larvae in response to insecticide

Spodoptera litura is a major insect with a cosmopolitan distribution and strong resistance to multiple insecticides. Determining the molecular basis and key candidate genes of the insecticide resistance of S. litura may help in managing this insect. In this study, fifth-instar S. litura larvae were subjected to transcriptome analysis at 6, 12, 24, 48, and 72 h after feeding on an LC20 dose of avermectin. The result showed that genes responding to avermectin changed dynamically with different gene counts and resistance mechanisms at the fifth instar based on a metabolic pathway map. These responses included degrading the insecticide by a series of P450 and glutathione-S-transferase enzymes starting at the 12 h time point, with subsequent increases in the number of genes involved and shifts to TOLL and immune deficiency (IMD) pathways at 48 h after feeding the insecticide. Weighted correlation network analysis (WGCNA) determined a co-expression module related to the avermectin response at 12 and 24 h (r = 0.403, p = 0.0371; r = 0.436, p = 0.023), in which a hub gene (LOC111358940) related to metalloproteinase activity was identified. In addition, Analysis of the genes in the co-expression module further revealed that eight genes encoding UDP-glucuronosyltransferases were directly associated with insecticide response in S. litura. These results provide better understanding of the avermectin response mechanism of S. litura and may be useful in developing improved control strategies for this species.

SUPPLEMENTARY INFORMATION

The online version of this article (10.1007/s13205-021-02651-9) contains supplementary material, which is available to authorized users.

Transcriptional analysis reveals sodium nitroprusside affects alfalfa in response to PEG-induced osmotic stress at germination stage

Drought is one of the most common environmental factors that affect alfalfa germination and development. Nitric oxide (NO) could mediate stress tolerance in plants. The goal of this study was to determine exogenous NO donor-mediated drought adaption molecular mechanisms during the alfalfa germination stage. In this study, physiological and transcriptome analyses were performed on 7 days of the growth period seedlings by sodium nitroprusside (SNP) and polyethylene glycol (PEG) treatment. The results showed that SNP supplementation alleviated malondialdehyde accumulation, increased levels of proline and soluble sugars, and enhanced antioxidant enzyme activity under osmotic stress conditions. RNA-Seq experiments identified 5828 genes exhibiting differential expression in seedlings treated with PEG, SNP, or SNP+PEG relative to seedlings treated with distilled water. Of these DEGs, 3235 were upregulated, and 2593 were downregulated relative to the controls. Fifteen DEGs were amplified by qRT-PCR to verify the changes in expression determined by RNA-Seq, revealing that PIF3, glnA, PLCG1, and RP-S11e exhibited enhanced expression under the SNP+PEG treatment. SNP was found to modulate redox homeostasis-related genes such as GSTs, SOD2, GPX, and RBOH, and triggered calcium signaling transduction. It also induced some key genes relating to the abscisic acid, ethylene, and auxin signaling transduction in response to PEG stress. Conversely, genes associated with secondary metabolite biosynthesis and the metabolism of starch and sucrose during osmotic stress were downregulated by SNP. These results provide new insights into SNP-mediated drought adaption mechanisms at transcriptome-wide in alfalfa and reveal key drought tolerance pathways in this species.

VviERF6Ls: an expanded clade in Vitis responds transcriptionally to abiotic and biotic stresses and berry development

Background

VviERF6Ls are an uncharacterized gene clade in Vitis with only distant Arabidopsis orthologs. Preliminary data indicated these transcription factors may play a role in berry development and extreme abiotic stress responses. To better understand this highly duplicated, conserved clade, additional members of the clade were identified in four Vitis genotypes. A meta-data analysis was performed on publicly available microarray and RNA-Seq data (confirmed and expanded with RT-qPCR), and Vitis VviERF6L1 overexpression lines were established and characterized with phenotyping and RNA-Seq.

Results

A total of 18 PN40024 VviERF6Ls were identified; additional VviERF6Ls were identified in Cabernet Sauvignon, Chardonnay, and Carménère. The amino acid sequences of VviERF6Ls were found to be highly conserved. VviERF6L transcripts were detected in numerous plant organs and were differentially expressed in response to numerous abiotic stresses including water deficit, salinity, and cold as well as biotic stresses such as red blotch virus, N. parvum, and E. necator. VviERF6Ls were differentially expressed across stages of berry development, peaking in the pre-veraison/veraison stage and retaining conserved expression patterns across different vineyards, years, and Vitis cultivars. Co-expression network analysis identified a scarecrow-like transcription factor and a calmodulin-like gene with highly similar expression profiles to the VviERF6L clade. Overexpression of VviERF6L1 in a Seyval Blanc background did not result in detectable morphological phenotypes. Genes differentially expressed in response to VviERF6L1 overexpression were associated with abiotic and biotic stress responses.

Conclusions

VviERF6Ls represent a large and distinct clade of ERF transcription factors in grapevine. The high conservation of protein sequence between these 18 transcription factors may indicate these genes originate from a duplication event in Vitis. Despite high sequence similarity and similar expression patterns, VviERF6Ls demonstrate unique levels of expression supported by similar but heterogeneous promoter sequences. VviERF6L gene expression differed between Vitis species, cultivars and organs including roots, leaves and berries. These genes respond to berry development and abiotic and biotic stresses. VviERF6L1 overexpression in Vitis vinifera results in differential expression of genes related to phytohormone and immune system signaling. Further investigation of this interesting gene family is warranted.

Identification of key genes and modules in response to Cadmium stress in different rice varieties and stem nodes by weighted gene co-expression network analysis

Soil cadmium (Cd) pollution threatens food safety. This study aimed to identify genes related to Cd accumulation in rice. Low- (Shennong 315, short for S315) and high- (Shendao 47, short for S47) Cd-accumulative rice cultivars were incubated with CdCl₂·2.5H₂O. RNA-seq and weighted gene co-expression network analysis (WGCNA) were performed to identify the modules and genes associated with Cd-accumulative traits of rice. After Cd stress treatment, the Cd content in various tissues of S315 was significantly higher than that of S47. In the stem nodes, the Cd distribution results of the two varieties indicated that the unelongated nodes near the root (short for node A) had a stronger ability to block Cd transfer upwards than the panicle node (short for node B). Cd stress induced huge changes in gene expression profiles. After analyzing the differentially expressed genes (DEGs) in significantly correlated WGCNA modules, we found that genes related to heavy metal transportation had higher expression levels in node A than that in node B, such as Copper transporter 6 (OS04G0415600), Zinc transporter 10 (OS06G0566300), and some heavy-metal associated proteins (OS11G0147500, OS03G0861400, and OS10G0506100). In the comparison results between S315 and S47, the expression of chitinase (OS03G0679700 and OS06G0726200) was increased by Cd treatment in S315. In addition, OsHSPs (OS05G0460000, OS08G0500700), OsHSFC2A (OS02G0232000), and OsDJA5 (OS03G0787300) were found differentially expressed after Cd treatment in S315, but changed less in S47. In summary, different rice varieties have different processes and intensities in response to Cd stress. The node A might function as the key tissue for blocking Cd upward transport into the panicle via vigorous processes, including of heavy metal transportation, response to stress, and cell wall.

Can metabolic tightening and expansion of co-expression network play a role in stress response and tolerance?

Plants respond and adapt to changes in their environment by employing a wide variety of genetic, molecular, and biochemical mechanisms. When so doing, they trigger large-scale rearrangements at the metabolic and transcriptional levels. The dynamics and patterns of these rearrangements and how they govern a stress response is not clear. In this opinion, we discuss a plant's response to stress from the perspective of the metabolic gene co-expression network and its rearrangement upon stress. As a case study, we use publicly available expression data of Arabidopsis thaliana plants exposed to heat and drought stress to evaluate and compare the co-expression networks of metabolic genes. The analysis highlights that stress conditions can lead to metabolic tightening and expansion of the co-expression network. We argue that this rearrangement could play a role in a plant's response to stress and thus may be an additional tool to assess and understand stress tolerance/sensitivity. Additional studies are needed to evaluate the metabolic network in response to multiple stresses at various intensities and across different genetic backgrounds (e.g., intra- and inter-species, sensitive and tolerant eco/genotypes).

Drought tolerance of the grapevine, Vitis champinii cv. Ramsey, is associated with higher photosynthesis and greater transcriptomic responsiveness of abscisic acid biosynthesis and signaling

BACKGROUND

Grapevine is an economically important crop for which yield and berry quality is strongly affected by climate change. Large variations in drought tolerance exist across Vitis species. Some of these species are used as rootstock to enhance abiotic and biotic stress tolerance. In this study, we investigated the physiological and transcriptomic responses to water deficit of four different genotypes that differ in drought tolerance: Ramsey (Vitis champinii), Riparia Gloire (Vitis riparia), Cabernet Sauvignon (Vitis vinifera), and SC2 (Vitis vinifera x Vitis girdiana).

RESULTS

Ramsey was particularly more drought tolerant than the other three genotypes. Ramsey maintained a higher stomatal conductance and photosynthesis at equivalent levels of moderate water deficit. We identified specific and common transcriptomic responses shared among the four different Vitis species using RNA sequencing analysis. A weighted gene co-expression analysis identified a water deficit core gene set with the ABA biosynthesis and signaling genes, NCED3, RD29B and ABI1 as potential hub genes. The transcript abundance of many abscisic acid metabolism and signaling genes was strongly increased by water deficit along with genes associated with lipid metabolism, galactinol synthases and MIP family proteins. This response occurred at smaller water deficits in Ramsey and with higher transcript abundance than the other genotypes. A number of aquaporin genes displayed differential and unique responses to water deficit in Ramsey leaves. Genes involved in cysteine biosynthesis and metabolism were constitutively higher in the roots of Ramsey; thus, linking the gene expression of a known factor that influences ABA biosynthesis to this genotype's increased NCED3 transcript abundance.

CONCLUSION

The drought tolerant Ramsey maintained higher photosynthesis at equivalent water deficit than the three other grapevine genotypes. Ramsey was more responsive to water deficit; its transcriptome responded at smaller water deficits, whereas the other genotypes did not respond until more severe water deficits were reached. There was a common core gene network responding to water deficit for all genotypes that included ABA metabolism and signaling. The gene clusters and sub-networks identified in this work represent interesting gene lists to explore and to better understand drought tolerance molecular mechanisms.

A sense of place: transcriptomics identifies environmental signatures in Cabernet Sauvignon berry skins in the late stages of ripening

BACKGROUND

Grape berry ripening is influenced by climate, the main component of the "terroir" of a place. Light and temperature are major factors in the vineyard that affect berry development and fruit metabolite composition.

RESULTS

To better understand the effect of "place" on transcript abundance during the late stages of berry ripening, Cabernet Sauvignon berries grown in Bordeaux and Reno were compared at similar sugar levels (19 to 26 °Brix (total soluble solids)). Day temperatures were warmer and night temperatures were cooler in Reno. °Brix was lower in Bordeaux berries compared to Reno at maturity levels considered optimum for harvest. RNA-Seq analysis identified 5528 differentially expressed genes between Bordeaux and Reno grape skins at 22°Brix. Weighted Gene Coexpression Network Analysis for all expressed transcripts for all four °Brix levels measured indicated that the majority (75%) of transcript expression differed significantly between the two locations. Top gene ontology categories for the common transcript sets were translation, photosynthesis, DNA metabolism and catabolism. Top gene ontology categories for the differentially expressed genes at 22°Brix involved response to stimulus, biosynthesis and response to stress. Some differentially expressed genes encoded terpene synthases, cell wall enzymes, kinases, transporters, transcription factors and photoreceptors. Most circadian clock genes had higher transcript abundance in Bordeaux. Bordeaux berries had higher transcript abundance with differentially expressed genes associated with seed dormancy, light, auxin, ethylene signaling, powdery mildew infection, phenylpropanoid, carotenoid and terpenoid metabolism, whereas Reno berries were enriched with differentially expressed genes involved in water deprivation, cold response, ABA signaling and iron homeostasis.

CONCLUSIONS

Transcript abundance profiles in the berry skins at maturity were highly dynamic. RNA-Seq analysis identified a smaller (25% of total) common core set of ripening genes that appear not to depend on rootstock, vineyard management, plant age, soil and climatic conditions. Much of the gene expression differed between the two locations and could be associated with multiple differences in environmental conditions that may have affected the berries in the two locations; some of these genes may be potentially controlled in different ways by the vinegrower to adjust final berry composition and reach a desired result.

Weighted Gene Co-expression Network Analysis (WGCNA) Reveals the Hub Role of Protein Ubiquitination in the Acquisition of Desiccation Tolerance in Boea hygrometrica

Boea hygrometrica can survive extreme drought conditions and has been used as a model to study desiccation tolerance. A genome-wide transcriptome analysis of B. hygrometrica showed that the plant can survive rapid air-drying after experiencing a slow soil-drying acclimation phase. In addition, a weighted gene co-expression network analysis was used to study the transcriptomic datasets. A network comprising 22 modules was constructed, and seven modules were found to be significantly related to desiccation response using an enrichment analysis. Protein ubiquitination was observed to be a common process linked to hub genes in all the seven modules. Ubiquitin-modified proteins with diversified functions were identified using immunoprecipitation coupled with mass spectrometry. The lowest level of ubiquitination was noted at the full soil drying priming stage, which coincided the accumulation of dehydration-responsive gene BhLEA2. The highly conserved RY motif (CATGCA) was identified from the promoters of ubiquitin-related genes that were downregulated in the desiccated samples. An in silico gene expression analysis showed that the negative regulation of ubiquitin-related genes is potentially mediated via a B3 domain-containing transcription repressor VAL1. This study suggests that priming may involve the transcriptional regulation of several major processes, and the transcriptional regulation of genes in protein ubiquitination may play a hub role to deliver acclimation signals to posttranslational level in the acquisition of desiccation tolerance in B. hygrometrica.

Integration of Metabolite Profiling and Transcriptome Analysis Reveals Genes Related to Volatile Terpenoid Metabolism in Finger Citron (C. medica var. sarcodactylis)

Finger citron (Citrus medica var. sarcodactylis) is a popular ornamental tree and an important source of essential oils rich in terpenoids, but the mechanisms behind volatile formation are poorly understood. We investigated gene expression changes combined with volatile profiling of ten samples from three developing organs: flower, leaf, and fruit. A total of 62 volatiles were identified with limonene and γ-terpinene being the most abundant ones. Six volatiles were identified using partial least squares discriminant analysis (PLS-DA) that could be used as markers for distinguishing finger citron from other citrus species. RNA-Seq revealed 1,611,966,118 high quality clean reads that were assembled into 32,579 unigenes. From these a total of 58 terpene synthase (TPS) gene family members were identified and the spatial and temporal distribution of their transcripts was measured in developing organs. Transcript levels of transcription factor genes AP2/ERF (251), bHLH (169), bZIP (76), MYB (155), NAC (184), and WRKY (66) during finger citron development were also analyzed. From extracted subnetworks of three modules constructed by weighted gene co-expression network analysis (WGCNA), thirteen TPS genes and fifteen transcription factors were suggested to be related to volatile terpenoid formation. These results provide a framework for future investigations into the identification and regulatory network of terpenoids in finger citron.

Distinct leaf transcriptomic response of water deficient Eucalyptus grandis submitted to potassium and sodium fertilization

While potassium fertilization increases growth yield in Brazilian eucalyptus plantations, it could also increase water requirements, making trees more vulnerable to drought. Sodium fertilization, which has been shown to promote eucalyptus growth compared to K-deficient trees, could partially mitigate this adverse effect of potassium. However, little is known about the influence of K and Na fertilization on the tree metabolic response to water deficit. The aim of the present study was thus to analyze the transcriptome of leaves sampled from Eucalyptus grandis trees subjected to 37% rainfall reduction, and fertilized with potassium (K), sodium (Na), compared to control trees (C). The multifactorial experiment was set up in a field with a throughfall exclusion system. Transcriptomic analysis was performed on leaves from two-year-old trees, and data analyzed using multifactorial statistical analysis and weighted gene co-expression network analysis (WGCNA). Significant sets of genes were seen to respond to rainfall reduction, in interaction with K or Na fertilization, or to fertilization only (regardless of the water supply regime). The genes were involved in stress signaling, primary and secondary metabolism, secondary cell wall formation and photosynthetic activity. Our focus on key genes related to cation transporters and aquaporins highlighted specific regulation of ion homeostasis, and plant adjustment to water deficit. While water availability significantly affects the transcriptomic response of eucalyptus species, this study points out that the transcriptomic response is highly dependent on the fertilization regime. Our study is based on the first large-scale field trial in a tropical region, specifically designed to study the interaction between water availability and nutrition in eucalyptus. To our knowledge, this is the first global transcriptomic analysis to compare the influence of K and Na fertilization on tree adaptive traits in water deficit conditions.

Expression profiling of CTR1-like and EIN2-like genes in buds and leaves of Populus tremula, and in vitro study of the interaction between their polypeptides

In Arabidopsis, the serine/threonine protein kinase Constitutive Triple Response 1 (CTR1) and Ethylene Insensitive 2 polypeptide (EIN2) functions are key negative and positive components, respectively, in the ethylene signalling route. Here, we report on an in silico study of members of the CTR1-like and EIN2-like polypeptide families from poplars. The expression of CTR1-like and EIN2-like genes such as Ptre-CTR1, Ptre-CTR3 and Ptre-EIN2a was studied in Populus tremula buds and leaves in response to dehydration, various light conditions and under senescence. In buds under dehydration, the maximal fold-change of the Ptre-CTR1, Ptre-CTR3 and Ptre-EIN2a expression level recorded almost identical values. This suggests that maintenance of a constant ratio between the transcript levels of genes encoding positive and negative ethylene signalling components is required under stress. The expression of the studied genes was 1.4-to 3-fold higher in response to darkness, but 4.5- to 51.2-fold and 21.6- to 51.2-fold higher under the early and moderate leaf senescence, respectively. It is worth noting that the senescence-related Ptre-EIN2a and Ptre-CTR3a expression profiles were very similar. Using in vitro assays, we demonstrated the ability of the catalytic domain of Ptre-CTR1 to phosphorylate the Ptre-EIN2a-like polypeptide, which is similar to that in Arabidopsis. The target substrate, the Ptre-CEND2a polypeptide (C-terminal part of Ptre-EIN2a), was only phosphorylated by the protein kinase Ptre-CTR1 and not by Ptre-CTR3. Moreover, the addition of Ptre-CTR3 polypeptides (-CTR3a or -CTR3b forms) to the reaction mixture had an inhibitory effect on Ptre-CTR1 auto- and trans-phosphorylation. In contrast to Ptre-CTR1, Ptre-CTR3 may act as a positive regulator in ethylene signalling in poplar; however, this hypothesis requires in vivo confirmation. Thus, the ethylene signalling route in poplar seems to be under the control of certain additional mechanisms which have not been reported in Arabidopsis.

Dynamic Transcriptome Changes Related to Oil Accumulation in Developing Soybean Seeds

Soybean is one of the most important oil crops in the world. Revealing the molecular basis and exploring key candidate genes for seed oil synthesis has great significance for soybean improvement. In this study, we found that oil accumulation rates and gene expression levels changed dynamically during soybean seed development. The expression levels of genes in metabolic pathways such as carbon fixation, photosynthesis, glycolysis, and fatty acid biosynthesis were significantly up-regulated during the rapid accumulation of oil in developing soybean seeds. Through weighted correlation network analysis, we identified six co-expression modules associated with soybean seed oil content and the pink module was the most positively correlated (r = 0.83, p = 7 × 10-4) network. Through the integration of differential expression and co-expression analysis, we predicted 124 candidate genes potentially affecting soybean seed oil content, including seven genes in lipid metabolism pathway, two genes involved in glycolysis, one gene in sucrose metabolism, and 12 genes belonged to transcription factors as well as other categories. Among these, three genes (GmABI3b, GmNFYA and GmFAD2-1B) have been shown to control oil and fatty acid content in soybean seeds, and other newly identified candidate genes would broaden our knowledge to understand the molecular basis for oil accumulation in soybean seeds.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Transcriptomic response is more sensitive to water deficit in shoots than roots of Vitis riparia (Michx.)

BACKGROUND

Drought is an important constraint on grapevine sustainability. Vitis riparia, widely used in rootstock and scion breeding, has been studied in isolated leaf drying response studies; however, it is essential to identify key root and shoot water deficit signaling traits in intact plants. This information will aid improved scion and rootstock selection and management practices in grapevine. RNAseq data were generated from V. riparia roots and shoots under water deficit and well-watered conditions to determine root signaling and shoot responses to water deficit.

RESULTS

Shoot elongation, photosynthetic rate, and stomatal conductance were significantly reduced in water deficit (WD) treated than in well-watered grapevines. RNAseq analysis indicated greater transcriptional differences in shoots than in roots under WD, with 6925 and 1395 genes differentially expressed, respectively (q-value < 0.05). There were 50 and 25 VitisNet pathways significantly enriched in WD relative to well-watered treatments in grapevine shoots and roots, respectively. The ABA biosynthesis genes beta-carotene hydroxylase, zeaxanthin epoxidase, and 9-cis-epoxycarotenoid dioxygenases were up-regulated in WD root and WD shoot. A positive enrichment of ABA biosynthesis genes and signaling pathways in WD grapevine roots indicated enhanced root signaling to the shoot. An increased frequency of differentially expressed reactive oxygen species scavenging (ROS) genes were found in the WD shoot. Analyses of hormone signaling genes indicated a strong ABA, auxin, and ethylene network and an ABA, cytokinin, and circadian rhythm network in both WD shoot and WD root.

CONCLUSIONS

This work supports previous findings in detached leaf studies suggesting ABA-responsive binding factor 2 (ABF2) is a central regulator in ABA signaling in the WD shoot. Likewise, ABF2 may have a key role in V. riparia WD shoot and WD root. A role for ABF3 was indicated only in WD root. WD shoot and WD root hormone expression analysis identified strong ABA, auxin, ethylene, cytokinin, and circadian rhythm signaling networks. These results present the first ABA, cytokinin, and circadian rhythm signaling network in roots under water deficit. These networks point to organ specific regulators that should be explored to further define the communication network from soil to shoot.

Swift metabolite changes and leaf shedding are milestones in the acclimation process of grapevine under prolonged water stress

BACKGROUND

Grape leaves provide the biochemical substrates for berry development. Thus, understanding the regulation of grapevine leaf metabolism can aid in discerning processes fundamental to fruit development and berry quality. Here, the temporal alterations in leaf metabolism in Merlot grapevine grown under sufficient irrigation and water deficit were monitored from veraison until harvest.

RESULTS

The vines mediated water stress gradually and involving multiple strategies: osmotic adjustment, transcript-metabolite alteration and leaf shedding. Initially stomatal conductance and leaf water potential showed a steep decrease together with the induction of stress related metabolism, e.g. up-regulation of proline and GABA metabolism and stress related sugars, and the down-regulation of developmental processes. Later, progressive soil drying was associated with an incremental contribution of Ca²⁺ and sucrose to the osmotic adjustment concomitant with the initiation of leaf shedding. Last, towards harvest under progressive stress conditions following leaf shedding, incremental changes in leaf water potential were measured, while the magnitude of perturbation in leaf metabolism lessened.

CONCLUSIONS

The data present evidence that over time grapevine acclimation to water stress diversifies in temporal responses encompassing the alteration of central metabolism and gene expression, osmotic adjustments and reduction in leaf area. Together these processes mitigate leaf water stress and aid in maintaining the berry-ripening program.

Dehydration-induced proteomic landscape of mitochondria in chickpea reveals large-scale coordination of key biological processes

Mitochondria play crucial roles in regulating multiple biological processes particularly electron transfer and energy metabolism in eukaryotic cells. Exposure to water-deficit or dehydration may affect mitochondrial function, and dehydration response may dictate cell fate decisions. iTRAQ-based quantitative proteome of a winter legume, chickpea, demonstrated the central metabolic alterations in mitochondria, presumably involved in dehydration adaptation. Three-week-old chickpea seedlings were subjected to progressive dehydration and the magnitude of dehydration-induced compensatory physiological responses was monitored in terms of physicochemical characteristics and mitochondrial architecture. The proteomics analysis led to the identification of 40 dehydration-responsive proteins whose expressions were significantly modulated by dehydration. The differentially expressed proteins were implicated in different metabolic processes, with obvious functional tendencies toward purine-thiamine metabolic network, pathways of carbon fixation and oxidative phosphorylation. The linearity of dehydration-induced proteome alteration was examined with transcript abundance of randomly selected candidates under multivariate stress conditions. The differentially regulated proteins were validated through sequence analysis. An extensive sequence based localization prediction revealed >62.5% proteins to be mitochondrial resident by, at least, one prediction algorithm. The results altogether provide intriguing insights into the dehydration-responsive metabolic pathways and useful clues to identify crucial proteins linked to stress tolerance. BIOLOGICAL SIGNIFICANCE: Investigation on plant mitochondrial proteome is of significance because it would allow a better understanding of mitochondrial function in plant adaptation to stress. Mitochondria are the unique organelles, which play a crucial role in energy metabolism and cellular homeostasis, particularly when exposed to stress conditions. Chickpea is one of the cultivated winter legumes, which enriches soil nitrogen and has very low water footprint and thus contributes to fortification of sustainable agriculture. We therefore examined the dehydration-responsive mitochondrial proteome landscape of chickpea and queried whether molecular interplay of mitochondrial proteins modulate dehydration tolerance. A total of 40 dehydration-induced mitochondrial proteins were identified, predicted to be involved in key metabolic processes. Our future efforts would focus on understanding both posttranslational modification and processing for comprehensive characterization of mitochondrial protein function. This approach will facilitate mining of more biomarkers linked to the tolerance trait and contribute to crop adaptation to climate change.

Patterns of Drought Response of 38 WRKY Transcription Factors of Zanthoxylum bungeanum Maxim

The WRKY family of transcription factors (TFs) includes a number of transcription-specific groupings that play important roles in plant growth and development and in plant responses to various stresses. To screen for WRKY transcription factors associated with drought stress in Zanthoxylum bungeanum, a total of 38 ZbWRKY were identified and these were then classified and identified with Arabidopsis WRKY. Using bioinformatics analyses based on the structural characteristics of the conservative domain, 38 WRKY transcription factors were identified and categorized into three groups: Groups I, II, and III. Of these, Group II can be divided into four subgroups: subgroups IIb, IIc, IId, and IIe. No ZbWRKY members of subgroup IIa were found in the sequencing data. In addition, 38 ZbWRKY were identified by real-time PCR to determine the behavior of this family of genes under drought stress. Twelve ZbWRKY transcription factors were found to be significantly upregulated under drought stress and these were identified by relative quantification. As predicted by the STRING website, the results show that the WRKYs are involved in four signaling pathways—the jasmonic acid (JA), the salicylic acid (SA), the mitogen-activated protein kinase (MAPK), and the ethylene signaling pathways. ZbWRKY33 is the most intense transcription factor in response to drought stress. We predict that WRKY33 binds directly to the ethylene synthesis precursor gene ACS6, to promote ethylene synthesis. Ethylene then binds to the ethylene activator release signal to activate a series of downstream genes for cold stress and osmotic responses. The roles of ZbWRKY transcription factors in drought stress rely on a regulatory network center on the JA signaling pathway.

Pathways and Network Based Analysis of Candidate Genes to Reveal Cross-Talk and Specificity in the Sorghum (Sorghum bicolor (L.) Moench) Responses to Drought and It's Co-occurring Stresses

Drought alone or in combination with other stresses forms the major crop production constraint worldwide. Sorghum, one of the most important cereal crops is affected by drought alone or in combination with co-occurring stresses; notwithstanding, sorghum has evolved adaptive responses to combined stresses. Furthermore, an impressive number of sorghum genes have been investigated for drought tolerance. However, the molecular mechanism underling drought response remains poorly understood. We employed a systems biology approach to elucidate regulatory and broad functional features of these genes. Their interaction network would provide insight into understanding the molecular mechanisms of drought tolerance and underpinning signal pathways. Functional analysis was undertaken to determine significantly enriched genesets for pathways involved in drought tolerance. Analysis of distinct pathway cross-talk network was performed and drought-specific subnetwork was extracted. Investigation of various data sources such as gene expression, regulatory pathways, sorghumCyc, sorghum protein-protein interaction (PPI) and Gene Ontology (GO) revealed 14 major drought stress related hub genes (DSRhub genes). Significantly enriched genesets have shown association with various biological processes underlying drought-related responses. Key metabolic pathways were significantly enriched in the drought-related genes. Systematic analysis of pathways cross-talk and gene interaction network revealed major cross-talk pathway modules associated with drought tolerance. Further investigation of the major DSRhub genes revealed distinct regulatory genes such as ZEP, NCED, AAO, and MCSU and CYP707A1. These were involved in the regulation of ABA biosynthesis and signal transduction. Other protein families, namely, aldehyde and alcohol dehydrogenases, mitogene activated protein kinases (MAPKs), and Ribulose-1,5-biphosphate carboxylase (RuBisCO) were shown to be involved in the drought-related responses. This shows a diversity of complex functional features in sorghum to respond to various abiotic stresses. Finally, we constructed a drought-specific subnetwork, characterized by unique candidate genes that were associated with DSRhub genes. According to our knowledge, this is the first in sorghum drought investigation that introduces pathway and network-based candidate gene approach for analysis of drought tolerance. We provide novel information about pathways cross-talk and signaling networks used in further systems level analysis for understanding the molecular mechanism behind drought tolerance and can, therefore, be adapted to other model and non-model crops.

Overexpression of 9-cis-Epoxycarotenoid Dioxygenase Cisgene in Grapevine Increases Drought Tolerance and Results in Pleiotropic Effects

9-cis-epoxycarotenoid dioxygenase (NCED) is a key enzyme involved in the biosynthesis of abscisic acid (ABA), which is associated with drought tolerance in plants. An osmotic-inducible VaNCED1 gene was isolated from a drought-resistant cultivar of Vitis amurensis and constitutively overexpressed in a drought-sensitive cultivar of Vitis vinifera. Transgenic plants showed significantly improved drought tolerance, including a higher growth rate and better drought resistant under drought conditions, compared to those of wild-type (WT) plants. After water was withheld for 50 days, the upper leaves of transgenic plants remained green, whereas most leaves of WT plants turned yellow and fell. Besides the increase in ABA content, overexpression of VaNCED1 induced the production of jasmonic acid (JA) and accumulation of JA biosynthesis-related genes, including allene oxide cyclase (AOC) and 12-oxophytodienoate reductase (OPR3). Moreover, transgenic plants possessed advantageous physiological indices, including lower leaf stomatal density, lower photosynthesis rate, and lower accumulation of proline and superoxide dismutase (SOD), compared to those of WT plants, indicating increased resistance to drought stress. Quantitative real time polymerase chain reaction (RT-qPCR) analysis revealed that overexpression of VaNCED1 enhanced the expression of drought-responsive genes, such as ABA-responsive element1 (ABRE1), ABRE binding factors 2 (ABF2), plasma membrane intrinsic proteins 2 (PIP2), C-repeat/DRE-Binding Factor 4 (VvCBF4) and ABA-insensitive 5 (ABI5). Although the development of transgenic plants was delayed by 4 months than WT plants, because of seed dormancy and abnormal seedlings, the surviving transgenic plants provided a solid method for protection of woody plants from drought stress.

Involvement of an ABI-like protein and a Ca2+-ATPase in drought tolerance as revealed by transcript profiling of a sweetpotato somatic hybrid and its parents Ipomoea batatas (L.) Lam. and I. triloba L

Previously, we obtained the sweetpotato somatic hybrid KT1 from a cross between sweetpotato (Ipomoea batatas (L.) Lam.) cv. Kokei No. 14 and its drought-tolerant wild relative I. triloba L. KT1 not only inherited the thick storage root characteristic of Kokei No. 14 but also the drought-tolerance trait of I. triloba L. The aim of this study was to explore the molecular mechanism of the drought tolerance of KT1. Four-week-old in vitro-grown plants of KT1, Kokei No. 14, and I. triloba L. were subjected to a simulated drought stress treatment (30% PEG6000) for 0, 6, 12 and 24 h. Total RNA was extracted from samples at each time point, and then used for transcriptome sequencing. The gene transcript profiles of KT1 and its parents were compared to identify differentially expressed genes, and drought-related modules were screened by a weighted gene co-expression network analysis. The functions of ABI-like protein and Ca2+-ATPase, two proteins screened from the cyan and light yellow modules, were analyzed in terms of their potential roles in drought tolerance in KT1 and its parents. These analyses of the drought responses of KT1 and its somatic donors at the transcriptional level provide new annotations for the molecular mechanism of drought tolerance in the somatic hybrid KT1 and its parents.

Divergence in the transcriptional landscape between low temperature and freeze shock in cultivated grapevine (Vitis vinifera)

Low-temperature stresses limit the sustainability and productivity of grapevines when early spring frosts damage young grapevine leaves. Spring conditions often expose grapevines to low, but not damaging, chilling temperatures and these temperatures have been shown to increase freeze resistance in other model systems. In this study, we examined whole-transcriptome gene expression patterns of young leaf tissue from cuttings of five different grapevine cultivars, exposed to chill and freeze shock, in order to understand the underlying transcriptional landscape associated with cold stress response. No visible damage was observed when grapevine leaves were exposed to chilling temperatures while freeze temperatures resulted in variable damage in all cultivars. Significant differences in gene expression were observed between warm control conditions and all types of cold stress. Exposure to chill stress (4 °C) versus freezing stress (-3 °C) resulted in very different patterns of gene expression and enriched pathway responses. Genes from the ethylene signaling, ABA signaling, the AP2/ERF, WRKY, and NAC transcription factor families, and starch/sucrose/galactose pathways were among the most commonly observed to be differentially regulated. Preconditioning leaves to chill temperatures prior to freezing temperatures resulted in slight buffering of gene expression responses, suggesting that differences between chill and freeze shock perception complicates identification of candidate genes for cold resistance in grapevine. Overall, the transcriptional landscape contrasts observed between low temperature and freezing stresses demonstrate very different activation of candidate pathways impacting grapevine cold response.

A Chemical Genetics Strategy that Identifies Small Molecules which Induce the Triple Response in Arabidopsis

To explore small molecules with ethylene-like biological activity, we conducted a triple response-based assay system for chemical library screening. Among 9600 compounds, we found N-[(1,3,5-trimethyl-1H-pyrazol-4-yl)methyl]-N-methyl-2-naphthalenesulfonamide (EH-1) displayed promising biological activity on inducing a triple response in Arabidopsis seedlings. Chemical synthesis and structure-activity relationship (SAR) analysis of EH-1 analogues with different substitution patterns on the phenyl ring structure of the sulfonamide group indicated that 3,4-dichloro-N-methyl-N-(1,3,5-trimethyl-1H-pyrazol-4-yl-methyl) benzenesulfonamide (8) exhibits the most potent biological activity. To determine the mechanism of action, we conducted RNA sequencing (RNA-Seq) analysis of the effect of EH-1 and 1-aminocyclopropane-1-carboxylate (ACC), the precursor of ethylene biosynthesis, following the quantitative real-time polymerase chain reaction (RT-PCR) confirmation. Data obtained from RNA-Seq analysis indicated that EH-1 and ACC significantly induced the expression of 39 and 48 genes, respectively (above 20 fold of control), among which five genes are up-regulated by EH-1 as well as by ACC. We also found 67 and 32 genes that are significantly down-regulated, respectively, among which seven genes are in common. For quantitative RT-PCR analysis. 12 up-regulated genes were selected from the data obtained from RNA-Seq analysis. We found a good correlation of quantitative RT-PCR analysis and RNA-Seq analysis. Based on these results, we conclude that the action mechanism of EH-1 on inducing triple response in Arabidopsis is different from that of ACC.

A Functional Genomic Perspective on Drought Signalling and its Crosstalk with Phytohormone-mediated Signalling Pathways in Plants

INTRODUCTION

Drought stress is one of the most important abiotic stresses that negatively influence crop performance and productivity. Plants acclimatize to drought stress conditions through altered molecular, biochemical and physiological responses. Gene and/or protein expression and regulation are thought to be modulated upon stress perception and signal transduction for providing requisite endurance to plants.Plant growth regulators or phytohormones are important molecules required for various biological processes in plants and are also central to stress signalling pathways. Among various phytohormones, Abscisic Acid (ABA) and Ethylene (ET) are considered to be the most vital growth regulators implicated in drought stress signalling and tolerance. Besides the above two known classical phytohormones, Salicylic Acid (SA) and Jasmonic Acid (JA) have also been found to potentially enhance abiotic stress tolerance particularly that of drought, salinity, and heat stress tolerance in plants. Apart from these several other growth regulators such as Cytokinins (CKs), Auxin (AUX), Gibberellic Acid (GA), Brassinosteroids (BRs) and Strigolactones (SLs) have also been reported to actively participate in abiotic stress responses and tolerance in plants. The abiotic stress signalling in plants regulated by these hormones further depends upon the nature, intensity, and duration of exposure to various environmental stresses. It has been reported that all these phytohormones are also involved in extensive crosstalk and signal transduction among themselves and/or with other factors.

CONCLUSION

This review thus summarizes the molecular mechanism of drought signalling and its crosstalk with various phytohormone signalling pathways implicated in abiotic stress response and tolerance.

Co-expression network analysis of the transcriptomes of rice roots exposed to various cadmium stresses reveals universal cadmium-responsive genes

BACKGROUND

The migration of cadmium (Cd) from contaminated soil to rice is a cause for concern. However, the molecular mechanism underlying the response of rice roots to various Cd stresses remains to be clarified from the viewpoint of the co-expression network at a system-wide scale.

RESULTS

We employed a comparative RNAseq-based approach to identify early Cd-responsive differentially expressed genes (DEGs) in rice 'Nipponbare' seedling roots after 1 h of high-Cd treatment. A multiplicity of the identified 1772 DEGs were implicated in hormone signaling and transcriptional regulation, particularly NACs and WRKYs were all upregulated under Cd stress. All of the 6 Cd-upregulated ABC transporters were pleiotropic drug resistance proteins (PDRs), whereas all of the 6 ZRT/IRT-like proteins (ZIPs) were consistently downregulated by Cd treatment. To further confirm our results of this early transcriptomic response to Cd exposure, we then conducted weighted gene co-expression network analysis (WGCNA) to re-analyze our RNAseq data in combination with other 11 previously published RNAseq datasets for rice roots exposed to diverse concentrations of Cd for extended treatment periods. This integrative approach identified 271 transcripts as universal Cd-regulated DEGs that are key components of the Cd treatment coupled co-expression module. A global view of the 164 transcripts with annotated functions in pathway networks revealed several Cd-upregulated key functional genes, including transporter ABCG36/OsPDR9, 12-oxo-phytodienoic acid reductases (OPRs) for JA synthesis, and ZIM domain proteins JAZs in JA signaling, as well as OsWRKY10, NAC, and ZFP transcription factors. More importantly, 104 of these, including ABCG36/OsPDR9, OsNAC3, as well as several orthologs in group metalloendoproteinase, plastocyanin-like domain containing proteins and pectin methylesterase inhibitor, may respond specifically to various Cd pressures, after subtracting the 60 general stress-responsive genes reported to be commonly upregulated following multiple stresses.

CONCLUSION

An integrative approach was implemented to identify DEGs and co-expression network modules in response to various Cd pressures, and 104 of the 164 annotatable universal Cd-responsive DEGs may specifically respond to various Cd pressures. These results provide insight into the universal molecular mechanisms beneath the Cd response in rice roots, and suggest many promising targets for improving the rice acclimation process against Cd toxicity.

A Nucleus-Localized Long Non-Coding RNA Enhances Drought and Salt Stress Tolerance

Long noncoding RNAs (lncRNAs) affect gene expression through a wide range of mechanisms and are considered as important regulators in many essential biological processes. A large number of lncRNA transcripts have been predicted or identified in plants in recent years. However, the biological functions for most of them are still unknown. In this study, we identified an Arabidopsis (Arabidopsis thaliana) lncRNA, DROUGHT INDUCED lncRNA (DRIR), as a novel positive regulator of the plant response to drought and salt stress. DRIR was expressed at a low level under nonstress conditions but can be significantly activated by drought and salt stress as well as by abscisic acid (ABA) treatment. We identified a T-DNA insertion mutant, drirD , which had higher expression of the DRIR gene than the wild-type plants. The drirD mutant exhibits increased tolerance to drought and salt stress. Overexpressing DRIR in Arabidopsis also increased tolerance to drought and salt stress of the transgenic plants. The drirD mutant and the overexpressing seedlings are more sensitive to ABA than the wild type in stomata closure and seedling growth. Genome-wide transcriptome analysis demonstrated that the expression of a large number of genes was altered in drirD and the overexpressing plants. These include genes involved in ABA signaling, water transport, and other stress-relief processes. Our study reveals a mechanism whereby DRIR regulates the plant response to abiotic stress by modulating the expression of a series of genes involved in the stress response.

Multi-Omics and Integrated Network Analyses Reveal New Insights into the Systems Relationships between Metabolites, Structural Genes, and Transcriptional Regulators in Developing Grape Berries (Vitis vinifera L.) Exposed to Water Deficit

Grapes are one of the major fruit crops and they are cultivated in many dry environments. This study comprehensively characterizes the metabolic response of grape berries exposed to water deficit at different developmental stages. Increases of proline, branched-chain amino acids, phenylpropanoids, anthocyanins, and free volatile organic compounds have been previously observed in grape berries exposed to water deficit. Integrating RNA-sequencing analysis of the transcriptome with large-scale analysis of central and specialized metabolites, we reveal that these increases occur via a coordinated regulation of key structural pathway genes. Water deficit-induced up-regulation of flavonoid genes is also coordinated with the down-regulation of many stilbene synthases and a consistent decrease in stilbenoid concentration. Water deficit activated both ABA-dependent and ABA-independent signal transduction pathways by modulating the expression of several transcription factors. Gene-gene and gene-metabolite network analyses showed that water deficit-responsive transcription factors such as bZIPs, AP2/ERFs, MYBs, and NACs are implicated in the regulation of stress-responsive metabolites. Enrichment of known and novel cis-regulatory elements in the promoters of several ripening-specific/water deficit-induced modules further affirms the involvement of a transcription factor cross-talk in the berry response to water deficit. Together, our integrated approaches show that water deficit-regulated gene modules are strongly linked to key fruit-quality metabolites and multiple signal transduction pathways may be critical to achieve a balance between the regulation of the stress-response and the berry ripening program. This study constitutes an invaluable resource for future discoveries and comparative studies, in grapes and other fruits, centered on reproductive tissue metabolism under abiotic stress.

Multi-Omics and Integrated Network Analyses Reveal New Insights into the Systems Relationships between Metabolites, Structural Genes, and Transcriptional Regulators in Developing Grape Berries (Vitis vinifera L.) Exposed to Water Deficit

Grapes are one of the major fruit crops and they are cultivated in many dry environments. This study comprehensively characterizes the metabolic response of grape berries exposed to water deficit at different developmental stages. Increases of proline, branched-chain amino acids, phenylpropanoids, anthocyanins, and free volatile organic compounds have been previously observed in grape berries exposed to water deficit. Integrating RNA-sequencing analysis of the transcriptome with large-scale analysis of central and specialized metabolites, we reveal that these increases occur via a coordinated regulation of key structural pathway genes. Water deficit-induced up-regulation of flavonoid genes is also coordinated with the down-regulation of many stilbene synthases and a consistent decrease in stilbenoid concentration. Water deficit activated both ABA-dependent and ABA-independent signal transduction pathways by modulating the expression of several transcription factors. Gene-gene and gene-metabolite network analyses showed that water deficit-responsive transcription factors such as bZIPs, AP2/ERFs, MYBs, and NACs are implicated in the regulation of stress-responsive metabolites. Enrichment of known and novel cis-regulatory elements in the promoters of several ripening-specific/water deficit-induced modules further affirms the involvement of a transcription factor cross-talk in the berry response to water deficit. Together, our integrated approaches show that water deficit-regulated gene modules are strongly linked to key fruit-quality metabolites and multiple signal transduction pathways may be critical to achieve a balance between the regulation of the stress-response and the berry ripening program. This study constitutes an invaluable resource for future discoveries and comparative studies, in grapes and other fruits, centered on reproductive tissue metabolism under abiotic stress.

Neofusicoccum parvum Colonization of the Grapevine Woody Stem Triggers Asynchronous Host Responses at the Site of Infection and in the Leaves

Grapevine trunk diseases cause important economic losses in vineyards worldwide. Neofusicoccum parvum, one of the most aggressive causal agents of the trunk disease Botryosphaeria dieback, colonizes cells and tissues of the grapevine wood, leading to the formation of an internal canker. Symptoms then extend to distal shoots, with wilting of leaves and bud mortality. Our aim was to characterize the transcriptional dynamics of grapevine genes in the woody stem and in the leaves during Neofusicoccum parvum colonization. Genome-wide transcriptional profiling at seven distinct time points (0, 3, and 24 hours; 2, 6, 8, and 12 weeks) showed that both stems and leaves undergo extensive transcriptomic reprogramming in response to infection of the stem. While most intense transcriptional responses were detected in the stems at 24 hours, strong responses were not detected in the leaves until the next sampling point at 2 weeks post-inoculation. Network co-expression analysis identified modules of co-expressed genes common to both organs and showed most of these genes were asynchronously modulated. The temporal shift between stem vs. leaf responses affected transcriptional modulation of genes involved in both signal perception and transduction, as well as downstream biological processes, including oxidative stress, cell wall rearrangement and cell death. Promoter analysis of the genes asynchronously modulated in stem and leaves during N. parvum colonization suggests that the temporal shift of transcriptional reprogramming between the two organs might be due to asynchronous co-regulation by common transcriptional regulators. Topology analysis of stem and leaf co-expression networks pointed to specific transcription factor-encoding genes, including WRKY and MYB, which may be associated with the observed transcriptional responses in the two organs.

The common transcriptional subnetworks of the grape berry skin in the late stages of ripening

BACKGROUND

Wine grapes are important economically in many countries around the world. Defining the optimum time for grape harvest is a major challenge to the grower and winemaker. Berry skins are an important source of flavor, color and other quality traits in the ripening stage. Senescent-like processes such as chloroplast disorganization and cell death characterize the late ripening stage.

RESULTS

To better understand the molecular and physiological processes involved in the late stages of berry ripening, RNA-seq analysis of the skins of seven wine grape cultivars (Cabernet Franc, Cabernet Sauvignon, Merlot, Pinot Noir, Chardonnay, Sauvignon Blanc and Semillon) was performed. RNA-seq analysis identified approximately 2000 common differentially expressed genes for all seven cultivars across four different berry sugar levels (20 to 26 °Brix). Network analyses, both a posteriori (standard) and a priori (gene co-expression network analysis), were used to elucidate transcriptional subnetworks and hub genes associated with traits in the berry skins of the late stages of berry ripening. These independent approaches revealed genes involved in photosynthesis, catabolism, and nucleotide metabolism. The transcript abundance of most photosynthetic genes declined with increasing sugar levels in the berries. The transcript abundance of other processes increased such as nucleic acid metabolism, chromosome organization and lipid catabolism. Weighted gene co-expression network analysis (WGCNA) identified 64 gene modules that were organized into 12 subnetworks of three modules or more and six higher order gene subnetworks. Some gene subnetworks were highly correlated with sugar levels and some subnetworks were highly enriched in the chloroplast and nucleus. The petal R package was utilized independently to construct a true small-world and scale-free complex gene co-expression network model. A subnetwork of 216 genes with the highest connectivity was elucidated, consistent with the module results from WGCNA. Hub genes in these subnetworks were identified including numerous members of the core circadian clock, RNA splicing, proteolysis and chromosome organization. An integrated model was constructed linking light sensing with alternative splicing, chromosome remodeling and the circadian clock.

CONCLUSIONS

A common set of differentially expressed genes and gene subnetworks from seven different cultivars were examined in the skin of the late stages of grapevine berry ripening. A densely connected gene subnetwork was elucidated involving a complex interaction of berry senescent processes (autophagy), catabolism, the circadian clock, RNA splicing, proteolysis and epigenetic regulation. Hypotheses were induced from these data sets involving sugar accumulation, light, autophagy, epigenetic regulation, and fruit development. This work provides a better understanding of berry development and the transcriptional processes involved in the late stages of ripening.

Design of a Custom RT-qPCR Array for Assignment of Abiotic Stress Tolerance in Traditional Portuguese Grapevine Varieties

Widespread agricultural losses attributed to drought, often combined with high temperatures, frequently occur in the field, namely in Mediterranean climate areas, where the existing scenarios for climate change indicate an increase in the frequency of heat waves and severe drought events in summer. Grapevine (Vitis vinifera L.) is the most cultivated fruit species in the world and the most valuable one and is a traditional Mediterranean species. Currently, viticulture must adjust to impending climate changes that are already pushing vine-growers toward the use of ancient and resilient varieties. Portugal is very rich in grapevine biodiversity, however, currently, 90% of the total producing area is planted with only 16 varieties. There is a pressing need to understand the existing genetic diversity and the physiological potential of the varieties/genotypes available to be able to respond to climate changes. With the above scenario in mind, an assembly of 65 differentially expresses genes (DEGs) previously identified as responsive to abiotic stresses in two well studied genotypes, 'Touriga Nacional' and 'Trincadeira,' was designed to scan the gene expression of leaf samples from 10 traditional Portuguese varieties growing in two regions with distinct environmental conditions. Forty-five of those DEGs proved to be associated to "abiotic stress" and were chosen to build a custom qPCR array to identify uncharacterized genotypes as sensitive or tolerant to abiotic stress. According to the experimental set-up behind the array design these DEGs can also be used as indicators of the main abiotic stress that the plant is subjected and responding to (drought, heat, or excess light).

Transcriptomic comparison between two Vitis vinifera L. varieties (Trincadeira and Touriga Nacional) in abiotic stress conditions

BACKGROUND

Predicted climate changes announce an increase of extreme environmental conditions including drought and excessive heat and light in classical viticultural regions. Thus, understanding how grapevine responds to these conditions and how different genotypes can adapt, is crucial for informed decisions on accurate viticultural actions. Global transcriptome analyses are useful for this purpose as the response to these abiotic stresses involves the interplay of complex and diverse cascades of physiological, cellular and molecular events. The main goal of the present work was to evaluate the response to diverse imposed abiotic stresses at the transcriptome level and to compare the response of two grapevine varieties with contrasting physiological trends, Trincadeira (TR) and Touriga Nacional (TN).

RESULTS

Leaf transcriptomic response upon heat, high light and drought treatments in growth room controlled conditions, as well as full irrigation and non-irrigation treatments in the field, was compared in TR and TN using GrapeGene GeneChips®. Breakdown of metabolism in response to all treatments was evidenced by the functional annotation of down-regulated genes. However, circa 30 % of the detected stress-responsive genes are still annotated as «Unknown» function. Selected differentially expressed genes from the GrapeGene GeneChip® were analysed by RT-qPCR in leaves of growth room plants under the combination of individual stresses and of field plants, in both varieties. The transcriptomic results correlated better with those obtained after each individual stress than with the results of plants from field conditions.

CONCLUSIONS

From the transcriptomic comparison between the two Portuguese grapevine varieties Trincadeira and Touriga Nacional under abiotic stress main conclusions can be drawn: 1. A different level of tolerance to stress is evidenced by a lower transcriptome reprogramming in TN than in TR. Interestingly, this lack of response in TN associates with its higher adaptation to extreme conditions including environmental conditions in a changing climate; 2. A complex interplay between stress transcriptional cascades is evidenced by antagonistic and, in lower frequency, synergistic effects on gene expression when several stresses are imposed together; 3. The grapevine responses to stress under controlled conditions are not fully extrapolated to the complex vineyard scenario and should be cautiously considered for agronomic management decision purposes.