Identification of 118 Arabidopsis transcription factor and 30 ubiquitin-ligase genes responding to chitin, a plant-defense elicitor

Transcriptomic analyses reveal complex and interconnected sucrose signaling cascades in developing seeds of castor bean

Seeds are highly specific organs that strongly sink sucrose resources from leaf and stem tissues to trigger seed metabolism and development. In particular, for heterotrophic non-green seeds, the potential molecular mechanism underlying sucrose-driven seed development remains unanswered. Castor bean (Ricinus communis L.), a typical non-green seed, has been considered as a model plant for seed biology study in dicotyledonous plants due to its heterotrophic seeds with persistent endosperms. In the present study, the fast-developing castor bean seeds were treated with exogenous sucrose and mannitol for four hours. The global transcriptomic data were obtained by high-throughput RNA-seq technique, resulting in 468 differentially expressed genes (DGEs). Further analyses revealed that sucrose functioned as both metabolic substrates and signal molecules. Specifically, 73 DGEs involved in carbohydrate and nitrogen metabolism, 42 differentially expressed transcription factors, and 35 DGEs involved in diverse signaling pathways such as auxin, brassinosteroid, ethelyene, cytokinin, gibberellin, and calcium signals, were identified, suggesting that the sucrose signaling pathway might have complex and multi-connected cross-talks with other signals to regulate castor bean seed development. Taken together, this study provides novel data to improve understanding of the potential molecular mechanisms of sucrose in regulating non-green seed development and storage reservoir accumulation during seed development.

Revisiting the Soybean GmNAC Superfamily

The NAC (NAM, ATAF, and CUC) genes encode transcription factors involved with the control of plant morph-physiology and stress responses. The release of the last soybean (Glycine max) genome assembly (Wm82.a2.v1) raised the possibility that new NAC genes would be present in the soybean genome. Here, we interrogated the last version of the soybean genome against a conserved NAC domain structure. Our analysis identified 32 putative novel NAC genes, updating the superfamily to 180 gene members. We also organized the genes in 15 phylogenetic subfamilies, which showed a perfect correlation among sequence conservation, expression profile, and function of orthologous Arabidopsis thaliana genes and NAC soybean genes. To validate our in silico analyses, we monitored the stress-mediated gene expression profiles of eight new NAC-genes by qRT-PCR and monitored the GmNAC senescence-associated genes by RNA-seq. Among ER stress, osmotic stress and salicylic acid treatment, all the novel tested GmNAC genes responded to at least one type of stress, displaying a complex expression profile under different kinetics and extension of the response. Furthermore, we showed that 40% of the GmNACs were differentially regulated by natural leaf senescence, including eight (8) newly identified GmNACs. The developmental and stress-responsive expression profiles of the novel NAC genes fitted perfectly with their phylogenetic subfamily. Finally, we examined two uncharacterized senescence-associated proteins, GmNAC065 and GmNAC085, and a novel, previously unidentified, NAC protein, GmNAC177, and showed that they are nuclear localized, and except for GmNAC065, they display transactivation activity in yeast. Consistent with a role in leaf senescence, transient expression of GmNAC065 and GmNAC085 induces the appearance of hallmarks of leaf senescence, including chlorophyll loss, leaf yellowing, lipid peroxidation and accumulation of H₂O₂. GmNAC177 was clustered to an uncharacterized subfamily but in close proximity to the TIP subfamily. Accordingly, it was rapidly induced by ER stress and by salicylic acid under late kinetic response and promoted cell death in planta. Collectively, our data further substantiated the notion that the GmNAC genes display functional and expression profiles consistent with their phylogenetic relatedness and established a complete framework of the soybean NAC superfamily as a foundation for future analyses.

MorphDB: Prioritizing Genes for Specialized Metabolism Pathways and Gene Ontology Categories in Plants

Recent times have seen an enormous growth of "omics" data, of which high-throughput gene expression data are arguably the most important from a functional perspective. Despite huge improvements in computational techniques for the functional classification of gene sequences, common similarity-based methods often fall short of providing full and reliable functional information. Recently, the combination of comparative genomics with approaches in functional genomics has received considerable interest for gene function analysis, leveraging both gene expression based guilt-by-association methods and annotation efforts in closely related model organisms. Besides the identification of missing genes in pathways, these methods also typically enable the discovery of biological regulators (i.e., transcription factors or signaling genes). A previously built guilt-by-association method is MORPH, which was proven to be an efficient algorithm that performs particularly well in identifying and prioritizing missing genes in plant metabolic pathways. Here, we present MorphDB, a resource where MORPH-based candidate genes for large-scale functional annotations (Gene Ontology, MapMan bins) are integrated across multiple plant species. Besides a gene centric query utility, we present a comparative network approach that enables researchers to efficiently browse MORPH predictions across functional gene sets and species, facilitating efficient gene discovery and candidate gene prioritization. MorphDB is available at http://bioinformatics.psb.ugent.be/webtools/morphdb/morphDB/index/. We also provide a toolkit, named "MORPH bulk" (https://github.com/arzwa/morph-bulk), for running MORPH in bulk mode on novel data sets, enabling researchers to apply MORPH to their own species of interest.

Comparative Transcriptome Profiles of Near-Isogenic Hexaploid Wheat Lines Differing for Effective Alleles at the 2DL FHB Resistance QTL

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum, represents one of the major wheat diseases worldwide, determining severe yield losses and reduction of grain quality due to the accumulation of mycotoxins. The molecular response associated with the wheat 2DL FHB resistance QTL was mined through a comprehensive transcriptomic analysis of the early response to F. graminearum infection, at 3 days post-inoculation, in spikelets and rachis. The analyses were conducted on two near isogenic lines (NILs) differing for the presence of the 2DL QTL (2-2618, resistant 2DL+ and 2-2890, susceptible null). The general response to fungal infection in terms of mRNAs accumulation trend was similar in both NILs, even though involving an higher number of DEGs in the susceptible NIL, and included down-regulation of the primary and energy metabolism, up-regulation of enzymes implicated in lignin and phenylpropanoid biosynthesis, activation of hormons biosynthesis and signal transduction pathways and genes involved in redox homeostasis and transcriptional regulation. The search for candidate genes with expression profiles associated with the 2DL QTL for FHB resistance led to the discovery of processes differentially modulated in the R and S NILs related to cell wall metabolism, sugar and JA signaling, signal reception and transduction, regulation of the redox status and transcription factors. Wheat FHB response-related miRNAs differentially regulated were also identified as putatively implicated in the superoxide dismutase activities and affecting genes regulating responses to biotic/abiotic stresses and auxin signaling. Altered gene expression was also observed for fungal non-codingRNAs. The putative targets of two of these were represented by the wheat gene WIR1A, involved in resistance response, and a gene encoding a jacalin-related lectin protein, which participate in biotic and abiotic stress response, supporting the presence of a cross-talk between the plant and the fungus.

Specificity versus redundancy in the RAP2.4 transcription factor family of Arabidopsis thaliana: transcriptional regulation of genes for chloroplast peroxidases

BACKGROUND

The Arabidopsis ERFIb / RAP2.4 transcription factor family consists of eight members with highly conserved DNA binding domains. Selected members have been characterized individually, but a systematic comparison is pending. The redox-sensitive transcription factor RAP2.4a mediates chloroplast-to-nucleus redox signaling and controls induction of the three most prominent chloroplast peroxidases, namely 2-Cys peroxiredoxin A (2CPA) and thylakoid- and stromal ascorbate peroxidase (tAPx and sAPx). To test the specificity and redundancy of RAP2.4 transcription factors in the regulation of genes for chloroplast peroxidases, we compared the DNA-binding sites of the transcription factors in tertiary structure models, analyzed transcription factor and target gene regulation by qRT-PCR in RAP2.4, 2-Cys peroxiredoxin and ascorbate peroxidase T-DNA insertion lines and RAP2.4 overexpressing lines of Arabidopsis thaliana and performed promoter binding studies.

RESULTS

All RAP2.4 proteins bound the tAPx promoter, but only the four RAP2.4 proteins with identical DNA contact sites, namely RAP2.4a, RAP2.4b, RAP2.4d and RAP2.4h, interacted stably with the redox-sensitive part of the 2CPA promoter. Gene expression analysis in RAP2.4 knockout lines revealed that RAP2.4a is the only one supporting 2CPA and chloroplast APx expression. Rap2.4h binds to the same promoter region as Rap2.4a and antagonizes 2CPA expression. Like the other six RAP2.4 proteins, Rap2.4 h promotes APx mRNA accumulation. Chloroplast ROS signals induced RAP2.4b and RAP2.4d expression, but these two transcription factor genes are (in contrast to RAP2.4a) insensitive to low 2CP availability, and their expression decreased in APx knockout lines. RAP2.4e and RAP2.4f gradually responded to chloroplast APx availability and activated specifically APx expression. These transcription factors bound, like RAP2.4c and RAP2.4g, the tAPx promoter, but hardly the 2CPA promoter.

CONCLUSIONS

The RAP2.4 transcription factors form an environmentally and developmentally regulated transcription factor network, in which the various members affect the expression intensity of the others. Within the transcription factor family, RAP2.4a has a unique function as a general transcriptional activator of chloroplast peroxidase activity. The other RAP2.4 proteins mediate the fine-control and adjust the relative availability of 2CPA, sAPx and tAPx.

Arabidopsis E3 ubiquitin ligase PLANT U-BOX13 (PUB13) regulates chitin receptor LYSIN MOTIF RECEPTOR KINASE5 (LYK5) protein abundance

Long-chain chitooligosaccharides are fungal microbe-associated molecular patterns (MAMPs) that are recognized by LYSIN MOTIF RECEPTOR KINASE5 (LYK5), inducing the formation of a complex with CHITIN ELICITOR RECEPTOR KINASE1 (CERK1). Formation of this complex leads to activation of the CERK1 intracellular kinase domain and induction of plant innate immunity in Arabidopsis. We found that addition of chitooctaose induced LYK5 protein accumulation as a result of de novo gene expression and the inhibition of LYK5 protein degradation. Screening the putative E3 ligases for interaction with LYK5 identified PLANT U-BOX13 (PUB13), which complexed with LYK5, but this complex dissociated upon addition of chitooctaose. Consistent with these results, LYK5 protein abundance was higher in pub13 mutants compared with the wild type without chitooctaose treatment, while similar abundance was detected with the addition of chitooctaose. The pub13 mutants showed hypersensitivity to chitooctaose-induced rapid responses, such as the production of reactive oxygen species (ROS) and mitogen-activated protein (MAP) kinase phosphorylation, but exhibited normal responses to subsequent long-term chitooctaose treatment, such as gene expression and callose deposition. In addition, PUB13 could ubiquitinate the LYK5 kinase domain in vitro. Taken together, our results suggest an important regulatory function for the turnover of LYK5 mediated by the E3 ligase PUB13.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses

BACKGROUND

Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt).

RESULTS

The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors.

CONCLUSIONS

Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Transgenic expression of fungal accessory hemicellulases in Arabidopsis thaliana triggers transcriptional patterns related to biotic stress and defense response

The plant cell wall is an abundant and renewable resource for lignocellulosic applications such as the production of biofuel. Due to structural and compositional complexities, the plant cell wall is, however, recalcitrant to hydrolysis and extraction of platform sugars. A cell wall engineering strategy to reduce this recalcitrance makes use of microbial cell wall modifying enzymes that are expressed directly in plants themselves. Previously, we constructed transgenic Arabidopsis thaliana constitutively expressing the fungal hemicellulases: Phanerochaete carnosa glucurnoyl esterase (PcGCE) and Aspergillus nidulans α-arabinofuranosidase (AnAF54). While the PcGCE lines demonstrated improved xylan extractability, they also displayed chlorotic leaves leading to the hypothesis that expression of such enzymes in planta resulted in plant stress. The objective of this study is to investigate the impact of transgenic expression of the aforementioned microbial hemicellulases in planta on the host arabidopsis. More specifically, we investigated transcriptome profiles by short read high throughput sequencing (RNAseq) from developmentally distinct parts of the plant stem. When compared to non-transformed wild-type plants, a subset of genes was identified that showed differential transcript abundance in all transgenic lines and tissues investigated. Intriguingly, this core set of genes was significantly enriched for those involved in plant defense and biotic stress responses. While stress and defense-related genes showed increased transcript abundance in the transgenic plants regardless of tissue or genotype, genes involved in photosynthesis (light harvesting) were decreased in their transcript abundance potentially reflecting wide-spread effects of heterologous microbial transgene expression and the maintenance of plant homeostasis. Additionally, an increase in transcript abundance for genes involved in salicylic acid signaling further substantiates our finding that transgenic expression of microbial cell wall modifying enzymes induces transcriptome responses similar to those observed in defense responses.

Short-Chain Chitin Oligomers: Promoters of Plant Growth

Chitin is the second most abundant biopolymer in nature after cellulose, and it forms an integral part of insect exoskeletons, crustacean shells, krill and the cell walls of fungal spores, where it is present as a high-molecular-weight molecule. In this study, we showed that a chitin oligosaccharide of lower molecular weight (tetramer) induced genes in Arabidopsis that are principally related to vegetative growth, development and carbon and nitrogen metabolism. Based on plant responses to this chitin tetramer, a low-molecular-weight chitin mix (CHL) enriched to 92% with dimers (2mer), trimers (3mer) and tetramers (4mer) was produced for potential use in biotechnological processes. Compared with untreated plants, CHL-treated plants had increased in vitro fresh weight (10%), radicle length (25%) and total carbon and nitrogen content (6% and 8%, respectively). Our data show that low-molecular-weight forms of chitin might play a role in nature as bio-stimulators of plant growth, and they are also a known direct source of carbon and nitrogen for soil biomass. The biochemical properties of the CHL mix might make it useful as a non-contaminating bio-stimulant of plant growth and a soil restorer for greenhouses and fields.

Metabolo-transcriptome profiling of barley reveals induction of chitin elicitor receptor kinase gene (HvCERK1) conferring resistance against Fusarium graminearum

We report plausible disease resistance mechanisms induced by barley resistant genotype CI89831 against Fusarium head blight (FHB) based on metabolo-transcriptomics approach. We identified HvCERK1 as a candidate gene for FHB resistance, which is functional in resistant genotype CI9831 but non-functional in susceptible cultivars H106-371 and Zhedar-2. For the first time, we were able to show a hierarchy of regulatory genes that regulated downstream biosynthetic genes that eventually produced resistance related metabolites that reinforce the cell walls to contain the pathogen progress in plant. The HvCERK1 can be used for replacing in susceptible commercial cultivars, if non-functional, based on genome editing. Fusarium head blight (FHB) management is a great challenge in barley and wheat production worldwide. Though barley genome sequence and advanced omics technologies are available, till date none of the resistance mechanisms has been clearly deciphered. Hence, this study was aimed at identifying candidate gene(s) and elucidating resistance mechanisms induced by barley resistant genotype CI9831 based on integrated metabolomics and transcriptomics approach. Following Fusarium graminearum infection, we identified accumulation of specific set of induced secondary metabolites, belonging to phenylpropanoid, hydroxycinnamic acid (HCAA) and jasmonic acid pathways, and their biosynthetic genes. In association with these, receptor kinases such as chitin elicitor receptor kinase (HvCERK1) and protein kinases such as MAP kinase 3 (HvMPK3) and MAPK substrate 1 (HvMKS1), and transcription factors such as HvERF1/5, HvNAC42, HvWRKY23 and HvWRKY70 were also found upregulated with high fold change. Polymorphism studies across three barley genotypes confirmed the presence of mutations in HvCERK1 gene in two susceptible genotypes, isolating this gene as a potential candidate for FHB resistance. Further, the silencing of functional HvCERK1 gene in the resistant genotype CI9831, followed by gene expression and metabolite analysis revealed its role as an elicitor recognition receptor that triggered downstream regulatory genes, which in turn, regulated downstream metabolic pathway genes to biosynthesize resistance related (RR) metabolites to contain the pathogen to spikelet infection. A putative model on metabolic pathway regulation is proposed.

Transcriptome response of Lolium arundinaceum to its fungal endophyte Epichloë coenophiala

Tall fescue (Lolium arundinaceum) is one of the primary forage and turf grasses in temperate regions of the world. A number of favourable characteristics of tall fescue are enhanced by its seed-transmissible fungal symbiont (endophyte) Epichloë coenophiala. Our approach was to assemble the tall fescue transcriptome, then identify differentially expressed genes (DEGs) for endophyte-symbiotic (E+) vs endophyte-free (E-) clones in leaf blades, pseudostems, crowns and roots. RNA-seq reads were used to construct a tall fescue reference transcriptome and compare gene expression profiles. Over all tissues examined, 478 DEGs were identified between the E+ and E- clones for at least one tissue (more than two-fold; P < 0.0001, 238 E+ > E- and 240 E- > E+), although no genes were differentially expressed in all four tissues. Gene ontology (GO) terms, GO:0010200 (response to chitin), GO:0002679 (respiratory burst during defence response) and GO:0035556 (intracellular signal transduction) were significantly overrepresented among 25 E- > E+ DEGs in leaf blade, and a number of other DEGs were associated with defence and abiotic response. In particular, endophyte effects on various WRKY transcription factors may have implications for symbiotic stability, endophyte distribution in the plant, or defence against pathogens.

Expression Pattern of ERF Gene Family under Multiple Abiotic Stresses in Populus simonii × P. nigra

Identification of gene expression patterns of key genes across multiple abiotic stresses is critical for mechanistic understanding of stress resistance in plant. In the present study, we identified differentially expressed genes (DEGs) in di-haploid Populus simonii × P. nigra under respective stresses of NaCl, KCl, CdCl_2, and PEG. On the basis of RNA-Seq, we detected 247 DEGs that are shared by the four stresses in wild type poplar, and mRNA abundance of the DEGs were validated in transgenic poplar overexpressing ERF76 gene by RNA-Seq and RT-qPCR. Results from gene ontology analysis indicated that these genes are enriched in significant pathways, such as phenylpropanoid biosynthesis, phenylalanine metabolism, starch and sucrose metabolism, and plant hormone signal transduction. Ethylene response factor (ERF) gene family plays significant role in plant abiotic stress responses. We also investigated expression pattern of ERF gene family under the four stresses. The ERFs and DEGs share similar expression pattern across the four stresses. The transgenic poplar is superior to WT in morphologic, physiological and biochemical traits, which demonstrated the ERF76 gene plays a significant role in stress resistance. These studies will give a rise in understanding the stress response mechanisms in poplar.

Transcriptomic analysis reveals key early events of narciclasine signaling in Arabidopsis root apex

Histochemical staining and RNA-seq data demonstrated that the ROS- and plant hormone-regulated stress responses are the key early events of narciclasine signaling in Arabidopsis root cells. Narciclasine, an amaryllidaceae alkaloid isolated from Narcissus tazetta bulbs, employs a broad range of functions on plant development and growth. However, its molecular interactions that modulate these roles in plants are not fully understood. To elucidate the global responses of Arabidopsis roots to short-term narciclasine exposure, we first measured the accumulation of H₂O₂ and O₂^- with histochemical staining, and then profiled the gene expression pattern in Arabidopsis root tips treated with 0.5 µM narciclasine across different exposure times by RNA-seq. Physiological measurements showed a significant increase in H₂O₂ began at 30-60 min of narciclasine treatment and O₂^- accumulated by 120 min. Compared with controls, 236 genes were upregulated and 54 genes were downregulated with 2 h of narciclasine treatment, while 968 genes were upregulated and 835 genes were downregulated with 12 h of treatment. The Gene Ontology analysis revealed that the differentially expressed genes were highly enriched during oxidative stress, including those involved in the "regulation of transcription", "response to oxidative stress", "plant-pathogen interaction", "ribonucleotide binding", "plant cell wall organization", and "ribosome biogenesis". Moreover, Kyoto Encyclopedia of Genes and Genomes pathway enrichment statistics suggested that carbohydrate metabolism, amino acid metabolism, amino sugar and nucleotide sugar metabolism, and biosynthesis of phenylpropanoid and secondary metabolites were significantly inhibited by 12 h of narciclasine exposure. Hence, our results demonstrate that hormones and H₂O₂ are important regulators of narciclasine signaling and help to uncover the factors involved in the molecular interplay between narciclasine and phytohormones in Arabidopsis root cells.

A novel exopolysaccharide elicitor from endophytic fungus Gilmaniella sp. AL12 on volatile oils accumulation in Atractylodes lancea

Endophytes and plants can establish specific long-term symbiosis through the accumulation of secondary metabolites. Previous studies have shown that the endophytic fungus Gilmaniella sp. AL12 can stimulate Atractylodes lancea to produce volatile oils. The purpose of this report is to investigate key factors involved in the stimulation of A. lancea by AL12 and reveal the mechanism. We identified the active component from AL12 as an extracellular mannan with a polymerization degree of 26–42. Differential membrane proteomics of A. lancea was performed by 2D electrophoresis. The results showed that there were significant differences in the expression of 83 proteins. Based on these results, we conclude that AL12 secreted mannan contributes to the antagonistic balance seen in interactions between AL12 and A. lancea. One portion of the mannan was degraded to mannose for hexokinase activation, promoting photosynthesis and energy metabolism, with a potential metabolic fluxes flowing towards terpenoid biosynthesis. The other portion of the mannan directly enhanced autoimmunity of A. lancea through G protein-mediated signal transduction and the mannan-binding lectin pathway. Volatile oil accumulation was ultimately promoted in subsequent defense reactions. This study provides a new perspective on the regulation of secondary metabolites by endophytic fungal elicitors in medicinal plants.

Exposure to Sound Vibrations Lead to Transcriptomic, Proteomic and Hormonal Changes in Arabidopsis

Sound vibration (SV) is considered as an external mechanical force that modulates plant growth and development like other mechanical stimuli (e.g., wind, rain, touch and vibration). A number of previous and recent studies reported developmental responses in plants tailored against SV of varied frequencies. This strongly suggests the existence of sophisticated molecular mechanisms for SV perception and signal transduction. Despite this there exists a huge gap in our understanding regarding the SV-mediated molecular alterations, which is a prerequisite to gain insight into SV-mediated plant development. Herein, we investigated the global gene expression changes in Arabidopsis thaliana upon treatment with five different single frequencies of SV at constant amplitude for 1 h. As a next step, we also studied the SV-mediated proteomic changes in Arabidopsis. Data suggested that like other stimuli, SV also activated signature cellular events, for example, scavenging of reactive oxygen species (ROS), alteration of primary metabolism, and hormonal signaling. Phytohormonal analysis indicated that SV-mediated responses were, in part, modulated by specific alterations in phytohormone levels; especially salicylic acid (SA). Notably, several touch regulated genes were also up-regulated by SV treatment suggesting a possible molecular crosstalk among the two mechanical stimuli, sound and touch. Overall, these results provide a molecular basis to SV triggered global transcriptomic, proteomic and hormonal changes in plant.

Insight into the AP2/ERF transcription factor superfamily in sesame and expression profiling of DREB subfamily under drought stress

BACKGROUND

Sesame is an important oilseed crop mainly grown in inclement areas with high temperatures and frequent drought. Thus, drought constitutes one of the major constraints of its production. The AP2/ERF is a large family of transcription factors known to play significant roles in various plant processes including biotic and abiotic stress responses. Despite their importance, little is known about sesame AP2/ERF genes. This constitutes a limitation for drought-tolerance candidate genes discovery and breeding for tolerance to water deficit.

RESULTS

One hundred thirty-two AP2/ERF genes were identified in the sesame genome. Based on the number of domains, conserved motifs, genes structure and phylogenetic analysis including 5 relatives species, they were classified into 24 AP2, 41 DREB, 61 ERF, 4 RAV and 2 Soloist. The number of sesame AP2/ERF genes was relatively few compared to that of other relatives, probably due to gene loss in ERF and DREB subfamilies during evolutionary process. In general, the AP2/ERF genes were expressed differently in different tissues but exhibited the highest expression levels in the root. Mostly all DREB genes were responsive to drought stress. Regulation by drought is not specific to one DREB group but depends on the genes and the group A6 and A1 appeared to be more actively expressed to cope with drought.

CONCLUSIONS

This study provides insights into the classification, evolution and basic functional analysis of AP2/ERF genes in sesame which revealed their putative involvement in multiple tissue-/developmental stages. Out of 20 genes which were significantly up- /down-regulated under drought stress, the gene AP2si16 may be considered as potential candidate gene for further functional validation as well for utilization in sesame improvement programs for drought stress tolerance.

Characterization and expression profiling of MYB transcription factors against stresses and during male organ development in Chinese cabbage (Brassica rapa ssp. pekinensis)

MYB proteins comprise a large family of plant transcription factors that play regulatory roles in different biological processes such as plant development, metabolism, and defense responses. To gain insight into this gene superfamily and to elucidate its roles in stress resistance, we performed a comprehensive genome-wide identification, characterization, and expression analysis of MYB genes in Chinese cabbage (Brassica rapa ssp. pekinensis). We identified 475 Chinese cabbage MYB genes, among which most were from R2R3-MYB (256 genes) and MYB-related (202) subfamilies. Analysis of sequence characteristics, phylogenetic classification, and protein motif structures confirmed the existence of several categories (1R, 2R, 3R, 4R, and 5R) of Chinese cabbage MYB genes, which is comparable with MYB genes of other crops. An extensive in silico functional analysis, based on established functional properties of MYB genes from different crop species, revealed 11 and four functional clades within the Chinese cabbage R2R3-MYB and MYB-related subfamilies, respectively. In this study, we reported a MYB-like group within the MYB-related subfamily contains 77 MYB genes. Expression analysis using low temperature-treated whole-genome microarray data revealed variable transcript abundance of 1R/2R/3R/4R/5R-MYB genes in 11 clusters between two inbred lines of Chinese cabbage, Chiifu and Kenshin, which differ in cold tolerance. In further validation tests, we used qRT-PCR to examine the cold-responsive expression patterns of 27 BrMYB genes; surprisingly, the MYB-related genes were induced more highly than the R2R3-MYB genes. In addition, we identified 10 genes with corresponsive expression patterns from a set of salt-, drought-, ABA-, JA-, and SA-induced R2R3-MYB genes. We identified 11 R2R3-MYBs functioning in resistance against biotic stress, including 10 against Fusarium oxysporum f.sp. conglutinans and one against Pectobacterium carotovoram subsp. caratovorum. Furthermore, based on organ-specific expression data, we identified nine R2R3-MYBs that were constitutively expressed in male reproductive tissue, which may provide an important key for studying male sterility in Chinese cabbage. The extensive annotation and transcriptome profiling reported in this study will be useful for understanding the involvement of MYB genes in stress resistance in addition to their growth regulatory functions, ultimately providing the basis for functional characterization and exploitation of the candidate MYB genes for genetic engineering of Chinese cabbage.

Tomato Plant Proteins Actively Responding to Fungal Applications and Their Role in Cell Physiology

The pattern of protein induction in tomato plants has been investigated after the applications of pathogenic and non-pathogenic fungal species. Moreover, particular roles of the most active protein against biological applications were also determined using chromatographic techniques. Alternaria alternata and Penicillium oxalicum were applied as a pathogenic and non-pathogenic fungal species, respectively. Protein profile analysis revealed that a five protein species (i.e., protein 1, 6, 10, 12, and 13) possessed completely coupled interaction with non-pathogenic inducer application (P. oxalicum). However, three protein species (i.e., 10, 12, and 14) recorded a strong positive interaction with both fungal species. Protein 14 exhibited the maximum interaction with fungal applications, and its role in plant metabolism was studied after its identification as protein Q9M1W6. It was determined that protein Q1M1W6 was involved in guaiacyl lignin biosynthesis, and its inhibition increased the coumarin contents in tomato plants. Moreover, it was also observed that the protein Q9M1W6 takes significant part in the biosynthesis of jasmonic acid and Indole acetic acid contents, which are defense and growth factors of tomato plants. The study will help investigators to design fundamental rules of plant proteins affecting cell physiology under the influence of external fungal applications.

An unbiased nuclear proteomics approach reveals novel nuclear protein components that participates in MAMP-triggered immunity

(MAMP)-triggered immunity (MTI) is the first layer of molecular defense encountered by pathogens. Genetic screens have contributed to our knowledge of MTI, but are limited to phenotype-causing mutations. Here we attempt to identify novel factors involved in the early event leading to plant MTI by comparing the nuclear proteomes of two Arabidopsis genotypes treated with chitosan. Our approach revealed that following chitosan treatment, cerk1 plants had many nuclear accumulating proteins in common, but also some unique ones, when compared with Col-0 plants. Analysis of the identified proteins revealed a nuclear accumulation of DNA-modifying enzymes, RNA-binding proteins and ribosomal proteins. Our results demonstrate that nuclear proteomic is a valid, phenotype-independent approach to uncover factor involved in cellular processes.

Transcriptome Profiling of Resistance to Fusarium oxysporum f. sp. conglutinans in Cabbage (Brassica oleracea) Roots

Fusarium wilt caused by Fusarium oxysporum f. sp. conglutinans (FOC) is a destructive disease of Brassica crops, which results in severe yield losses. There is little information available about the mechanism of disease resistance. To obtain an overview of the transcriptome profiles in roots of R4P1, a Brassica oleracea variety that is highly resistant to fusarium wilt, we compared the transcriptomes of samples inoculated with FOC and samples inoculated with distilled water. RNA-seq analysis generated more than 136 million 100-bp clean reads, which were assembled into 62,506 unigenes (mean size = 741 bp). Among them, 49,959 (79.92%) genes were identified based on sequence similarity searches, including SwissProt (29,050, 46.47%), Gene Ontology (GO) (33,767, 54.02%), Clusters of Orthologous Groups (KOG) (14,721, 23.55%) and Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG) (12,974, 20.76%) searches; digital gene expression analysis revealed 885 differentially expressed genes (DEGs) between infected and control samples at 4, 12, 24 and 48 hours after inoculation. The DEGs were assigned to 31 KEGG pathways. Early defense systems, including the MAPK signaling pathway, calcium signaling and salicylic acid-mediated hypersensitive response (SA-mediated HR) were activated after pathogen infection. SA-dependent systemic acquired resistance (SAR), ethylene (ET)- and jasmonic (JA)-mediated pathways and the lignin biosynthesis pathway play important roles in plant resistance. We also analyzed the expression of defense-related genes, such as genes encoding pathogenesis-related (PR) proteins, UDP-glycosyltransferase (UDPG), pleiotropic drug resistance, ATP-binding cassette transporters (PDR-ABC transporters), myrosinase, transcription factors and kinases, which were differentially expressed. The results of this study may contribute to efforts to identify and clone candidate genes associated with disease resistance and to uncover the molecular mechanism underlying FOC resistance in cabbage.

Identification of Powdery Mildew Responsive Genes in Hevea brasiliensis through mRNA Differential Display

Powdery mildew is an important disease of rubber trees caused by Oidium heveae B. A. Steinmann. As far as we know, none of the resistance genes related to powdery mildew have been isolated from the rubber tree. There is little information available at the molecular level regarding how a rubber tree develops defense mechanisms against this pathogen. We have studied rubber tree mRNA transcripts from the resistant RRIC52 cultivar by differential display analysis. Leaves inoculated with the spores of O. heveae were collected from 0 to 120 hpi in order to identify pathogen-regulated genes at different infection stages. We identified 78 rubber tree genes that were differentially expressed during the plant–pathogen interaction. BLAST analysis for these 78 ESTs classified them into seven functional groups: cell wall and membrane pathways, transcription factor and regulatory proteins, transporters, signal transduction, phytoalexin biosynthesis, other metabolism functions, and unknown functions. The gene expression for eight of these genes was validated by qRT-PCR in both RRIC52 and the partially susceptible Reyan 7-33-97 cultivars, revealing the similar or differential changes of gene expressions between these two cultivars. This study has improved our overall understanding of the molecular mechanisms of rubber tree resistance to powdery mildew.

Plant resistance against the parasitic nematode Heterodera schachtii is mediated by MPK3 and MPK6 kinases, which are controlled by the MAPK phosphatase AP2C1 in Arabidopsis

Plant-parasitic cyst nematodes infect plants and form highly sophisticated feeding sites in roots. It is not known which plant cell signalling mechanisms trigger plant defence during the early stages of nematode parasitism. Mitogen-activated protein kinases (MAPKs) are central components of protein phosphorylation cascades transducing extracellular signals to plant defence responses. MAPK phosphatases control kinase activities and the signalling outcome. The involvement and the role of MPK3 and MPK6, as well as the MAPK phosphatase AP2C1, is demonstrated during parasitism of the beet cyst nematode Heterodera schachtii in Arabidopsis. Our data reveal notable activation patterns of plant MAPKs and the induction of AP2C1 suggesting the attenuation of defence signalling in plant cells during early nematode infection. It is demonstrated that the ap2c1 mutant that is lacking AP2C1 is more attractive but less susceptible to nematodes compared with the AP2C1-overexpressing line. This implies that the function of AP2C1 is a negative regulator of nematode-induced defence. By contrast, the enhanced susceptibility of mpk3 and mpk6 plants indicates a positive role of stress-activated MAPKs in plant immunity against nematodes. Evidence is provided that phosphatase AP2C1, as well as AP2C1-targeted MPK3 and MPK6, are important regulators of plant-nematode interaction, where the co-ordinated action of these signalling components ensures the timely activation of plant defence.

Transcriptional regulatory networks in Arabidopsis thaliana during single and combined stresses

Differentially evolved responses to various stress conditions in plants are controlled by complex regulatory circuits of transcriptional activators, and repressors, such as transcription factors (TFs). To understand the general and condition-specific activities of the TFs and their regulatory relationships with the target genes (TGs), we have used a homogeneous stress gene expression dataset generated on ten natural ecotypes of the model plant Arabidopsis thaliana, during five single and six combined stress conditions. Knowledge-based profiles of binding sites for 25 stress-responsive TF families (187 TFs) were generated and tested for their enrichment in the regulatory regions of the associated TGs. Condition-dependent regulatory sub-networks have shed light on the differential utilization of the underlying network topology, by stress-specific regulators and multifunctional regulators. The multifunctional regulators maintain the core stress response processes while the transient regulators confer the specificity to certain conditions. Clustering patterns of transcription factor binding sites (TFBS) have reflected the combinatorial nature of transcriptional regulation, and suggested the putative role of the homotypic clusters of TFBS towards maintaining transcriptional robustness against cis-regulatory mutations to facilitate the preservation of stress response processes. The Gene Ontology enrichment analysis of the TGs reflected sequential regulation of stress response mechanisms in plants.

Transcriptome-wide high-throughput deep m(6)A-seq reveals unique differential m(6)A methylation patterns between three organs in Arabidopsis thaliana

Background

m⁶A is a ubiquitous RNA modification in eukaryotes. Transcriptome-wide m⁶A patterns in Arabidopsis have been assayed recently. However, differential m⁶A patterns between organs have not been well characterized.

Results

Over two-third of the transcripts in Arabidopsis are modified by m⁶A. In contrast to a recent observation of m⁶A enrichment in 5′ mRNA, we find that m⁶A is distributed predominantly near stop codons. Interestingly, 85 % of the modified transcripts show high m⁶A methylation extent compared to their transcript level. The 290 highly methylated transcripts are mainly associated with transporters, stress responses, redox, regulation factors, and some non-coding RNAs. On average, the proportion of transcripts showing differential methylation between two plant organs is higher than that showing differential transcript levels. The transcripts with extensively higher m⁶A methylation in an organ are associated with the unique biological processes of this organ, suggesting that m⁶A may be another important contributor to organ differentiation in Arabidopsis. Highly expressed genes are relatively less methylated and vice versa, and different RNAs have distinct m⁶A patterns, which hint at mRNA fate. Intriguingly, most of the transposable element transcripts maintained a fragmented form with a relatively low transcript level and high m⁶A methylation in the cells.

Conclusions

This is the first study to comprehensively analyze m⁶A patterns in a variety of RNAs, the relationship between transcript level and m⁶A methylation extent, and differential m⁶A patterns across organs in Arabidopsis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0839-2) contains supplementary material, which is available to authorized users.

Sustained exposure to abscisic acid enhances the colonization potential of the mutualist fungus Piriformospora indica on Arabidopsis thaliana roots

Root colonization by the beneficial fungus Piriformospora indica is controlled by plant innate immunity, but factors that channel this interaction into a mutualistic relationship are not known. We have explored the impact of abscisic acid (ABA) and osmotic stress on the P. indica interaction with Arabidopsis thaliana. The activation of plant innate immunity in roots was determined by measuring the concentration of the phytoalexin camalexin and expression of transcription factors regulating the biosynthesis of tryptophan-related defence metabolites. Furthermore, the impact of the fungus on the content of ABA, salicylic acid, jasmonic acid (JA) and JA-related metabolites was examined. We demonstrated that treatment with exogenous ABA or the ABA analogue pyrabactin increased fungal colonization efficiency without impairment of plant fitness. Concomitantly, ABA-deficient mutants of A. thaliana (aba1-6 and aba2-1) were less colonized, while plants exposed to moderate stress were more colonized than corresponding controls. Sustained exposure to ABA attenuated expression of transcription factors MYB51, MYB122 and WRKY33 in roots upon P. indica challenge or chitin treatment, and prevented an increase in camalexin content. The results indicate that ABA can strengthen the interaction with P. indica as a consequence of its impact on plant innate immunity. Consequently, ABA will be relevant for the establishment and outcome of the symbiosis under stress conditions.

WRKY40 and WRKY6 act downstream of the green leaf volatile E-2-hexenal in Arabidopsis

Plants are known to be responsive to volatiles, but knowledge about the molecular players involved in transducing their perception remains scarce. We study the response of Arabidopsis thaliana to E-2-hexenal, one of the green leaf volatiles (GLV) that is produced upon wounding, herbivory or infection with pathogens. We have taken a transcriptomics approach to identify genes that are induced by E-2-hexenal, but not by defence hormones or other GLVs. Furthermore, by studying the promoters of early E-2-hexenal-induced genes we determined that the only statistically enriched cis-element was the W-box motif. Since members of the plant-specific family of WRKY transcription factors act in trans on this cis-element, we focused on WRKY6, 40 and 53 that were most strongly induced by E-2-hexenal. Root elongation of Arabidopsis seedlings of the wrky40 wrky6 double mutant was much less inhibited than in wt plants, similar to the E-2-hexenal-responsive mutant her1, which is perturbed in γ-amino butyric acid (GABA) metabolism. The induction of several of the E-2-hexenal-specific genes was much higher in the wrky40, wrky6 or wrky40 wrky6 mutants, including GAD4, a glutamate decarboxylase that catalyzes the formation of GABA from glutamate. In conclusion, WRKY6 and 40 seem to act as important players transducing E-2-hexenal perception.

Transcriptional profiling of Medicago truncatula during Erysiphe pisi infection

Resistance to powdery mildew has been studied in a number of plant species, yet the molecular mechanisms remain largely unknown. Transcription factors (TFs) play a critical role in the plant defense response by regulating the transcriptional machinery which coordinates the expression of a large group of genes involved in plant defense. Using high-throughput quantitative real-time PCR (qPCR) technology more than 1000 Medicago truncatula TFs were screened in a pair of susceptible and resistant genotypes of M. truncatula after 4 h of Erysiphe pisi infection. Seventy nine TF genes, belonging to 33 families showed a significant transcriptional change in response to E. pisi infection. Forty eight TF genes were differentially expressed in the resistant genotypes compared to the susceptible one in response to E. pisi infection, including pathogenesis-related transcriptional factors, AP2/EREBP (APETALA2/ETHYLENE-RESPONSIVE ELEMENT BINDING FACTORS), WRKY (highly conserved WRKYGQK amino-acid sequence), MYB (Myeloblastoma), homeodomain (HD) and zinc finger C2C2 (CYS2-CYS2), C2H2, (CYS2-HIS2), LIM (Lin-11, Isl-1, Mec-3) gene families, which are involved in known defense responses. Our results suggest that these TF genes are among the E. pisi responsive genes in resistant M. truncatula that may constitute a regulatory network which controls the transcriptional changes in defense genes involved in resistance to E. pisi.

The Mutualist Laccaria bicolor Expresses a Core Gene Regulon During the Colonization of Diverse Host Plants and a Variable Regulon to Counteract Host-Specific Defenses

The coordinated transcriptomic responses of both mutualistic ectomycorrhizal (ECM) fungi and their hosts during the establishment of symbiosis are not well-understood. This study characterizes the transcriptomic alterations of the ECM fungus Laccaria bicolor during different colonization stages on two hosts (Populus trichocarpa and Pseudotsuga menziesii) and compares this to the transcriptomic variations of P. trichocarpa across the same time-points. A large number of L. bicolor genes (≥ 8,000) were significantly regulated at the transcriptional level in at least one stage of colonization. From our data, we identify 1,249 genes that we hypothesize is the 'core' gene regulon necessary for the mutualistic interaction between L. bicolor and its host plants. We further identify a group of 1,210 genes that are regulated in a host-specific manner. This variable regulon encodes a number of genes coding for proteases and xenobiotic efflux transporters that we hypothesize act to counter chemical-based defenses simultaneously activated at the transcriptomic level in P. trichocarpa. The transcriptional response of the host plant P. trichocarpa consisted of differential waves of gene regulation related to signaling perception and transduction, defense response, and the induction of nutrient transfer in P. trichocarpa tissues. This study, therefore, gives fresh insight into the shifting transcriptomic landscape in both the colonizing fungus and its host and the different strategies employed by both partners in orchestrating a mutualistic interaction.

The interaction of Arabidopsis with Piriformospora indica shifts from initial transient stress induced by fungus-released chemical mediators to a mutualistic interaction after physical contact of the two symbionts

BACKGROUND

Piriformospora indica, an endophytic fungus of Sebacinales, colonizes the roots of many plant species including Arabidopsis thaliana. The symbiotic interaction promotes plant performance, growth and resistance/tolerance against abiotic and biotic stress.

RESULTS

We demonstrate that exudated compounds from the fungus activate stress and defense responses in the Arabidopsis roots and shoots before the two partners are in physical contact. They induce stomata closure, stimulate reactive oxygen species (ROS) production, stress-related phytohormone accumulation and activate defense and stress genes in the roots and/or shoots. Once a physical contact is established, the stomata re-open, ROS and phytohormone levels decline, and the number and expression level of defense/stress-related genes decreases.

CONCLUSIONS

We propose that exudated compounds from P. indica induce stress and defense responses in the host. Root colonization results in the down-regulation of defense responses and the activation of genes involved in promoting plant growth, metabolism and performance.

Interplays between nitric oxide and reactive oxygen species in cryptogein signalling

Nitric oxide (NO) has many functions in plants. Here, we investigated its interplays with reactive oxygen species (ROS) in the defence responses triggered by the elicitin cryptogein. The production of NO induced by cryptogein in tobacco cells was partly regulated through a ROS-dependent pathway involving the NADPH oxidase NtRBOHD. In turn, NO down-regulated the level of H2O2. Both NO and ROS synthesis appeared to be under the control of type-2 histone deacetylases acting as negative regulators of cell death. Occurrence of an interplay between NO and ROS was further supported by the finding that cryptogein triggered a production of peroxynitrite (ONOO(-)). Next, we showed that ROS, but not NO, negatively regulate the intensity of activity of the cryptogein-induced protein kinase NtOSAK. Furthermore, using a DNA microarray approach, we identified 15 genes early induced by cryptogein via NO. A part of these genes was also modulated by ROS and encoded proteins showing sequence identity to ubiquitin ligases. Their expression appeared to be negatively regulated by ONOO(-), suggesting that ONOO(-) mitigates the effects of NO and ROS. Finally, we provided evidence that NO required NtRBOHD activity for inducing cell death, thus confirming previous assumption that ROS channel NO through cell death pathways.

Loss of Arabidopsis thaliana Dynamin-Related Protein 2B reveals separation of innate immune signaling pathways

Vesicular trafficking has emerged as an important means by which eukaryotes modulate responses to microbial pathogens, likely by contributing to the correct localization and levels of host components necessary for effective immunity. However, considering the complexity of membrane trafficking in plants, relatively few vesicular trafficking components with functions in plant immunity are known. Here we demonstrate that Arabidopsis thaliana Dynamin-Related Protein 2B (DRP2B), which has been previously implicated in constitutive clathrin-mediated endocytosis (CME), functions in responses to flg22 (the active peptide derivative of bacterial flagellin) and immunity against flagellated bacteria Pseudomonas syringae pv. tomato (Pto) DC3000. Consistent with a role of DRP2B in Pattern-Triggered Immunity (PTI), drp2b null mutant plants also showed increased susceptibility to Pto DC3000 hrcC-, which lacks a functional Type 3 Secretion System, thus is unable to deliver effectors into host cells to suppress PTI. Importantly, analysis of drp2b mutant plants revealed three distinct branches of the flg22-signaling network that differed in their requirement for RESPIRATORY BURST OXIDASE HOMOLOGUE D (RBOHD), the NADPH oxidase responsible for flg22-induced apoplastic reactive oxygen species production. Furthermore, in drp2b, normal MAPK signaling and increased immune responses via the RbohD/Ca2+-branch were not sufficient for promoting robust PR1 mRNA expression nor immunity against Pto DC3000 and Pto DC3000 hrcC-. Based on live-cell imaging studies, flg22-elicited internalization of the plant flagellin-receptor, FLAGELLIN SENSING 2 (FLS2), was found to be partially dependent on DRP2B, but not the closely related protein DRP2A, thus providing genetic evidence for a component, implicated in CME, in ligand-induced endocytosis of FLS2. Reduced trafficking of FLS2 in response to flg22 may contribute in part to the non-canonical combination of immune signaling defects observed in drp2b. In conclusion, this study adds DRP2B to the relatively short list of known vesicular trafficking proteins with roles in flg22-signaling and PTI in plants.

Genomic and transcriptomic characterization of the transcription factor family R2R3-MYB in soybean and its involvement in the resistance responses to Phakopsora pachyrhizi

Myb genes constitute one of the largest transcription factor families in the plant kingdom. Soybean MYB transcription factors have been related to the plant response to biotic stresses. Their involvement in response to Phakopsora pachyrhizi infection has been reported by several transcriptional studies. Due to their apparently highly diverse functions, these genes are promising targets for developing crop varieties resistant to diseases. In the present study, the identification and phylogenetic analysis of the soybean R2R3-MYB (GmMYB) transcription factor family was performed and the expression profiles of these genes under biotic stress were determined. GmMYBs were identified from the soybean genome using bioinformatic tools, and their putative functions were determined based on the phylogenetic tree and classified into subfamilies using guides AtMYBs describing known functions. The transcriptional profiles of GmMYBs upon infection with different pathogen were revealed by in vivo and in silico analyses. Selected target genes potentially involved in disease responses were assessed by RT-qPCR after different times of inoculation with P. pachyrhizi using different genetic backgrounds related to resistance genes (Rpp2 and Rpp5). R2R3-MYB transcription factors related to lignin synthesis and genes responsive to chitin were significantly induced in the resistant genotypes.

Transcriptome analysis during berry development provides insights into co-regulated and altered gene expression between a seeded wine grape variety and its seedless somatic variant

BACKGROUND

Seedless grapes are greatly appreciated for fresh and dry fruit consumption. Parthenocarpy and stenospermocarpy have been described as the main phenomena responsible for seedlessness in Vitis vinifera. However, the key genes underpinning molecular and cellular processes that play a significant role in seed development are not well characterized. To identify important regulators and mechanisms that may be altered in the seedless phenotype, we performed a comprehensive transcriptional analysis to compare the transcriptomes of a popular seeded wine cultivar (wild-type) and its seedless somatic variant (mutant) at three key developmental stages.

RESULTS

The transcriptomes revealed by Illumina mRNA-Seq technology had approximately 98% of grapevine annotated transcripts and about 80% of them were commonly expressed in the two lines. Differential gene expression analysis revealed a total of 1075 differentially expressed genes (DE) in the pairwise comparison of developmental stages, which included DE genes specific to the wild-type background, DE genes specific to the mutant background and DE genes commonly shared in both backgrounds. The analysis of differential expression patterns and functional category enrichment of wild-type and mutant DE genes highlighted significant coordination and enrichment of pollen and ovule developmental pathways. The expression of some selected DE genes was further confirmed by real-time RT-PCR analysis.

CONCLUSIONS

This study represents the most comprehensive attempt to characterize the genetic bases of seed formation in grapevine. With a high throughput method, we have shown that a seeded wine grape and its seedless somatic variant are similar in several biological processes. Nevertheless, we could identify an inventory of genes with altered expression in the mutant compared to the wild-type, which may be responsible for the seedless phenotype. The genes located within known genomic regions regulating seed content may be used for the development of molecular tools to assist table grape breeding. Therefore the data reported here have provided a rich genomic resource for practical use and functional characterization of the genes that potentially underpin seedlessness in grapevine.

Molecular effects of resistance elicitors from biological origin and their potential for crop protection

Plants contain a sophisticated innate immune network to prevent pathogenic microbes from gaining access to nutrients and from colonizing internal structures. The first layer of inducible response is governed by the plant following the perception of microbe- or modified plant-derived molecules. As the perception of these molecules results in a plant response that can provide efficient resistance toward non-adapted pathogens they can also be described as "defense elicitors." In compatible plant/microbe interactions, adapted microorganisms have means to avoid or disable this resistance response and promote virulence. However, this requires a detailed spatial and temporal response from the invading pathogens. In agricultural practice, treating plants with isolated defense elicitors in the absence of pathogens can promote plant resistance by uncoupling defense activation from the effects of pathogen virulence determinants. The plant responses to plant, bacterial, oomycete, or fungal-derived elicitors are not, in all cases, universal and need elucidating prior to the application in agriculture. This review provides an overview of currently known elicitors of biological rather than synthetic origin and places their activity into a molecular context.

Altered expression of Arabidopsis genes in response to a multifunctional geminivirus pathogenicity protein

BACKGROUND

Geminivirus AC2 is a multifunctional protein that acts as a pathogenicity factor. Transcriptional regulation by AC2 appears to be mediated through interaction with a plant specific DNA binding protein, PEAPOD2 (PPD2), that specifically binds to sequences known to mediate activation of the CP promoter of Cabbage leaf curl virus (CaLCuV) and Tomato golden mosaic virus (TGMV). Suppression of both basal and innate immune responses by AC2 in plants is mediated through inactivation of SnRK1.2, an Arabidopsis SNF1 related protein kinase, and adenosine kinase (ADK). An indirect promoter targeting strategy, via AC2-host dsDNA binding protein interactions, and inactivation of SnRK1.2-mediated defense responses could provide the opportunity for geminiviruses to alter host gene expression and in turn, reprogram the host to support virus infection. The goal of this study was to identify changes in the transcriptome of Arabidopsis induced by the transcription activation function of AC2 and the inactivation of SnRK1.2.

RESULTS

Using full-length and truncated AC2 proteins, microarray analyses identified 834 genes differentially expressed in response to the transcriptional regulatory function of the AC2 protein at one and two days post treatment. We also identified 499 genes differentially expressed in response to inactivation of SnRK1.2 by the AC2 protein at one and two days post treatment. Network analysis of these two sets of differentially regulated genes identified several networks consisting of between four and eight highly connected genes. Quantitative real-time PCR analysis validated the microarray expression results for 10 out of 11 genes tested.

CONCLUSIONS

It is becoming increasingly apparent that geminiviruses manipulate the host in several ways to facilitate an environment conducive to infection, predominantly through the use of multifunctional proteins. Our approach of identifying networks of highly connected genes that are potentially co-regulated by geminiviruses during infection will allow us to identify novel pathways of co-regulated genes that are stimulated in response to pathogen infection in general, and virus infection in particular.

The ubiquitination machinery of the ubiquitin system

The protein ubiquitin is a covalent modifier of proteins, including itself. The ubiquitin system encompasses the enzymes required for catalysing attachment of ubiquitin to substrates as well as proteins that bind to ubiquitinated proteins leading them to their final fate. Also included are activities that remove ubiquitin independent of, or in concert with, proteolysis of the substrate, either by the proteasome or proteases in the vacuole. In addition to ubiquitin encoded by a family of fusion proteins, there are proteins with ubiquitin-like domains, likely forming ubiquitin's β-grasp fold, but incapable of covalent modification. However, they serve as protein-protein interaction platforms within the ubiquitin system. Multi-gene families encode all of these types of activities. Within the ubiquitination machinery "half" of the ubiquitin system are redundant, partially redundant, and unique components affecting diverse developmental and environmental responses in plants. Notably, multiple aspects of biotic and abiotic stress responses require, or are modulated by, ubiquitination. Finally, aspects of the ubiquitin system have broad utility: as components to enhance gene expression or to regulate protein abundance. This review focuses on the ubiquitination machinery: ubiquitin, unique aspects about the synthesis of ubiquitin and organization of its gene family, ubiquitin activating enzymes (E1), ubiquitin conjugating enzymes (E2) and ubiquitin ligases, or E3s. Given the large number of E3s in Arabidopsis this review covers the U box, HECT and RING type E3s, with the exception of the cullin-based E3s.

Transcriptome analysis of Gossypium hirsutum flower buds infested by cotton boll weevil (Anthonomus grandis) larvae

Background

Cotton is a major fibre crop grown worldwide that suffers extensive damage from chewing insects, including the cotton boll weevil larvae (Anthonomus grandis). Transcriptome analysis was performed to understand the molecular interactions between Gossypium hirsutum L. and cotton boll weevil larvae. The Illumina HiSeq 2000 platform was used to sequence the transcriptome of cotton flower buds infested with boll weevil larvae.

Results

The analysis generated a total of 327,489,418 sequence reads that were aligned to the G. hirsutum reference transcriptome. The total number of expressed genes was over 21,697 per sample with an average length of 1,063 bp. The DEGseq analysis identified 443 differentially expressed genes (DEG) in cotton flower buds infected with boll weevil larvae. Among them, 402 (90.7%) were up-regulated, 41 (9.3%) were down-regulated and 432 (97.5%) were identified as orthologues of A. thaliana genes using Blastx. Mapman analysis of DEG indicated that many genes were involved in the biotic stress response spanning a range of functions, from a gene encoding a receptor-like kinase to genes involved in triggering defensive responses such as MAPK, transcription factors (WRKY and ERF) and signalling by ethylene (ET) and jasmonic acid (JA) hormones. Furthermore, the spatial expression pattern of 32 of the genes responsive to boll weevil larvae feeding was determined by “in situ” qPCR analysis from RNA isolated from two flower structures, the stamen and the carpel, by laser microdissection (LMD).

Conclusion

A large number of cotton transcripts were significantly altered upon infestation by larvae. Among the changes in gene expression, we highlighted the transcription of receptors/sensors that recognise chitin or insect oral secretions; the altered regulation of transcripts encoding enzymes related to kinase cascades, transcription factors, Ca²⁺ influxes, and reactive oxygen species; and the modulation of transcripts encoding enzymes from phytohormone signalling pathways. These data will aid in the selection of target genes to genetically engineer cotton to control the cotton boll weevil.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-854) contains supplementary material, which is available to authorized users.

OsNAC111, a blast disease-responsive transcription factor in rice, positively regulates the expression of defense-related genes

Plants respond to pathogen attack by transcriptionally regulating defense-related genes via various types of transcription factors. We identified a transcription factor in rice, OsNAC111, belonging to the TERN subgroup of the NAC family that was transcriptionally upregulated after rice blast fungus (Magnaporthe oryzae) inoculation. OsNAC111 was localized in the nucleus of rice cells and had transcriptional activation activity in yeast and rice cells. Transgenic rice plants overexpressing OsNAC111 showed increased resistance to the rice blast fungus. In OsNAC111-overexpressing plants, the expression of several defense-related genes, including pathogenesis-related (PR) genes, was constitutively high compared with the control. These genes all showed blast disease-responsive expression in leaves. Among them, two chitinase genes and one β-1,3-glucanase gene showed reduced expression in transgenic rice plants in which OsNAC111 function was suppressed by a chimeric repressor (OsNAC111-SRDX). OsNAC111 activated transcription from the promoters of the chitinase and β-1,3-glucanase genes in rice cells. In addition, brown pigmentation at the infection sites, a defense response of rice cells to the blast fungus, was lowered in OsNAC111-SRDX plants at the early infection stage. These results indicate that OsNAC111 positively regulates the expression of a specific set of PR genes in the disease response and contributes to disease resistance.

Roles for jasmonate- and ethylene-induced transcription factors in the ability of Arabidopsis to respond differentially to damage caused by two insect herbivores

Plant responses to insects and wounding involve substantial transcriptional reprogramming that integrates hormonal, metabolic, and physiological events. The ability to respond differentially to various stresses, including wounding, generally involves hormone signaling and trans-acting regulatory factors. Evidence of the importance of transcription factors (TFs) in responses to insects is also accumulating. However, the relationships among hormone signaling, TF activity, and ability to respond specifically to different insects are uncertain. We examined transcriptional and hormonal changes in Arabidopsis thaliana after herbivory by larvae of two lepidopteran species, Spodoptera exigua (Hübner) and Pieris rapae L. over a 24-h time course. Transcriptional responses to the two insects differed and were frequently weaker or absent in response to the specialist P. rapae. Using microarray analysis and qRT-PCR, we found 141 TFs, including many AP2/ERFs (Ethylene Response Factors) and selected defense-related genes, to be differentially regulated in response to the two insect species or wounding. Jasmonic Acid (JA), JA-isoleucine (JA-IL), and ethylene production by Arabidopsis plants increased after attack by both insect species. However, the amounts and timing of ethylene production differed between the two herbivory treatments. Our results support the hypothesis that the different responses to these two insects involve modifications of JA-signaling events and activation of different subsets of ERF TFs, resulting in different degrees of divergence from responses to wounding alone.

Photoreceptor partner FHY1 has an independent role in gene modulation and plant development under far-red light

To incorporate the far-red light (FR) signal into a strategy for optimizing plant growth, FAR-RED ELONGATED HYPOCOTYL1 (FHY1) mediates the nuclear translocation of the FR photoreceptor phytochrome A (phyA) and facilitates the association of phyA with the promoters of numerous associated genes crucial for the response to environmental stimuli. However, whether FHY1 plays additional roles after FR irradiation remains elusive. Here, through the global identification of FHY1 chromatin association sites through ChIP-seq analysis and by the comparison of FHY1-associated sites with phyA-associated sites, we demonstrated that nuclear FHY1 can either act independently of phyA or act in association with phyA to activate the expression of distinct target genes. We also determined that phyA can act independently of FHY1 in regulating phyA-specific target genes. Furthermore, we determined that the independent FHY1 nuclear pathway is involved in crucial aspects of plant development, as in the case of inhibited seed germination under FR during salt stress. Notably, the differential presence of cis-elements and transcription factors in common and unique FHY1- and/or phyA-associated genes are indicative of the complexity of the independent and coordinated FHY1 and phyA pathways. Our study uncovers previously unidentified aspects of FHY1 function beyond its currently recognized role in phyA-dependent photomorphogenesis.

LIK1, a CERK1-interacting kinase, regulates plant immune responses in Arabidopsis

Chitin, an integral component of the fungal cell wall, is one of the best-studied microbe-associated molecular patterns. Previous work identified a LysM receptor-like kinase (LysM-RLK1/CERK1) as the primary chitin receptor in Arabidopsis. In order to identify proteins that interact with CERK1, we conducted a yeast two-hybrid screen using the intracellular kinase domain of CERK1 as the bait. This screen identified 54 putative CERK1-interactors. Screening mutants defective in 43 of these interacting proteins identified only two, a calmodulin like protein (At3g10190) and a leucine-rich repeat receptor like kinase (At3g14840), which differed in their response to pathogen challenge. In the present work, we focused on characterizing the LRR-RLK gene where mutations altered responses to chitin elicitation. This LRR-RLK was named LysM RLK1-interacting kinase 1 (LIK1). The interaction between CERK1 and LIK1 was confirmed by co-immunoprecipitation using protoplasts and transgenic plants. In vitro experiments showed that LIK1 was directly phosphorylated by CERK1. In vivo phosphorylation assays showed that Col-0 wild-type plants have more phosphorylated LIK1 than cerk1 mutant plants, suggesting that LIK1 may be directly phosphorylated by CERK1. Lik1 mutant plants showed an enhanced response to both chitin and flagellin elicitors. In comparison to the wild-type plants, lik1 mutant plants were more resistant to the hemibiotrophic pathogen Pseudomonas syringae, but more susceptible to the necrotrophic pathogen Sclerotinia sclerotiorum. Consistent with the enhanced susceptibility to necrotrophs, lik1 mutants showed reduced expression of genes involved in jasmonic acid and ethylene signaling pathways. These data suggest that LIK1 directly interacts with CERK1 and regulates MAMP-triggered innate immunity.

The BBX family of plant transcription factors

The B-box (BBX) proteins are a class of zinc-finger transcription factors containing a B-box domain with one or two B-box motifs, and sometimes also feature a CCT (CONSTANS, CO-like, and TOC1) domain. BBX proteins are key factors in regulatory networks controlling growth and developmental processes that include seedling photomorphogenesis, photoperiodic regulation of flowering, shade avoidance, and responses to biotic and abiotic stresses. In this review we discuss the functions of BBX proteins and the role of B-box motif in mediating transcriptional regulation and protein-protein interaction in plant signaling. In addition, we provide novel insights into the molecular mechanisms of their action and the evolutionary significance of their functional divergence.

Synchronization of developmental processes and defense signaling by growth regulating transcription factors

Growth regulating factors (GRFs) are a conserved class of transcription factor in seed plants. GRFs are involved in various aspects of tissue differentiation and organ development. The implication of GRFs in biotic stress response has also been recently reported, suggesting a role of these transcription factors in coordinating the interaction between developmental processes and defense dynamics. However, the molecular mechanisms by which GRFs mediate the overlaps between defense signaling and developmental pathways are elusive. Here, we report large scale identification of putative target candidates of Arabidopsis GRF1 and GRF3 by comparing mRNA profiles of the grf1/grf2/grf3 triple mutant and those of the transgenic plants overexpressing miR396-resistant version of GRF1 or GRF3. We identified 1,098 and 600 genes as putative targets of GRF1 and GRF3, respectively. Functional classification of the potential target candidates revealed that GRF1 and GRF3 contribute to the regulation of various biological processes associated with defense response and disease resistance. GRF1 and GRF3 participate specifically in the regulation of defense-related transcription factors, cell-wall modifications, cytokinin biosynthesis and signaling, and secondary metabolites accumulation. GRF1 and GRF3 seem to fine-tune the crosstalk between miRNA signaling networks by regulating the expression of several miRNA target genes. In addition, our data suggest that GRF1 and GRF3 may function as negative regulators of gene expression through their association with other transcription factors. Collectively, our data provide new insights into how GRF1 and GRF3 might coordinate the interactions between defense signaling and plant growth and developmental pathways.

Validation of potential reference genes for qPCR in maize across abiotic stresses, hormone treatments, and tissue types

The reverse transcription quantitative polymerase chain reaction (RT-qPCR) is a powerful and widely used technique for the measurement of gene expression. Reference genes, which serve as endogenous controls ensure that the results are accurate and reproducible, are vital for data normalization. To bolster the literature on reference gene selection in maize, ten candidate reference genes, including eight traditionally used internal control genes and two potential candidate genes from our microarray datasets, were evaluated for expression level in maize across abiotic stresses (cold, heat, salinity, and PEG), phytohormone treatments (abscisic acid, salicylic acid, jasmonic acid, ethylene, and gibberellins), and different tissue types. Three analytical software packages, geNorm, NormFinder, and Bestkeeper, were used to assess the stability of reference gene expression. The results revealed that elongation factor 1 alpha (EF1α), tubulin beta (β-TUB), cyclophilin (CYP), and eukaryotic initiation factor 4A (EIF4A) were the most reliable reference genes for overall gene expression normalization in maize, while GRP (Glycine-rich RNA-binding protein), GLU1(beta-glucosidase), and UBQ9 (ubiquitin 9) were the least stable and most unsuitable genes. In addition, the suitability of EF1α, β-TUB, and their combination as reference genes was confirmed by validating the expression of WRKY50 in various samples. The current study indicates the appropriate reference genes for the urgent requirement of gene expression normalization in maize across certain abiotic stresses, hormones, and tissue types.

Accumulation of N-acetylglucosamine oligomers in the plant cell wall affects plant architecture in a dose-dependent and conditional manner

To study the effect of short N-acetylglucosamine (GlcNAc) oligosaccharides on the physiology of plants, N-ACETYLGLUCOSAMINYLTRANSFERASE (NodC) of Azorhizobium caulinodans was expressed in Arabidopsis (Arabidopsis thaliana). The corresponding enzyme catalyzes the polymerization of GlcNAc and, accordingly, β-1,4-GlcNAc oligomers accumulated in the plant. A phenotype characterized by difficulties in developing an inflorescence stem was visible when plants were grown for several weeks under short-day conditions before transfer to long-day conditions. In addition, a positive correlation between the oligomer concentration and the penetrance of the phenotype was demonstrated. Although NodC overexpression lines produced less cell wall compared with wild-type plants under nonpermissive conditions, no indications were found for changes in the amount of the major cell wall polymers. The effect on the cell wall was reflected at the transcriptome level. In addition to genes encoding cell wall-modifying enzymes, a whole set of genes encoding membrane-coupled receptor-like kinases were differentially expressed upon GlcNAc accumulation, many of which encoded proteins with an extracellular Domain of Unknown Function26. Although stress-related genes were also differentially expressed, the observed response differed from that of a classical chitin response. This is in line with the fact that the produced chitin oligomers were too small to activate the chitin receptor-mediated signal cascade. Based on our observations, we propose a model in which the oligosaccharides modify the architecture of the cell wall by acting as competitors in carbohydrate-carbohydrate or carbohydrate-protein interactions, thereby affecting noncovalent interactions in the cell wall or at the interface between the cell wall and the plasma membrane.

The heat shock factor A4A confers salt tolerance and is regulated by oxidative stress and the mitogen-activated protein kinases MPK3 and MPK6

Heat shock factors (HSFs) are principal regulators of plant responses to several abiotic stresses. Here, we show that estradiol-dependent induction of HSFA4A confers enhanced tolerance to salt and oxidative agents, whereas inactivation of HSFA4A results in hypersensitivity to salt stress in Arabidopsis (Arabidopsis thaliana). Estradiol induction of HSFA4A in transgenic plants decreases, while the knockout hsfa4a mutation elevates hydrogen peroxide accumulation and lipid peroxidation. Overexpression of HSFA4A alters the transcription of a large set of genes regulated by oxidative stress. In yeast (Saccharomyces cerevisiae) two-hybrid and bimolecular fluorescence complementation assays, HSFA4A shows homomeric interaction, which is reduced by alanine replacement of three conserved cysteine residues. HSFA4A interacts with mitogen-activated protein kinases MPK3 and MPK6 in yeast and plant cells. MPK3 and MPK6 phosphorylate HSFA4A in vitro on three distinct sites, serine-309 being the major phosphorylation site. Activation of the MPK3 and MPK6 mitogen-activated protein kinase pathway led to the transcriptional activation of the HEAT SHOCK PROTEIN17.6A gene. In agreement that mutation of serine-309 to alanine strongly diminished phosphorylation of HSFA4A, it also strongly reduced the transcriptional activation of HEAT SHOCK PROTEIN17.6A. These data suggest that HSFA4A is a substrate of the MPK3/MPK6 signaling and that it regulates stress responses in Arabidopsis.

Trichoderma asperelloides suppresses nitric oxide generation elicited by Fusarium oxysporum in Arabidopsis roots

Inoculations with saprophytic fungus Trichoderma spp. are now extensively used both to promote plant growth and to suppress disease development. The underlying mechanisms for both roles have yet to be fully described so that the use of Trichoderma spp. could be optimized. Here, we show that Trichoderma asperelloides effects include the manipulation of host nitric oxide (NO) production. NO was rapidly formed in Arabidopsis roots in response to the soil-borne necrotrophic pathogen Fusarium oxysporum and persisted for about 1 h but is only transiently produced (approximately 10 min) when roots interact with T. asperelloides (T203). However, inoculation of F. oxysporum-infected roots with T. asperelloides suppressed F. oxysporum-initiated NO production. A transcriptional study of 78 NO-modulated genes indicated most genes were suppressed by single and combinational challenge with F. oxysporum or T. asperelloides. Only two F. oxysporum-induced genes were suppressed by T. asperelloides inoculation undertaken either 10 min prior to or after pathogen infection: a concanavlin A-like lectin protein kinase (At4g28350) and the receptor-like protein RLP30. Thus, T. asperelloides can actively suppress NO production elicited by F. oxysporum and impacts on the expression of some genes reported to be NO-responsive. Of particular interest was the reduced expression of receptor-like genes that may be required for F. oxysporum-dependent necrotrophic disease development.