Regulation of specialized metabolism by WRKY transcription factors

De novo transcriptome analysis unravels tissue-specific expression of candidate genes involved in major secondary metabolite biosynthetic pathways of Plumbago zeylanica: implication for pharmacological potential

KEY MESSAGE

The present study provides comparative transcriptome analysis, besides identifying functional secondary metabolite genes of Plumbago zeylanica with pharmacological potential for future functional genomics, and metabolomic engineering of secondary metabolites from this plant towards diversified biomedical applications.

ABSTRACT

Plumbago zeylanica is a widely used medicinal plant of the traditional Indian system of medicine with wide pharmacological potential to treat several disorders. The present study aimed to carry out comparative transcriptome analysis in leaf and root tissue of P. zeylanica using Illumina paired end sequencing to identify tissue-specific functional genes involved in the biosynthesis of secondary metabolites, contributing to its therapeutic efficacy. De novo sequencing assembly resulted in the identification of 62,321 "Unigenes" transcripts with an average size of 1325 bp. Functional annotation using BLAST2GO resulted in the identification of 50,301 annotated transcripts (80.71%) and GO assigned to 18,814 transcripts. KEGG pathway annotation of the "Unigenes" revealed that 2465 transcripts could be assigned to 242 KEGG pathway maps wherein the number of transcripts involved in secondary metabolism was distinct in root and leaf transcriptome. Among the secondary metabolite biosynthesis pathways, the cluster of "Unigenes" encoding enzymes of 'Phenylpropanoid biosynthesis pathway' represents the largest group (84 transcripts) followed by 'Terpenoid Backbone biosynthesis' (48 transcripts). The transcript levels of the candidate unigenes encoding key enzymes of phenylpropanoid (PAL, TAL) and flavanoid biosynthesis (CHS, ANS, FLS) pathways were up-regulated in root, while the expression levels of candidate "Unigenes" transcript for monoterpenoid (DXS, ISPF), diterpenoid biosynthesis (SPS, SDS) and indole alkaloid pathways (STR) were significantly higher in leaf of P. zeylanica. Interestingly, validation of differential gene expression profile by qRT-PCR also confirmed that candidate "Unigenes" enzymes of phenylpropanoid and flavonoid biosynthesis were highly expressed in the root, while the key regulatory enzymes of terpenoid and indole alkaloid compounds were up-regulated in the leaf, suggesting that (differences in) the levels of these functional genes could be attributed to the (differential) pharmacological activity (between root and leaf) in tissues of P. zeylanica.

Phylogenetic and Transcription Analysis of Hibiscus hamabo Sieb. et Zucc. WRKY Transcription Factors

WRKY transcription factors are known to play important roles in the regulation of various aspects of plant growth and development, including germination, stress resistance, and senescence. Nevertheless, there is little information about the WRKY genes in Hibiscus hamabo Sieb. et Zucc., an important semimangrove plant. In this study, HhWRKY genes in H. hamabo were identificated based on Illumina RNA-sequencing and isoform sequencing from salt-treated roots. Then phylogenetic analysis and conserved motif analysis of the WRKY family in H. hamabo and Arabidopsis thaliana were used to classify WRKY genes. Sixteen HhWRKY genes were selected from different groups to detect their expression patterns using real-time quantitative PCR in different organ (root, old leaf, tender leaf, receptacle, petal, or stamen) from 10-year-old H. hamabo plants grown in their natural environment and in seedlings with 8 to 10 true leaves challenged by phytohormone (salicylic acid, methyl jasmonate, or abscisic acid) and abiotic stress (drought, salt, or high temperature). As a result, the identified 78 HhWRKY genes were divided into two major groups and several subgroups based on their structural and phylogenetic features. Most transcripts of the selected 16 HhWRKY genes were more abundant in old than in tender leaves of H. hamabo. HhWRKY genes were regulated in reaction to abiotic stresses and phytohormone treatments and may participate in signaling networks to improve plant stress resistance. Some of HhWRKY genes behaved as would be predicted based on their homology with A. thaliana WRKY genes, but others showed divergent behavior. This systematic analysis lays the foundation for further identification of WRKY gene functions, with the aim of improving woody plants.

AsHSP26.8a, a creeping bentgrass small heat shock protein integrates different signaling pathways to modulate plant abiotic stress response

BACKGROUND

Small heat shock proteins (sHSPs) are critical for plant response to biotic and abiotic stresses, especially heat stress. They have also been implicated in various aspects of plant development. However, the acting mechanisms of the sHSPs in plants, especially in perennial grass species, remain largely elusive.

RESULTS

In this study, AsHSP26.8a, a novel chloroplast-localized sHSP gene from creeping bentgrass (Agrostis stolonifera L.) was cloned and its role in plant response to environmental stress was studied. AsHSP26.8a encodes a protein of 26.8 kDa. Its expression was strongly induced in both leaf and root tissues by heat stress. Transgenic Arabidopsis plants overexpressing AsHSP26.8a displayed reduced tolerance to heat stress. Furthermore, overexpression of AsHSP26.8a resulted in hypersensitivity to hormone ABA and salinity stress. Global gene expression analysis revealed AsHSP26.8a-modulated expression of heat-shock transcription factor gene, and the involvement of AsHSP26.8a in ABA-dependent and -independent as well as other stress signaling pathways.

CONCLUSIONS

Our results suggest that AsHSP26.8a may negatively regulate plant response to various abiotic stresses through modulating ABA and other stress signaling pathways.

Genome-wide investigation of WRKY transcription factors in sweet osmanthus and their potential regulation of aroma synthesis

WRKY transcription factors, one of the largest transcription factor families, play important roles in regulating the synthesis of secondary metabolites. In sweet osmanthus (Osmanthus fragrans), the monoterpenes have been demonstrated as the most important volatile compounds, and the W-box, which is the cognate binding site of WRKY transcription factors, could be identified in most of the terpene-synthesis-related genes' promoters. However, the role of the WRKY family in terpene synthesis in sweet osmanthus has rarely been examined. In this study, 154 WRKY genes with conserved WRKY domain were identified and classified into three groups. The group II was further divided into five subgroups, and almost all members of IId contained a plant zinc cluster domain. Eight OfWRKYs (OfWRKY7/19/36/38/42/84/95/139) were screened from 20 OfWRKYs for their flower-specific expression patterns in different tissues. Simultaneously, the expression patterns of OfWRKYs and emission patterns of volatile compounds during the flowering process were determined and gas chromatography-mass spectrometry results showed that monoterpenes, such as linalool and ocimene, accounted for the highest proportion, contributing to the floral scent of sweet osmanthus in two cultivars. In addition, correlation analysis revealed the expression patterns of OfWRKYs (OfWRKY7/19/36/139) were each correlated with distinct monoterpenes (linalool, linalool derivatives, ocimene and ocimene derivatives). Subcellular localization analysis showed that p35S::GFP-OfWRKY7/38/95/139 were localized in the nucleus and OfWRKY139 had very strong transactivation activity. Collectively, the results indicated potential roles of OfWRKY139 and OfWRKYs with plant zinc cluster domain in regulating synthesis of aromatic compounds in sweet osmanthus, laying the foundation for use of OfWRKYs to improve the aroma of ornamental plants.

Physiological, Metabolic, and Transcriptomic Analyses Reveal the Responses of Arabidopsis Seedlings to Carbon Nanohorns

Carbon-based nanomaterials have potential applications in nanoenabled agriculture. However, the physiological and molecular mechanisms underlying single-walled carbon nanohorn (SWCNH)-mediated plant growth remain unclear. Here, we investigated the effects of SWCNHs on Arabidopsis grown in ¹/₄-strength Murashige and Skoog medium via physiological, genetic, and molecular analyses. Treatment with 0.1 mg/L SWCNHs promoted primary root (PR) growth and lateral root (LR) formation; 50 and 100 mg/L SWCNHs inhibited PR growth. Treatment with 0.1 mg/L SWCNHs increased the lengths of the meristematic and elongation zones, and transcriptomic and genetic analyses confirmed the positive effects of SWCNHs on root tip stem cell niche activity and meristematic cell division potential. Increased expression of YUC3 and YUC5 and increased PIN2 abundance improved PR growth and LR development in 0.1 mg/L SWCNH-treated seedlings. Metabolomic analyses revealed that SWCNHs altered the levels of sugars, amino acids, and organic acids, suggesting that SWCNHs reprogrammed carbon/nitrogen metabolism in plants. SWCNHs also regulate plant growth and development by increasing the levels of several secondary metabolites; transcriptomic analyses further supported these results. The present results are valuable for continued use of SWCNHs in agri-nanotechnology, and these molecular approaches could serve as examples for studies on the effects of nanomaterials in plants.

Metabolic Profiling of Varronia curassavica Jacq. Terpenoids by Flow Modulated Two-Dimensional Gas Chromatography Coupled to Mass Spectrometry

In this study, a metabolomic approach was used to investigate the effect of seasonality on the chemical composition and yield of anti-inflammatory active principle, α-humulene, in the essential oil of three genotypes of Varronia curassavica Jacq. (Syn. Cordia verbenaceae). The essential oils were extracted by hydrodistillation and analyzed by comprehensive two-dimensional gas chromatography coupled to mass spectrometry (GC×GC-MS). The GC×GC approach a three-fold improvement in qualitative analysis (48 compounds were identified by GC-MS versus 135 by GC×GC-MS). The improved resolving power of GC×GC resolved important coelutions and enabled the detection of unusual substances in V. curassavica essential oil. The chromatographic data was analyzed by using peak table-based chemometrics, namely, principal component analysis (PCA) and hierarchical cluster analysis (HCA). The metabolic study showed that seasonality has a significant effect on the chemical composition. The α-humulene content was affected by genotype and season. Spring and summer were the best harvest seasons for the yield of the active ingredient, found in higher concentrations in the VC2 genotype. The proposed metabolomic workflow was successfully applied to terpene analysis found in V. curassavica essential oil, and such results have broadened our understanding of the influence of seasonal factors on the specialized metabolism of the species.

Signal transduction networks during stress combination

Episodes of heat waves combined with drought can have a devastating impact on agricultural production worldwide. These conditions, as well as many other types of stress combinations, impose unique physiological and developmental demands on plants and require the activation of dedicated pathways. Here, we review recent RNA sequencing studies of stress combination in plants, and conduct a meta-analysis of the transcriptome response of plants to different types of stress combination. Our analysis reveals that each different stress combination is accompanied by its own set of stress combination-specific transcripts, and that the response of different transcription factor families is unique to each stress combination. The alarming rate of increase in global temperatures, coupled with the predicted increase in future episodes of extreme weather, highlight an urgent need to develop crop plants with enhanced tolerance to stress combination. The uniqueness and complexity of the physiological and molecular response of plants to each different stress combination, highlighted here, demonstrate the daunting challenge we face in accomplishing this goal. Dedicated efforts combining field experimentation, omics, and network analyses, coupled with advanced phenotyping and breeding methods, will be needed to address specific crops and particular stress combinations relevant to maintaining our future food chain secured.

Genome-wide analysis of WRKY transcription factors in Aquilaria sinensis (Lour.) Gilg

The WRKY proteins are a superfamily of transcription factor that regulate diverse developmental and physiological processes in plants. Completion of the whole-genome sequencing of Aquilaria sinensis allowed us to perform a genome-wide investigation for WRKY proteins. Here, we predicted 70 WRKY genes from the A. sinensis genome and undertaken a comprehensive bioinformatic analysis. Due to their diverse structural features, the 70 AsWRKY genes are classified into three main groups (group I-III), with five subgroups (IIa-IIe) in group II, except two belong to none of them. Distinct expression profiles of AsWRKYs with RNA sequencing data revealed their diverse expression patterns among different tissues and in the process of whole-tree-inducing agarwood formation. Based on the expression characteristics, we predict some AsWRKYs are pseudogenes, and some may be involved in the biosynthesis of agarwood sesquiterpenes as activators or repressors. Among the tested genes treated with MeJA and H₂O₂, most of them are induced by H₂O₂, but downregulated by MeJA, implying the complexity of their involvement in signal transduction regulation. Our results not only provide a basic platform for functional identification of WRKYs in A. sinensis but important clues for further analysis their regulation role in agarwood formation.

Enhancement of secondary metabolites in Bacopa monnieri (L.) Pennell plants treated with copper-based nanoparticles in vivo

The study aims to document the effect of starch-stabilised copper-based nanoparticles (CuNPs) on the biosynthesis of pharmaceutically valuable secondary metabolites, especially saponins, of the reputed nootropic herb Bacopa monnieri (L.) Pennell. CuNPs were synthesised chemically by the reduction of cupric sulphate pentahydrate with ascorbic acid using starch as the capping agent. They were characterised by UV-visible spectrophotometry, Fourier-transform infra-red spectroscopy, X-ray diffraction, high-resolution transmission electron microscopy and zeta potential. The nanoparticles consisted of cuprous oxide and metallic copper, were approximately spherical, polydispersed with diameter <20 nm. Hydroponically grown B. monnieri plants were treated in vivo with the CuNPs between the concentrations of 0-100 mg l-1. Spectrophotometric estimation of the total contents of saponins, alkaloids, phenolics, flavonoids and DPPH radical scavenging capacity from the methanolic extracts of the whole plants showed a hormetic increase in the content of secondary metabolites in a concentration-dependent manner from 5 mg l-1 until it declined at toxic metabolic concentration. This was accompanied by an increase in ROS markers hydrogen peroxide and malondialdehyde as well as a hormetic effect on activities of phenylalanine ammonia lyase and antioxidant enzymes catalase, ascorbate peroxidase and superoxide dismutase. CuNPs at sub-toxic concentrations were found to enhance secondary metabolism and antioxidant capacity in Bacopa monnieri through ROS-mediated defence response.

Hierarchical and Dynamic Regulation of Defense-Responsive Specialized Metabolism by WRKY and MYB Transcription Factors

The plant kingdom produces hundreds of thousands of specialized bioactive metabolites, some with pharmaceutical and biotechnological importance. Their biosynthesis and function have been studied for decades, but comparatively less is known about how transcription factors with overlapping functions and contrasting regulatory activities coordinately control the dynamics and output of plant specialized metabolism. Here, we performed temporal studies on pathogen-infected intact host plants with perturbed transcription factors. We identified WRKY33 as the condition-dependent master regulator and MYB51 as the dual functional regulator in a hierarchical gene network likely responsible for the gene expression dynamics and metabolic fluxes in the camalexin and 4-hydroxy-indole-3-carbonylnitrile (4OH-ICN) pathways. This network may have also facilitated the regulatory capture of the newly evolved 4OH-ICN pathway in Arabidopsis thaliana by the more-conserved transcription factor MYB51. It has long been held that the plasticity of plant specialized metabolism and the canalization of development should be differently regulated; our findings imply a common hierarchical regulatory architecture orchestrated by transcription factors for specialized metabolism and development, making it an attractive target for metabolic engineering.

A Comprehensive Transcriptome-Wide Identification and Screening of WRKY Gene Family Engaged in Abiotic Stress in Glycyrrhiza glabra

The study reports 147 full-length WRKY genes based on the transcriptome analysis of Glycyrrhiza genus (G. glabra and G. uralensis). Additional motifs in G. glabra included DivIVA (GgWRKY20) and SerS Superfamily (GgWRKY21) at the C-terminal, and Coat family motifs (GgWRKY55) at the N-terminal of the proteins, while Exo70 exo cyst complex subunit of 338 amino acid (GuWRKY9) was present at the N-terminal of G. uralensis only. Plant Zn cluster super-family domain (17 WRKYs) and bZIP domain (2 WRKYs) were common between the two species. Based on the number of WRKY domains, sequence alignment and phylogenesis, the study identified GuWRKY27 comprising of 3 WRKY domains in G. uralensis and a new subgroup-IIf (10 members), having novel zinc finger pattern (C-X4-C-X22-HXH) in G. glabra. Multiple WRKY binding domains (1-11) were identified in the promoter regions of the GgWRKY genes indicating strong interacting network between the WRKY proteins. Tissue-specific expression of 25 GgWRKYs, under normal and treated conditions, revealed 11 of the 18 induction factor triggered response corroborating to response observed in AtWRKYs. The study identified auxin-responsive GgWRKY 55 & GgWRKY38; GA3 responsive GgWRKYs15&59 in roots and GgWRKYs8, 20, 38, 57 &58 in the shoots of the treated plant. GgWRKYs induced under various stresses included GgWRKY33 (cold), GgWRKY4 (senescence), GgWRKYs2, 28 & 33 (salinity) and GgWRKY40 (wounding). Overall, 23 GgWRKYs responded to abiotic stress, and 17 WRKYs were induced by hormonal signals. Of them 13 WRKYs responded to both suggesting inter-connection between hormone signalling and stress response. The present study will help in understanding the transcriptional reprogramming, protein-protein interaction and cross-regulation during stress and other physiological processes in the plant.

Genomic Variance and Transcriptional Comparisons Reveal the Mechanisms of Leaf Color Affecting Palatability and Stressed Defense in Tea Plant

Leaves are one of the most important organs of plants, and yet, the association between leaf color and consumable traits remains largely unclear. Tea leaves are an ideal study system with which to investigate the mechanism of how leaf coloration affects palatability, since tea is made from the leaves of the crop Camellia sinensis. Our genomic resequencing analysis of a tea cultivar ZiJuan (ZJ) with purple leaves and altered flavor revealed genetic variants when compared with the green-leaf, wild type cultivar YunKang(YK). RNA-Seq based transcriptomic comparisons of the bud and two youngest leaves in ZJ and YK identified 93%, 9% and 5% expressed genes that were shared in YK- and ZJ-specific cultivars, respectively. A comparison of both transcript abundance and particular metabolites revealed that the high expression of gene UFGT for anthocyanin biosynthesis is responsible for purple coloration, which competes with the intermediates for catechin-like flavanol biosynthesis. Genes with differential expression are enriched in response to stress, heat and defense, and are casually correlated with the environmental stress of ZJ plant origin in the Himalayas. In addition, the highly expressed C4H and LDOX genes for synthesizing flavanol precursors, ZJ-specific CLH1 for degrading chlorophyll, alternatively spliced C4H and FDR and low photosynthesis also contributed to the altered color and flavor of ZJ. Thus, our study provides a better molecular understanding of the effect of purple coloration on leaf flavor, and helps to guide future engineering improvement of palatability.

Understanding the Effect of Structural Diversity in WRKY Transcription Factors on DNA Binding Efficiency through Molecular Dynamics Simulation

A precise understanding of the molecular mechanism involved in stress conditions has great importance for crop improvement. Biomolecules, such as WRKY proteins, which are the largest transcription factor family that is widely distributed in higher plants, plays a significant role in plant defense response against various biotic and abiotic stressors. In the present study, an extensive homology-based three-dimensional model construction and subsequent interaction study of WRKY DNA-binding domain (DBD) in CcWRKY1 (Type I), CcWRKY51 (Type II), and CcWRKY70 (Type III) belonging to pigeonpea, a highly tolerant crop species, was performed. Evaluation of the generated protein models was done to check their reliability and accuracy based on the quantitative and qualitative parameters. The final model was subjected to investigate the comparative binding analysis of different types of WRKY–DBD with DNA-W-box (a cis-acting element) by protein–DNA docking and molecular dynamics (MD) simulation. The DNA binding specificity with WRKY variants was scrutinized through protein–DNA interaction using the HADDOCK server. The stability, as well as conformational changes of protein–DNA complex, was investigated through molecular dynamics (MD) simulations for 100 ns using GROMACS. Additionally, the comparative stability and dynamic behavior of each residue of the WRKY–DBD type were analyzed in terms of root mean square deviation (RMSD), root mean square fluctuation (RMSF)values of the backbone atoms for each frame taking the minimized structure as a reference. The details of DNA binding activity of three different types of WRKY–DBD provided here will be helpful to better understand the regulation of WRKY gene family members in plants.

De novo transcriptomic analysis of light-induced flavonoid pathway, transcription factors in the flower buds of Lonicera japonica

KEY MESSAGE

Transcriptomic analysis of the relationship between gene expression patterns and flavonoid contents in the flower buds of Lonicera japonica under light-induced conditions, especially the flavonoid pathway genes and transcription factors.

ABSTRACT

Flos Lonicerae Japonicae (FLJ), the flower buds of Lonicera japonica Thunb., has been used to treat some human diseases including severe respiratory syndromes and hand-foot-and-mouth diseases owing to its putative antibacterial, and antiviral effects. Luteoloside is a flavonoid that is used by the Chinese Pharmacopoeia to evaluate the quality of FLJ. Light is an important environmental factor that affects flavonoid biosynthesis in the flower buds of L. japonica. However, how light triggers increases in flavonoid production remains unclear. To enhance our understanding of the mechanism involved in light-regulated flavonoid biosynthesis, we sequenced the transcriptomes of L. japonica exposed to three different light conditions: 100% light intensity (CK), 50% light intensity (LI50), and 25% light intensity (LI25) using an Illumina HiSeq 4000 System. A total of 77,297 unigenes with an average length of 809 bp were obtained. Among them, 43,334 unigenes (56.06%) could be matched to at least one biomolecular database. Additionally, 4188, 1545 and 1023 differentially expressed genes (DEGs) were identified by comparative transcriptomics LI25-vs-CK, LI50-vs-CK, and LI25-vs-LI50, respectively. Of note, genes known to be involved in flavonoid biosynthesis, such as 4-coumarate coenzyme A ligase (4CL), and chalcone synthase (CHS) were up-regulated. In addition, a total of 1649 transcription factors (TFs) were identified and divided into 58 TF families; 98 TFs exhibited highly dynamic changes in response to light intensity. Quantitative real-time PCR (qRT-PCR) was used to test the expression profiles of the RNA sequencing (RNA-Seq) data. This study offers insight into how transcriptional expression pattern is influenced by light in the flower buds of L. japonica, and will enhance the understanding of molecular mechanisms of flavonoid biosynthesis in response to light in L. japonica.

A WRKY transcription factor PbrWRKY53 from Pyrus betulaefolia is involved in drought tolerance and AsA accumulation

WRKY comprises a large family of transcription factors in plants, but most WRKY members are still poorly understood. In this study, we report the identification and functional characterization of PbrWRKY53 isolated from Pyrus betulaefolia. PbrWRKY53 was greatly up-regulated by drought and abscisic acid, but slightly induced by salt and cold. Subcellar localization analyses showed that PbrWRKY53 was located in the nucleus. Ectopic expression of PbrWRKY53 in tobacco and Pyrus ussuriensis conferred enhanced tolerance to drought stress. The transgenic plants exhibited better water status, less reactive oxygen species generation and higher levels of antioxidant enzyme activities and metabolites than the wild type. In addition, overexpression of PbrWRKY53 in transgenic tobacco resulted in enhanced expression level of PbrNCED1, and led to the increase in larger amount of vitamin C accumulation in comparison to WT. Knock-down of PbrWRKY53 in P. ussuriensis down-regulated PbrNCED1 abundance, accompanied by compromised drought tolerance. Yeast one-hybrid assay, EMSA and transient expression analysis demonstrated that PbrWRKY53 could bind to the W-box element in the promoter region of PbrNCED1. Taken together, these results demonstrated that PbrWRKY53 plays a positive role in drought tolerance, which might be, at least in part, promoting production of vitamin C via regulating PbrNCED1 expression.

Expansion of a core regulon by transposable elements promotes Arabidopsis chemical diversity and pathogen defense

Plants synthesize numerous ecologically specialized, lineage-specific metabolites through biosynthetic gene duplication and functional specialization. However, it remains unclear how duplicated genes are wired into existing regulatory networks. We show that the duplicated gene CYP82C2 has been recruited into the WRKY33 regulon and indole-3-carbonylnitrile (ICN) biosynthetic pathway through exaptation of a retroduplicated LINE retrotransposon (EPCOT3) into an enhancer. The stepwise development of a chromatin-accessible WRKY33-binding site on EPCOT3 has potentiated the regulatory neofunctionalization of CYP82C2 and the evolution of inducible defense metabolite 4-hydroxy-ICN in Arabidopsis thaliana. Although transposable elements (TEs) have long been recognized to have the potential to rewire regulatory networks, these results establish a more complete understanding of how duplicated genes and TEs contribute in concert to chemical diversity and pathogen defense.

Arabidopsis plants can produce 4-hydroxyindole-3-carbonitrile (4OH-ICN) upon pathogen infection. Here, the authors show that EPCOT3, a retrotransposonderived enhancer, mediates WRKY33-binding, pathogen-responsive transcription of CYP82C2, and synthesis of 4OH-ICN.

Differential physiology and expression of phenylalanine ammonia lyase (PAL) and universal stress protein (USP) in the endangered species Astragalus fridae following seed priming with cold plasma and manipulation of culture medium with silica nanoparticles

KEY MESSAGE

Seed priming with cold plasma in combination with manipulation of culture medium with silica nanoparticle provokes anatomical, physiological and molecular changes, thereby reinforcing the plant growth and protection.

ABSTRACT

This study addressed responses of Astragalus fridae to seed priming with cold plasma (0.84 W/cm²; 0, 30, 60, and 90 s) and applications of SiO₂ nanoparticle (nSi; 0, 5, 40, and 80 mgl^-1) in culture medium (an in vitro study). FE-SEM confirmed nSi uptake and translocation. Bulk Si at high concentrations reduced biomass accumulation (mean = 45%), while nSi did not make significant differences. The growth-enhancing effects of plasma by 41.5% were promoted by the nSi supplementation and reached 71%. Plasma did not make significant changes in Chla, while led to the slightly higher (mean = 14%) Chlb. The presence of nSi at high doses caused slight reductions in Chlb (mean = 25%) which were mitigated by plasma. The plasma and/or nSi treatments modified activities of phenylalanine ammonia lyase (PAL) in both roots (mean = 32%) and leaves (mean = 44%). With a similar trend, both individual and combined treatments of plasma and nSi provoked inductions in peroxidase activities in roots and leaves. The simultaneous treatments of plasma and nSi had the highest expression rates of PAL gene. The individual treatments of plasma did not make a significant difference in the expression of universal stress protein (USP) gene, whereas the nSi-treated seedlings exhibited the higher expression rates of USP. Leaf thicknesses and development of the vascular system (xylem and phloem) were reinforced in response to plasma and nSi. The findings provide evidence on potential benefits and phytotoxicity of nSi and plasma which may be employed as a theoretical basis for possible exploitation.

Abscisic acid and ethylene biosynthesis-related genes are associated with anthocyanin accumulation in purple ornamental cabbage (Brassica oleracea var. acephala)

Purple ornamental cabbage (Brassica oleracea var. acephala) is a popular decorative plant, cultivated for its colorful leaf rosettes that persist in cool weather. It is characterized by green outer leaves and purple inner leaves, whose purple pigmentation is due to the accumulation of anthocyanin pigments. Phytohormones play important roles in anthocyanin biosynthesis in other species. Here, we identified 14 and 19 candidate genes putatively involved in abscisic acid (ABA) and ethylene (ET) biosynthesis, respectively, in B. oleracea. We determined the expression patterns of these candidate genes by reverse-transcription quantitative PCR (RT-qPCR). Among candidate ABA biosynthesis-related genes, the expressions of BoNCED2.1, BoNCED2.2, BoNCED6, BoNCED9.1, and BoAAO3.2 were significantly higher in purple compared to green leaves. Likewise, most of the ET biosynthetic genes (BoACS6, BoACS9.1, BoACS11, BoACO1.1, BoACO1.2, BoACO3.1, BoACO4, and BoACO5) had significantly higher expression in purple compared to green leaves. Among these genes, BoNCED2.1, BoNCED2.2, BoACS11, and BoACO4 showed particularly strong associations with total anthocyanin content of the purple inner leaves. Our results suggest that ABA and ET might promote the intense purple pigmentation of the inner leaves of purple ornamental cabbage.

Transcriptomic profiles of 33 opium poppy samples in different tissues, growth phases, and cultivars

Opium poppy is one of the most important medicinal plants and remains the only commercial resource of morphinan-based painkillers. However, little is known about the regulatory mechanisms involved in benzylisoquinoline alkaloids (BIAs) biosynthesis in opium poppy. Herein, the full-length transcriptome dataset of opium poppy was constructed for the first time in accompanied with the 33 samples of Illumina transcriptome data from different tissues, growth phases and cultivars. The long-read sequencing produced 902,140 raw reads with 55,114 high-quality transcripts, and short-read sequencing produced 1,923,679,864 clean reads with an average Q30 rate of 93%. The high-quality transcripts were subsequently quantified using the short reads, and the expression of each unigene among different samples was calculated as reads per kilobase per million mapped reads (RPKM). These data provide a foundation for opium poppy transcriptomic analysis, which may aid in capturing splice variants and some non-coding RNAs involved in the regulation of BIAs biosynthesis. It can also be used for genome assembly and annotation which will favor in new transcript identification.

Pitaya HpWRKY3 Is Associated with Fruit Sugar Accumulation by Transcriptionally Modulating Sucrose Metabolic Genes HpINV2 and HpSuSy1

Sugar level is an important determinant of fruit taste and consumer preferences. However, upstream regulators that control sugar accumulation during fruit maturation are poorly understood. In the present work, we found that glucose is the main sugar in mature pitaya (Hylocereus) fruit, followed by fructose and sucrose. Expression levels of two sucrose-hydrolyzing enzyme genes HpINV2 and HpSuSy1 obviously increased during fruit maturation, which were correlated well with the elevated accumulation of glucose and fructose. A WRKY transcription factor HpWRKY3 was further identified as the putative binding protein of the HpINV2 and HpSuSy1 promoters by yeast one-hybrid and gel mobility shift assays. HpWRKY3 was localized exclusively in the nucleus and possessed trans-activation ability. HpWRKY3 exhibited the similar expression pattern with HpINV2 and HpSuSy1. Finally, transient expression assays in tobacco leaves showed that HpWRKY3 activated the expressions of HpINV2 and HpSuSy1. Taken together, we propose that HpWRKY3 is associated with pitaya fruit sugar accumulation by activating the transcriptions of sucrose metabolic genes. Our findings thus shed light on the transcriptional mechanism that regulates the sugar accumulation during pitaya fruit quality formation.

The genome of the medicinal plant Andrographis paniculata provides insight into the biosynthesis of the bioactive diterpenoid neoandrographolide

Andrographis paniculata is a herbaceous dicot plant widely used for its anti-inflammatory and anti-viral properties across its distribution in China, India and other Southeast Asian countries. A. paniculata was used as a crucial therapeutic treatment during the influenza epidemic of 1919 in India, and is still used for the treatment of infectious disease in China. A. paniculata produces large quantities of the anti-inflammatory diterpenoid lactones andrographolide and neoandrographolide, and their analogs, which are touted to be the next generation of natural anti-inflammatory medicines for lung diseases, hepatitis, neurodegenerative disorders, autoimmune disorders and inflammatory skin diseases. Here, we report a chromosome-scale A. paniculata genome sequence of 269 Mb that was assembled by Illumina short reads, PacBio long reads and high-confidence (Hi-C) data. Gene annotation predicted 25 428 protein-coding genes. In order to decipher the genetic underpinning of diterpenoid biosynthesis, transcriptome data from seedlings elicited with methyl jasmonate were also obtained, which enabled the identification of genes encoding diterpenoid synthases, cytochrome P450 monooxygenases, 2-oxoglutarate-dependent dioxygenases and UDP-dependent glycosyltransferases potentially involved in diterpenoid lactone biosynthesis. We further carried out functional characterization of pairs of class-I and -II diterpene synthases, revealing the ability to produce diversified labdane-related diterpene scaffolds. In addition, a glycosyltransferase able to catalyze O-linked glucosylation of andrograpanin, yielding the major active product neoandrographolide, was also identified. Thus, our results demonstrate the utility of the combined genomic and transcriptomic data set generated here for the investigation of the production of the bioactive diterpenoid lactone constituents of the important medicinal herb A. paniculata.

Physiological and transcriptional responses to low-temperature stress in rice genotypes at the reproductive stage

Rice growth and productivity is adversely affected by low-temperature stress. From a previous screen of diverse rice genotypes for cold tolerance parameters at the vegetative stage, we selected the tolerant Nipponbare and M202 genotypes and sensitive Cypress and Secano do Brazil genotypes for further analysis at the reproductive stage for physiological and yield parameters. Cold stress severely affected grain yield as estimated by the number of grain per panicle, panicle length, and 100 seed weight. Analysis of gene expression of 21 genes involved in physiological responses to low temperature tested, in the flag leaf and inflorescence tissue of these genotypes, showed an increased expression of the Lipid Transfer Protein genes LTP7 and LTP10 in flag leaf tissue of the tolerant Nipponbare and M202, along with a significant increase in the relative expression of stress-responsive transcription factors (TFs) and cold-inducible genes. In flag leaf tissue OsNAC9, OsNAC10 and OsNAP genes showed high correlation with photosynthesis, stomatal conductance, transpiration and Quantum Efficiency of PSII. In consequence of the foregoing results, we conclude that Nipponbare and M202 are cold tolerant genotypes and that LTP7, LTP10, OsNAC9, OsNAC10 and OsNAP genes can be used as markers in screening for cold tolerance at the reproductive stage. Furthermore based on the results we propose a model of low-temperature tolerance mechanism of how stress is perceived, and how the signal cascade acts to promote tolerance at below-ideal temperatures.

Multiple molecular defense strategies in Brachypodium distachyon surmount Hessian fly (Mayetiola destructor) larvae-induced susceptibility for plant survival

The Hessian fly is a destructive pest of wheat causing severe economic damage. Numerous genes and associated biological pathways have been implicated in defense against Hessian fly. However, due to limited genetic resources, compounded with genome complexity, functional analysis of the candidate genes are challenging in wheat. Physically, Brachypodium distachyon (Bd) exhibits nonhost resistance to Hessian fly, and with a small genome size, short life cycle, vast genetic resources and amenability to transformation, it offers an alternate functional genomic model for deciphering plant-Hessian fly interactions. RNA-sequencing was used to reveal thousands of Hessian fly-responsive genes in Bd one, three, and five days after egg hatch. Genes encoding defense proteins, stress-regulating transcription factors, signaling kinases, and secondary metabolites were strongly up-regulated within the first 24 hours of larval feeding indicating an early defense, similar to resistant wheat. Defense was mediated by a hypersensitive response that included necrotic lesions, up-regulated ROS-generating and -scavenging enzymes, and H₂O₂ production. Suppression of cell wall-associated proteins and increased cell permeability in Bd resembled susceptible wheat. Thus, Bd molecular responses shared similarities to both resistant and susceptible wheat, validating its suitability as a model genome for undertaking functional studies of candidate Hessian fly-responsive genes.

GmWRKY16 Enhances Drought and Salt Tolerance Through an ABA-Mediated Pathway in Arabidopsis thaliana

The WRKY transcription factors (TFs) are one of the largest families of TFs in plants and play multiple roles in plant development and stress response. In the present study, GmWRKY16 encoding a WRKY transcription factor in soybean was functionally characterized in Arabidopsis. GmWRKY16 is a nuclear protein that contains a highly conserved WRKY domain and a C2H2 zinc-finger structure, and has the characteristics of transcriptional activation ability, presenting a constitutive expression pattern with relative expression levels of over fourfold in the old leaves, flowers, seeds and roots of soybean. The results of quantitative real time polymerase chain reaction (qRT-PCR) showed that GmWRKY16 could be induced by salt, alkali, ABA, drought and PEG-6000. As compared with the control, overexpression of GmWRKY16 in Arabidopsis increased the seed germination rate and root growth of seedlings in transgenic lines under higher concentrations of mannitol, NaCl and ABA. In the meantime, GmWRKY16 transgenic lines showed over 75% survival rate after rehydration and enhanced Arabidopsis tolerance to salt and drought with higher proline and lower MDA accumulation, less water loss of the detached leaves, and accumulated more endogenous ABA than the control under stress conditions. Further studies showed that AtWRKY8, KIN1, and RD29A were induced in GmWRKY16 transgenic plants under NaCl treatment. The expressions of the ABA biosynthesis gene (NCED3), signaling genes (ABI1, ABI2, ABI4, and ABI5), responsive genes (RD29A, COR15A, COR15B, and RD22) and stress-related marker genes (KIN1, LEA14, LEA76, and CER3) were regulated in transgenic lines under drought stress. In summary, these results suggest that GmWRKY16 as a WRKY TF may promote tolerance to drought and salt stresses through an ABA-mediated pathway.

Genome-wide identification and expression analyses of WRKY transcription factor family members from chickpea (Cicer arietinum L.) reveal their role in abiotic stress-responses

BACKGROUND

WRKY proteins play a vital role in the regulation of several imperative plant metabolic processes and pathways, especially under biotic and abiotic stresses. Although WRKY genes have been characterized in various major crop plants, their identification and characterization in pulse legumes is still in its infancy. Chickpea (Cicer arietinum L.) is the most important pulse legume grown in arid and semi-arid tropics.

OBJECTIVE

In silico identification and characterization of WRKY transcription factor-encoding genes in chickpea genome.

METHODS

For this purpose, a systematic genome-wide analysis was carried out to identify the non-redundant WRKY transcription factors in the chickpea genome.

RESULTS

We have computationally identified 70 WRKY-encoding non-redundant genes which were randomly distributed on all the chickpea chromosomes except chromosome 8. The evolutionary phylogenetic analysis classified the WRKY proteins into three major groups (I, II and III) and seven sub-groups (IN, IC, IIa, IIb, IIc, IId and IIe). The gene structure analysis revealed the presence of 2-7 introns among the family members. Along with the presence of absolutely conserved signatory WRKY domain, 19 different domains were also found to be conserved in a group-specific manner. Insights of gene duplication analysis revealed the predominant role of segmental duplications for the expansion of WRKY genes in chickpea. Purifying selection seems to be operated during the evolution and expansion of paralogous WRKY genes. The transcriptome data-based in silico expression analysis revealed the differential expression of CarWRKY genes in root and shoot tissues under salt, drought, and cold stress conditions. Moreover, some of these genes showed identical expression pattern under these stresses, revealing the possibility of involvement of these genes in conserved abiotic stress-response pathways.

CONCLUSION

This genome-wide computational analysis will serve as a base to accelerate the functional characterization of WRKY TFs especially under biotic and abiotic stresses.

VvWRKY8 represses stilbene synthase genes through direct interaction with VvMYB14 to control resveratrol biosynthesis in grapevine

Resveratrol (Res) is a stilbenoid, a group of plant phenolic metabolites derived from stilbene that possess activities against pests, pathogens, and abiotic stresses. Only a few species, including grapevine (Vitis), synthesize and accumulate Res. Although stilbene synthases (STSs) have been isolated and characterized in several species, the gene regulatory mechanisms underlying stilbene biosynthesis are still largely unknown. Here, we characterize a grapevine WRKY transcription factor, VvWRKY8, that regulates the Res biosynthetic pathway. Transient and stable overexpression of VvWRKY8 in grapevine results in decreased expression of VvSTS15/21 and VvMYB14, as well as in a reduction of Res accumulation. VvWRKY8 does not bind to or activate the promoters of VvMYB14 and VvSTS15/21; however, it physically interacts with VvMYB14 proteins through their N-terminal domains to prevent them from binding to the VvSTS15/21 promoter. Application of exogenous Res results in the stimulation of VvWRKY8 expression and in a decrease of VvMYB14 and VvSTS15/21 expression in grapevine suspension cells, and in the activation of the VvWRKY8 promoter in tobacco leaves. These results demonstrate that VvWRKY8 represses VvSTS15/21 expression and Res biosynthesis through interaction with VvMYB14. In this context, the VvMYB14-VvSTS15/21-Res-VvWRKY8 regulatory loop may be an important mechanism for the fine-tuning of Res biosynthesis in grapevine.

Regulation of resveratrol biosynthesis in grapevine: new approaches for disease resistance?

This article comments on: Jiang J, Xi H, Dai Z, Lecourieux F, Yuan L, Liu X, Patra B, Wei Y, Li S, Wang L. 2019. VvWRKY8 negatively regulates VvSTS through direct interaction with VvMYB14 to balance resveratrol biosynthesis in grapevine. Journal of Experimental Botany 70, 715–729.

Transcriptome driven characterization of curly- and smooth-leafed endives reveals molecular differences in the sesquiterpenoid pathway

Endives (Cichorium endivia L.) are popular vegetables, diversified into curly/frisée- and smooth/broad-leafed (escaroles) cultivar types (cultigroups), and consumed as fresh and bagged salads. They are rich in sesquiterpene lactones (STL) that exert proven function on bitter taste and human health. The assembly of a reference transcriptome of 77,022 unigenes and RNA-sequencing experiments were carried out to characterize the differences between endives and escaroles at the gene structural and expression levels. A set of 3177 SNPs distinguished smooth from curly cultivars, and an SNP-supported phylogenetic tree separated the cultigroups into two distinct clades, consistently with the botanical varieties of origin (crispum and latifolium, respectively). A pool of 699 genes maintained differential expression pattern (core-DEGs) in pairwise comparisons between curly vs smooth cultivars grown in the same environment. Accurate annotation allowed the identification of 26 genes in the sesquiterpenoid biosynthesis pathway, which included several g ermacrene A s ynthase, g ermacrene A o xidase and co stunolide s ynthase members (GAS/GAO/COS module), required for the synthesis of costunolide, a key precursor of lactucopicrin- and lactucin-like sesquiterpene lactones. The core-DEGs contained a GAS gene (contig83192) that was positively correlated with STL levels and recurrently more expressed in curly than smooth endives, suggesting a cultigroup-specific behavior. The significant positive correlation of GAS/GAO/COS transcription and STL abundance (2.4-fold higher in frisée endives) suggested that sesquiterpenoid pathway control occurs at the transcriptional level. Based on correlation analyses, five transcription factors (MYB, MYB-related and WRKY) were inferred to act on contig83192/GAS and specific STL, suggesting the occurrence of two distinct routes in STL biosynthesis.

Tomato SlWRKY3 acts as a positive regulator for resistance against the root-knot nematode Meloidogyne javanica by activating lipids and hormone-mediated defense-signaling pathways

Diseases caused by plant-parasitic nematodes in vegetables, among them Meloidogyne spp. root-knot nematodes (RKNs), lead to extensive yield decline. A molecular understanding of the mechanisms underlying plants' innate resistance may enable developing safe alternatives to harmful chemical nematicides in controlling RKNs. A tight relationship has been revealed between the WRKY transcription factors and RKN parasitism on tomato roots. We investigated the function role of tomato SlWRK3 and SlWRKY35 in regulating nematode disease development. Using promoter-GUS reporter gene fusions, we show that both SlWRKY3 and SlWRKY35 are induced within 5 days of infection and through feeding-site development and gall maturation, with a much stronger response of the former vs. the latter to nematode infection. Histological analysis of nematode-feeding sites indicated a high expression of SlWRKY3 in developing and mature feeding cells and associated vasculature cells, whereas SlWRKY35 expression was only observed in mature feeding sites. Both SlWRKY3 and SlWRKY35 promoters were induced by the defense phytohormones salicylic acid and indole-3-butyric acid, with no response to either jasmonic acid or methyl jasmonate. SlWRKY3 overexpression resulted in lower infection of the RKN Meloidogyne javanica, whereas knocking down SlWRKY3 resulted in increased infection. Phytohormone and oxylipin profiles determined by LC-MS/MS showed that the enhanced resistance in the former is coupled with an increased accumulation of defense molecules from the shikimate and oxylipin pathways. Our results pinpoint SlWRKY3 as a positive regulator of induced resistance in response to nematode invasion and infection, mostly during the early stages of nematode infection.

Comparative Transcriptome Analysis Identifies Genes Involved in Diosgenin Biosynthesis in Trigonella foenum-graecum L

Trigonella foenum-graecum L. (fenugreek) is a valuable resource of producing diosgenin which serves as a substrate for synthesizing more than two hundred kinds of steroidal drugs. Phytochemical analysis indicated that methyl jasmonate (MeJA) efficiently induced diosgenin biosynthesis in fenugreek seedlings. Though early steps up to cholesterol have recently been elucidated in plants, cytochrome P450 (CYP)- and glycosyltransferase (GT)-encoding genes involved in the late steps from cholesterol to diosgenin remain unknown. This study established comparative fenugreek transcriptome datasets from the MeJA-treated seedlings and the corresponding control lines. Differential gene expression analysis identified a number of MeJA-induced CYP and GT candidate genes. Further gene expression pattern analysis across a different MeJA-treating time points, together with a phylogenetic analysis, suggested specific family members of CYPs and GTs that may participate in the late steps during diosgenin biosynthesis. MeJA-induced transcription factors (TFs) that may play regulatory roles in diosgenin biosynthesis were also discussed. This study provided a valuable genetic resource to functionally characterize the genes involved in diosgenin biosynthesis, which will push forward the production of diosgenin in microbial organisms using a promising synthetic biology strategy.

Molecular Characterization of WRKY Transcription Factors That Act as Negative Regulators of O-Methylated Catechin Biosynthesis in Tea Plants ( Camellia sinensis L.)

Tea O-methylated catechins, especially (-)-epigallocatechin 3- O-(3- O-methyl)gallate (EGCG3″Me), have been attracting much attention as a result of their positive health effects. The transcription regulators of O-methylated catechin biosynthesis remain elusive. In this study, the expression pattern of genes related to O-methylated catechin biosynthesis, including CsLAR, CsANS, CsDFR, CsANR, and CCoAOMT, in three tea cultivars with different contents of EGCG3″Me was investigated. Two WRKY transcription factors (TFs), designated as CsWRKY31 and CsWRKY48, belonging to groups IIb and IIc of the WRKY family, respectively, were further identified. CsWRKY31 and CsWRKY48 were nuclear-localized proteins and possessed transcriptional repression ability. Furthermore, expression of CsWRKY31 and CsWRKY48 showed negative correlation with CsLAR, CsDFR, and CCoAOMT during EGCG3″Me accumulation in tea leaves. More importantly, W-box (C/T)TGAC(T/C) elements were located in the promoter of CsLAR, CsDFR, and CCoAOMT, and further assays revealed that CsWRKY31 and CsWRKY48 were capable of repressing the transcription of CsLAR, CsDFR, and CCoAOMT via the attachment of their promoters to the W-box elements. Collectively, our findings identify two novel negative regulators of O-methylated catechin biosynthesis in tea plants, which might provide a potential strategy to breed high-quality tea cultivar.

Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus

3-O-caffeoylquinic acid, also known as chlorogenic acid (CGA), functions as an intermediate in lignin biosynthesis in the phenylpropanoid pathway. It is widely distributed among numerous plant species and acts as an antioxidant in both plants and animals. Using GC-MS, we discovered consistent and extreme variation in CGA content across a population of 739 4-yr-old Populus trichocarpa accessions. We performed genome-wide association studies (GWAS) from 917 P. trichocarpa accessions and expression-based quantitative trait loci (eQTL) analyses to identify key regulators. The GWAS and eQTL analyses resolved an overlapped interval encompassing a hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyl transferase 2 (PtHCT2) that was significantly associated with CGA and partially characterized metabolite abundances. PtHCT2 leaf expression was significantly correlated with CGA abundance and it was regulated by cis-eQTLs containing W-box for WRKY binding. Among all nine PtHCT homologs, PtHCT2 is the only one that responds to infection by the fungal pathogen Sphaerulina musiva (a Populus pathogen). Validation using protoplast-based transient expression system suggests that PtHCT2 is regulated by the defense-responsive WRKY. These results are consistent with reports of CGA functioning as an antioxidant in response to biotic stress. This study provides insights into data-driven and omics-based inference of gene function in woody species.

Activation of rice WRKY transcription factors: an army of stress fighting soldiers?

Rice WRKYs comprise a large family of transcription factors and present remarkable structure features and a unique DNA binding site. Their importance in plants goes beyond the response to stressful stimuli, since they participate in hormonal pathways and developmental processes. Indeed, the majority of WRKYs present an independent activation since they are able to perform self-transcriptional regulation. However, some WRKY activation depends on epigenetic and transcript regulation by micro RNAs. Their protein function depends, almost always, on the posttranslational changes. Taking to account its properties of auto-activation, all these regulators process are extremely important for complete WRKY regulation. In this sense, here we provide an overview of transcriptional activation and posttranscriptional and posttranslational regulation of rice WRKY genes under stresses.

Coregulation of Biosynthetic Genes and Transcription Factors for Aporphine-Type Alkaloid Production in Wounded Lotus Provides Insight into the Biosynthetic Pathway of Nuciferine

Lotus (Nelumbo nucifera Gaertn.) contains various bioactive compounds, with benzylisoquinoline alkaloids (BIAs) as one of the major groups. The biosynthetic pathways of two major bioactive BIAs in this plant, nuciferine and N-nornuciferine, are still not clear. Therefore, several genes related to BIA biosynthesis were searched from the lotus database to study the role of key genes in regulating these pathways. In this study, the expression profiles of NCS, CNMT, 6OMT, CYP80G2, and WRKY TFs were investigated in mechanically wounded lotus leaves. It was found that the accumulation of nuciferine and N-nornuciferine significantly increased in the mechanically wounded lotus leaves in accordance with the relative expression of putative CYP80G2 and one WRKY transcription factor (NNU_24385), with the coregulation of CNMT. Furthermore, the role of methyltransferase-related genes in this study suggested that methylation of the isoquinoline nucleus to yield a methylated-BIA structure may occur at the N position before the O position. Altogether, this study provides improved understanding of the genes regulating BIA biosynthesis under stressed conditions, which could lead to improvements in BIA production from the commercial lotus.

Coactivation of MEP-biosynthetic genes and accumulation of abietane diterpenes in Salvia sclarea by heterologous expression of WRKY and MYC2 transcription factors

Plant abietane diterpenoids (e.g. aethiopinone, 1- oxoaethiopinone, salvipisone and ferruginol), synthesized in the roots of several Salvia spp, have antibacterial, antifungal, sedative and anti-proliferative properties. Recently we have reported that content of these compounds in S. sclarea hairy roots is strongly depending on transcriptional regulation of genes belonging to the plastidial MEP-dependent terpenoid pathway, from which they mostly derive. To boost the synthesis of this interesting class of compounds, heterologous AtWRKY18, AtWRKY40, and AtMYC2 TFs were overexpressed in S. sclarea hairy roots and proved to regulate in a coordinated manner the expression of several genes encoding enzymes of the MEP-dependent pathway, especially DXS, DXR, GGPPS and CPPS. The content of total abietane diterpenes was enhanced in all overexpressing lines, although in a variable manner due to a negative pleiotropic effect on HR growth. Interestingly, in the best performing HR lines overexpressing the AtWRKY40 TF induced a significant 4-fold increase in the final yield of aethiopinone, for which we have reported an interesting anti-proliferative activity against resistant melanoma cells. The present results are also informative and instrumental to enhance the synthesis of abietane diterpenes derived from the plastidial MEP-derived terpenoid pathway in other Salvia species.

De novo sequencing, assembly and characterisation of Aloe vera transcriptome and analysis of expression profiles of genes related to saponin and anthraquinone metabolism

BACKGROUND

Aloe vera is a perennial, succulent, drought-resistant plant that exhibits many pharmacological characteristics such as wound healing ability against skin burns, anti-ulcer, anti-inflammatory, anti-tumor, anti-viral, anti-hypercholesterolemic, anti-hyperglycemic, anti-asthmatic and much more. Despite great medicinal worth, little genomic information is available on Aloe vera. This study is an initiative to explore the full-scale functional genomics of Aloe vera by generating whole transcriptome sequence database, using Illumina HiSeq technology and its progressive annotation specifically with respect to the metabolic specificity of the plant.

RESULTS

Transcriptome sequencing of root and leaf tissue of Aloe vera was performed using Illumina paired-end sequencing technology. De novo assembly of high quality paired-end reads, resulted into 1,61,733 and 2,21,792 transcripts with mean length of 709 and 714 nucleotides for root and leaf respectively. The non-redundant transcripts were clustered using CD-HIT-EST, yielding a total of 1,13,063 and 1,41,310 unigenes for root and leaf respectively. A total of 6114 and 6527 CDS for root and leaf tissue were enriched into 24 different biological pathway categories using KEGG pathway database. DGE profile prepared by calculating FPKM values was analyzed for differential expression of specific gene encoding enzymes involved in secondary metabolite biosynthesis. Sixteen putative genes related to saponin, lignin, anthraquinone, and carotenoid biosynthesis were selected for quantitative expression by real-time PCR. DGE as well as qRT PCR expression analysis represented up-regulation of secondary metabolic genes in root as compared to leaf. Furthermore maximum number of genes was found to be up-regulated after the induction of methyl jasmonate, which stipulates the association of secondary metabolite synthesis with the plant's defense mechanism during stress. Various transcription factors including bHLH, NAC, MYB were identified by searching predicted CDS against PlantTFdb.

CONCLUSIONS

This is the first transcriptome database of Aloe vera and can be potentially utilized to characterize the genes involved in the biosynthesis of important secondary metabolites, metabolic regulation, signal transduction mechanism, understanding function of a particular gene in the biology and physiology of plant of this species as well as other species of Aloe genus.

Transcriptional Profiles of SmWRKY Family Genes and Their Putative Roles in the Biosynthesis of Tanshinone and Phenolic Acids in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge is a Chinese traditional herb for treating cardiovascular and cerebrovascular diseases, and tanshinones and phenolic acids are the dominated medicinal and secondary metabolism constituents of this plant. WRKY transcription factors (TFs) can function as regulators of secondary metabolites biosynthesis in many plants. However, studies on the WRKY that regulate tanshinones and phenolics biosynthesis are limited. In this study, 69 SmWRKYs were identified in the transcriptome database of S. miltiorrhiza, and phylogenetic analysis indicated that some SmWRKYs had closer genetic relationships with other plant WRKYs, which were involved in secondary metabolism. Hairy roots of S. miltiorrhiza were treated by methyl jasmonate (MeJA) to detect the dynamic change trend of SmWRKY, biosynthetic genes, and medicinal ingredients accumulation. Base on those date, a correlation analysis using Pearson’s correlation coefficient was performed to construct gene-to-metabolite network and identify 9 SmWRKYs (SmWRKY1, 7, 19, 29, 45, 52, 56, 58, and 68), which were most likely to be involved in tanshinones and phenolic acids biosynthesis. Taken together, this study has provided a significant resource that could be used for further research on SmWRKY in S. miltiorrhiza and especially could be used as a cue for further investigating SmWRKY functions in secondary metabolite accumulation.

Comparative analysis of the root transcriptomes of cultivated and wild rice varieties in response to Magnaporthe oryzae infection revealed both common and species-specific pathogen responses

BACKGROUND

Magnaporthe oryzae, the causal fungus of rice blast disease, negatively impacts global rice production. Wild rice (Oryza rufipogon), a relative of cultivated rice (O. sativa), possesses unique attributes that enable it to resist pathogen invasion. Although wild rice represents a major resource for disease resistance, relative to current cultivated rice varieties, no prior studies have compared the immune and transcriptional responses in the roots of wild and cultivated rice to M. oryzae.

RESULTS

In this study, we showed that M. oryzae could act as a typical root-infecting pathogen in rice, in addition to its common infection of leaves, and wild rice roots were more resistant to M. oryzae than cultivated rice roots. Next, we compared the differential responses of wild and cultivated rice roots to M. oryzae using RNA-sequencing (RNA-seq) to unravel the molecular mechanisms underlying the enhanced resistance of the wild rice roots. Results indicated that both common and genotype-specific mechanisms exist in both wild and cultivated rice that are associated with resistance to M. oryzae. In wild rice, resistance mechanisms were associated with lipid metabolism, WRKY transcription factors, chitinase activities, jasmonic acid, ethylene, lignin, and phenylpropanoid and diterpenoid metabolism; while the pathogen responses in cultivated rice were mainly associated with phenylpropanoid, flavone and wax metabolism. Although modulations in primary metabolism and phenylpropanoid synthesis were common to both cultivated and wild rice, the modulation of secondary metabolism related to phenylpropanoid synthesis was associated with lignin synthesis in wild rice and flavone synthesis in cultivated rice. Interestingly, while the expression of fatty acid and starch metabolism-related genes was altered in both wild and cultivated rice in response to the pathogen, changes in lipid acid synthesis and lipid acid degradation were dominant in cultivated and wild rice, respectively.

CONCLUSIONS

The response mechanisms to M. oryzae were more complex in wild rice than what was observed in cultivated rice. Therefore, this study may have practical implications for controlling M. oryzae in rice plantings and will provide useful information for incorporating and assessing disease resistance to M. oryzae in rice breeding programs.

Overexpression of the alfalfa WRKY11 gene enhances salt tolerance in soybean

The WRKY transcription factors play an important role in the regulation of transcriptional reprogramming associated with plant abiotic stress responses. In this study, the WRKY transcription factor MsWRKY11, containing the plant-specific WRKY zinc finger DNA-binding motif, was isolated from alfalfa. The MsWRKY11 gene was detected in all plant tissues (root, stem, leaf, flower, and fruit), with high expression in root and leaf tissues. MsWRKY11 was upregulated in response to a variety of abiotic stresses, including salinity, alkalinity, cold, abscisic acid, and drought. Overexpression of MsWRKY11 in soybean enhanced the salt tolerance at the seedling stage. Transgenic soybean had a better salt-tolerant phenotype, and the hypocotyls were significantly longer than those of wild-type seeds after salt treatment. Furthermore, MsWRKY11 overexpression increased the contents of chlorophyll, proline, soluble sugar, superoxide dismutase, and catalase, but reduced the relative electrical conductivity and the contents of malonaldehyde, H2O2, and O2-. Plant height, pods per plant, seeds per plant, and 100-seed weight of transgenic MsWRKY11 soybean were higher than those of wild-type soybean, especially OX2. Results of the salt experiment showed that MsWRKY11 is involved in salt stress responses, and its overexpression improves salt tolerance in soybean.

WRKY23 is a component of the transcriptional network mediating auxin feedback on PIN polarity

Auxin is unique among plant hormones due to its directional transport that is mediated by the polarly distributed PIN auxin transporters at the plasma membrane. The canalization hypothesis proposes that the auxin feedback on its polar flow is a crucial, plant-specific mechanism mediating multiple self-organizing developmental processes. Here, we used the auxin effect on the PIN polar localization in Arabidopsis thaliana roots as a proxy for the auxin feedback on the PIN polarity during canalization. We performed microarray experiments to find regulators of this process that act downstream of auxin. We identified genes that were transcriptionally regulated by auxin in an AXR3/IAA17- and ARF7/ARF19-dependent manner. Besides the known components of the PIN polarity, such as PID and PIP5K kinases, a number of potential new regulators were detected, among which the WRKY23 transcription factor, which was characterized in more detail. Gain- and loss-of-function mutants confirmed a role for WRKY23 in mediating the auxin effect on the PIN polarity. Accordingly, processes requiring auxin-mediated PIN polarity rearrangements, such as vascular tissue development during leaf venation, showed a higher WRKY23 expression and required the WRKY23 activity. Our results provide initial insights into the auxin transcriptional network acting upstream of PIN polarization and, potentially, canalization-mediated plant development.

The plant hormone auxin belongs to the major plant-specific developmental regulators. It mediates or modifies almost all aspects of plant life. One of the fascinating features of the auxin action is its directional movement between cells, whose direction can be regulated by auxin signaling itself. This plant-specific feedback regulation has been proposed decades ago and allows for the self-organizing formation of distinct auxin channels shown to be crucial for processes, such as the regular pattern formation of leaf venation, organ formation, and regeneration of plant tissues. Despite the prominent importance of this so called auxin canalization process, the insight into the underlying molecular mechanism is very limited. Here, we identified a number of genes that are transcriptionally regulated and act downstream of the auxin signaling to mediate the auxin feedback on the polarized auxin transport. One of them is the WRKY23 transcription factor that has previously been unsuspected to play a role in this process. Our work provides the first insights into the transcriptional regulation of the auxin canalization and opens multiple avenues to further study this crucial process.

Identification, Expression, and Functional Analysis of the Group IId WRKY Subfamily in Upland Cotton (Gossypium hirsutum L.)

WRKY transcription factors have diverse functions in regulating stress response, leaf senescence, and plant growth and development. However, knowledge of the group IId WRKY subfamily in cotton is largely absent. This study identified 34 group IId WRKY genes in the Gossypium hirsutum genome, and their genomic loci were investigated. Members clustered together in the phylogenetic tree had similar motif compositions and gene structural features, revealing similarity and conservation within group IId WRKY genes. During the evolutionary process, 14 duplicated genes appeared to undergo purification selection. Public RNA-seq data were used to examine the expression patterns of group IId WRKY genes in various tissues and under drought and salt stress conditions. Ten highly expressed genes were identified, and the ten candidate genes revealed distinct expression patterns under drought and salt treatments by qRT-PCR analysis. Among them, Gh_A11G1801 was used for functional characterization. GUS activity was differentially induced by various stresses in Gh_A11G1801p::GUS transgenic Arabidopsis plants. The virus-induced gene silencing (VIGS) of Gh_A11G1801 resulted in drought sensitivity in cotton plants, which was accompanied by elevated malondialdehyde (MDA) content and reduced catalase (CAT) content. Taken together, these findings obtained in this study provide valuable resources for further studying group IId WRKY genes in cotton. Our results also enrich the gene resources for the genetic improvements of cotton varieties that are suitable for growth in stressful conditions.

A Novel WRKY Transcription Factor, MuWRKY3 (Macrotyloma uniflorum Lam. Verdc.) Enhances Drought Stress Tolerance in Transgenic Groundnut (Arachis hypogaea L.) Plants

Drought stress has adverse effects on growth, water relations, photosynthesis and yield of groundnut. WRKY transcription factors (TFs) are the plant-specific TFs which regulate several down-stream stress-responsive genes and play an essential role in plant biotic and abiotic stress responses. We found that WRKY3 gene is highly up-regulated under drought stress conditions and therefore isolated a new WRKY3TF gene from a drought-adapted horsegram (Macrotyloma uniflorum Lam. Verdc.). Conserved domain studies revealed that protein encoded by this gene contains highly conserved regions of two WRKY domains and two C2H2 zinc-finger motifs. The fusion protein localization studies of transient MuWRKY3-YFP revealed its nuclear localization. Overexpression of MuWRKY3 TF gene in groundnut (Arachis hypogaea L.) showed increased tolerance to drought stress compared to wild-type (WT) plants. MuWRKY3 groundnut transgenics displayed lesser and delayed wilting symptoms than WT plants after 10-days of drought stress imposition. The transgenic groundnut plants expressing MuWRKY3 showed less accumulation of malondialdehyde, hydrogen peroxide (H₂O₂), and superoxide anion (O₂^∙-), accompanied by more free proline, total soluble sugar content, and activities of antioxidant enzymes than WT plants under drought stress. Moreover, a series of stress-related LEA, HSP, MIPS, APX, SOD, and CAT genes found up-regulated in the transgenic groundnut plants. The study demonstrates that nuclear-localized MuWRKY3 TF regulates the expression of stress-responsive genes and the activity of ROS scavenging enzymes which results in improved drought tolerance in groundnut. We conclude that MuWRKY3 may serve as a new putative candidate gene for the improvement of stress resistance in plants.

Characterization of Soybean WRKY Gene Family and Identification of Soybean WRKY Genes that Promote Resistance to Soybean Cyst Nematode

WRKY proteins are a superfamily of plant transcription factors with important roles in plants. WRKY proteins have been extensively analyzed in plant species including Arabidopsis and rice. Here we report characterization of soybean WRKY gene family and their functional analysis in resistance to soybean cyst nematode (SCN), the most important soybean pathogen. Through search of the soybean genome, we identified 174 genes encoding WRKY proteins that can be classified into seven groups as established in other plants. WRKY variants including a WRKY-related protein unique to legumes have also been identified. Expression analysis reveals both diverse expression patterns in different soybean tissues and preferential expression of specific WRKY groups in certain tissues. Furthermore, a large number of soybean WRKY genes were responsive to salicylic acid. To identify soybean WRKY genes that promote soybean resistance to SCN, we first screened soybean WRKY genes for enhancing SCN resistance when over-expressed in transgenic soybean hairy roots. To confirm the results, we transformed five WRKY genes into a SCN-susceptible soybean cultivar and generated transgenic soybean lines. Transgenic soybean lines overexpressing three WRKY transgenes displayed increased resistance to SCN. Thus, WRKY genes could be explored to develop new soybean cultivars with enhanced resistance to SCN.

Transcriptome analysis of WRKY gene family in Oryza officinalis Wall ex Watt and WRKY genes involved in responses to Xanthomonas oryzae pv. oryzae stress

Oryza officinalis Wall ex Watt, a very important and special wild rice species, shows abundant genetic diversity and disease resistance features, especially high resistance to bacterial blight. The molecular mechanisms of bacterial blight resistance in O. officinalis have not yet been elucidated. The WRKY transcription factor family is one of the largest gene families involved in plant growth, development and stress response. However, little is known about the numbers, structure, molecular phylogenetics, and expression of the WRKY genes under Xanthomonas oryzae pv. oryzae (Xoo) stress in O. officinalis due to lacking of O. officinalis genome. Therefore, based on the RNA-sequencing data of O. officinalis, we performed a comprehensive study of WRKY genes in O. officinalis and identified 89 OoWRKY genes. Then 89 OoWRKY genes were classified into three groups based on the WRKY domains and zinc finger motifs. Phylogenetic analysis strongly supported that the evolution of OoWRKY genes were consistent with previous studies of WRKYs, and subgroup IIc OoWRKY genes were the original ancestors of some group II and group III OoWRKYs. Among the 89 OoWRKY genes, eight OoWRKYs displayed significantly different expression (>2-fold, p<0.01) in the O. officinalis transcriptome under Xoo strains PXO99 and C5 stress 48 h, suggesting these genes might play important role in PXO99 and C5 stress responses in O. officinalis. QRT-PCR analysis and confirmation of eight OoWRKYs expression patterns revealed that they responded strongly to PXO99 and C5 stress 24 h, 48 h, and 72 h, and the trends of these genes displaying marked changes were consistent with the 48 h RNA-sequencing data, demonstrated these genes played important roles in response to biotic stress and might even involved in the bacterial blight resistance. Tissue expression profiles of eight OoWRKY genes revealed that they were highly expressed in root, stem, leaf, and flower, especially in leaf (except OoWRKY71), suggesting these genes might be also important for plant growth and organ development. In this study, we analyzed the WRKY family of transcription factors in O.officinalis. Insight was gained into the classification, evolution, and function of the OoWRKY genes, revealing the putative roles of eight significantly different expression OoWRKYs in Xoo strains PXO99 and C5 stress responses in O.officinalis. This study provided a better understanding of the evolution and functions of O. officinalis WRKY genes, and suggested that manipulating eight significantly different expression OoWRKYs would enhance resistance to bacterial blight.

RNA sequencing of Brassica napus reveals cellular redox control of Sclerotinia infection

Brassica napus is one of the world's most valuable oilseeds and is under constant pressure by the necrotrophic fungal pathogen, Sclerotinia sclerotiorum, the causal agent of white stem rot. Despite our growing understanding of host pathogen interactions at the molecular level, we have yet to fully understand the biological processes and underlying gene regulatory networks responsible for determining disease outcomes. Using global RNA sequencing, we profiled gene activity at the first point of infection on the leaf surface 24 hours after pathogen exposure in susceptible (B. napus cv. Westar) and tolerant (B. napus cv. Zhongyou 821) plants. We identified a family of ethylene response factors that may contribute to host tolerance to S. sclerotiorum by activating genes associated with fungal recognition, subcellular organization, and redox homeostasis. Physiological investigation of redox homeostasis was further studied by quantifying cellular levels of the glutathione and ascorbate redox pathway and the cycling enzymes associated with host tolerance to S. sclerotiorum. Functional characterization of an Arabidopsis redox mutant challenged with the fungus provides compelling evidence into the role of the ascorbate-glutathione redox hub in the maintenance and enhancement of plant tolerance against fungal pathogens.

Genome-wide characterization of the WRKY gene family in radish (Raphanus sativus L.) reveals its critical functions under different abiotic stresses

The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. The WRKY is among the largest transcription factors (TFs) associated with multiple biological activities for plant survival, including control response mechanisms against abiotic stresses such as heat, salinity, and heavy metals. Radish is an important root vegetable crop and therefore characterization and expression pattern investigation of WRKY transcription factors in radish is imperative. In the present study, 126 putative WRKY genes were retrieved from radish genome database. Protein sequence and annotation scrutiny confirmed that RsWRKY proteins possessed highly conserved domains and zinc finger motif. Based on phylogenetic analysis results, RsWRKYs candidate genes were divided into three groups (Group I, II and III) with the number 31, 74, and 20, respectively. Additionally, gene structure analysis revealed that intron-exon patterns of the WRKY genes are highly conserved in radish. Linkage map analysis indicated that RsWRKY genes were distributed with varying densities over nine linkage groups. Further, RT-qPCR analysis illustrated the significant variation of 36 RsWRKY genes under one or more abiotic stress treatments, implicating that they might be stress-responsive genes. In total, 126 WRKY TFs were identified from the R. sativus genome wherein, 35 of them showed abiotic stress-induced expression patterns. These results provide a genome-wide characterization of RsWRKY TFs and baseline for further functional dissection and molecular evolution investigation, specifically for improving abiotic stress resistances with an ultimate goal of increasing yield and quality of radish.

Transcriptome analysis of terpene chemotypes of Melaleuca alternifolia across different tissues

Plant chemotypes or chemical polymorphisms are defined by discrete variation in secondary metabolites within a species. This variation can have consequences for ecological interactions or the human use of plants. Understanding the molecular basis of chemotypic variation can help to explain how variation of plant secondary metabolites is controlled. We explored the transcriptomes of the 3 cardinal terpene chemotypes of Melaleuca alternifolia in young leaves, mature leaves, and stem and compared transcript abundance to variation in the constitutive profile of terpenes. Leaves from chemotype 1 plants (dominated by terpinen-4-ol) show a similar pattern of gene expression when compared to chemotype 5 plants (dominated by 1,8-cineole). Only terpene synthases in young leaves were differentially expressed between these chemotypes, supporting the idea that terpenes are mainly synthetized in young tissue. Chemotype 2 plants (dominated by terpinolene) show a greater degree of differential gene expression compared to the other chemotypes, which might be related to the isolation of plant populations that exhibit this chemotype and the possibility that the terpinolene synthase gene in M. alternifolia was derived by introgression from a closely related species, Melaleuca trichostachya. By using multivariate analyses, we were able to associate terpenes with candidate terpene synthases.