Characterization of WRKY co-regulatory networks in rice and Arabidopsis

Ecophysiological and transcriptional landscapes of arbuscular mycorrhiza fungi enhancing yield, quality, and stalk rot resistance in Anoectochilus roxburghii

Anoectochilus roxburghii (Wall.) Lindl. (A. roxburghii) is an increasingly popular medicinal herb. Arbuscular mycorrhiza (AM) fungi, known for their symbiotic relationships with plant roots, enhance nutrient uptake and disease resistance in host plants. However, their specific regulatory mechanisms in A. roxburghii are not fully understood. In this study, Fujian A. roxburghii was inoculated with the AM fungus Glomus intraradices, and successful root colonization was observed. Following AM fungal colonization, there was a significant upregulation of photosynthesis-related genes in the stems, accompanied by improved canopy phenotypes and root architecture. Consequently, AM-inoculated plants exhibited increased fresh and dry biomass, as well as elevated levels of polysaccharides and flavonoids. Additionally, the incidence of Fusarium oxysporum-induced stalk rot was reduced in AM-inoculated plants. Analysis of defense-related enzymes indicated that AM-inoculated plants exhibited a rapid and robust response to pathogen infection, mitigating oxidative stress. Transcriptomic analysis revealed significant upregulation of genes associated "Fatty acid degradation", "MAPK signaling pathway-plant", and "Plant-pathogen interaction", suggesting their involvement in enhanced disease resistance. A regulatory network centered on ACX1 and calmodulin, involving multiple transcription factors such as WRKY, bHLH, ERF, NAC, and HSF, was implicated in defense responses. These findings demonstrated the beneficial effects of AM fungi on yield, quality, and disease resistance in A. roxburghii, providing a theoretical foundation for its cultivation and genetic improvement.

Genome-wide characterization of the WRKY gene family and the role of LsfWRKY29 in regulating somatic embryogenesis in hybrid sweetgum (Liquidambar styraciflua × L. formosana)

Hybrid sweetgum (Liquidambar styraciflua × L. formosana) is a globally significant forest tree resource, exhibiting significant economic, ornamental, ecological and medicinal values. Somatic embryogenesis (SE) is an effective reproductive strategy, having great application potential and economic value in large-scale propagation, artificial seed production, genetic transformation, germplasm preservation and biotechnology. It's an essential technical approach for the industrial application of hybrid sweetgum. The WRKY gene family, one of the largest groups in plants, is crucial for regulating responses related to plant growth, while its function and molecular mechanism in SE are unknown. Here, we identified 61 members of the WRKY gene family from the genome of hybrid sweetgum and characterized their amino acid sequence characteristics, phylogenetic relationships, conserved motifs, cis-acting elements and gene structures at the genome-wide level. LsfWRKY29, identified as a key gene in Weighted Gene Co-expression Network Analysis (WGCNA), is significantly expressed during the morphogenesis stage of SE. Subcellular localization studies indicate that WRKY29 is located in the nucleus. Functional assays indicate that overexpression of LsfWRKY29 leads to more abnormal embryos and subsequent developmental arrest, suggesting its negative role during SE. Yeast two-hybrid (Y2H) and luciferase complementation imaging (LCI) assays validate the interaction between LsfWRKY29 and the growth regulator LsfGRF2, a member of the growth regulators that has been extensively proven to boost plant regeneration. This interaction implies that the WRKY29-GRF2 might be an important regulatory module in SE of hybrid sweetgum. These findings provide new insights into the role of WRKY genes in regulating SE in plants.

Genome-wide identification and expression profiling of WRKY gene family in grain Amaranth (Amaranthus hypochondriacus L.) under salinity and drought stresses

BACKGROUND

The WRKY gene family plays a significant role in plant growth, development, and responses to biotic and abiotic stresses. However, the role of the WRKY gene family has not been reported in Amaranthus hypochondriacus. This study presents a comprehensive genome-wide analysis of the WRKY gene family in grain amaranth (A. hypochondriacus L.), a resilient crop known for its high nutritional value and adaptability to challenging environments.

RESULTS

In this study, 55 WRKY genes (AhyWRKY1-55) were identified in A. hypochondriacus and distributed unevenly across 16 scaffolds. Of these, 50 contained conserved WRKY domains and were classified into three main groups. Group II was further divided into five subgroups (IIa-IIe) based on phylogenetic analysis, with each clade being well supported by conserved motifs. Additionally, the gene structure analysis revealed variations in exon-intron organization. In contrast, motif analysis showed the presence of conserved domains that were similar within the group but differed between groups, suggesting their functional diversity. Cis-acting elements related to plant growth and development and light, hormones, and stress responses were identified. Synteny analysis revealed that 34 (61.8%) of the genes originated from tandem duplication, indicating the role of tandem duplication in the expansion of the A. hypochondriacus WRKY gene family. Protein-protein interaction analysis suggested that AhyWRKY3, AhyWRKY27, AhyWRKY28, AhyWRKY36, and AhyWRKY52 were hub genes involved in the complex protein interaction network. Using in silico and real-time quantitative PCR, expression analysis revealed tissue- and condition-specific expression patterns of AhyWRKY genes. Notably, under drought stress, AhyWRKY39, AhyWRKY40, AhyWRKY54, and AhyWRKY01 showed increased expression, while under salt stress, AhyWRKY40, AhyWRKY54, AhyWRKY39, AhyWRKY49, and AhyWRKY8 were upregulated at 30 days, suggesting that these genes may play key role in response to salinity stress.

CONCLUSIONS

The present study provides valuable insights into the organization and evolutionary patterns of the WRKY gene family in amaranth. It also identifies putative candidate WRKY genes that may play a role in conferring drought and salt tolerance. Overall, this study lays a foundation for further functional validation of these WRKY candidate genes, facilitating their exploitation in the amaranth genetic improvement programs to develop stress-resilient varieties.

Genome-wide identification and characterization of WRKYs family involved in responses to Cylindrocarpon destructans in Panax notoginseng

BACKGROUND

WRKY transcription factors (TFs) are key regulators of plant responses to biotic and abiotic stresses. Previous studies demonstrated that the role of WRKY TFs play in the disease resistance of Panax notoginseng, the causal agent of root rot disease. However, comprehensive genome-wide analyses of WRKY genes in this species remain scarce.

RESULT

We identified 79 WRKY genes in the P. notoginseng genome, classifing them into three groups based on structural features and phylogenetic relationships: Class I (14 genes), Class II (55 genes), and Class III (10 genes). Of these, 58 PnWRKY genes were mapped to the P. notoginseng chromosomes and showed collinearity with Arabidopsis thaliana, Daucus carota, and three Solanaceae species. Expression analysis revealed that 53 PnWRKY genes were actively transcribed across various tissues, including roots, flowers, stems, rhizomes, and different root parts. Furthermore, PnWRKY genes responded to Cylindrocarpon destructans infection and were induced by jasmonic acid (JA) and salicylic acid (SA). Notably, the ectopic expression of PnWRKY35 in tobacco enhanced resistance to C. destructans, accompanied by increased levels of gibberellins (GA24, GA3, GA8), JA, SA, and brassinolide. KEGG enrichment analysis of 547 differentially expressed genes (DEGs) indicated their involvement pathways related to disease resistance, including the biosynthesis of sesquiterpenes and triterpenes, JA biosynthesis, plant-pathogen interactions and the MAPK signaling pathway.

CONCLUSION

This study demonstrated that the WRKY family in P. notoginseng plays a significant role in resistance to root rot disease. PnWRKY genes are responsive to MeJA and SA induction as well as C. destructans infection. Moreover, ectopic expression of PnWRKY35 activates multiple plant disease resistance pathways, increases phytohormone levels, and enhance resistance to C. destructans. These findings provide a foundation for future exploration of the mechanism underlying P. notoginseng resistance to root rot disease.

Genome-wide identification of the WRKY gene family in blueberry (Vaccinium spp.) and expression analysis under abiotic stress

Introduction

The WRKY transcription factor (TF) family is one of the largest TF families in plants and is widely involved in responses to both biotic and abiotic stresses.

Methods

To clarify the function of the WRKY family in blueberries, this study identified the WRKY genes within the blueberry genome and systematically analyzed gene characteristics, phylogenetic evolution, promoter cis-elements, expression patterns, and subcellular localization of the encoded products.

Results

In this study, 57 VcWRKY genes were identified, and all encoding products had a complete WRKY heptapeptide structure and zinc-finger motif. The VcWRKY genes were divided into three subgroups (I-III) by phylogenetic analysis. Group II was divided into five subgroups: IIa, IIb, IIc, IId, and IIe. 57 VcWRKY genes were distributed unevenly across 32 chromosomes. The amino acids ranged from 172 to 841, and molecular weights varied from 19.75 to 92.28 kD. Intra-group syntenic analysis identified 12 pairs of duplicate segments. Furthermore, 34 cis-element recognition sites were identified in the promoter regions of VcWRKY genes, primarily comprising phytohormone-responsive and light-responsive elements. Comparative syntenic maps were generated to investigate the evolutionary relationships of VcWRKY genes, revealing the closest homology to dicotyledonous WRKY gene families. VcWRKY genes were predominantly expressed in the fruit flesh and roots of blueberries. Gene expression analysis showed that the responses of VcWRKY genes to stress treatments were more strongly in leaves than in roots. Notably, VcWRKY13 and VcWRKY25 exhibited significant upregulation under salt stress, alkali stress, and saline-alkali stress, and VcWRKY1 and VcWRKY13 showed notable induction under drought stress. Subcellular localization analysis confirmed that VcWRKY13 and VcWRKY25 function within the nucleus.

Conclusion

These findings establish a foundation for further investigation into the functions and regulatory mechanisms of VcWRKY genes and provide guidance for selecting stress-tolerant genes in the development of blueberry cultivars.

Comparative transcriptome profiling and co-expression network analysis uncover the key genes associated with pear petal defense responses against Monilinia laxa infection

Pear brown rot and blossom blight caused by Monilinia laxa seriously affect pear production worldwide. Here, we compared the transcriptomic profiles of petals after inoculation with M. laxa using two pear cultivars with different levels of sensitivity to disease (Sissy, a relatively tolerant cultivar, and Kristalli, a highly susceptible cultivar). Physiological indexes were also monitored in the petals of both cultivars at 2 h and 48 h after infection (2 HAI and 48 HAI). RNA-seq data and weighted gene co-expression network analysis (WGCNA) allowed the identification of key genes and pathways involved in immune- and defense-related responses that were specific for each cultivar in a time-dependent manner. In particular, in the Kristalli cultivar, a significant transcriptome reprogramming occurred early at 2 HAI and was accompanied either by suppression of key differentially expressed genes (DEGs) involved in the modulation of any defense responses or by activation of DEGs acting as sensitivity factors promoting susceptibility. In contrast to the considerably high number of DEGs induced early in the Kristalli cultivar, upregulation of specific DEGs involved in pathogen perception and signal transduction, biosynthesis of secondary and primary metabolism, and other defense-related responses was delayed in the Sissy cultivar, occurring at 48 HAI. The WGCNA highlighted one module that was significantly and highly correlated to the relatively tolerant cultivar. Six hub genes were identified within this module, including three WRKY transcription factor-encoding genes: WRKY 65 (pycom05g27470), WRKY 71 (pycom10g22220), and WRKY28 (pycom17g13130), which may play a crucial role in enhancing the tolerance of pear petals to M. laxa. Our results will provide insights into the interplay of the molecular mechanisms underlying immune responses of petals at the pear-M. laxa pathosystem.

Comparative Transcriptomics and Metabolomics Analyses of Avicennia marina and Kandelia obovata under Chilling Stress during Seedling Stage

One of the most productive ecosystems in the world, mangroves are susceptible to cold stress. However, there is currently insufficient knowledge of the adaptation mechanisms of mangrove plants in response to chilling stress. This study conducted a comparative analysis of transcriptomics and metabolomics to investigate the adaptive responses of Kandelia obovata (chilling-tolerant) and Avicennia marina (chilling-sensitive) to 5 °C. The transcriptomics results revealed that differentially expressed genes (DEGs) were mostly enriched in signal transduction, photosynthesis-related pathways, and phenylpropanoid biosynthesis. The expression pattern of genes involved in photosynthesis-related pathways in A. marina presented a downregulation of most DEGs, which correlated with the decrease in total chlorophyll content. In the susceptible A. marina, all DEGs encoding mitogen-activated protein kinase were upregulated. Phenylpropanoid-related genes were observed to be highly induced in K. obovata. Additionally, several metabolites, such as 4-aminobutyric acid, exhibited higher levels in K. obovata than in A. marina, suggesting that chilling-tolerant varieties regulated more metabolites in response to chilling. The investigation defined the inherent distinctions between K. obovata and A. marina in terms of signal transduction gene expression, as well as phenylpropanoid and flavonoid biosynthesis, during exposure to low temperatures.

QTLs and Candidate Loci Associated with Drought Tolerance Traits of Kaybonnet x ZHE733 Recombinant Inbred Lines Rice Population

Rice is the most important staple crop for the sustenance of the world's population, and drought is a major factor limiting rice production. Quantitative trait locus (QTL) analysis of drought-resistance-related traits was conducted on a recombinant inbred line (RIL) population derived from the self-fed progeny of a cross between the drought-resistant tropical japonica U.S. adapted cultivar Kaybonnet and the drought-sensitive indica cultivar ZHE733. K/Z RIL population of 198 lines was screened in the field at Fayetteville (AR) for three consecutive years under controlled drought stress (DS) and well-watered (WW) treatment during the reproductive stage. The effects of DS were quantified by measuring morphological traits, grain yield components, and root architectural traits. A QTL analysis using a set of 4133 single nucleotide polymorphism (SNP) markers and the QTL IciMapping identified 41 QTLs and 184 candidate genes for drought-related traits within the DR-QTL regions. RT-qPCR in parental lines was used to confirm the putative candidate genes. The comparison between the drought-resistant parent (Kaybonnet) and the drought-sensitive parent (ZHE733) under DS conditions revealed that the gene expression of 15 candidate DR genes with known annotations and two candidate DR genes with unknown annotations within the DR-QTL regions was up-regulated in the drought-resistant parent (Kaybonnet). The outcomes of this research provide essential information that can be utilized in developing drought-resistant rice cultivars that have higher productivity when DS conditions are prevalent.

Genome-Wide Identification and Analysis of the WRKY Gene Family in Asparagus officinalis

In recent years, the related research of the WRKY gene family has been gradually promoted, which is mainly reflected in the aspects of environmental stress and hormone response. However, to make the study of the WRKY gene family more complete, we also need to focus on the whole-genome analysis and identification of the family. In previous studies, the whole WRKY gene family of Arabidopsis, legumes and other plants has been thoroughly studied. However, since the publication of Asparagus officinalis genome-wide data, there has never been an analysis of the whole WRKY gene family. To understand more broadly the function of the WRKY gene family, the whole genome and salt stress transcriptome data of asparagus were used for comprehensive analysis in this study, including WRKY gene family identification, phylogenetic tree construction, analysis of conserved mods and gene domains, extraction of cis-acting elements, intron/exon analysis, species collinearity analysis, and WRKY expression analysis under salt stress. The results showed that a total of 70 genes were selected and randomly distributed on 10 chromosomes and one undefined chromosome. According to the functional classification of Arabidopsis thaliana, the WRKY family of asparagus was divided into 11 subgroups (C1-C9, U1, U2). It is worth considering that the distribution rules of gene-conserved motifs, gene domains and introns/exons in the same subfamily are similar, which suggests that genes in the same subfamily may regulate similar physiological processes. In this study, 11 cis-acting elements of WRKY family were selected, among which auxin, gibberellin, abscisic acid, salicylic acid and other hormone-regulated induction elements were involved. In addition, environmental stress (such as drought stress and low-temperature response) also accounted for a large proportion. Interestingly, we analyzed a total of two tandem duplicate genes and 13 segmental duplication genes, suggesting that this is related to the amplification of the WRKY gene family. Transcriptome data analysis showed that WRKY family genes could regulate plant growth and development by up-regulating and down-regulating gene expression under salt stress. Volcanic maps showed that 3 and 15 AoWRKY genes were significantly up-regulated or down-regulated in NI&NI+S and AMF&AMF+S, respectively. These results provide a new way to analyze the evolution and function of the WRKY gene family, and can provide a reference for the production and research of asparagus.

Genome-wide identification and analysis of the WRKY gene family and low-temperature stress response in Prunus sibirica

BACKGROUND

WRKY transcription factors are a prominent gene family in plants, playing a crucial role in various biological processes including development, metabolism, defense, differentiation, and stress response. Although the WRKY gene family has been extensively studied and analysed in numerous plant species, research on Prunus sibirica's WRKY genes (PsWRKY) remains lacking.

RESULTS

This study analysed the basic physicochemical properties, phylogeny, gene structure, cis-acting elements, and Gene ontology (GO) annotation of PsWRKY gene family members using bioinformatics methods based on the whole-genome data of P. sibirica. In total, 55 WRKYs were identified in P. sibirica and were heterogeneously distributed on eight chromosomes. Based on the phylogenetic analysis, these WRKYs were classified into three major groups: Group I, Group II (II-a, II-b, II-c, II-d, II-e), and Group III. Members of different subfamilies have different cis-acting elements, conserved motifs, and intron-exon structures, indicating functional heterogeneity of the WRKY family. Prediction of subcellular localisation indicated that PsWRKYs were mainly located in the nucleus. Twenty pairs of duplicated genes were identified, and segmental duplication events may play an important role in PsWRKY gene family expansion. Analysis of the Ka/Ks ratio showed that the PsWRKY family's homologous genes were primarily purified by selection. Additionally, GO annotation analysis showed that the WRKY gene family was mainly involved in responses to stimuli, immune system processes, and reproductive processes. Furthermore, quantitative real-time PCR (qRT-PCR) analysis showed that 23 PsWRKYs were highly expressed in one or more tissues (pistils and roots) and PsWRKYs showed specific expression patterns under different low-temperature stress conditions.

CONCLUSIONS

Our results provide a scientific basis for the further exploration and functional validation of WRKYs in P. sibirica.

Comparative transcriptomic analysis of germinating rice seedlings to individual and combined anaerobic and cold stress

BACKGROUND

Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination.

RESULTS

The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein).

CONCLUSION

A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.

OsLPXC negatively regulates tolerance to cold stress via modulating oxidative stress, antioxidant defense and JA accumulation in rice

Exposure of crops to low temperature (LT) during emerging and reproductive stages influences their growth and development. In this study, we have isolated a cold induced, nucleus-localized lipid A gene from rice named OsLPXC, which encodes a protein of 321 amino acids. Knockout of OsLPXC resulted in enhance sensitivity to LT stress in rice, with increased accumulation of reactive oxygen species (ROS), malondialdehyde and electrolyte leakage, while expression and activities of antioxidant enzymes were significantly suppressed. The accumulation of chlorophyll content and net photosynthetic rate of knockout plants were also decreased compared with WT under LT stress. The functional analysis of differentially expressed genes (DEGs), showed that numerous genes associated with antioxidant defense, photosynthesis, cold signaling were solely expressed and downregulated in oslpxc plants compared with WT under LT. The accumulation of methyl jasmonate (MeJA) in leave and several DEGs related to the jasmonate biosynthesis pathway were significantly downregulated in OsLPXC knockout plants, which showed differential levels of MeJA regulation in WT and knockout plants in response to cold stress. These results indicated that OsLPXC positively regulates cold tolerance in rice via stabilizing the expression and activities of ROS scavenging enzymes, photosynthetic apparatus, cold signaling genes, and jasmonate biosynthesis.

Overexpression of GmWRKY172 enhances cadmium tolerance in plants and reduces cadmium accumulation in soybean seeds

INTRODUCTION

Cadmium (Cd) stress is a significant threat to soybean production, and enhancing Cd tolerance in soybean is the focus of this study. The WRKY transcription factor family is associated with abiotic stress response processes. In this study, we aimed to identify a Cd-responsive WRKY transcription factor GmWRKY172 from soybean and investigate its potential for enhancing Cd tolerance in soybean.

METHODS

The characterization of GmWRKY172 involved analyzing its expression pattern, subcellular localization, and transcriptional activity. To assess the impact of GmWRKY172, transgenic Arabidopsis and soybean plants were generated and examined for their tolerance to Cd and Cd content in shoots. Additionally, transgenic soybean plants were evaluated for Cd translocation and various physiological stress indicators. RNA sequencing was performed to identify the potential biological pathways regulated by GmWRKY172.

RESULTS

GmWRKY172 was significantly upregulated by Cd stress, highly expressed in leaves and flowers, and localized to the nucleus with transcriptional activity. Transgenic plants overexpressing GmWRKY172 showed enhanced Cd tolerance and reduced Cd content in shoots compared to WT. Lower Cd translocation from roots to shoots and seeds was also observed in transgenic soybean. Under Cd stress, transgenic soybean accumulated less malondialdehyde (MDA) and hydrogen peroxide (H₂O₂) than WT plants, with higher flavonoid and lignin contents, and peroxidase (POD) activity. RNA sequencing analysis revealed that many stress-related pathways were regulated by GmWRKY172 in transgenic soybean, including flavonoid biosynthesis, cell wall synthesis, and peroxidase activity.

DISCUSSION

Our findings demonstrated that GmWRKY172 enhances Cd tolerance and reduces seed Cd accumulation in soybean by regulating multiple stress-related pathways, and could be a promising candidate for breeding Cd-tolerant and low Cd soybean varieties.

Genome-Wide Identification, Evolutionary and Functional Analyses of WRKY Family Members in Ginkgo biloba

WRKY transcription factors (TFs) are one of the largest families in plants which play essential roles in plant growth and stress response. Ginkgo biloba is a living fossil that has remained essentially unchanged for more than 200 million years, and now has become widespread worldwide due to the medicinal active ingredients in its leaves. Here, 37 WRKY genes were identified, which were distributed randomly in nine chromosomes of G. biloba. Results of the phylogenetic analysis indicated that the GbWRKY could be divided into three groups. Furthermore, the expression patterns of GbWRKY genes were analyzed. Gene expression profiling and qRT-PCR revealed that different members of GbWRKY have different spatiotemporal expression patterns in different abiotic stresses. Most of the GbWRKY genes can respond to UV-B radiation, drought, high temperature and salt treatment. Meanwhile, all GbWRKY members performed phylogenetic tree analyses with the WRKY proteins of other species which were known to be associated with abiotic stress. The result suggested that GbWRKY may play a crucial role in regulating multiple stress tolerances. Additionally, GbWRKY13 and GbWRKY37 were all located in the nucleus, while GbWRKY15 was located in the nucleus and cytomembrane.

A Network-Based Approach for Improving Annotation of Transcription Factor Functions and Binding Sites in Arabidopsis thaliana

Transcription factors are an integral component of the cellular machinery responsible for regulating many biological processes, and they recognize distinct DNA sequence patterns as well as internal/external signals to mediate target gene expression. The functional roles of an individual transcription factor can be traced back to the functions of its target genes. While such functional associations can be inferred through the use of binding evidence from high-throughput sequencing technologies available today, including chromatin immunoprecipitation sequencing, such experiments can be resource-consuming. On the other hand, exploratory analysis driven by computational techniques can alleviate this burden by narrowing the search scope, but the results are often deemed low-quality or non-specific by biologists. In this paper, we introduce a data-driven, statistics-based strategy to predict novel functional associations for transcription factors in the model plant Arabidopsis thaliana. To achieve this, we leverage one of the largest available gene expression compendia to build a genome-wide transcriptional regulatory network and infer regulatory relationships among transcription factors and their targets. We then use this network to build a pool of likely downstream targets for each transcription factor and query each target pool for functionally enriched gene ontology terms. The results exhibited sufficient statistical significance to annotate most of the transcription factors in Arabidopsis with highly specific biological processes. We also perform DNA binding motif discovery for transcription factors based on their target pool. We show that the predicted functions and motifs strongly agree with curated databases constructed from experimental evidence. In addition, statistical analysis of the network revealed interesting patterns and connections between network topology and system-level transcriptional regulation properties. We believe that the methods demonstrated in this work can be extended to other species to improve the annotation of transcription factors and understand transcriptional regulation on a system level.

Genome and Transcriptome-Wide Analysis of OsWRKY and OsNAC Gene Families in Oryza sativa and Their Response to White-Backed Planthopper Infestation

Plants are threatened by a wide variety of herbivorous insect assaults, and display a variety of inherent and induced defenses that shield them against herbivore attacks. Looking at the massive damage caused by the white-backed planthopper (WBPH), Sogatella furcifera, we undertook a study to identify and functionally annotate OsWRKY and OsNAC transcription factors (TFs) in rice, especially their involvement in WBPH stress. OsWRKY and OsNAC TFs are involved in various developmental processes and responses to biotic and abiotic stresses. However, no comprehensive reports are available on the specific phycological functions of most of the OsWRKY and OsNAC genes in rice during WBPH infestation. The current study aimed to comprehensively explore the OsWRKY and OsNAC genes by analyzing their phylogenetic relationships, subcellular localizations, exon-intron arrangements, conserved motif identities, chromosomal allocations, interaction networks and differential gene expressions during stress conditions. Comparative phylogenetic trees of 101 OsWRKY with 72 AtWRKY genes, and 121 OsNAC with 110 AtNAC genes were constructed to study relationships among these TFs across species. Phylogenetic relationships classified OsWRKY and OsNAC into eight and nine clades, respectively. Most TFs in the same clade had similar genomic features that represented similar functions, and had a high degree of co-expression. Some OsWRKYs (Os09g0417800 (OsWRKY62), Os11g0117600 (OsWRKY50), Os11g0117400 (OsWRKY104) and OsNACs (Os05g0442700, Os12g0630800, Os01g0862800 and Os12g0156100)) showed significantly higher expressions under WBPH infestation, based on transcriptome datasets. This study provides valuable information and clues about predicting the potential roles of OsWRKYs and OsNACs in rice, by combining their genome-wide characterization, expression profiling, protein-protein interactions and gene expressions under WBPH stress. These findings may require additional investigation to understand their metabolic and expression processes, and to develop rice cultivars that are resistant to WBPH.

Genome-Wide Identification and Characterization of the Oat (Avena sativa L.) WRKY Transcription Factor Family

The WRKY family is widely involved in the regulation of plant growth and stress response and is one of the largest gene families related to plant environmental adaptation. However, no systematic studies on the WRKY family in oat (Avena sativa L.) have been conducted to date. The recently published complete genome sequence of oat enables the systematic analysis of the AsWRKYs. Based on a genome-wide study of oat, we identified 162 AsWRKYs that were unevenly distributed across 21 chromosomes; a phylogenetic tree of WRKY domains divided these genes into three groups (I, II, and III). We also analyzed the gene duplication events and identified a total of 111 gene pairs that showed strong purifying selection during the evolutionary process. Surprisingly, almost all genes evolved after the completion of subgenomic differentiation of hexaploid oat. Further studies on the functional analysis indicated that AsWRKYs were widely involved in various biological processes. Notably, expression patterns of 16 AsWRKY genes revealed that the response of AsWRKYs were affected by stress level and time. In conclusion, this study provides a reference for further analysis of the role of WRKY transcription factors in species evolution and functional differentiation.

Overexpression of Myrothamnus flabellifolia MfWRKY41 confers drought and salinity tolerance by enhancing root system and antioxidation ability in Arabidopsis

Myrothamnus flabellifolia is the only woody resurrection plant discovered so far and could recover from extreme desiccation condition. However, few genes related to its strong drought tolerance have been characterized, and the underlying molecular mechanisms remains mysterious. Members of WRKY transcription factor family are effective in regulating abiotic stress responses or tolerance in various plants. An early dehydration-induced gene encoding a WRKY transcription factor namely MfWRKY41 was isolated from M. flabellifolia, which is homologous to AtWRKY41 of Arabidopsis. It contains a typical WRKY domain and zinc finger motif, and is located in the nucleus. Comparing to wild type, the four transgenic lines overexpressing MfWRKY41 showed better growth performance under drought and salt treatments, and exhibited higher chlorophyll content, lower water loss rate and stomatal aperture and better osmotic adjustment capacity. These results indicated that MfWRKY41 of M. flabellifolia positively regulates drought as well as salinity responses. Interestingly, the root system architecture, including lateral root number and primary root length, of the transgenic lines was enhanced by MfWRKY41 under both normal and stressful conditions, and the antioxidation ability was also significantly improved. Therefore, MfWRKY41 may have potential application values in genetic improvement of plant tolerance to drought and salinity stresses. The molecular mechanism involving in the regulatory roles of MfWRKY41 is worthy being explored in the future.

Assessment of Rice Sheath Blight Resistance Including Associations with Plant Architecture, as Revealed by Genome-Wide Association Studies

BACKGROUND

Sheath blight (ShB) disease caused by Rhizoctonia solani Kühn, is one of the most economically damaging rice (Oryza sativa L.) diseases worldwide. There are no known major resistance genes, leaving only partial resistance from small-effect QTL to deploy for cultivar improvement. Many ShB-QTL are associated with plant architectural traits detrimental to yield, including tall plants, late maturity, or open canopy from few or procumbent tillers, which confound detection of physiological resistance.

RESULTS

To identify QTL for ShB resistance, 417 accessions from the Rice Diversity Panel 1 (RDP1), developed for association mapping studies, were evaluated for ShB resistance, plant height and days to heading in inoculated field plots in Arkansas, USA (AR) and Nanning, China (NC). Inoculated greenhouse-grown plants were used to evaluate ShB using a seedling-stage method to eliminate effects from height or maturity, and tiller (TN) and panicle number (PN) per plant. Potted plants were used to evaluate the RDP1 for TN and PN. Genome-wide association (GWA) mapping with over 3.4 million SNPs identified 21 targeted SNP markers associated with ShB which tagged 18 ShB-QTL not associated with undesirable plant architecture traits. Ten SNPs were associated with ShB among accessions of the Indica subspecies, ten among Japonica subspecies accessions, and one among all RDP1 accessions. Across the 18 ShB QTL, only qShB4-1 was not previously reported in biparental mapping studies and qShB9 was not reported in the GWA ShB studies. All 14 PN QTL overlapped with TN QTL, with 15 total TN QTL identified. Allele effects at the five TN QTL co-located with ShB QTL indicated that increased TN does not inevitably increase disease development; in fact, for four ShB QTL that overlapped TN QTL, the alleles increasing resistance were associated with increased TN and PN, suggesting a desirable coupling of alleles at linked genes.

CONCLUSIONS

Nineteen accessions identified as containing the most SNP alleles associated with ShB resistance for each subpopulation were resistant in both AR and NC field trials. Rice breeders can utilize these accessions and SNPs to develop cultivars with enhanced ShB resistance along with increased TN and PN for improved yield potential.

Genome-Wide Identification of the SAMS Gene Family in Upland Cotton (Gossypium hirsutum L.) and Expression Analysis in Drought Stress Treatments

Cotton is an important commercial crop whose growth and yield are severely affected by drought. S-adenosylmethionine (SAM) is widely involved in the plant stress response and growth regulation; however, the role of the S-adenosylmethionine synthase (SAMS) gene family in this process is poorly understood. Here, we systematically analyzed the expression of SAMS genes in Upland Cotton (Gossypium hirsutum L.). A total of 16 SAMS genes were identified, each with a similar predicted structure. A large number of cis-acting elements involved in the response to abiotic stress were predicted based on promoter analysis, indicating a likely important role in abiotic stress responses. The results of qRT-PCR validation showed that GhSAMS genes had different expression patterns after drought stress and in response to drought stress. Analysis of a selected subset of GhSAMS genes showed increased expression in cultivar Xinluzhong 39 (drought resistant) when compared to cultivar Xinluzao 26 (drought-sensitive) upland cotton. This study provides important relevant information for further study of SAMS genes in drought resistance research of upland cotton, which is helpful for drought-resistance improvement of upland cotton.

The transcription factor WRKY75 regulates the development of adventitious roots, lateral buds and callus by modulating hydrogen peroxide content in poplar

Hydrogen peroxide (H2O2) plays important roles in plant development. Adventitious roots (AR), lateral buds (LB) and callus formation are important traits for plants. Here, a gene encoding RESPIRATORY BURST OXIDASE HOMOLOG B (PdeRBOHB) from poplar line 'NL895' (Populus. deltoides × P. euramericana) was predicted to be involved in H2O2 accumulation, and lines with reduced expression were generated. H2O2 content was decreased, and the development of adventitious roots, lateral buds, and callus was inhibited in reduced expression PdeRBOHB lines. A gene encoding PdeWRKY75 was identified as the upstream transcription factor positively regulating PdeRBOHB. This regulation was confirmed by dual luciferase reporter assay, GUS transient expression analysis and electrophoretic mobility shift assay. In the reduced expression PdeWRKY75 lines, H2O2 content was decreased and the development of adventitious roots, lateral buds, and callus development was inhibited, while in the overexpression lines, H2O2 content was increased and the development of adventitious roots and lateral buds was inhibited, but callus formation was enhanced. Additionally, reduced expression PdeRBOHB lines showed lowered expression of PdeWRKY75, while exogenous application of H2O2 showed the opposite effect. Together, these results suggest that PdeWRKY75 and PdeRBOHB are part of a regulatory module in H2O2 accumulation, which is involved in the regulation of multiple biological processes.

Genome-Wide Expression Analysis of Root Tips in Contrasting Rice Genotypes Revealed Novel Candidate Genes for Water Stress Adaptation

Root system architecture (RSA) is an important agronomic trait with vital roles in plant productivity under water stress conditions. A deep and branched root system may help plants to avoid water stress by enabling them to acquire more water and nutrient resources. Nevertheless, our knowledge of the genetics and molecular control mechanisms of RSA is still relatively limited. In this study, we analyzed the transcriptome response of root tips to water stress in two well-known genotypes of rice: IR64, a high-yielding lowland genotype, which represents a drought-susceptible and shallow-rooting genotype; and Azucena, a traditional, upland, drought-tolerant and deep-rooting genotype. We collected samples from three zones (Z) of root tip: two consecutive 5 mm sections (Z1 and Z2) and the following next 10 mm section (Z3), which mainly includes meristematic and maturation regions. Our results showed that Z1 of Azucena was enriched for genes involved in cell cycle and division and root growth and development whereas in IR64 root, responses to oxidative stress were strongly enriched. While the expansion of the lateral root system was used as a strategy by both genotypes when facing water shortage, it was more pronounced in Azucena. Our results also suggested that by enhancing meristematic cell wall thickening for insulation purposes as a means of confronting stress, the sensitive IR64 genotype may have reduced its capacity for root elongation to extract water from deeper layers of the soil. Furthermore, several members of gene families such as NAC, AP2/ERF, AUX/IAA, EXPANSIN, WRKY, and MYB emerged as main players in RSA and drought adaptation. We also found that HSP and HSF gene families participated in oxidative stress inhibition in IR64 root tip. Meta-quantitative trait loci (QTL) analysis revealed that 288 differentially expressed genes were colocalized with RSA QTLs previously reported under drought and normal conditions. This finding warrants further research into their possible roles in drought adaptation. Overall, our analyses presented several major molecular differences between Azucena and IR64, which may partly explain their differential root growth responses to water stress. It appears that Azucena avoided water stress through enhancing growth and root exploration to access water, whereas IR64 might mainly rely on cell insulation to maintain water and antioxidant system to withstand stress. We identified a large number of novel RSA and drought associated candidate genes, which should encourage further exploration of their potential to enhance drought adaptation in rice.

Network organization of the plant immune system: from pathogen perception to robust defense induction

The plant immune system has been explored essentially through the study of qualitative resistance, a simple form of immunity, and from a reductionist point of view. The recent identification of genes conferring quantitative disease resistance revealed a large array of functions, suggesting more complex mechanisms. In addition, thanks to the advent of high-throughput analyses and system approaches, our view of the immune system has become more integrative, revealing that plant immunity should rather be seen as a distributed and highly connected molecular network including diverse functions to optimize expression of plant defenses to pathogens. Here, we review the recent progress made to understand the network complexity of regulatory pathways leading to plant immunity, from pathogen perception, through signaling pathways and finally to immune responses. We also analyze the topological organization of these networks and their emergent properties, crucial to predict novel immune functions and test them experimentally. Finally, we report how these networks might be regulated by environmental clues. Although system approaches remain extremely scarce in this area of research, a growing body of evidence indicates that the plant response to combined biotic and abiotic stresses cannot be inferred from responses to individual stresses. A view of possible research avenues in this nascent biology domain is finally proposed.

WRKY Transcription Factors in Cassava Contribute to Regulation of Tolerance and Susceptibility to Cassava Mosaic Disease through Stress Responses

Among the numerous biological constraints that hinder cassava (Manihot esculenta Crantz) production, foremost is cassava mosaic disease (CMD) caused by virus members of the family Geminiviridae, genus Begomovirus. The mechanisms of CMD tolerance and susceptibility are not fully understood; however, CMD susceptible T200 and tolerant TME3 cassava landraces have been shown to exhibit different large-scale transcriptional reprogramming in response to South African cassava mosaic virus (SACMV). Recent identification of 85 MeWRKY transcription factors in cassava demonstrated high orthology with those in Arabidopsis, however, little is known about their roles in virus responses in this non-model crop. Significant differences in MeWRKY expression and regulatory networks between the T200 and TME3 landraces were demonstrated. Overall, WRKY expression and associated hormone and enriched biological processes in both landraces reflect oxidative and other biotic stress responses to SACMV. Notably, MeWRKY11 and MeWRKY81 were uniquely up and downregulated at 12 and 67 days post infection (dpi) respectively in TME3, implicating a role in tolerance and symptom recovery. AtWRKY28 and AtWRKY40 homologs of MeWRKY81 and MeWRKY11, respectively, have been shown to be involved in regulation of jasmonic and salicylic acid signaling in Arabidopsis. AtWRKY28 is an interactor in the RPW8-NBS resistance (R) protein network and downregulation of its homolog MeWRKY81 at 67 dpi in TME3 suggests a negative role for this WRKY in SACMV tolerance. In contrast, in T200, nine MeWRKYs were differentially expressed from early (12 dpi), middle (32 dpi) to late (67 dpi) infection. MeWRKY27 (homolog AtWRKY33) and MeWRKY55 (homolog AtWRKY53) were uniquely up-regulated at 12, 32 and 67 dpi in T200. AtWRKY33 and AtWRKY53 are positive regulators of leaf senescence and oxidative stress in Arabidopsis, suggesting MeWRKY55 and 27 contribute to susceptibility in T200.

When rice gets the chills: comparative transcriptome profiling at germination shows WRKY transcription factor responses

Rice is vital for food security. Due to its tropical origin, rice suffers from cold temperatures that affect its entire life cycle. Key genes have been identified involved in cold tolerance. WRKYs are generally downstream of the MAPK cascade and can act together with VQ proteins to regulate stress-responsive genes. Chilling treatment was applied at germination to two rice genotypes (tolerant and sensitive). Shoots at S3 stage were collected for RNA-seq to identify OsWRKY, OsMAPKs and OsVQs expression. Relationships among MAPKs, WRKYs and VQs were predicted through correlation analysis. OsWRKY transcriptional regulation was predicted by in silico analysis of cis-regulatory elements. A total of 39 OsWRKYs were differentially expressed. OsWRKY21, OsWRK24 and OsWRKY69 are potential positive regulators, while OsWRKY10, OsWRK47, OsWRKY62, OsWRKY72 and OsWRKY77 are potential negative regulators, of chilling tolerance. 12 OsMAPKs were differentially expressed. OsMAPKs were downregulated and negatively correlated with the upregulated OsWRKYs in the tolerant genotype. 19 OsVQs were differentially expressed, three and six OsVQs were positively correlated with OsWRKYs in the tolerant and sensitive genotypes, respectively. Seven differentially expressed OsWRKYs have cold-responsive elements in their promoters and five upregulated OsWRKYs in the tolerant genotype contained the W-box motif. Chilling causes changes in OsWRKY, OsMAPK and OsVQ gene expression at germination. OsWRKYs may not act downstream of the MAPK cascade to coordinate chilling tolerance, but OsWRKYs may act with VQs to regulate chilling tolerance. Candidate OsWRKYs are correlated and have a W-box in the promoter, suggesting an auto-regulation mechanism.

Comparative transcriptome and metabolome profiling reveal molecular mechanisms underlying OsDRAP1-mediated salt tolerance in rice

Integration of transcriptomics and metabolomics data can provide detailed information for better understanding the molecular mechanisms underlying salt tolerance in rice. In the present study, we report a comprehensive analysis of the transcriptome and metabolome of rice overexpressing the OsDRAP1 gene, which encodes an ERF transcription factor and was previously identified to be conferring drought tolerance. Phenotypic analysis showed that OsDRAP1 overexpression (OE) improved salt tolerance by increasing the survival rate under salt stress. OsDRAP1 affected the physiological indices such as superoxide dismutase (SOD), catalase (CAT) and malondialdehyde (MDA) to enhance redox homeostasis and membrane stability in response to salt stress. Higher basal expression of OsDRAP1 resulted in differential expression of genes that potentially function in intrinsic salt tolerance. A core set of genes with distinct functions in transcriptional regulation, organelle gene expression and ion transport were substantially up-regulated in the OE line in response to salt stress, implying their important role in OsDRAP1-mediated salt tolerance. Correspondingly, metabolome profiling detected a number of differentially metabolites in the OE line relative to the wild type under salt stress. These metabolites, including amino acids (proline, valine), organic acids (glyceric acid, phosphoenolpyruvic acid and ascorbic acid) and many secondary metabolites, accumulated to higher levels in the OE line, demonstrating their role in salt tolerance. Integration of transcriptome and metabolome analysis highlights the crucial role of amino acids and carbohydrate metabolism pathways in OsDRAP1-mediated salt tolerance.

Grain Inorganic Arsenic Content in Rice Managed Through Targeted Introgressions and Irrigation Management

Arsenic (As) accumulation in rice grain is a significant public health concern. Inorganic As (iAs) is of particular concern because it has increased toxicity as compared to organic As. Irrigation management practices, such as alternate wetting and drying (AWD), as well as genotypic differences between cultivars, have been shown to influence As accumulation in rice grain. A 2 year field study using a Lemont × TeQing backcross introgression line (TIL) mapping population examined the impact of genotype and AWD severity on iAs grain concentrations. The "Safe"-AWD [35-40% soil volumetric water content (VWC)] treatment did not reduce grain iAs levels, whereas the more severe AWD30 (25-30% VWC) consistently reduced iAs concentrations across all genotypes. The TILs displayed a range of iAs concentrations by genotype, from less than 10 to up to 46 μg kg-1 under AWD30 and from 28 to 104 μg kg-1 under Safe-AWD. TIL grain iAs concentrations for flood treatments across both years ranged from 26 to 127 μg kg-1. Additionally, seven quantitative trait loci (QTLs) were identified in the mapping population associated with grain iAs. A subset of eight TILs and their parents were grown to confirm field-identified grain iAs QTLs in a controlled greenhouse environment. Greenhouse results confirmed the genotypic grain iAs patterns observed in the field; however, iAs concentrations were higher under greenhouse conditions as compared to the field. In the greenhouse, the number of days under AWD was negatively correlated with grain iAs concentrations. Thus, longer drying periods to meet the same soil VWC resulted in lower grain iAs levels. Both the number and combinations of iAs-affecting QTLs significantly impacted grain iAs concentrations. Therefore, identifying more grain iAs-affecting QTLs could be important to inform future breeding efforts for low iAs rice varieties. Our study suggests that coupling AWD practices targeting a soil VWC of less than or equal to 30% coupled with the use of cultivars developed to possess multiple QTLs that negatively regulate grain iAs concentrations will be helpful in mitigating exposure of iAs from rice consumption.

Function and Mechanism of WRKY Transcription Factors in Abiotic Stress Responses of Plants

The WRKY gene family is a plant-specific transcription factor (TF) group, playing important roles in many different response pathways of diverse abiotic stresses (drought, saline, alkali, temperature, and ultraviolet radiation, and so forth). In recent years, many studies have explored the role and mechanism of WRKY family members from model plants to agricultural crops and other species. Abiotic stress adversely affects the growth and development of plants. Thus, a review of WRKY with stress responses is important to increase our understanding of abiotic stress responses in plants. Here, we summarize the structural characteristics and regulatory mechanism of WRKY transcription factors and their responses to abiotic stress. We also discuss current issues and future perspectives of WRKY transcription factor research.

The overexpression of OsACBP5 protects transgenic rice against necrotrophic, hemibiotrophic and biotrophic pathogens

The most devastating diseases in rice (Oryza sativa) are sheath blight caused by the fungal necrotroph Rhizoctonia solani, rice blast by hemibiotrophic fungus Magnaporthe oryzae, and leaf blight by bacterial biotroph Xanthomonas oryzae (Xoo). It has been reported that the Class III acyl-CoA-binding proteins (ACBPs) such as those from dicots (Arabidopsis and grapevine) play a role in defence against biotrophic pathogens. Of the six Arabidopsis (Arabidopsis thaliana) ACBPs, AtACBP3 conferred protection in transgenic Arabidopsis against Pseudomonas syringae, but not the necrotrophic fungus, Botrytis cinerea. Similar to Arabidopsis, rice possesses six ACBPs, designated OsACBPs. The aims of this study were to test whether OsACBP5, the homologue of AtACBP3, can confer resistance against representative necrotrophic, hemibiotrophic and biotrophic phytopathogens and to understand the mechanisms in protection. Herein, when OsACBP5 was overexpressed in rice, the OsACBP5-overexpressing (OsACBP5-OE) lines exhibited enhanced disease resistance against representative necrotrophic (R. solani & Cercospora oryzae), hemibiotrophic (M. oryzae & Fusarium graminearum) and biotrophic (Xoo) phytopathogens. Progeny from a cross between OsACBP5-OE9 and the jasmonate (JA)-signalling deficient mutant were more susceptible than the wild type to infection by the necrotroph R. solani. In contrast, progeny from a cross between OsACBP5-OE9 and the salicylic acid (SA)-signalling deficient mutant was more susceptible to infection by the hemibiotroph M. oryzae and biotroph Xoo. Hence, enhanced resistance of OsACBP5-OEs against representative necrotrophs appears to be JA-dependent whilst that to (hemi)biotrophs is SA-mediated.

The transcriptional reprograming and functional identification of WRKY family members in pepper's response to Phytophthora capsici infection

BACKGROUND

Plant transcription factors (TFs) are key transcriptional regulators to manipulate the regulatory network of host immunity. However, the globally transcriptional reprogramming of plant TF families in response to pathogens, especially between the resistant and susceptible host plants, remains largely unknown.

RESULTS

Here, we performed time-series RNA-seq from a resistant pepper line CM334 and a susceptible pepper line EC01 upon challenged with Phytophthora capsici, and enrichment analysis indicated that WRKY family most significantly enriched in both CM334 and EC01. Interestingly, we found that nearly half of the WRKY family members were significantly up-regulated, whereas none of them were down-regulated in the two lines. These induced WRKY genes were greatly overlapped between CM334 and EC01. More strikingly, most of these induced WRKY genes were expressed in time-order patterns, and could be mainly divided into three subgroups: early response (3 h-up), mid response (24 h-up) and mid-late response (ML-up) genes. Moreover, it was found that the responses of these ML-up genes were several hours delayed in EC01. Furthermore, a total of 19 induced WRKY genes were selected for functional identification by virus-induced gene silencing. The result revealed that silencing of CaWRKY03-6, CaWRKY03-7, CaWRKY06-5 or CaWRKY10-4 significantly increase the susceptibility to P. capsici both in CM334 and EC01, indicating that they might contribute to pepper's basal defense against P. capsici; while silencing of CaWRKY08-4 and CaWRKY01-10 significantly impaired the disease resistance in CM334 but not in EC01, suggesting that these two WRKY genes are prominent modulators specifically in the resistant pepper plants.

CONCLUSIONS

These results considerably extend our understanding of WRKY gene family in pepper's resistance against P. capsici and provide potential applications for genetic improvement against phytophthora blight.

WRKY33-PIF4 loop is required for the regulation of H2O2 homeostasis

The reactive oxygen species (ROS) are continuously produced and are essential for mediating the growth and development of plants. However too much accumulation of ROS can result in the oxidative damage to cells, especially under the adverse environmental conditions. Plants have evolved sophisticated strategies to regulate the homeostasis of H₂O₂. In this study, we generated transgenic Arabidopsis plants in the Ws ecotype (Ws) background in which WRKY33 is co-suppressed (csWRKY33/Ws). Compared with Ws, csWRKY33/Ws plants accumulate more H₂O₂. RNA-seq analysis indicated that in csWRKY33/Ws plants, expression of oxidative stress related genes such as ascorbate peroxidase 2 (APX2) is affected. Over-expression of APX2 can rescue the phenotype of csWRKY33/Ws, suggesting that the changes in the growth of csWRKY33/Ws is duo to the higher accumulation of H₂O₂. Analysis of the CHIP-seq data suggested that WRKY33 can directly regulate the expression of PIF4, vice versa. qPCR analysis also confirmed that the mutual regulation between WRKY33 and PIF4. Similar to that of csWRKY33/Ws, and the accumulation of H₂O₂ in pif4 also increased. Taken together, our results reveal a WRKY33-PIF4 regulatory loop that appears to play an important role in regulating the growth and development of seedlings by mediating H₂O₂ homeostasis.

Comparative Transcriptome-Based Mining and Expression Profiling of Transcription Factors Related to Cold Tolerance in Peanut

Plants tolerate cold stress by regulating gene networks controlling cellular and physiological traits to modify growth and development. Transcription factor (TF)-directed regulation of transcription within these gene networks is key to eliciting appropriate responses. Identifying TFs related to cold tolerance contributes to cold-tolerant crop breeding. In this study, a comparative transcriptome analysis was carried out to investigate global gene expression of entire TFs in two peanut varieties with different cold-tolerant abilities. A total of 87 TF families including 2328 TF genes were identified. Among them, 445 TF genes were significantly differentially expressed in two peanut varieties under cold stress. The TF families represented by the largest numbers of differentially expressed members were bHLH (basic helix—loop—helix protein), C2H2 (Cys2/His2 zinc finger protein), ERF (ethylene-responsive factor), MYB (v-myb avian myeloblastosis viral oncogene homolog), NAC (NAM, ATAF1/2, CUC2) and WRKY TFs. Phylogenetic evolutionary analysis, temporal expression profiling, protein–protein interaction (PPI) network, and functional enrichment of differentially expressed TFs revealed the importance of plant hormone signal transduction and plant-pathogen interaction pathways and their possible mechanism in peanut cold tolerance. This study contributes to a better understanding of the complex mechanism of TFs in response to cold stress in peanut and provides valuable resources for the investigation of evolutionary history and biological functions of peanut TFs genes involved in cold tolerance.

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.

The Chalcone Isomerase Family in Cotton: Whole-Genome Bioinformatic and Expression Analyses of the Gossypium barbadense L. Response to Fusarium Wilt Infection

Chalcone isomerase (CHI) is a key component of phenylalanine metabolism that can produce a variety of flavonoids. However, little information and no systematic analysis of CHI genes is available for cotton. Here, we identified 33 CHI genes in the complete genome sequences of four cotton species (Gossypium arboretum L., Gossypium raimondii L., Gossypium hirsutum L., and Gossypium barbadense L.). Cotton CHI proteins were classified into two main groups, and whole-genome/segmental and dispersed duplication events were important in CHI gene family expansion. qRT-PCR and semiquantitative RT-PCR results suggest that CHI genes exhibit temporal and spatial variation and respond to infection with Fusarium wilt race 7. A preliminary model of CHI gene involvement in cotton evolution was established. Pairwise comparison revealed that seven CHI genes showed higher expression in cultivar 06-146 than in cultivar Xinhai 14. Overall, this whole-genome identification unlocks a new approach to the comprehensive functional analysis of the CHI gene family, which may be involved in adaptation to plant pathogen stress.

IbMPK3/IbMPK6-mediated IbSPF1 phosphorylation promotes tolerance to bacterial pathogen in sweetpotato

KEY MESSAGE

IbSPF1, a novel target of IbMPK3/IbMPK6, regulates biotic stress response in sweetpotato. Environmental stresses due to biotic and abiotic factors negatively affect crop quality and productivity. To minimize the damage caused by these factors, numerous stress signaling pathways are activated in plants. Among these, the mitogen-activated protein kinase (MAPK) signaling cascade plays a pivotal role in diverse plant stress responses. MPK3 and MPK6 function in several cellular signaling pathways by phosphorylating downstream partner proteins in response to environmental stresses. However, little is known about the MPK3/MPK6 signaling pathway in sweetpotato [Ipomoea batatas (L.) Lam]. We recently confirmed that IbMPK3 and IbMPK6, two pathogen-responsive MAPKs, play essential roles in defense gene activation in sweetpotato. In this study, we show that sweetpotato SP8-binding factor (IbSPF1), a substrate of IbMPK3/IbMPK6, functions as a transcriptional regulator of biotic stress signaling in sweetpotato. IbSPF1 specifically interacts with IbMPK3 and IbMPK6, which phosphorylate Ser75 and Ser110 residues of IbSPF1. This increases the affinity of IbSPF1 for the W-box element in target gene promoters. Additionally, the expression of IbSPF1 was up-regulated under various stress conditions and different hormone treatments involved in plant defense responses. Interestingly, the phospho-mimicking mutant of IbSPF1 showed enhanced resistance to Pseudomonas syringae pv. tabaci, and transient expression of mutant IbSPF1 induced the expression of pathogenesis-related genes. These results indicate that the phosphorylation of IbSPF1 by IbMPK3/IbMPK6 plays a critical role in plant immunity by up-regulating the expression of downstream genes.

Gene networks underlying the early regulation of Paraburkholderia phytofirmans PsJN induced systemic resistance in Arabidopsis

Plant defense responses to biotic stresses are complex biological processes, all governed by sophisticated molecular regulations. Induced systemic resistance (ISR) is one of these defense mechanisms where beneficial bacteria or fungi prime plants to resist pathogens or pest attacks. In ISR, the defense arsenal in plants remains dormant and it is only triggered by an infection, allowing a better allocation of plant resources. Our group recently described that the well-known beneficial bacterium Paraburkholderia phytofirmans PsJN is able to induce Arabidopsis thaliana resistance to Pseudomonas syringae pv. tomato (Pst) DC3000 through ISR, and that ethylene, jasmonate and salicylic acid are involved in this protection. Nevertheless, the molecular networks governing this beneficial interaction remain unknown. To tackle this issue, we analyzed the temporal changes in the transcriptome of PsJN-inoculated plants before and after being infected with Pst DC3000. These data were used to perform a gene network analysis to identify highly connected transcription factors. Before the pathogen challenge, the strain PsJN regulated 405 genes (corresponding to 1.8% of the analyzed genome). PsJN-inoculated plants presented a faster and stronger transcriptional response at 1-hour post infection (hpi) compared with the non-inoculated plants, which presented the highest transcriptional changes at 24 hpi. A principal component analysis showed that PsJN-induced plant responses to the pathogen could be differentiated from those induced by the pathogen itself. Forty-eight transcription factors were regulated by PsJN at 1 hpi, and a system biology analysis revealed a network with four clusters. Within these clusters LHY, WRKY28, MYB31 and RRTF1 are highly connected transcription factors, which could act as hub regulators in this interaction. Concordantly with our previous results, these clusters are related to jasmonate, ethylene, salicylic, acid and ROS pathways. These results indicate that a rapid and specific response of PsJN-inoculated plants to the virulent DC3000 strain could be the pivotal element in the protection mechanism.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

Genes, pathways and transcription factors involved in seedling stage chilling stress tolerance in indica rice through RNA-Seq analysis

BACKGROUND

Rice plants show yellowing, stunting, withering, reduced tillering and utimately low productivity in susceptible varieties under low temperature stress. Comparative transcriptome analysis was performed to identify novel transcripts, gain new insights into different gene expression and pathways involved in cold tolerance in rice.

RESULTS

Comparative transcriptome analyses of 5 treatments based on chilling stress exposure revealed more down regulated genes in susceptible and higher up regulated genes in tolerant genotypes. A total of 13930 and 10599 differentially expressed genes (DEGs) were detected in cold susceptible variety (CSV) and cold tolerant variety (CTV), respectively. A continuous increase in DEGs at 6, 12, 24 and 48 h exposure of cold stress was detected in both the genotypes. Gene ontology (GO) analysis revealed 18 CSV and 28 CTV term significantly involved in molecular function, cellular component and biological process. GO classification showed a significant role of transcription regulation, oxygen, lipid binding, catalytic and hydrolase activity for tolerance response. Absence of photosynthesis related genes, storage products like starch and synthesis of other classes of molecules like fatty acids and terpenes during the stress were noticed in susceptible genotype. However, biological regulations, generation of precursor metabolites, signal transduction, photosynthesis, regulation of cellular process, energy and carbohydrate metabolism were seen in tolerant genotype during the stress. KEGG pathway annotation revealed more number of genes regulating different pathways resulting in more tolerant. During early response phase, 24 and 11 DEGs were enriched in CTV and CSV, respectively in energy metabolism pathways. Among the 1583 DEG transcription factors (TF) genes, 69 WRKY, 46 bZIP, 41 NAC, 40 ERF, 31/14 MYB/MYB-related, 22 bHLH, 17 Nin-like 7 HSF and 4C3H were involved during early response phase. Late response phase showed 30 bHLH, 65 NAC, 30 ERF, 26/20 MYB/MYB-related, 11 C3H, 12 HSF, 86 Nin-like, 41 AP2/ERF, 55 bZIP and 98 WRKY members TF genes. The recovery phase included 18 bHLH, 50 NAC, 31 ERF, 24/13 MYB/MYB-related, 4 C3H, 4 HSF, 14 Nin-like, 31 bZIP and 114 WRKY TF genes.

CONCLUSIONS

Transcriptome analysis of contrasting genotypes for cold tolerance detected the genes, pathways and transcription factors involved in the stress tolerance.

QTL-BSA: A Bulked Segregant Analysis and Visualization Pipeline for QTL-seq

In recent years, the application of Whole Genome Sequencing (WGS) on plants has generated sufficient data for the identification of trait-associated genomic loci or genes. A high-throughput genome-assisted QTL-seq strategy, combined with bulked-segregant analysis and WGS of two bulked populations from a segregating progeny with opposite phenotypic trait values, has gained increasing popularities in research community. However, there is no publicly available user friendly software for the identification and visualization. Hence, we developed a tool named QTL-BSA (QTL-bulked segregant analysis and visualization pipeline), which could facilitate the rapid identification and visualization of candidate QTLs from QTL-seq. As a proof-of-concept study, we have applied the tool for the rapid discovery and the identification of genes related with the partial blast resistance in rice. Genomic region of the major QTL identified on chromosome 6, is located between 1.52 and 4.32 Mb, which is consistent with previous studies (2.39-4.39 Mb). We also derived the gene and QTLs functional annotation of this region. QTL-BSA offers a comprehensive solution to facilitate a wide range of programming and visualization tasks in QTL-seq analysis, is expected to be used widely by the research community.

The Involvement of Jasmonic Acid, Ethylene, and Salicylic Acid in the Signaling Pathway of Clonostachys rosea-Induced Resistance to Gray Mold Disease in Tomato

Tomato gray mold disease caused by Botrytis cinerea is a serious disease that threatens tomato production around the world. Clonostachys rosea has been used successfully as a biocontrol agent against divergent plant pathogens, including B. cinerea. To understand the signal transduction pathway of C. rosea-induced resistance to tomato gray mold disease, the effects of C. rosea on gray mold tomato leaves along with changes in the activities of three defense enzymes (phenylalanine ammonialyase [PAL], polyphenol oxidase [PPO], and catalase [CAT]), second messengers (nitric oxide [NO], hydrogen peroxide [H₂O₂], and superoxide anion radical [O₂^-]), and stress-related genes (mitogen-activated protein kinase [MAPK], WRKY, Lexyl2, and atpA) in four different hormone-deficient (jasmonic acid [JA], ethylene [ET], salicylic acid [SA], and gibberellin) tomato mutants were investigated. The results revealed that C. rosea significantly inhibited the growth of mycelia and spore germination of B. cinerea. Furthermore, it reduced the incidence of gray mold disease, induced higher levels of PAL and PPO, and induced lower levels of CAT activities in tomato leaves. Moreover, it also increased NO, H₂O₂, and O₂^- levels and the gene expression levels of WRKY, MAPK, atpA, and Lexyl2. The incidence of gray mold disease in four hormone-deficient mutants was higher than that in the corresponding wild-type tomato plants. Among all of these hormone-deficient tomato mutants, JA had the most significant effect in regulating the different signal molecules. Additional study suggested that JA upregulated the expression of Lexyl2, MAPK, and WRKY but downregulated atpA. Furthermore, JA also enhanced the activity of PAL, PPO, and CAT and the production of NO and H₂O₂. SA downregulated CAT and PAL, whereas ET upregulated PAL but downregulated CAT. This study is of significance in understanding the regulatory pathways and biocontrol mechanism of C. rosea against B. cinerea.

CaMPK9 increases the stability of CaWRKY40 transcription factor which triggers defense response in chickpea upon Fusarium oxysporum f. sp. ciceri Race1 infection

Physical interaction and phosphorylation by CaMPK9 protects the degradation of CaWRKY40 that induces resistance response in chickpea to Fusarium wilt disease by modulating the transcription of defense responsive genes. WRKY transcription factors (TFs) are the global regulators of plant defense signaling that modulate immune responses in host plants by regulating transcription of downstream target genes upon challenged by pathogens. However, very little is known about immune responsive role of Cicer arietinum L. (Ca) WRKY TFs particularly. Using two contrasting chickpea genotypes with respect to resistance against Fusarium oxysporum f. sp. ciceri Race1 (Foc1), we demonstrate transcript accumulation of different CaWRKYs under multiple stresses and establish that CaWRKY40 triggers defense. CaWRKY40 overexpressing chickpea mounts resistance to Foc1 by positively modulating the defense related gene expression. EMSA, ChIP assay and real-time PCR analyses suggest CaWRKY40 binds at the promoters and positively regulates transcription of CaDefensin and CaWRKY33. Further studies revealed that mitogen Activated Protein Kinase9 (CaMPK9) phosphorylates CaWRKY40 by directly interacting with its two canonical serine residues. Interestingly, CaMPK9 is unable to interact with CaWRKY40 when the relevant two serine residues were replaced by alanine. Overexpression of serine mutated WRKY40 isoform in chickpea fails to provide resistance against Foc1. Mutated WRKY40^{Ser.224/225 to AA} overexpressing chickpea resumes its ability to confer resistance against Foc1 after application of 26S proteasomal inhibitor MG132, suggests that phosphorylation is essential to protect CaWRKY40 from proteasomal degradation. CaMPK9 silencing also led to susceptibility in chickpea to Foc1. Altogether, our results elucidate positive regulatory roles of CaMPK9 and CaWRKY40 in modulating defense response in chickpea upon Foc1 infection.

Variations in the Structure and Evolution of Rice WRKY Genes in Indica and Japonica Genotypes and their Co-expression Network in Mediating Disease Resistance

WRKY transcription factor (TF) family regulates many functions in plant growth and development and also during biotic and abiotic stress. In this study, 101 WRKY TF gene models in indica and japonica rice were used to conduct evolutionary analysis, gene structure analysis, and motif composition. Co-expression analysis was carried out first by selecting the differentially expressing genes that showed a significant change in response to the pathogens from Rice Oligonucleotide Array Database (ROAD). About 82 genes showed responses to infection by Magnaporthe oryzae or Xanthomonas oryzae pv. oryzae. Co-expression gene network was constructed using direct neighborhood and context associated inbuilt mode in RiceNetv2 tool. Only 41 genes showed interaction with 2299 non-WRKY genes. Variations exist in the structure and evolution of WRKY genes among indica and japonica genotypes which have important implications in their differential roles including disease resistance. WRKY genes mediate a complex networking and co-express along with other WRKY and non-WRKY genes to mediate resistance against fungal and bacterial pathogens in rice.

Comparing time series transcriptome data between plants using a network module finding algorithm

BACKGROUND

Comparative transcriptome analysis is the comparison of expression patterns between homologous genes in different species. Since most molecular mechanistic studies in plants have been performed in model species, including Arabidopsis and rice, comparative transcriptome analysis is particularly important for functional annotation of genes in diverse plant species. Many biological processes, such as embryo development, are highly conserved between different plant species. The challenge is to establish one-to-one mapping of the developmental stages between two species.

RESULTS

In this manuscript, we solve this problem by converting the gene expression patterns into co-expression networks and then apply network module finding algorithms to the cross-species co-expression network. We describe how such analyses are carried out using bash scripts for preliminary data processing followed by using the R programming language for module finding with a simulated annealing method. We also provide instructions on how to visualize the resulting co-expression networks across species.

CONCLUSIONS

We provide a comprehensive pipeline from installing software and downloading raw transcriptome data to predicting homologous genes and finding orthologous co-expression networks. From the example provided, we demonstrate the application of our method to reveal functional conservation and divergence of genes in two plant species.

The potential application of genome editing by using CRISPR/Cas9, and its engineered and ortholog variants for studying the transcription factors involved in the maintenance of phosphate homeostasis in model plants

Phosphorus (P), an essential macronutrient, is pivotal for growth and development of plants. Availability of phosphate (Pi), the only assimilable P, is often suboptimal in rhizospheres. Pi deficiency triggers an array of spatiotemporal adaptive responses including the differential regulation of several transcription factors (TFs). Studies on MYB TF PHR1 in Arabidopsis thaliana (Arabidopsis) and its orthologs OsPHRs in Oryza sativa (rice) have provided empirical evidence of their significant roles in the maintenance of Pi homeostasis. Since the functional characterization of PHR1 in 2001, several other TFs have now been identified in these model plants. This raised a pertinent question whether there are any likely interactions across these TFs. Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) system has provided an attractive paradigm for editing genome in plants. Here, we review the applications and challenges of this technique for genome editing of the TFs for deciphering the function and plausible interactions across them. This technology could thus provide a much-needed fillip towards engineering TFs for generating Pi use efficient plants for sustainable agriculture. Furthermore, we contemplate whether this technology could be a viable alternative to the controversial genetically modified (GM) rice or it may also eventually embroil into a limbo.

Water impacts nutrient dose responses genome-wide to affect crop production

Changes in nutrient dose have dramatic effects on gene expression and development. One outstanding question is whether organisms respond to changes in absolute nutrient amount (moles) vs. its concentration in water (molarity). This question is particularly relevant to plants, as soil drying can alter nutrient concentration, without changing its absolute amount. To compare the effects of amount vs. concentration, we expose rice to a factorial matrix varying the dose of nitrogen (N) and water (W) over a range of combinations, and quantify transcriptome and phenotype responses. Using linear models, we identify distinct dose responses to either N-moles, W-volume, N-molarity (N/W), or their synergistic interaction (N×W). Importantly, genes whose expression patterns are best explained by N-dose and W interactions (N/W or N×W) in seedlings are associated with crop outcomes in replicated field trials. Such N-by-W responsive genes may assist future efforts to develop crops resilient to increasingly arid, low nutrient soils.

Transcriptome analysis of a rice cultivar reveals the differentially expressed genes in response to wild and mutant strains of Xanthomonas oryzae pv. oryzae

Bacterial blight (BB), caused by Xanthomonas oryzae pv. oryzae (Xoo), is a devastating disease in most of the rice growing regions worldwide. Among the 42 BB resistance (R) genes, Xa23 is an executor R gene, conferring broad-spectrum disease resistance to all naturally occurring biotypes of Xoo. In this study, CBB23, a rice line carrying Xa23 gene, was inoculated with wild PXO99^A and its mutant, P99M2, to retrieve the differentially expressed genes (DEGs). RNA-Seq analysis retrieved 1,235 DEGs (p-value ≤ 0.05) at 12, 24, 36, and 48 hours of post inoculation (hpi). Gene ontology (GO) analysis classified the DEGs functionally into biological process, cellular component and molecular function. KEGG pathway analysis categorized the DEGs into 11 different pathways, and the ribosome is a prominent pathway followed by biosynthesis of phenylpropanoids. Gene co-expression network analysis identified the clusters of transcription factors (TFs) which may be involved in PXO99^A resistance. Additionally, we retrieved 67 differentially expressed TFs and 26 peroxidase responsive genes which may be involved in disease resistance mechanism. DEGs involved in the host-pathogen interaction, e.g., signaling mechanism, cell wall and plant hormones were identified. This data would be a valuable resource for researchers to identify the candidate genes associated with Xoo resistance.

Gene network mediated by WRKY13 to regulate resistance against sheath infecting fungi in rice (Oryza sativa L.)

OsWRKY13 TF gene is known to play a regulatory role of signaling in physiological pathways related to either development or disease resistance in rice plants. Rice cultivars IR 50 and TRY 3, resistant and susceptible respectively to sheath blight, TRY 3 and CO 43 resistant and susceptible respectively to sheath rot were challenged with fungal pathogens and disease scoring was carried out. Percent Disease Index (PDI) was significantly higher in susceptible varieties than resistant varieties. RT-PCR and qPCR analyses of WRKY13 using RNA extracted from the plant tissues revealed higher WRKY13 expression in resistant varieties (both diseases) upon pathogen challenge compared to uninfected control and also the susceptible varieties. To compute and evaluate the possible molecular mechanism for observed resistance correlated to WRKY13 gene expression, rice gene expression profiles against bacterial leaf blight and leaf blast disease from ROAD database were used to prioritize the locus IDs that were used as input in RiceNet v2 tool. The expression of WRKY13-regulated TIFY9 gene was predicted and validated using RT-PCR and qRT-PCR along with WRKY12 and PR2. All three genes showed induced expression in R. solani challenged sheath blight resistant variety. WRKY12 and PR2 expression in S. oryzae challenged sheath rot resistant variety was higher. Agrobacterium mediated transformation was carried out in rice plants using overexpression construct of WRKY13 (agroinfection in seeds of varieties susceptible to sheath blight and sheath rot, followed by selection in antibiotic media, germinating and hardening of putative transgenic lines). Based on qPCR analysis, the expression level of WRKY13 and the co-expression levels of WRKY12, TIFY9 and PR2 were found higher in PCR-positive T₁ plants compared to wild-type. Infection bioassays in the transgenic plants of both varieties revealed enhanced resistance to the pathogens. A mechanism by which WRKY13 would influence the MAPK cascade with TIFY9 acting as a mediator, is proposed.

The WRKY Transcription Factor GmWRKY12 Confers Drought and Salt Tolerance in Soybean

WRKYs are important regulators in plant development and stress responses. However, knowledge of this superfamily in soybean is limited. In this study, we characterized the drought- and salt-induced gene GmWRKY12 based on RNA-Seq and qRT-PCR. GmWRKY12, which is 714 bp in length, encoded 237 amino acids and grouped into WRKY II. The promoter region of GmWRKY12 included ABER4, MYB, MYC, GT-1, W-box and DPBF cis-elements, which possibly participate in abscisic acid (ABA), drought and salt stress responses. GmWRKY12 was minimally expressed in different tissues under normal conditions but highly expressed under drought and salt treatments. As a nucleus protein, GmWRKY12 was responsive to drought, salt, ABA and salicylic acid (SA) stresses. Using a transgenic hairy root assay, we further characterized the roles of GmWRKY12 in abiotic stress tolerance. Compared with control (Williams 82), overexpression of GmWRKY12 enhanced drought and salt tolerance, increased proline (Pro) content and decreased malondialdehyde (MDA) content under drought and salt treatment in transgenic soybean seedlings. These results may provide a basis to understand the functions of GmWRKY12 in abiotic stress responses in soybean.

Expression Characteristics and Functional Analysis of the ScWRKY3 Gene from Sugarcane

The plant-specific WRKY transcriptional regulatory factors have been proven to play vital roles in plant growth, development, and responses to biotic and abiotic stresses. However, there are few studies on the WRKY gene family in sugarcane (Saccharum spp.). In the present study, the characterization of a new subgroup, IIc WRKY protein ScWRKY3, from a Saccharum hybrid cultivar is reported. The ScWRKY3 protein was localized in the nucleus of Nicotiana benthamiana leaves and showed no transcriptional activation activity and no toxic effects on the yeast strain Y2HGold. An interaction between ScWRKY3 and a reported sugarcane protein ScWRKY4, was confirmed in the nucleus. The ScWRKY3 gene had the highest expression level in sugarcane stem pith. The transcript of ScWRKY3 was stable in the smut-resistant Saccharum hybrid cultivar Yacheng05-179, while it was down-regulated in the smut-susceptible Saccharum hybrid cultivar ROC22 during inoculation with the smut pathogen (Sporisorium scitamineum) at 0⁻72 h. ScWRKY3 was remarkably up-regulated by sodium chloride (NaCl), polyethylene glycol (PEG), and plant hormone abscisic acid (ABA), but it was down-regulated by salicylic acid (SA) and methyl jasmonate (MeJA). Moreover, transient overexpression of the ScWRKY3 gene in N. benthamiana indicated a negative regulation during challenges with the fungal pathogen Fusarium solani var. coeruleum or the bacterial pathogen Ralstonia solanacearum in N. benthamiana. The findings of the present study should accelerate future research on the identification and functional characterization of the WRKY family in sugarcane.

Principles and characteristics of the Arabidopsis WRKY regulatory network during early MAMP-triggered immunity

During microbe-associated molecular pattern-triggered immunity more than 5000 Arabidopsis genes are significantly altered in their expression, and the question arises, how such an enormous reprogramming of the transcriptome can be regulated in a safe and robust manner? For the WRKY transcription factors (TFs), which are important regulators of numerous defense responses, it appears that they act in a complex regulatory sub-network rather than in a linear fashion, which would be much more vulnerable to gene function loss either by pathogen-derived effectors or by mutations. In this study we employed RNA-seq, mass spectrometry and chromatin immunoprecipitation-seq to find evidence for and uncover principles and characteristics of this network. Upon flg22-treatment, one can distinguish between two sets of WRKY genes: constitutively expressed and induced WRKY genes. Prior to elicitation the induced WRKY genes appear to be maintained in a repressed state mainly by the constitutively expressed WRKY factors, which themselves appear to be regulated by non-WRKY TFs. Upon elicitation, induced WRKYs rapidly bind to induced WRKY gene promoters and by auto- and cross-regulation build up the regulatory network. Maintenance of this flg22-induced network appears highly robust as removal of three key WRKY factors can be physically and functionally compensated for by other WRKY family members.