WRKY transcription factors in plant responses to stresses

The Wheat GENIE3 Network Provides Biologically-Relevant Information in Polyploid Wheat

Gene regulatory networks are powerful tools which facilitate hypothesis generation and candidate gene discovery. However, the extent to which the network predictions are biologically relevant is often unclear. Recently a GENIE3 network which predicted targets of wheat transcription factors was produced. Here we used an independent RNA-Seq dataset to test the predictions of the wheat GENIE3 network for the senescence-regulating transcription factor NAM-A1 (TraesCS6A02G108300). We re-analyzed the RNA-Seq data against the RefSeqv1.0 genome and identified a set of differentially expressed genes (DEGs) between the wild-type and nam-a1 mutant which recapitulated the known role of NAM-A1 in senescence and nutrient remobilisation. We found that the GENIE3-predicted target genes of NAM-A1 overlap significantly with the DEGs, more than would be expected by chance. Based on high levels of overlap between GENIE3-predicted target genes and the DEGs, we identified candidate senescence regulators. We then explored genome-wide trends in the network related to polyploidy and found that only homeologous transcription factors are likely to share predicted targets in common. However, homeologs which vary in expression levels across tissues are less likely to share predicted targets than those that do not, suggesting that they may be more likely to act in distinct pathways. This work demonstrates that the wheat GENIE3 network can provide biologically-relevant predictions of transcription factor targets, which can be used for candidate gene prediction and for global analyses of transcription factor function. The GENIE3 network has now been integrated into the KnetMiner web application, facilitating its use in future studies.

MdWRKY30, a group IIa WRKY gene from apple, confers tolerance to salinity and osmotic stresses in transgenic apple callus and Arabidopsis seedlings

Abiotic stresses threaten the productivity and quality of economically important perennial fruit crops such as apple (Malus × domestica Borkh.). WRKY transcription factors play various roles in plant responses to abiotic stress, but little is known regarding WRKY genes in apple. Here, we carried out functional characterization of an apple Group IIa WRKY gene (MdWRKY30). qRT-PCR analysis found that MdWRKY30 expression was induced by salt and drought stress. A subcellular localization assay showed that MdWRKY30 is localized to the nucleus. A transactivation assay found that MdWRKY30 has no transcriptional activation activity. A Y2H assay indicated that MdWRKY26, MdWRKY28, and MdWRKY30 interact with each other to form heterodimers and homodimers. Transgenic analysis revealed that the overexpression of MdWRKY30 in Arabidopsis enhanced salt and osmotic tolerance in the seedling stage, as well as during the seed germination and greening cotyledon stages. MdWRKY30 overexpression enhanced tolerance to salt and osmotic stresses in transgenic apple callus through transcriptional regulation of stress-related genes. Together, our results demonstrate that MdWRKY30 is an important regulator of salinity and osmotic stress tolerance in apple.

Molecular characterization and functional analysis of a novel WRKY transcription factor HbWRKY83 possibly involved in rubber production of Hevea brasiliensis

WRKY transcription factors play important roles in plant growth and developmental processes and various stress responses, and are also associated with jasmonic acid (JA) signaling in the regulation of secondary metabolite biosynthesis in plants. The regulatory networks mediated by WRKY proteins in the latex production of Hevea brasiliensis (the Pará rubber tree) are poorly understood. In this study, one novel WRKY gene (designated HbWRKY83) was identified from the latex of H. brasiliensis, and its functions were characterized via gene expression analysis in both the latex and HbWRKY83-overexpressing transgenic Arabidopsis. HbWRKY83 gene contains an open reading frame (ORF) of 921 bp encoding a 306-amino-acid protein which is clustered with group IIc WRKY TF. HbWRKY83 is a nuclear-localized protein with transcriptional activity. Real-time quantitative PCR analysis demonstrated that the transcription level of HbWRKY83 was up-regulated by exogenous methyl jasmonate, Ethrel (ethylene releaser) stimulation, and bark tapping (mechanical wounding). Compared with the wild-type plants, overexpression of HbWRKY83 improved the tolerance of transgenic Arabidopsis lines to drought and salt stresses by enhancing the expression levels of ethylene-insensitive3 transcription factors (EIN3s) and several stress-responsive genes, including Cu/Zn superoxide dismutases CSD1 (Cu/Zn-SOD1) and CSD2 (Cu/Zn-SOD2), related to reactive oxygen species scavenging. Additionally, these genes were also significantly up-regulated by bark tapping. In combination, these results suggest that HbWRKY83 might act as a positive regulator of rubber production by activating the expression of JA-, ethylene-, and wound-responsive genes in the laticiferous cells of rubber trees.

Investigation of an Antioxidative System for Salinity Tolerance in Oenanthe javanica

Abiotic stress, such as drought and salinity, severely affect the growth and yield of many plants. Oenanthe javanica (commonly known as water dropwort) is an important vegetable that is grown in the saline-alkali soils of East Asia, where salinity is the limiting environmental factor. To study the defense mechanism of salt stress responses in water dropwort, we studied two water dropwort cultivars, V11E0022 and V11E0135, based on phenotypic and physiological indexes. We found that V11E0022 were tolerant to salt stress, as a result of good antioxidant defense system in the form of osmolyte (proline), antioxidants (polyphenols and flavonoids), and antioxidant enzymes (APX and CAT), which provided novel insights for salt-tolerant mechanisms. Then, a comparative transcriptomic analysis was conducted, and Gene Ontology (GO) analysis revealed that differentially expressed genes (DEGs) involved in the carbohydrate metabolic process could reduce oxidative stress and enhance energy production that can help in adaptation against salt stress. Similarly, lipid metabolic processes can also enhance tolerance against salt stress by reducing the transpiration rate, H₂O₂, and oxidative stress. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs involved in hormone signals transduction pathway promoted the activities of antioxidant enzymes and reduced oxidative stress; likewise, arginine and proline metabolism, and flavonoid pathways also stimulated the biosynthesis of proline and flavonoids, respectively, in response to salt stress. Moreover, transcription factors (TFs) were also identified, which play an important role in salt stress tolerance of water dropwort. The finding of this study will be helpful for crop improvement under salt stress.

Suppression of a leucine-rich repeat receptor-like kinase enhances host plant resistance to a specialist herbivore

The mechanisms by which herbivores induce plant defenses are well studied. However, how specialized herbivores suppress plant resistance is still poorly understood. Here, we discovered a rice (Oryza sativa) leucine-rich repeat receptor-like kinase, OsLRR-RLK2, which is induced upon attack by gravid females of a specialist piercing-sucking herbivore, the brown planthopper (BPH, Nilaparvata lugens). Silencing OsLRR-RLK2 decreases the constitutive activity of mitogen-activated protein kinase (OsMPK6) and alters BPH-induced transcript levels of several defense-related WRKY transcription factors. Moreover, silencing OsLRR-RLK2 reduces BPH-induction of jasmonic acid and ethylene but promotes the biosynthesis of both elicited salicylic acid and H₂ O₂ ; silencing also enhances the production of volatiles emitted from rice plants infested with gravid BPH females. These changes decrease BPH preference and performance in the glasshouse and the field. These findings suggest that OsLRR-RLK2, by regulating the plant's defense-related signaling profile, increases the susceptibility of rice to BPH, and that BPH infestation influences the expression of OsLRR-RLK2, suppressing the resistance of rice to BPH.

Identifying key regulatory genes of maize root growth and development by RNA sequencing

Root system architecture (RSA), the spatio-temporal configuration of roots, plays vital roles in maize (Zea mays L.) development and productivity. We sequenced the maize root transcriptome of four key growth and development stages: the 6th leaf stage, the 12th leaf stage, the tasseling stage and the milk-ripe stage. Differentially expressed genes (DEGs) were detected. 81 DEGs involved in plant hormone signal transduction pathway and 26 transcription factor (TF) genes were screened. These DEGs and TFs were predicted to be potential candidate genes during maize root growth and development. Several of these genes are homologous to well-known genes regulating root architecture or development in Arabidopsis or rice, such as, Zm00001d005892 (AtERF109), Zm00001d027925 (AtERF73/HRE1), Zm00001d047017 (AtMYC2, OsMYC2), Zm00001d039245 (AtWRKY6). Identification of these key genes will provide a further understanding of the molecular mechanisms responsible for maize root growth and development, it will be beneficial to increase maize production and improve stress resistance by altering RSA traits in modern breeding.

Co-expression network analysis uncovers key candidate genes related to the regulation of volatile esters accumulation in Woodland strawberry

FveERF (FvH4_5g04470.1), FveAP2 (FvH4_1g16370.1) and FveWRKY (FvH4_6g42870.1) might be involved in fruit maturation of strawberry. Overexpression of FveERF could activate the expression of AAT gene and ester accumulation. Volatile esters play an important role in the aroma of strawberry fruits, whose flavor is the result of a complex mixture of various esters. The accumulation of these volatiles is closely tied to changes in metabolism during fruit ripening. Acyltransferase (AAT) is recognized as having a significant effect in ester formation. However, there is little knowledge about the regulation network of AAT. Here, we collected the data of RNA-seq and headspace GC-MS at five time points during fruit maturation of Hawaii4 and Ruegen strawberry varieties. A total of 106 volatile compounds were identified in the fruit of woodland strawberries, including 58 esters, which occupied 41.09% (Hawaii4) or 33.40% (Ruegen) of total volatile concentration. Transcriptome analysis revealed eight transcription factors highly associated with AAT genes. Through the changes in esters and the weight co-expression network analysis (WGCNA), a detailed gene network was established. This demonstrated that ERF gene (FvH4_5g04470.1), AP2 gene (FvH4_1g16370.1) and one WRKY gene (FvH4_6g42870.1) might be involved in expression of AAT genes, especially ERF genes. Overexpression of FveERF (FvH4_5g04470.1) does activate expression of AAT genes and ester accumulation in fruits of strawberry. Our findings provide valuable clues to gain better insight into the ester formation process of numerous fruits.

Identification of QTLs for grain yield and other traits in tropical maize under Striga infestation

Striga is an important biotic factor limiting maize production in sub-Saharan Africa and can cause yield losses as high as 100%. Marker-assisted selection (MAS) approaches hold a great potential for improving Striga resistance but requires identification and use of markers associated with Striga resistance for adequate genetic gains from selection. However, there is no report on the discovery of quantitative trait loci (QTL) for resistance to Striga in maize under artificial field infestation. In the present study, 198 BC₁S₁ families obtained from a cross involving TZEEI 29 (Striga resistant inbred line) and TZEEI 23 (Striga susceptible inbred line) plus the two parental lines were screened under artificial Striga-infested conditions at two Striga-endemic locations in Nigeria in 2018, to identify QTL associated with Striga resistance indicator traits, including grain yield, ears per plant, Striga damage and number of emerged Striga plants. Genetic map was constructed using 1,386 DArTseq markers distributed across the 10 maize chromosomes, covering 2076 cM of the total genome with a mean spacing of 0.11 cM between the markers. Using composite interval mapping (CIM), fourteen QTL were identified for key Striga resistance/tolerance indicator traits: 3 QTL for grain yield, 4 for ears per plant and 7 for Striga damage at 10 weeks after planting (WAP), across environments. Putative candidate genes which encode major transcription factor families WRKY, bHLH, AP2-EREBPs, MYB, and bZIP involved in plant defense signaling were detected for Striga resistance/tolerance indicator traits. The QTL detected in the present study would be useful for rapid transfer of Striga resistance/tolerance genes into Striga susceptible but high yielding maize genotypes using MAS approaches after validation. Further studies on validation of the QTL in different genetic backgrounds and in different environments would help verify their reproducibility and effective use in breeding for Striga resistance/tolerance.

Overexpression of transcription factor SlWRKY28 improved the tolerance of Populus davidiana × P. bolleana to alkaline salt stress

BACKGROUND

WRKY transcription factors (TFs) have been suggested to play crucial roles in the response to biotic and abiotic stresses. This study is the first to report the alkaline salt regulation of the WRKY gene.

RESULTS

In this study, we cloned a WRKY gene (SlWRKY28) from the Salix linearistipularis and then transferred to the Populus davidiana × P. bolleana for expression. Sequence analysis on the transcriptome of Salix linearistipular showed the significant up-regulation of WRKY gene expression in response to salt-alkali stress in seedlings. Our data showed that SlWRKY28 localized to the nucleus. Furthermore, the expression of the SlWRKY28 from female plants increased with saline-alkali stress according to the northern blot analysis results. The results of 3,3'-Diaminobenzidine (DAB) staining showed that hydrogen peroxide (H2O2) concentration was lower under stress, but ascorbate peroxidase (APX) enzyme activity was significantly higher in the overexpressed plants than that in non-transgenic (NT) plants.

CONCLUSIONS

We found out the SlWRKY28 induced regulation of the enzyme gene in the reactive oxygen species (ROS) scavenging pathway is a potential mechanism for transgenic lines to improve their resistance to alkaline salt. This study shows theoretical and practical significance in determining SlWRKY28 transcription factors involved in the regulation of alkaline salt tolerance.

Origin and evolution of jasmonate signaling

Jasmonate (JA) signaling is a key mediator of plant development and defense which arose during plants transition from an aqueous to terrestrial environment. Elucidating the evolution of JA signaling is important for understanding plant development, defense, and production of specialized metabolites. The lineage of key protein domains characterizing JA signaling factors was traced to identify the origins of CORONITINE INSENSITIVE 1 (COI1), JASMONATE ZIM-DOMAIN (JAZ), NOVEL INTERACTOR OF JAZ, MYC2, TOPLESS, and MEDIATOR SUBUNIT 25. Charophytes do not possess genes encoding key JA signaling components, including COI1, JAZ, MYC2, and the JAZ-interacting bHLH factors, yet their orthologs are present in bryophytes. TIFY family genes were found in charophyta and chlorophya algae. JAZs evolved from ZIM genes of the TIFY family through changes to several key amino acids. Dating placed the origin of JA signaling 515 to 473 million years ago during the middle Cambrian to early Ordovician periods. This time is known for rapid biodiversification and mass extinction events. An increased predation from the diversifying and changing fauna may have driven evolution of JA signaling and plant defense.

Integrated transcriptomic and metabolomic analyses of yellow horn (Xanthoceras sorbifolia) in response to cold stress

Xanthoceras sorbifolia, a medicinal and oil-rich woody plant, has great potential for biodiesel production. However, little study explores the link between gene expression level and metabolite accumulation of X. sorbifolia in response to cold stress. Herein, we performed both transcriptomic and metabolomic analyses of X. sorbifolia seedlings to investigate the regulatory mechanism of resistance to low temperature (4 °C) based on physiological profile analyses. Cold stress resulted in a significant increase in the malondialdehyde content, electrolyte leakage and activity of antioxidant enzymes. A total of 1,527 common differentially expressed genes (DEGs) were identified, of which 895 were upregulated and 632 were downregulated. Annotation of DEGs revealed that amino acid metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, galactose metabolism, fructose and mannose metabolism, and the citrate cycle (TCA) were strongly affected by cold stress. In addition, DEGs within the plant mitogen-activated protein kinase (MAPK) signaling pathway and TF families of ERF, WRKY, NAC, MYB, and bHLH were transcriptionally activated. Through metabolomic analysis, we found 51 significantly changed metabolites, particularly with the analysis of primary metabolites, such as sugars, amino acids, and organic acids. Moreover, there is an overlap between transcript and metabolite profiles. Association analysis between key genes and altered metabolites indicated that amino acid metabolism and sugar metabolism were enhanced. A large number of specific cold-responsive genes and metabolites highlight a comprehensive regulatory mechanism, which will contribute to a deeper understanding of the highly complex regulatory program under cold stress in X. sorbifolia.

VcMYB4a, an R2R3-MYB transcription factor from Vaccinium corymbosum, negatively regulates salt, drought, and temperature stress

MYB transcription factors (TFs) play important roles in the plant's response to abiotic stress. In this study, we cloned a novel MYB TF gene from Vaccinium corymbosum (blueberry) using rapid amplification of cDNA ends (RACE). The cDNA contained a 798-bp open reading frame that encodes a 265-amino acid protein. VcMYB4a possessed a C2/EAR-repressor motif domain and phylogenetic analysis showed that it clustered into a subgroup 4 with six Arabidopsis thaliana MYBs. Quantitative RT-PCR analysis demonstrated that VcMYB4a expression was downregulated by salt, drought, and cold treatment, but was induced by freezing and heat. Overexpression of VcMYB4a in blueberry callus enhanced sensitivity to salt, drought, cold, freezing, and heat stress. These results indicate that VcMYB4a may be an important repressor of abiotic stress in blueberry.

Screening of abiotic stress-responsive cotton genes using a cotton full-length cDNA overexpressing Arabidopsis library

Cotton (Gossypium hirsutum L.) is a major crop and the main source of natural fiber worldwide. Because various abiotic and biotic stresses strongly influence cotton fiber yield and quality, improved stress resistance of this crop plant is urgently needed. In this study, we used Gateway technology to construct a normalized full-length cDNA overexpressing (FOX) library from upland cotton cultivar ZM12 under various stress conditions. The library was transformed into Arabidopsis to produce a cotton-FOX-Arabidopsis library. Screening of this library yielded 6,830 transgenic Arabidopsis lines, of which 757 were selected for sequencing to ultimately obtain 659 cotton ESTs. GO and KEGG analyses mapped most of the cotton ESTs to plant biological process, cellular component, and molecular function categories. Next, 156 potential stress-responsive cotton genes were identified from the cotton-FOX-Arabidopsis library under drought, salt, ABA, and other stress conditions. Four stress-related genes identified from the library, designated as GhCAS, GhAPX, GhSDH, and GhPOD, were cloned from cotton complementary DNA, and their expression patterns under stress were analyzed. Phenotypic experiments indicated that overexpression of these cotton genes in Arabidopsis affected the response to abiotic stress. The method developed in this study lays a foundation for high-throughput cloning and rapid identification of cotton functional genes.

The specific W-boxes of GAPC5 promoter bound by TaWRKY are involved in drought stress response in wheat

Drought is one of the most common abiotic stresses, and can limit wheat yield, crops and productivity. GAPCs play vital roles under drought stress conditions in multiple species. The aim of this experiment was to determine the regulatory mechanism of TaGAPC5 under drought stress. In this study, the genes and promoters of TaGAPC5 in diverse drought-tolerant cultivars were cloned. The amino acid sequences were conserved, while the promoter fragments were not identical. Under abiotic stress, the expression level of TaGAPC5 was substantially different among the diverse drought-tolerant cultivars and the promoter activities were significantly improved. The yeast one-hybrid system and Electrophoretic mobility shift assay (EMSA) demonstrated that TaWRKYs bound to specific W-boxes: TaWRKY28, TaWRKY33, TaWRKY40 and TaWRKY47 bind to G/ATGACG/C/A, C/G/ATGACG, C/ATGACC and C/ATGACC/G, respectively. By analyzing different 5' deletion mutants of these promoters, it was determined that these W-boxes in CW-TaGAPC5 promoter (-1262, -1202, -904, -880 and -207) and ZY-TaGAPC5 promoter (-697 and -220) bound by these four TaWRKYs and were functional under drought stress. The deletion or addition of specific W-boxes in the promoter fragments significantly restrained or advanced the promoter activity under drought stress, and these results further confirmed that these W-boxes play vital roles in improving transcription levels under drought stress. The W-boxes in CW-TaGAPC5P (-1262, -1202, -904, -880 and -207) and ZY-TaGAPC5P (-697 and -220) were identified as the key cis-elements for responding to drought stress and were bound by the transcription factor TaWRKY.

Crop reproductive meristems in the genomic era: a brief overview

Modulation of traits beneficial for cultivation and yield is one of the main goals of crop improvement. One of the targets for enhancing productivity is changing the architecture of inflorescences since in many species it determines fruit and seed yield. Inflorescence shape and organization is genetically established during the early stages of reproductive development and depends on the number, arrangement, activities, and duration of meristems during the reproductive phase of the plant life cycle. Despite the variety of inflorescence architectures observable in nature, many key aspects of inflorescence development are conserved among different species. For instance, the genetic network in charge of specifying the identity of the different reproductive meristems, which can be indeterminate or determinate, seems to be similar among distantly related species. The availability of a large number of published transcriptomic datasets for plants with different inflorescence architectures, allowed us to identify transcription factor gene families that are differentially expressed in determinate and indeterminate reproductive meristems. The data that we review here for Arabidopsis, rice, barley, wheat, and maize, particularly deepens our knowledge of their involvement in meristem identity specification.

Transcriptome Analysis Reveals Potential Roles of Abscisic Acid and Polyphenols in Adaptation of Onobrychis viciifolia to Extreme Environmental Conditions in the Qinghai-Tibetan Plateau

A detailed understanding of the molecular mechanisms of plant stress resistance in the face of ever-changing environmental stimuli will be helpful for promoting the growth and production of crop and forage plants. Investigations of plant responses to various single abiotic or biotic factors, or combined stresses, have been extensively reported. However, the molecular mechanisms of plants in responses to environmental stresses under natural conditions are not clearly understood. In this study, we carried out a transcriptome analysis using RNA-sequencing to decipher the underlying molecular mechanisms of Onobrychis viciifolia responding and adapting to the extreme natural environment in the Qinghai-Tibetan Plateau (QTP). The transcriptome data of plant samples collected from two different altitudes revealed a total of 8212 differentially expressed genes (DEGs), including 5387 up-regulated and 2825 down-regulated genes. Detailed analysis of the identified DEGs uncovered that up-regulation of genes potentially leading to changes in hormone homeostasis and signaling, particularly abscisic acid-related ones, and enhanced biosynthesis of polyphenols play vital roles in the adaptive processes of O. viciifolia. Interestingly, several DEGs encoding uridine diphosphate glycosyltransferases, which putatively regulate phytohormone homeostasis to resist environmental stresses, were also discovered. Furthermore, numerous DEGs encoding transcriptional factors, such as members of the myeloblastosis (MYB), homeodomain-leucine zipper (HD-ZIP), WRKY, and nam-ataf1,2-cuc2 (NAC) families, might be involved in the adaptive responses of O. viciifolia to the extreme natural environmental conditions. The DEGs identified in this study represent candidate targets for improving environmental stress resistance of O. viciifolia grown in higher altitudes of the QTP, and can provide deep insights into the molecular mechanisms underlying the responses of this plant species to the extreme natural environmental conditions of the QTP.

Genome-Wide Identification of the Physic Nut WUSCHEL-Related Homeobox Gene Family and Functional Analysis of the Abiotic Stress Responsive Gene JcWOX5

Plant-specific WOX transcription factors have important regulatory functions in plant development and response to abiotic stress. However, the identification and functional analysis of members of the WOX family have rarely been reported in the physic nut plant until now. Our research identified 12 WOX genes (JcWOXs) in physic nut, and these genes were divided into three groups corresponding to the ancient clade, WUS clade, and intermediate clade. Expression analysis based on RNA-seq and qRT-PCR showed that most of the JcWOX genes were expressed in at least one of the tissues tested, whereas five genes were identified as being highly responsive to drought and salt stresses. Subcellular localization analysis in Arabidopsis protoplast cells showed that JcWOX5 encoded a nuclear-localized protein. JcWOX5-overexpression plants increased sensitivity to drought stress, and transgenic plants suggested a lower proline content and CAT activity, higher relative electrolyte leakage, higher MDA content, and higher rate of water loss under drought conditions. Expression of some stress-related genes was obviously lower in the transformed rice lines as compared to their expression in wild-type rice lines under drought stress. Further data on JcWOX5-overexpressing plants reducing drought tolerance verified the potential role of JcWOX genes in responsive to abiotic stress. Collectively, the study provides a foundation for further functional analysis of JcWOX genes and the improvement of physic nut crops.

Regional differences in the abiotic environment contribute to genomic divergence within a wild tomato species

The wild currant tomato Solanum pimpinellifolium inhabits a wide range of abiotic habitats across its native range of Ecuador and Peru. Although it has served as a key genetic resource for the improvement of domestic cultivars, little is known about the genetic basis of traits underlying local adaptation in this species, nor what abiotic variables are most important for driving differentiation. Here we use redundancy analysis (RDA) and other multivariate statistical methods (structural equation modelling [SEM] and generalized dissimilarity modelling [GDM]) to quantify the relationship of genomic variation (6,830 single nucleotide polymorphisms [SNPs]) with climate and geography, among 140 wild accessions. RDA, SEM and GDM each identified environment as explaining more genomic variation than geography, suggesting that local adaptation to heterogeneous abiotic habitats may be an important source of genetic diversity in this species. Environmental factors describing temporal variation in precipitation and evaporative demand explained the most SNP variation among accessions, indicating that these forces may represent key selective agents. Lastly, by studying how SNP-environment associations vary throughout the genome (44,064 SNPs), we mapped the location and investigated the functions of loci putatively contributing to climatic adaptations. Together, our findings indicate an important role for selection imposed by the abiotic environment in driving genomic differentiation between populations.

Identification and expression analysis of the DREB transcription factor family in pineapple (Ananas comosus (L.) Merr.)

Background

Dehydration responsive element-binding (DREB) transcription factors play a crucial role in plant growth, development and stress responses. Although DREB genes have been characterized in many plant species, genome-wide identification of the DREB gene family has not yet been reported in pineapple (Ananas comosus (L.) Merr.).

Results

Using comprehensive genome-wide screening, we identified 20 AcoDREB genes on 14 chromosomes. These were categorized into five subgroups. AcoDREBs within a group had similar gene structures and domain compositions. Using gene structure analysis, we showed that most AcoDREB genes (75%) lacked introns, and that the promoter regions of all 20 AcoDREB genes had at least one stress response-related cis-element. We identified four genes with high expression levels and six genes with low expression levels in all analyzed tissues. We detected expression changes under abiotic stress for eight selected AcoDREB genes.

Conclusions

This report presents the first genome-wide analysis of the DREB transcription factor family in pineapple. Our results provide preliminary data for future functional analysis of AcoDREB genes in pineapple, and useful information for developing new pineapple varieties with key agronomic traits such as stress tolerance.

WRKY genes family study reveals tissue-specific and stress-responsive TFs in wild potato species

Wild potatoes, as dynamic resource adapted to various environmental conditions, represent a powerful and informative reservoir of genes useful for breeding efforts. WRKY transcription factors (TFs) are encoded by one of the largest families in plants and are involved in several biological processes such as growth and development, signal transduction, and plant defence against stress. In this study, 79 and 84 genes encoding putative WRKY TFs have been identified in two wild potato relatives, Solanum commersonii and S. chacoense. Phylogenetic analysis of WRKY proteins divided ScWRKYs and SchWRKYs into three Groups and seven subGroups. Structural and phylogenetic comparative analyses suggested an interspecific variability of WRKYs. Analysis of gene expression profiles in different tissues and under various stresses allowed to select ScWRKY045 as a good candidate in wounding-response, ScWRKY055 as a bacterial infection triggered WRKY and ScWRKY023 as a multiple stress-responsive WRKY gene. Those WRKYs were further studied through interactome analysis allowing the identification of potential co-expression relationships between ScWRKYs/SchWRKYs and genes of various pathways. Overall, this study enabled the discrimination of WRKY genes that could be considered as potential candidates in both breeding programs and functional studies.

Molecular adaptation to salinity fluctuation in tropical intertidal environments of a mangrove tree Sonneratia alba

BACKGROUND

Mangroves have adapted to intertidal zones - the interface between terrestrial and marine ecosystems. Various studies have shown adaptive evolution in mangroves at physiological, ecological, and genomic levels. However, these studies paid little attention to gene regulation of salt adaptation by transcriptome profiles.

RESULTS

We sequenced the transcriptomes of Sonneratia alba under low (fresh water), medium (half the seawater salinity), and high salt (seawater salinity) conditions and investigated the underlying transcriptional regulation of salt adaptation. In leaf tissue, 64% potential salinity-related genes were not differentially expressed when salinity increased from freshwater to medium levels, but became up- or down-regulated when salt concentrations further increased to levels found in sea water, indicating that these genes are well adapted to the medium saline condition. We inferred that both maintenance and regulation of cellular environmental homeostasis are important adaptive processes in S. alba. i) The sulfur metabolism as well as flavone and flavonol biosynthesis KEGG pathways were significantly enriched among up-regulated genes in leaves. They are both involved in scavenging ROS or synthesis and accumulation of osmosis-related metabolites in plants. ii) There was a significantly increased percentage of transcription factor-encoding genes among up-regulated transcripts. High expressions of salt tolerance-related TF families were found under high salt conditions. iii) Some genes up-regulated in response to salt treatment showed signs of adaptive evolution at the amino acid level and might contribute to adaptation to fluctuating intertidal environments.

CONCLUSIONS

This study first elucidates the mechanism of high-salt adaptation in mangroves at the whole-transcriptome level by salt gradient experimental treatments. It reveals that several candidate genes (including salt-related genes, TF-encoding genes, and PSGs) and major pathways are involved in adaptation to high-salt environments. Our study also provides a valuable resource for future investigation of adaptive evolution in extreme environments.

Transcriptomic and Metabolic Analyses Provide New Insights into the Apple Fruit Quality Decline during Long-Term Cold Storage

Long-term low-temperature conditioning (LT-LTC) decreases apple fruit quality, but the underlying physiological and molecular basis is relatively uncharacterized. We identified 12 clusters of differentially expressed genes (DEGs) involved in multiple biological processes (i.e., sugar, malic acid, fatty acid, lipid, complex phytohormone, and stress-response pathways). The expression levels of genes in sugar pathways were correlated with decreasing starch levels during LT-LTC. Specifically, starch-synthesis-related genes (e.g., BE, SBE, and GBSS genes) exhibited downregulated expression, whereas sucrose-metabolism-related gene expression levels were up- or downregulated. The expression levels of genes in the malic acid pathway (ALMT9, AATP1, and AHA2) were upregulated, as well as the content of malic acid in apple fruit during LT-LTC. A total of 151 metabolites, mainly related to amino acids and their isoforms, amines, organic acids, fatty acids, sugars, and polyols, were identified during LT-LTC. Additionally, 35 organic-acid-related metabolites grouped into three clusters, I (3), II (22), and III (10), increased in abundance during LT-LTC. Multiple phytohormones regulated the apple fruit chilling injury response. The ethylene (ET) and abscisic acid (ABA) levels increased at CS2 and CS3, and jasmonate (JA) levels also increased during LT-LTC. Furthermore, the expression levels of genes involved in ET, ABA, and JA synthesis and response pathways were upregulated. Finally, some key transcription factor genes (MYB, bHLH, ERF, NAC, and bZIP genes) related to the apple fruit cold acclimation response were differentially expressed. Our results suggest that the multilayered mechanism underlying apple fruit deterioration during LT-LTC is a complex, transcriptionally regulated process involving cell structures, sugars, lipids, hormones, and transcription factors.

Role of MPK4 in pathogen-associated molecular pattern-triggered alternative splicing in Arabidopsis

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied alternative splicing induced by the PAMP flagellin in Arabidopsis. We identified 506 PAMP-induced differentially alternatively spliced (DAS) genes. Importantly, of the 506 PAMP-induced DAS genes, only 89 overlap with the set of 1950 PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DAS events. Global DAS analysis of mpk3, mpk4, and mpk6 mutants in the absence of PAMP treatment showed no major splicing changes. However, in contrast to MPK3 and MPK6, MPK4 was found to be a key regulator of PAMP-induced DAS events as the AS of a number of splicing factors and immunity-related protein kinases is affected, such as the calcium-dependent protein kinase CPK28, the cysteine-rich receptor like kinases CRK13 and CRK29 or the FLS2 co-receptor SERK4/BKK1. Although MPK4 is guarded by SUMM2 and consequently, the mpk4 dwarf and DEG phenotypes are suppressed in mpk4 summ2 mutants, MPK4-dependent DAS is not suppressed by SUMM2, supporting the notion that PAMP-triggered MPK4 activation mediates regulation of alternative splicing.

Alternative splicing (AS) of pre-mRNAs in plants is an important mechanism of gene regulation in environmental stress tolerance but plant signals involved are essentially unknown. Pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) is mediated by mitogen-activated protein kinases and the majority of PTI defense genes are regulated by MPK3, MPK4 and MPK6. These responses have been mainly analyzed at the transcriptional level, however many splicing factors are direct targets of MAPKs. Here, we studied PAMP-induced alternative splicing in Arabidopsis and identified several hundred differentially alternatively spliced (DAS) genes. Importantly, of these PAMP-induced DAS genes, only 18% overlap with the set of PAMP-induced differentially expressed genes (DEG), indicating that transcriptome analysis does not identify most DAS events. Global DAS analysis of MAPK mutants identified MPK4 as a key regulator of PAMP-induced DAS events. Although MPK4 is guarded by SUMM2 and consequently, the mpk4 dwarf and DEG phenotypes are suppressed in mpk4 summ2 mutants, MPK4-dependent DAS is not suppressed by SUMM2, showing that PAMP-triggered MPK4 activation mediates regulation of alternative splicing.

Genomics, molecular and evolutionary perspective of NAC transcription factors

NAC (NAM, ATAF1,2, and CUC2) transcription factors are one of the largest transcription factor families found in the plants and are involved in diverse developmental and signalling events. Despite the availability of comprehensive genomic information from diverse plant species, the basic genomic, biochemical, and evolutionary details of NAC TFs have not been established. Therefore, NAC TFs family proteins from 160 plant species were analyzed in the current study. Study revealed, Brassica napus (410) encodes highest number and Klebsormidium flaccidum (3) encodes the lowest number of TFs. The study further revealed the presence of NAC TF in the Charophyte algae K. flaccidum. On average, the monocot plants encode higher number (141.20) of NAC TFs compared to the eudicots (125.04), gymnosperm (75), and bryophytes (22.66). Furthermore, our analysis revealed that several NAC TFs are membrane bound and contain monopartite, bipartite, and multipartite nuclear localization signals. NAC TFs were also found to encode several novel chimeric proteins and regulate a complex interactome network. In addition to the presence of NAC domain, several NAC proteins were found to encode other functional signature motifs as well. Relative expression analysis of NAC TFs in A. thaliana revealed root tissue treated with urea and ammonia showed higher level of expression and leaf tissues treated with urea showed lower level of expression. The synonymous codon usage is absent in the NAC TFs and it appears that they have evolved from orthologous ancestors and undergone vivid duplications to give rise to paralogous NAC TFs. The presence of novel chimeric NAC TFs are of particular interest and the presence of chimeric NAC domain with other functional signature motifs in the NAC TF might encode novel functional properties in the plants.

Identification and expression analysis of WRKY gene family under drought stress in peanut (Arachis hypogaea L.)

WRKY transcription factors play crucial roles in regulation mechanism leading to the adaption of plants to the complex environment. In this study, AhWRKY family was comprehensively analyzed using bioinformatic approaches in combination with transcriptome sequencing data of the drought-tolerant peanut variety ‘L422’. A total of 158 AhWRKY genes were identified and named according to their distribution on the chromosomes. Based on the structural features and phylogenetic analysis of AhWRKY proteins, the AhWRKY family members were classified into three (3) groups, of which group II included five (5) subgroups. Results of structure and conserved motifs analysis for the AhWRKY genes confirmed the accuracy of the clustering analysis. In addition, 12 tandem and 136 segmental duplication genes were identified. The results indicated that segmental duplication events were the main driving force in the evolution of AhWRKY family. Collinearity analysis found that 32 gene pairs existed between Arachis hypogaea and two diploid wild ancestors (Arachis duranensis and Arachis ipaensis), which provided valuable clues for phylogenetic characteristics of AhWRKY family. Furthermore, 19 stress-related cis-acting elements were found in the promoter regions. During the study of gene expression level of AhWRKY family members in response to drought stress, 73 differentially expressed AhWRKY genes were obtained to have been influenced by drought stress. These results provide fundamental insights for further study of WRKY genes in peanut drought resistance.

The barley stripe mosaic virus expression system reveals the wheat C2H2 zinc finger protein TaZFP1B as a key regulator of drought tolerance

BACKGROUND

Drought stress is one of the major factors limiting wheat production globally. Improving drought tolerance is important for agriculture sustainability. Although various morphological, physiological and biochemical responses associated with drought tolerance have been documented, the molecular mechanisms and regulatory genes that are needed to improve drought tolerance in crops require further investigation. We have used a novel 4-component version (for overexpression) and a 3-component version (for underexpression) of a barley stripe mosaic virus-based (BSMV) system for functional characterization of the C2H2-type zinc finger protein TaZFP1B in wheat. These expression systems avoid the need to produce transgenic plant lines and greatly speed up functional gene characterization.

RESULTS

We show that overexpression of TaZFP1B stimulates plant growth and up-regulates different oxidative stress-responsive genes under well-watered conditions. Plants that overexpress TaZFP1B are more drought tolerant at critical periods of the plant's life cycle. Furthermore, RNA-Seq analysis revealed that plants overexpressing TaZFP1B reprogram their transcriptome, resulting in physiological and physical modifications that help wheat to grow and survive under drought stress. In contrast, plants transformed to underexpress TaZFP1B are significantly less tolerant to drought and growth is negatively affected.

CONCLUSIONS

This study clearly shows that the two versions of the BSMV system can be used for fast and efficient functional characterization of genes in crops. The extent of transcriptome reprogramming in plants that overexpress TaZFP1B indicates that the encoded transcription factor is a key regulator of drought tolerance in wheat.

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

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

Comparative transcriptome analysis reveals gene expression differences between two peach cultivars under saline-alkaline stress

BACKGROUND

Saline-alkaline stress is a major abiotic stress that is harmful to plant growth worldwide. Two peach cultivars (GF677 and Maotao) display distinct phenotypes under saline-alkaline stress. The molecular mechanism explaining the differences between the two cultivars is still unclear.

RESULTS

In the present study, we systematically analysed the changes in GF677 and Maotao leaves upon saline-alkaline stress by using cytological and biochemical technologies as well as comparative transcriptome analysis. Transmission electron microscopy (TEM) observations showed that the structure of granum was dispersive in Maotao chloroplasts. The biochemical analysis revealed that POD activity and the contents of chlorophyll a and chlorophyll b, as well as iron, were notably decreased in Maotao. Comparative transcriptome analysis detected 881 genes with differential expression (including 294 upregulated and 587 downregulated) under the criteria of |log2 Ratio| ≥ 1 and FDR ≤0.01. Gene ontology (GO) analysis showed that all differentially expressed genes (DEGs) were grouped into 30 groups. MapMan annotation of DEGs showed that photosynthesis, antioxidation, ion metabolism, and WRKY TF were activated in GF677, while cell wall degradation, secondary metabolism, starch degradation, MYB TF, and bHLH TF were activated in Maotao. Several iron and stress-related TFs (ppa024966m, ppa010295m, ppa0271826m, ppa002645m, ppa010846m, ppa009439m, ppa008846m, and ppa007708m) were further discussed from a functional perspective based on the phylogenetic tree integration of other species homologues.

CONCLUSIONS

According to the cytological and molecular differences between the two cultivars, we suggest that the integrity of chloroplast structure and the activation of photosynthesis as well as stress-related genes are crucial for saline-alkaline resistance in GF677. The results presented in this report provide a theoretical basis for cloning saline-alkaline tolerance genes and molecular breeding to improve saline-alkaline tolerance in peach.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Knockout of the OsNAC006 Transcription Factor Causes Drought and Heat Sensitivity in Rice

Rice (Oryza sativa) responds to various abiotic stresses during growth. Plant-specific NAM, ATAF1/2, and CUC2 (NAC) transcription factors (TFs) play an important role in controlling numerous vital growth and developmental processes. To date, 170 NAC TFs have been reported in rice, but their roles remain largely unknown. Herein, we discovered that the TF OsNAC006 is constitutively expressed in rice, and regulated by H₂O₂, cold, heat, abscisic acid (ABA), indole-3-acetic acid (IAA), gibberellin (GA), NaCl, and polyethylene glycol (PEG) 6000 treatments. Furthermore, knockout of OsNAC006 using the CRISPR-Cas9 system resulted in drought and heat sensitivity. RNA sequencing (RNA-seq) transcriptome analysis revealed that OsNAC006 regulates the expression of genes mainly involved in response to stimuli, oxidoreductase activity, cofactor binding, and membrane-related pathways. Our findings elucidate the important role of OsNAC006 in drought responses, and provide valuable information for genetic manipulation to enhance stress tolerance in future plant breeding programs.

De novo transcriptome assembly and analysis of Phragmites karka, an invasive halophyte, to study the mechanism of salinity stress tolerance

With the rapidly deteriorating environmental conditions, the development of stress tolerant plants has become a priority for sustaining agricultural productivity. Therefore, studying the process of stress tolerance in naturally tolerant species hold significant promise. Phragmites karka is an invasive plant species found abundantly in tropical and sub tropical regions, fresh water regions and brackish marshy areas, such as river banks and lake shores. The plant possesses the ability to adapt and survive under conditions of high salinity. We subjected P. karka seedlings to salt stress and carried out whole transcriptome profiling of leaf and root tissues. Assessing the global transcriptome changes under salt stress resulted in the identification of several genes that are differentially regulated under stress conditions in root and leaf tissue. A total of 161,403 unigenes were assembled and used as a reference for digital gene expression analysis. A number of key metabolic pathways were found to be over-represented. Digital gene expression analysis was validated using qRT-PCR. In addition, a number of different transcription factor families including WRKY, MYB, CCCH, NAC etc. were differentially expressed under salinity stress. Our data will facilitate further characterisation of genes involved in salinity stress tolerance in P. karka. The DEGs from our results are potential candidates for understanding and engineering abiotic stress tolerance in plants.

Expression Analyses of Soybean VOZ Transcription Factors and the Role of GmVOZ1G in Drought and Salt Stress Tolerance

Vascular plant one-zinc-finger (VOZ) transcription factor, a plant specific one-zinc-finger-type transcriptional activator, is involved in regulating numerous biological processes such as floral induction and development, defense against pathogens, and response to multiple types of abiotic stress. Six VOZ transcription factor-encoding genes (GmVOZs) have been reported to exist in the soybean (Glycine max) genome. In spite of this, little information is currently available regarding GmVOZs. In this study, GmVOZs were cloned and characterized. GmVOZ genes encode proteins possessing transcriptional activation activity in yeast cells. GmVOZ1E, GmVOZ2B, and GmVOZ2D gene products were widely dispersed in the cytosol, while GmVOZ1G was primarily located in the nucleus. GmVOZs displayed a differential expression profile under dehydration, salt, and salicylic acid (SA) stress conditions. Among them, GmVOZ1G showed a significantly induced expression in response to all stress treatments. Overexpression of GmVOZ1G in soybean hairy roots resulted in a greater tolerance to drought and salt stress. In contrast, RNA interference (RNAi) soybean hairy roots suppressing GmVOZ1G were more sensitive to both of these stresses. Under drought treatment, soybean composite plants with an overexpression of hairy roots had higher relative water content (RWC). In response to drought and salt stress, lower malondialdehyde (MDA) accumulation and higher peroxidase (POD) and superoxide dismutase (SOD) activities were observed in soybean composite seedlings with an overexpression of hairy roots. The opposite results for each physiological parameter were obtained in RNAi lines. In conclusion, GmVOZ1G positively regulates drought and salt stress tolerance in soybean hairy roots. Our results will be valuable for the functional characterization of soybean VOZ transcription factors under abiotic stress.

Genome-wide identification and expression analysis of WRKY transcription factors in pearl millet (Pennisetum glaucum) under dehydration and salinity stress

Background

Plants have developed various sophisticated mechanisms to cope up with climate extremes and different stress conditions, especially by involving specific transcription factors (TFs). The members of the WRKY TF family are well known for their role in plant development, phytohormone signaling and developing resistance against biotic or abiotic stresses. In this study, we performed a genome-wide screening to identify and analyze the WRKY TFs in pearl millet (Pennisetum glaucum; PgWRKY), which is one of the most widely grown cereal crops in the semi-arid regions.

Results

A total number of 97 putative PgWRKY proteins were identified and classified into three major Groups (I-III) based on the presence of WRKY DNA binding domain and zinc-finger motif structures. Members of Group II have been further subdivided into five subgroups (IIa-IIe) based on the phylogenetic analysis. In-silico analysis of PgWRKYs revealed the presence of various cis-regulatory elements in their promoter region like ABRE, DRE, ERE, EIRE, Dof, AUXRR, G-box, etc., suggesting their probable involvement in growth, development and stress responses of pearl millet. Chromosomal mapping evidenced uneven distribution of identified 97 PgWRKY genes across all the seven chromosomes of pearl millet. Synteny analysis of PgWRKYs established their orthologous and paralogous relationship among the WRKY gene family of Arabidopsis thaliana, Oryza sativa and Setaria italica. Gene ontology (GO) annotation functionally categorized these PgWRKYs under cellular components, molecular functions and biological processes. Further, the differential expression pattern of PgWRKYs was noticed in different tissues (leaf, stem, root) and under both drought and salt stress conditions. The expression pattern of PgWRKY33, PgWRKY62 and PgWRKY65 indicates their probable involvement in both dehydration and salinity stress responses in pearl millet.

Conclusion

Functional characterization of identified PgWRKYs can be useful in delineating their role behind the natural stress tolerance of pearl millet against harsh environmental conditions. Further, these PgWRKYs can be employed in genome editing for millet crop improvement.

Trait analysis reveals DOG1 determines initial depth of seed dormancy, but not changes during dormancy cycling that result in seedling emergence timing

Seedling emergence timing is crucial in competitive plant communities and so contributes to species fitness. To understand the mechanistic basis of variation in seedling emergence timing, we exploited the contrasting behaviour of two Arabidopsis thaliana ecotypes: Cape Verde Islands (Cvi) and Burren (Bur-0). We used RNA-Seq analysis of RNA from exhumed seeds and quantitative trait loci (QTL) analyses on a mapping population from crossing the Cvi and Bur-0 ecotypes. We determined genome-wide expression patterns over an annual dormancy cycle in both ecotypes, identifying nine major clusters based on the seasonal timing of gene expression, and variation in behaviour between them. QTL were identified for depth of seed dormancy and seedling emergence timing (SET). Both analyses showed a key role for DOG1 in determining depth of dormancy, but did not support a direct role for DOG1 in generating altered seasonal patterns of seedling emergence. The principle QTL determining SET (SET1: dormancy cycling) is physically close on chromosome 5, but is distinct from DOG1. We show that SET1 and two other SET QTLs each contain a candidate gene (AHG1, ANAC060, PDF1 respectively) closely associated with DOG1 and abscisic acid signalling and suggest a model for the control of SET in the field.

Transcriptome -wide modulation combined with morpho-physiological analyses of Typha orientalis roots in response to lead challenge

Lead (Pb) is a common pollutant in many environments, including in the soil, water, and/or air. Typha orientalis Presl, a large emergent aquatic plant, has been reported to function as a Pb-tolerant and Pb-accumulating plant; however, very little molecular information regarding the tolerance of T. orientalis towards Pb is known. In this study, Pb accumulation and key factors involved in the Pb stress response at different Pb concentrations were investigated. Pb was primarily accumulated in the roots and was mainly located in the cell wall and membrane systems. Differentially expressed genes (DEGs) were identified in T. orientalis roots after Pb exposure via RNA-seq analyses. In the 0.10 mM and 0.25 mM Pb²⁺-treated groups, a total of 3275 DEGs were detected relative to the control. Many of these genes were associated with oxidation-reduction processes, metal transport, protein kinase/phosphorylation, and DNA binding transcription factors, which were shown to be Pb-responsive DEGs. Mapping Kyoto Encyclopedia of Genes and Genomes (KEGG) database, "phenylpropanoid biosynthesis" was analyzed as the major pathway of the important modules of overlapping DEGs of 0.10 mM and 0.25 mM Pb²⁺ treatments. Furthermore, a lead response gene named ToLR1 with unknown function was of particular interest. The full-length of ToLR1 sequence was cloned using rapid amplification of cDNA ends (RACE) and overexpressed in Arabidopsis thaliana, which resulted in enhanced resistance to Pb stress. This is the first report providing genomic information detailing Pb responsive genes in T. orientalis. Moreover, this study provides novel insights into the molecular mechanisms underlying the response of T. orientalis and other accumulators towards Pb stress. The key genes identified in this study may serve as potential targets for genetic engineering targeting phytoremediation.

Comparative transcriptome analysis uncovers regulatory roles of long non-coding RNAs involved in resistance to powdery mildew in melon

Background

Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs with more than 200 nucleotides in length, which play vital roles in a wide range of biological processes. Powdery mildew disease (PM) has become a major threat to the production of melon. To investigate the potential roles of lncRNAs in resisting to PM in melon, it is necessary to identify lncRNAs and uncover their molecular functions. In this study, we compared the lncRNAs between a resistant and a susceptible melon in response to PM infection.

Results

It is reported that 11,612 lncRNAs were discovered, which were distributed across all 12 melon chromosomes, and > 85% were from intergenic regions. The melon lncRNAs have shorter transcript lengths and fewer exon numbers than protein-coding genes. In addition, a total of 407 and 611 lncRNAs were found to be differentially expressed after PM infection in PM-susceptible and PM-resistant melons, respectively. Furthermore, 1232 putative targets of differently expressed lncRNAs (DELs) were discovered and gene ontology enrichment (GO) analysis showed that these target genes were mainly enriched in stress-related terms. Consequently, co-expression patterns between LNC_018800 and CmWRKY21, LNC_018062 and MELO3C015771 (glutathione reductase coding gene), LNC_014937 and CmMLO5 were confirmed by qRT-PCR. Moreover, we also identified 24 lncRNAs that act as microRNA (miRNA) precursors, 43 lncRNAs as potential targets of 22 miRNA families and 13 lncRNAs as endogenous target mimics (eTMs) for 11 miRNAs.

Conclusion

This study shows the first characterization of lncRNAs involved in PM resistance in melon and provides a starting point for further investigation into the functions and regulatory mechanisms of lncRNAs in the resistance to PM.

Transcriptome-Wide Identification, Evolutionary Analysis, and GA Stress Response of the GRAS Gene Family in Panax ginseng C. A. Meyer

GRAS transcription factors are a kind of plant-specific transcription factor that have been found in a variety of plants. According to previous studies, GRAS proteins are widely involved in the physiological processes of plant signal transduction, stress, growth and development. The Jilin ginseng (Panax ginseng C.A. Meyer) is a heterogeneous tetraploid perennial herb of the Araliaceae family, ginseng genus. Important information regarding the GRAS transcription factors has not been reported in ginseng. In this study, 59 Panax ginseng GRAS (PgGRAS) genes were obtained from the Jilin ginseng transcriptome data and divided into 13 sub-families according to the classification of Arabidopsis thaliana. Through systematic evolution, structural variation, function and gene expression analysis, we further reveal GRAS’s potential function in plant growth processes and its stress response. The expression of PgGRAS genes responding to gibberellin acids (GAs) suggests that these genes could be activated after application concentration of GA. The qPCR analysis result shows that four PgGRAS genes belonging to the DELLA sub-family potentially have important roles in the GA stress response of ginseng hairy roots. This study provides not only a preliminary exploration of the potential functions of the GRAS genes in ginseng, but also valuable data for further exploration of the candidate PgGRAS genes of GA signaling in Jilin ginseng, especially their roles in ginseng hairy root development and GA stress response.

ABA signaling is negatively regulated by GbWRKY1 through JAZ1 and ABI1 to affect salt and drought tolerance

GbWRKY1 can function as a negative regulator of ABA signaling via JAZ1 and ABI1, with effects on salt and drought tolerance. WRKY transcription factors play important roles in plant development and stress responses. GbWRKY1 was initially identified as a defense-related gene in cotton and negatively regulates the response to fungal pathogens by activating the expression of JAZ1. Here, we characterized the role of GbWRKY1, an orthologue of the Arabidopsis gene AtWRKY75, in abiotic stress (salt and drought) and established novel connection between JAZ1 and ABA signaling in Arabidopsis. GbWRKY1 is nucleus localized and its expression is significantly induced by treatment with ABA and osmotic stresses NaCl and PEG. Increased levels of expression of GbWRKY1 in transgenic Arabidopsis enhance sensitivity to salt and drought as revealed by seed germination tests and soil stress experiments. Similarly, GbWRKY1 overexpression cotton plants also display increased sensitivity to PEG treatment and drought. Expression analysis shows that the induction of two ABA responsive genes, RAB18 and RD29A by NaCl, mannitol, and ABA treatment is significantly impaired in GbWRKY1 overexpression Arabidopsis lines. GbWRKY1 overexpression Arabidopsis displays a strong ABA-insensitive phenotype at both germination and early stages of seedling development. Further genetic evidence suggested that the ABA-insensitive phenotype of GbWRKY1 overexpression Arabidopsis was dependent on JAZ1, and overexpression of JAZ1 also displayed an ABA-insensitive phenotype. In addition, yeast two hybrid and bimolecular fluorescence complementation assays showed that JAZ1 directly interacts with ABI1, a key negative regulator of ABA signaling. We, therefore, demonstrate that GbWRKY1 acts as a negative regulator of ABA signaling, through an interaction network involving JAZ1 and ABI1, to regulate salt and drought tolerance.

Salicylic acid and nitric oxide signaling in plant heat stress

In agriculture, heat stress (HS) has become one of the eminent abiotic threats to crop growth, productivity and nutritional security because of the continuous increase in global mean temperature. Studies have annotated that the heat stress response (HSR) in plants is highly conserved, involving complex regulatory networks of various signaling and sensor molecules. In this context, the ubiquitous-signaling molecules salicylic acid (SA) and nitric oxide (NO) have diverted the attention of the plant science community because of their putative roles in plant abiotic and biotic stress tolerance. However, their involvement in the transcriptional regulatory networks in plant HS tolerance is still poorly understood. In this review, we have conceptualized current knowledge concerning how SA and NO sense HS in plants and how they trigger the HSR leading to the activation of transcriptional-signaling cascades. Fundamentals of functional components and signaling networks associated with molecular mechanisms involved in SA/NO-mediated HSR in plants have also been discussed. Increasing evidences have suggested the involvement of epigenetic modifications in the development of a 'stress memory', thereby provoking the role of epigenetic mechanisms in the regulation of plant's innate immunity under HS. Thus, we have also explored the recent advancements regarding the biological mechanisms and the underlying significance of epigenetic regulations involved in the activation of HS responsive genes and transcription factors by providing conceptual frameworks for understanding molecular mechanisms behind the 'transcriptional stress memory' as potential memory tools in the regulation of plant HSR.

The draft genome of a wild barley genotype reveals its enrichment in genes related to biotic and abiotic stresses compared to cultivated barley

Wild barley (Hordeum spontaneum) is the progenitor of cultivated barley (Hordeum vulgare) and provides a rich source of genetic variations for barley improvement. Currently, the genome sequences of wild barley and its differences with cultivated barley remain unclear. In this study, we report a high-quality draft assembly of wild barley accession (AWCS276; henceforth named as WB1), which consists of 4.28 Gb genome and 36 395 high-confidence protein-coding genes. BUSCO analysis revealed that the assembly included full lengths of 95.3% of the 956 single-copy plant genes, illustrating that the gene-containing regions have been well assembled. By comparing with the genome of the cultivated genotype Morex, it is inferred that the WB1 genome contains more genes involved in resistance and tolerance to biotic and abiotic stresses. The presence of the numerous WB1-specific genes indicates that, in addition to enhance allele diversity for genes already existing in the cultigen, exploiting the wild barley taxon in breeding should also allow the incorporation of novel genes. Furthermore, high levels of genetic variation in the pericentromeric regions were detected in chromosomes 3H and 5H between the wild and cultivated genotypes, which may be the results of domestication. This H. spontaneum draft genome assembly will help to accelerate wild barley research and be an invaluable resource for barley improvement and comparative genomics research.

Transcriptome analysis of potato shoots, roots and stolons under nitrogen stress

Potato crop requires high dose of nitrogen (N) to produce high tuber yield. Excessive application of N causes environmental pollution and increases cost of production. Hence, knowledge about genes and regulatory elements is essential to strengthen research on N metabolism in this crop. In this study, we analysed transcriptomes (RNA-seq) in potato tissues (shoot, root and stolon) collected from plants grown in aeroponic culture under controlled conditions with varied N supplies i.e. low N (0.2 milli molar N) and high N (4 milli molar N). High quality data ranging between 3.25 to 4.93 Gb per sample were generated using Illumina NextSeq500 that resulted in 83.60-86.50% mapping of the reads to the reference potato genome. Differentially expressed genes (DEGs) were observed in the tissues based on statistically significance (p ≤ 0.05) and up-regulation with ≥ 2 log2 fold change (FC) and down-regulation with ≤ -2 log2 FC values. In shoots, of total 19730 DEGs, 761 up-regulated and 280 down-regulated significant DEGs were identified. Of total 20736 DEGs in roots, 572 (up-regulated) and 292 (down-regulated) were significant DEGs. In stolons, of total 21494 DEG, 688 and 230 DEGs were significantly up-regulated and down-regulated, respectively. Venn diagram analysis showed tissue specific and common genes. The DEGs were functionally assigned with the GO terms, in which molecular function domain was predominant in all the tissues. Further, DEGs were classified into 24 KEGG pathways, in which 5385, 5572 and 5594 DEGs were annotated in shoots, roots and stolons, respectively. The RT-qPCR analysis validated gene expression of RNA-seq data for selected genes. We identified a few potential DEGs responsive to N deficiency in potato such as glutaredoxin, Myb-like DNA-binding protein, WRKY transcription factor 16 and FLOWERING LOCUS T in shoots; high-affinity nitrate transporter, protein phosphatase-2c, glutaredoxin family protein, malate synthase, CLE7, 2-oxoglutarate-dependent dioxygenase and transcription factor in roots; and glucose-6-phosphate/phosphate translocator 2, BTB/POZ domain-containing protein, F-box family protein and aquaporin TIP1;3 in stolons, and many genes of unknown function. Our study highlights that these potential genes play very crucial roles in N stress tolerance, which could be useful in augmenting research on N metabolism in potato.

Arabidopsis thaliana WRKY25 Transcription Factor Mediates Oxidative Stress Tolerance and Regulates Senescence in a Redox-Dependent Manner

Senescence is the last developmental step in plant life and is accompanied by a massive change in gene expression implying a strong participation of transcriptional regulators. In the past decade, the WRKY53 transcription factor was disclosed to be a central node of a complex regulatory network of leaf senescence and to underlie a tight multi-layer control of expression, activity and protein stability. Here, we identify WRKY25 as a redox-sensitive up-stream regulatory factor of WRKY53 expression. Under non-oxidizing conditions, WRKY25 binds to a specific W-box in the WRKY53 promoter and acts as a positive regulator of WRKY53 expression in a transient expression system using Arabidopsis protoplasts, whereas oxidizing conditions dampened the action of WRKY25. However, overexpression of WRKY25 did not accelerate senescence but increased lifespan of Arabidopsis plants, whereas the knock-out of the gene resulted in the opposite phenotype, indicating a more complex regulatory function of WRKY25 within the WRKY subnetwork of senescence regulation. In addition, overexpression of WRKY25 mediated higher tolerance to oxidative stress and the intracellular H2O2 level is lower in WRKY25 overexpressing plants and higher in wrky25 mutants compared to wildtype plants suggesting that WRKY25 is also involved in controlling intracellular redox conditions. Consistently, WRKY25 overexpressers had higher and wrky mutants lower H2O2 scavenging capacity. Like already shown for WRKY53, MEKK1 positively influenced the activation potential of WRKY25 on the WRKY53 promoter. Taken together, WRKY53, WRKY25, MEKK1 and H2O2 interplay with each other in a complex network. As H2O2 signaling molecule participates in many stress responses, WRKK25 acts most likely as integrators of environmental signals into senescence regulation.

Comparative Analysis between Wild and Cultivated Cucumbers Reveals Transcriptional Changes during Domestication Process

The cultivated cucumber (Cucumis sativus L.) was reported to have been developed from a wild cucumber (Cucumis hystrix Chakrav.), nevertheless, these two organisms exhibit noteworthy differences. For example, the wild cucumber is known for its high resistance to different biotic and abiotic stresses. Moreover, the leaves and fruits of the wild cucumber have a bitter taste compared to the cultivated cucumber. These differences could be attributed mainly to the differences in gene expression levels. In the present investigation, we analyzed the RNA-sequencing data to show the differentially expressed genes (DEGs) between the wild and cultivated cucumbers. The identified DEGs were further utilized for Gene Ontology (GO) and pathway enrichment analysis and for identification of transcription factors and regulators. In the results, several enriched GO terms in the biological process, cellular component, and molecular functions categories were identified and various enriched pathways, especially the biosynthesis pathways of secondary products were recognized. Plant-specific transcription factor families were differentially expressed between the wild and cultivated cucumbers. The results obtained provide preliminary evidence for the transcriptional differences between the wild and cultivated cucumbers which developed during the domestication process as a result of natural and/or artificial selection, and they formulate the basis for future genetic research and improvement of the cultivated cucumber.

Sugar Transporter ZjSWEET2.2 Mediates Sugar Loading in Leaves of Ziziphus jujuba Mill

In plants, sugar transporters play an important role in the allocation of sugars from cells in source organs to cells in sink organs. Hence, an understanding of the molecular basis and regulation of assimilate partitioning by sugar transporters is essential. Leaves are the main source of photosynthetic products. In jujube (Ziziphus jujuba Mill.), the mechanisms regulating initial sugar unloading in leaves are still unclear. In this study, an expression profiling analysis showed that ZjSWEET2.2, encoding a sugar transporter in the SWEET family, is highly expressed in leaves. Over-expression of ZjSWEET2.2 increased carbon fixation in photosynthetic organs. Our analyses showed that ZjSWEET2.2 encodes a plasma membrane-localized sugar transporter protein. Its expression levels were found to be suppressed under drought stress and by high concentrations of exogenous sugars, but increased by low concentrations of exogenous sugars. Finally, DNA sequence analyses revealed several cis-elements related to sugar signaling in the promoter of ZjSWEET2.2. Together, these results suggest that ZjSWEET2.2 functions to mediate photosynthesis by exporting sugars from photosynthetic cells in the leaves, and its gene expression is regulated by sugar signals.

Salicylic acid application alleviates the adverse effects of triclosan stress in tobacco plants through the improvement of plant photosynthesis and enhancing antioxidant system

Triclosan (TCS) is a chlorophenol which is highly bacteriostatic and used in a wide array of consumer products. TCS is now one of the most commonly detected organic pollutants in the sewage sludges. The sludge utilization for fertilizers on agricultural land would pose the risk of causing adverse effects on plant growth and yield by TCS. However, the toxicity of TCS toward plants is comparatively less understood. In this study, we assessed the effects of TCS on tobacco plants which were grown in MS medium or soils containing various concentrations of TCS. Our results indicated that TCS at the concentration of 2 mg/L could strongly inhibit the tobacco seed germination. TCS could suppress tobacco plant growth in soil with different concentrations (10, 20, and 50 mg/kg) of TCS through the downregulation of chlorophyll contents, restricting photosynthesis and increasing generation of reactive oxygen species (ROS). Salicylic acid (SA) plays important roles in the stress response of plants. The role of exogenous SA application in protecting tobacco plants from TCS stress was also investigated in this study. SA application could significantly increase net photosynthesis, enhance antioxidant enzyme activity, and thereby enhancing tobacco plant tolerance to TCS. Moreover, the activation of MPK3 and MPK6 induced by TCS was downregulated in plants with the treatment of SA. It was thus referred that mitogen-activated protein kinases (MAPKs) might play a key role in the signal transduction of TCS stress, and this process might be regulated by SA signaling. Overall, our results demonstrated that TCS had negative impacts on tobacco plants and SA played a protective role on tobacco plants against TCS stress.

The WRKY transcription factor WRKY8 promotes resistance to pathogen infection and mediates drought and salt stress tolerance in Solanum lycopersicum

WRKY transcription factors play a key role in the tolerance of biotic and abiotic stresses across various crop species, but the function of some WRKY genes, particularly in tomato, remains unexplored. Here, we characterize the roles of a previously unstudied WRKY gene, SlWRKY8, in the resistance to pathogen infection and the tolerance to drought and salt stresses. Expression of SlWRKY8 was up-regulated upon Pseudomonas syringae pv. tomato DC3000 (Pst. DC3000), abiotic stresses such as drought, salt and cold, as well as ABA and SA treatments. The SlWRKY8 protein was localized to the nucleus with no transcription activation in yeast, but it could activate W-box-dependent transcription in plants. The overexpression of SlWRKY8 in tomato conferred a greater resistance to the pathogen Pst. DC3000 and resulted in the increased transcription levels of two pathogen-related genes SlPR1a1 and SlPR7. Moreover, transgenic plants displayed the alleviated wilting or chlorosis phenotype under drought and salt stresses, with higher levels of stress-induced osmotic substances like proline and higher transcript levels of the stress-responsive genes SlAREB, SlDREB2A and SlRD29. Stomatal aperature was smaller under drought stress in transgenic plants, maintaining higher water content in leaves compared with wild-type plants. The oxidative pressure, indicated by the concentration of hydrogen peroxide (H₂ O₂ ) and malondialdehyde (MDA), was also reduced in transgenic plants, where we also observed higher levels of antioxidant enzyme activities under stress. Overall, our results suggest that SlWRKY8 functions as a positive regulator in plant immunity against pathogen infection as well as in plant responses to drought and salt stresses.

GhWRKY70D13 Regulates Resistance to Verticillium dahliae in Cotton Through the Ethylene and Jasmonic Acid Signaling Pathways

Verticillium wilt caused by Verticillium dahliae is a destructive cotton disease causing severe yield and quality losses worldwide. WRKY transcription factors play important roles in plant defense against pathogen infection. However, little has been reported on the functions of WRKYs in cotton's resistance to V. dahliae. Here, we identified 5, 5, and 10 WRKY70 genes in Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum, respectively, and investigated the expression profiles of all GhWRKY70 genes in various cotton tissues and in response to hormone treatment or V. dahliae infection. Reverse transcription-quantitative PCR analysis showed that GhWRKY70D13 was expressed higher in roots and stems than in other tissues, and up-regulated after V. dahliae inoculation. Knock-down of GhWRKY70D13 improved resistance to V. dahliae in both resistant and susceptible cotton cultivars. Comparative analysis of transcriptomes generated from wild-type and stable RNAi (RNA interference) plant with down-regulated GhWRKY70D13 showed that genes involved in ethylene (ET) and jasmonic acid (JA) biosynthesis and signaling were significantly upregulated in the GhWRKY70D13 RNAi plants. Consistently, the contents of 1-aminocyclopropane-1-carboxylic (ACC), JA, and JA-isoleucine levels were significantly higher in the GhWRKY70D13 RNAi plants than in wild-type. Following V. dahliae infection, the levels of ACC and JA decreased in the GhWRKY70D13 RNAi plants but still significantly higher (for ACC) than that in wild-type or at the same level (for JA) as in non-infected wild-type plants. Collectively, our results suggested that GhWRKY70D13 negatively regulates cotton's resistance to V. dahliae mainly through its effect on ET and JA biosynthesis and signaling pathways.

The Soybean GmNAC019 Transcription Factor Mediates Drought Tolerance in Arabidopsis in an Abscisic Acid-Dependent Manner

Being master regulators of gene expression, transcription factors (TFs) play important roles in determining plant growth, development and reproduction. To date, many TFs have been shown to positively mediate plant responses to environmental stresses. In the current study, the biological functions of a stress-responsive NAC [NAM (No Apical Meristem), ATAF1/2 (Arabidopsis Transcription Activation Factor1/2), CUC2 (Cup-shaped Cotyledon2)]-TF encoding gene isolated from soybean (GmNAC019) in relation to plant drought tolerance and abscisic acid (ABA) responses were investigated. By using a heterologous transgenic system, we revealed that transgenic Arabidopsis plants constitutively expressing the GmNAC019 gene exhibited higher survival rates in a soil-drying assay, which was associated with lower water loss rate in detached leaves, lower cellular hydrogen peroxide content and stronger antioxidant defense under water-stressed conditions. Additionally, the exogenous treatment of transgenic plants with ABA showed their hypersensitivity to this phytohormone, exhibiting lower rates of seed germination and green cotyledons. Taken together, these findings demonstrated that GmNAC019 functions as a positive regulator of ABA-mediated plant response to drought, and thus, it has potential utility for improving plant tolerance through molecular biotechnology.

Arabidopsis WRKY53, a Node of Multi-Layer Regulation in the Network of Senescence

Leaf senescence is an integral part of plant development aiming at the remobilization of nutrients and minerals out of the senescing tissue into developing parts of the plant. Sequential as well as monocarpic senescence maximize the usage of nitrogen, mineral, and carbon resources for plant growth and the sake of the next generation. However, stress-induced premature senescence functions as an exit strategy to guarantee offspring under long-lasting unfavorable conditions. In order to coordinate this complex developmental program with all kinds of environmental input signals, complex regulatory cues have to be in place. Major changes in the transcriptome imply important roles for transcription factors. Among all transcription factor families in plants, the NAC and WRKY factors appear to play central roles in senescence regulation. In this review, we summarize the current knowledge on the role of WRKY factors with a special focus on WRKY53. In contrast to a holistic multi-omics view we want to exemplify the complexity of the network structure by summarizing the multilayer regulation of WRKY53 of Arabidopsis.

Genome-wide identification, characterisation and functional evaluation of WRKY genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.) G. Don. under abiotic stresses

BACKGROUND

WRKY DNA-binding protein (WRKY) is a large gene family involved in plant responses and adaptation to salt, drought, cold and heat stresses. Sweet potato from the genus Ipomoea is a staple food crop, but the WRKY genes in Ipomoea species remain unknown to date. Hence, we carried out a genome-wide analysis of WRKYs in Ipomoea trifida (H.B.K.) G. Don., the wild ancestor of sweet potato.

RESULTS

A total of 83 WRKY genes encoding 96 proteins were identified in I. trifida, and their gene distribution, duplication, structure, phylogeny and expression patterns were studied. ItfWRKYs were distributed on 15 chromosomes of I. trifida. Gene duplication analysis showed that segmental duplication played an important role in the WRKY gene family expansion in I. trifida. Gene structure analysis showed that the intron-exon model of the ItfWRKY gene was highly conserved. Meanwhile, the ItfWRKYs were divided into five groups (I, IIa + IIb, IIc, IId + IIe and III) on the basis of the phylogenetic analysis on I. trifida and Arabidopsis thaliana WRKY proteins. In addition, gene expression profiles confirmed by quantitative polymerase chain reaction showed that ItfWRKYs were highly up-regulated or down-regulated under salt, drought, cold and heat stress conditions, implying that these genes play important roles in response and adaptation to abiotic stresses.

CONCLUSIONS

In summary, genome-wide identification, gene structure, phylogeny and expression analysis of WRKY gene in I. trifida provide basic information for further functional studies of ItfWRKYs and for the molecular breeding of sweet potato.