Functional studies of transcription factors involved in plant defenses in the genomics era

Interactive transcriptome analyses of Northern Wild Rice (Zizania palustris L.) and Bipolaris oryzae show convoluted communications during the early stages of fungal brown spot development

Fungal diseases, caused mainly by Bipolaris spp., are past and current threats to Northern Wild Rice (NWR) grain production and germplasm preservation in both natural and cultivated settings. Genetic resistance against the pathogen is scarce. Toward expanding our understanding of the global gene communications of NWR and Bipolaris oryzae interaction, we designed an RNA sequencing study encompassing the first 12 h and 48 h of their encounter. NWR activated numerous plant recognition receptors after pathogen infection, followed by active transcriptional reprogramming of signaling mechanisms driven by Ca2+ and its sensors, mitogen-activated protein kinase cascades, activation of an oxidative burst, and phytohormone signaling-bound mechanisms. Several transcription factors associated with plant defense were found to be expressed. Importantly, evidence of diterpenoid phytoalexins, especially phytocassane biosynthesis, among expression of other defense genes was found. In B. oryzae, predicted genes associated with pathogenicity including secreted effectors that could target plant defense mechanisms were expressed. This study uncovered the early molecular communication between the NWR-B. oryzae pathosystem, which could guide selection for allele-specific genes to boost NWR defenses, and overall aid in the development of more efficient selection methods in NWR breeding through the use of the most virulent fungal isolates.

Functional characterization of plant specific Indeterminate Domain (IDD) transcription factors in tomato (Solanum lycopersicum L.)

Plant-specific transcription factors (TFs) are responsible for regulating the genes involved in the development of plant-specific organs and response systems for adaptation to terrestrial environments. This includes the development of efficient water transport systems, efficient reproductive organs, and the ability to withstand the effects of terrestrial factors, such as UV radiation, temperature fluctuations, and soil-related stress factors, and evolutionary advantages over land predators. In rice and Arabidopsis, INDETERMINATE DOMAIN (IDD) TFs are plant-specific TFs with crucial functions, such as development, reproduction, and stress response. However, in tomatoes, IDD TFs remain uncharacterized. Here, we examined the presence, distribution, structure, characteristics, and expression patterns of SlIDDs. Database searches, multiple alignments, and motif alignments suggested that 24 TFs were related to Arabidopsis IDDs. 18 IDDs had two characteristic C2H2 domains and two C2HC domains in their coding regions. Expression analyses suggest that some IDDs exhibit multi-stress responsive properties and can respond to specific stress conditions, while others can respond to multiple stress conditions in shoots and roots, either in a tissue-specific or universal manner. Moreover, co-expression database analyses suggested potential interaction partners within IDD family and other proteins. This study functionally characterized SlIDDs, which can be studied using molecular and bioinformatics methods for crop improvement.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

Transcriptome analysis reveals key genes and pathways for prickle development in Zanthoxylumarmatum

Zanthoxylum armatum is an economically important tree species. However, well-developed prickles on its stems and leaves pose serious challenges in terms of management and harvesting. To investigate the molecular mechanism underlying prickle development, we sequenced different stages of prickle morphological development and transcriptomes of different tissues in the root tips (Gen), leaf buds (Ya), and fruits of Z. armatum. The results revealed that proteins related to cell division and genes related to the growth hormone signaling pathway were highly expressed in the prickle just protrusion (PC1). In addition, a high expression of lignin biosynthesis genes was observed during the developmental onset of lignification (PC2) and prickle lignification (PC3). These findings indicate that phenylpropanoid biosynthesis and plant hormone signal transduction are key pathways for the completion of lignification development in the prickle. During prickle development, ZaMYB2 and ZaWRKY3 were significantly upregulated in PC2 and PC3, suggesting their possible involvement in prickle development. Transcriptome and qRT-PCR analyses revealed differential gene expression of zaPAL3, za4CLL1, zaCOMT1, ZaWRKY3, and ZaCCD31 in the Gen, Ya, newly formed fruit (ZaF1), newly oil-spotted fruits (ZaF2), PC1, PC2, and PC3 of Zarmatum. zaCCD31 was highly expressed in leaf buds, whereas Za4CLL1 was highly expressed in root tips. During the lignification of prickles, the relative expression of genes including zaMYB2 increased gradually; however, the relative expression of zaCCD31 decreased during this process. Therefore, we inferred that these genes might be closely related to prickle development. Notably, zaMYB2 was expressed at higher levels in PC2 and PC3 than in PC1 and was not expressed in Gen, Ya, ZaF1, and ZaF2. Therefore, zaMYB2 is a key gene involved in prickle development of Z. armatum that exhibited tissue-specific expression. This study establishes a foundation for future analyses of the molecular mechanism underlying prickle development in Z. armatum.

Bacillus velezensis WB induces systemic resistance in watermelon against Fusarium wilt

BACKGROUND

Our previous findings indicated that Bacillus velezensis WB could control Fusarium wilt by changing the structure of the microbial community in the watermelon rhizosphere. However, there are few studies on its mechanism in the pathogen resistance of watermelon. Therefore, in this study, we determined the mechanism of B. velezensis WB-induced systemic resistance in watermelon against Fusarium wilt through glasshouse pot experiments.

RESULTS

The results showed that B. velezensis WB significantly reduced the incidence and disease index of Fusarium wilt in watermelon. B. velezensis WB can enhance the basal immunity of watermelon plants by: increasing the activity of phenylalanine ammonia-lyase (PAL), peroxidase (POD), superoxide dismutase (SOD) and β-1,3-glucanase; accumulating lignin, salicylic acid (SA) and jasmonic acid (JA); reducing malondialdehyde (MDA) concentrations; and inducing callus deposition in watermelon plant cells. RNA-seq analysis showed that 846 watermelon genes were upregulated and 612 watermelon genes were downregulated in the WF treatment. This process led to the activation of watermelon genes associated with auxin, gibberellin, SA, ethylene and JA, and the expression of genes in the phenylalanine biosynthetic pathway was upregulated. In addition, transcription factors involved in plant resistance to pathogens, such as MYB, NAC and WRKY, were induced. Gene correlation analysis showed that Cla97C10G195840 and Cla97C02G049930 in the phenylalanine biosynthetic pathway, and Cla97C02G041360 and Cla97C10G197290 in the plant hormone signal transduction pathway showed strong correlations with other genes.

CONCLUSION

Our results indicated that B. velezensis WB is capable of inducing systemic resistance in watermelon against Fusarium wilt. © 2023 Society of Chemical Industry.

In Silico Comparison of WRKY Transcription Factors in Wild and Cultivated Soybean and Their Co-expression Network Arbitrating Disease Resistance

WRKY Transcription factors (TFs) play critical roles in plant defence mechanisms that are activated in response to biotic and abiotic stresses. However, information on the Glycine soja WRKYs (GsoWRKYs) is scarce. Owing to its importance in soybean breeding, here we identified putative WRKY TFs in wild soybean, and compared the results with Glycine max WRKYs (GmaWRKYs) by phylogenetic, conserved motif, and duplication analyses. Moreover, we explored the expression trends of WRKYs in G. max (oomycete, fungi, virus, bacteria, and soybean cyst nematode) and G. soja (soybean cyst nematode), and identified commonly expressed WRKYs and their co-expressed genes. We identified, 181 and 180 putative WRKYs in G. max and G. soja, respectively. Though the number of WRKYs in both studied species is almost the same, they differ in many ways, i.e., the number of WRKYs on corresponding chromosomes, conserved domain structures, WRKYGQK motif variants, and zinc-finger motifs. WRKYs in both species grouped in three major clads, i.e., I-III, where group-II had sub-clads IIa-IIe. We found that GsoWRKYs expanded mostly through segmental duplication. A large number of WRKYs were expressed in response to biotic stresses, i.e., Phakospora pachyrhizi, Phytoplasma, Heterodera glycines, Macrophomina phaseolina, and Soybean mosaic virus; 56 GmaWRKYs were commonly expressed in soybean plants infected with these diseases. Finally, 30 and 63 GmaWRKYs and GsoWRKYs co-expressed with 205 and 123 non-WRKY genes, respectively, indicating that WRKYs play essential roles in biotic stress tolerance in Glycine species.

Comparative transcriptome profiling reveals differential defense responses among Alternaria brassicicola resistant Sinapis alba and susceptible Brassica rapa

Alternaria blight is a devastating disease that causes significant crop losses in oilseed Brassicas every year. Adoption of conventional breeding to generate disease-resistant varieties has so far been unsuccessful due to the lack of suitable resistant source germplasms of cultivated Brassica spp. A thorough understanding of the molecular basis of resistance, as well as the identification of defense-related genes involved in resistance responses in closely related wild germplasms, would substantially aid in disease management. In the current study, a comparative transcriptome profiling was performed using Illumina based RNA-seq to detect differentially expressed genes (DEGs) specifically modulated in response to Alternaria brassicicola infection in resistant Sinapis alba, a close relative of Brassicas, and the highly susceptible Brassica rapa. The analysis revealed that, at 48 hpi (hours post inoculation), 3396 genes were upregulated and 23239 were downregulated, whereas at 72 hpi, 4023 genes were upregulated and 21116 were downregulated. Furthermore, a large number of defense response genes were detected to be specifically regulated as a result of Alternaria infection. The transcriptome data was validated using qPCR-based expression profiling for selected defense-related DEGs, that revealed significantly higher fold change in gene expression in S. alba when compared to B. rapa. Expression of most of the selected genes was elevated across all the time points under study with significantly higher expression towards the later time point of 72 hpi in the resistant germplasm. S. alba activates a stronger defense response reaction against the disease by deploying an array of genes and transcription factors involved in a wide range of biological processes such as pathogen recognition, signal transduction, cell wall modification, antioxidation, transcription regulation, etc. Overall, the study provides new insights on resistance of S. alba against A. brassicicola, which will aid in devising strategies for breeding resistant varieties of oilseed Brassica.

Multifaceted roles of transcription factors during plant embryogenesis

Transcription factors (TFs) are diverse groups of regulatory proteins. Through their specific binding domains, TFs bind to their target genes and regulate their expression, therefore TFs play important roles in various growth and developmental processes. Plant embryogenesis is a highly regulated and intricate process during which embryos arise from various sources and undergo development; it can be further divided into zygotic embryogenesis (ZE) and somatic embryogenesis (SE). TFs play a crucial role in the process of plant embryogenesis with a number of them acting as master regulators in both ZE and SE. In this review, we focus on the master TFs involved in embryogenesis such as BABY BOOM (BBM) from the APETALA2/Ethylene-Responsive Factor (AP2/ERF) family, WUSCHEL and WUSCHEL-related homeobox (WOX) from the homeobox family, LEAFY COTYLEDON 2 (LEC2) from the B3 family, AGAMOUS-Like 15 (AGL15) from the MADS family and LEAFY COTYLEDON 1 (LEC1) from the Nuclear Factor Y (NF-Y) family. We aim to present the recent progress pertaining to the diverse roles these master TFs play in both ZE and SE in Arabidopsis, as well as other plant species including crops. We also discuss future perspectives in this context.

Transcriptome-based analysis of candidate gene markers associated with resistance mechanism to Phytophthora melonis that causes root and crown rot in pumpkin

Root and crown rot incited by an oomycete, Phytophthora melonis , causes significant yield losses in commercial pumpkin (Cucurbita pepo ) production worldwide. Currently, resistant cultivars and knowledge of molecular mechanism of C. pepo against P. melonis are scarce. Here, we analysed the quantitative gene expression changes of 10 candidate gene markers (bHLH87, ERF014, HSF, MYB, PR-1, WRKY21, CPI, POD, PSK, SGT ) in pumpkin roots and leaves at three time points (h post-inoculation, hpi) following inoculation with P. melonis in two resistant (Ghelyani and Tanbal), and two susceptible (Marmari and Khoreshti) varieties of pumpkin. Gene expression using quantitative real time PCR along a time course revealed the strongest transcriptomic response at 48 and 72hpi in resistant genotypes, 1.1-2.7-fold in roots and leaves, respectively, with a high significant correlation (r =0.857**-0.974**). We also found that CPI , PSK, SGT1 and POD act as a dual regulator that similarly modulate immunity not only against P. melonis , but also against other diseases such as early blight (Alternaria cucumerina) , powdery mildew (Podosphaera xanthii ), downy mildews (Pseudoperonospora cubensis ), and pathogenic plant nematodes (Meloidogyne javanica ). Furthermore, significantly higher activities of the ROS scavenging defence enzymes, catalase (1.6-fold increase) and peroxidase (6-fold increase) were observed in the roots of resistant cultivars at different hpi compared with non-inoculated controls. In addition, the biomass growth parameters including leaf and root length, stem and root diameter, root fresh weight and volume were significantly different among studied genotypes. Cumulatively, the transcriptome data provide novel insights into the response of pumpkins for improving pumpkin breeding to P. melonis .

The dynamic arms race during the early invasion of woodland strawberry by Botrytis cinerea revealed by dual dense high-resolution RNA-seq analyses

Necrotrophic pathogens replicate massively upon colonizing plants, causing large-scale wilting and death of plant tissues. Understanding both mechanisms of pathogen invasion and host response processes prior to symptom appearance and their key regulatory networks is therefore important for defense against pathogen attack. Here, we investigated the mechanisms of interaction between woodland strawberry (Fragaria vesca) leaves and gray mold pathogen (Botrytis cinerea) at 14 infection time points during the first 12 hours of the infection period using a dense, high-resolution time series dual transcriptomic analysis, characterizing the arms race between strawberry F. vesca and B. cinerea before the appearance of localized lesions. Strawberry leaves rapidly initiated strong systemic defenses at the first sign of external stimulation and showed lower levels of transcriptomic change later in the infection process. Unlike the host plants, B. cinerea showed larger-scale transcriptomic changes that persisted throughout the infection process. Weighted gene co-expression network analysis identified highly correlated genes in 32 gene expression modules between B. cinerea and strawberry. Yeast two-hybrid and bimolecular fluorescence complementation assays revealed that the disease response protein FvRLP2 from woodland strawberry interacted with the cell death inducing proteins BcXYG1 and BcPG3 from B. cinerea. Overexpression of FvRLP2 in both strawberry and Arabidopsis inhibited B. cinerea infection, confirming these genes' respective functions. These findings shed light on the arms race process by which B. cinerea invades host plants and strawberry to defend against pathogen infection.

SlWRKY16 and SlWRKY31 of tomato, negative regulators of plant defense, involved in susceptibility activation following root-knot nematode Meloidogyne javanica infection

The involvement of WRKY transcription factors in plant-nematode interactions, and in particular, how these WRKYs participate in regulating the complex morphological and physiological changes occurring after nematode infection, are the topic of active research. We characterized the functional role of the unstudied tomato WRKY genes SlWRKY16 and SlWRKY31 in regulating tomato roots' response to infection by the root-knot nematode Meloidogyne javanica. Using promoter-GUS reporter gene fusions and qRT-PCR, we show that both SlWRKYs are predominantly expressed during the first half of the parasitic life stages, when feeding-site induction and construction occur. Expression of SlWRKY16 increased sharply 15 days after inoculation, whereas SlWRKY31 was already induced earlier, but reached its maximum expression at this time. Both genes were downregulated at the mature female stage. To determine biological function, we produced transgenic lines overexpressing SlWRKY16 and SlWRKY31 in tomato hairy roots. Overexpression of both genes resulted in enhanced M. javanica infection, reflected by increased galling occurrence and reproduction. Expression profiling of marker genes responsive to defense-associated phytohormones indicated reductions in salicylic acid defense-related PR-1 and jasmonic acid defense-related PI in inoculated roots overexpressing SlWRK16 and SlWRKY31, respectively. Our results suggest that SlWRKY16 and SlWRKY31 function as negative regulators of plant immunity induced upon nematode infection.

Comparative transcriptome profiling of potato cultivars infected by late blight pathogen Phytophthora infestans: Diversity of quantitative and qualitative responses

The Estonia potato cultivar Ando has shown elevated field resistance to Phytophthora infestans, even after being widely grown for over 40 years. A comprehensive transcriptional analysis was performed using RNA-seq from plant leaf tissues to gain insight into the mechanisms activated for the defense after infection. Pathogen infection in Ando resulted in about 5927 differentially expressed genes (DEGs) compared to 1161 DEGs in the susceptible cultivar Arielle. The expression levels of genes related to plant disease resistance such as serine/threonine kinase activity, signal transduction, plant-pathogen interaction, endocytosis, autophagy, mitogen-activated protein kinase (MAPK), and others were significantly enriched in the upregulated DEGs in Ando, whereas in the susceptible cultivar, only the pathway related to phenylpropanoid biosynthesis was enriched in the upregulated DEGs. However, in response to infection, photosynthesis was deregulated in Ando. Multi-signaling pathways of the salicylic-jasmonic-ethylene biosynthesis pathway were also activated in response to Phytophthora infestans infection.

Binding of stress-responsive OsWRKY proteins through WRKYGQK heptapeptide residue with the promoter region of two rice blast disease resistance genes Pi2 and Pi54 is important for development of blast resistance

Molecular docking was done to investigate the interactions between five differentially expressed rice WRKY proteins when challenged with the rice blast disease caused by Magnaporthe oryzae and drought stresses applied either individually or overlapped, with the promoter region of two blast resistance genes (Pi2 and Pi54). Molecular docking was performed using the HDOCK server. Initially, the homology models for each of the five rice WRKY proteins were prepared using I-TASSER server, and then the secondary structure as well as the DNA-binding pockets were predicted using PSIPRED and BindUP servers, respectively. The molecular docking study revealed a differential binding pattern of the rice WRKYs with the two blast resistance genes. The WRKY proteins (OsWRKY88 and OsWRKY102), whose transcript levels decrease when drought and blast stresses are overlapped, interact with the two resistance genes mostly involving the residues of the zinc finger structure. On the other hand, the WRKY proteins (OsWRKY53-1 and OsWRKY113), whose transcript levels did not reduce significantly when challenged by drought and blast overlapped condition compared to individual treatment of blast, interact mostly involving the residues of the conserved WRKYGQK heptapeptide sequence. Interestingly, the protein OsWRKY74 whose transcript levels are unaffected in both individual and overlapped stresses, interacts with both the blast resistance genes involving few residues of both WRKYGQK heptapeptide and the zinc finger structure. The findings thus indicate that the interaction of OsWRKY proteins involving the conserved WRKYGQK heptapeptide sequence with the blast resistance genes Pi2 and Pi54 is important to mitigate the blast challenge in rice even during overlapping challenges of drought.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03711-y.

RAV1 family members function as transcriptional regulators and play a positive role in plant disease resistance

Phytopathogens pose a severe threat to agriculture and strengthening the plant defense response is an important strategy for disease control. Here, we report that AtRAV1, an AP2 and B3 domain-containing transcription factor, is required for basal plant defense in Arabidopsis thaliana. The atrav1 mutant lines demonstrate hyper-susceptibility against fungal pathogens (Rhizoctonia solani and Botrytis cinerea), whereas AtRAV1 overexpressing lines exhibit disease resistance against them. Enhanced expression of various defense genes and activation of mitogen-activated protein kinases (AtMPK3 and AtMPK6) are observed in the R. solani infected overexpressing lines, but not in the atrav1 mutant plants. An in vitro phosphorylation assay suggests AtRAV1 to be a novel phosphorylation target of AtMPK3. Bimolecular fluorescence complementation and yeast two-hybrid assays support physical interactions between AtRAV1 and AtMPK3. Overexpression of the native as well as phospho-mimic but not the phospho-defective variant of AtRAV1 imparts disease resistance in the atrav1 mutant A. thaliana lines. On the other hand, overexpression of AtRAV1 fails to impart disease resistance in the atmpk3 mutant. These analyses emphasize that AtMPK3-mediated phosphorylation of AtRAV1 is important for the elaboration of the defense response in A. thaliana. Considering that RAV1 homologs are conserved in diverse plant species, we propose that they can be gainfully deployed to impart disease resistance in agriculturally important crop plants. Indeed, overexpression of SlRAV1 (a member of the RAV1 family) imparts disease tolerance against not only fungal (R. solani and B. cinerea), but also against bacterial (Ralstonia solanacearum) pathogens in tomato, whereas silencing of the gene enhances disease susceptibility.

Genome-Wide Characterization and Expression Analysis of Transcription Factor Families in Desert Moss Syntrichia caninervis under Abiotic Stresses

Transcription factor (TF) families play important roles in plant stress responses. S. caninervis is a new model moss for plant desiccation tolerance studies. Here, we report a high-confidence identification and characterization of 591 TFs representing 52 families that covered all chromosomes in S. caninervis. GO term and KEGG pathway analysis showed that TFs were involved in the regulation of transcription, DNA-templated, gene expression, binding activities, plant hormone signal transduction, and circadian rhythm. A number of TF promoter regions have a mixture of various hormones-related cis-regulatory elements. AP2/ERF, bHLH, MYB, and C2H2-zinc finger TFs were the overrepresented TF families in S. caninervis, and the detailed classification of each family is performed based on structural features. Transcriptome analysis revealed the transcript abundances of some ScAP2/ERF, bHLH, MYB, and C2H2 genes were accumulated in the treated S. caninervis under cold, dehydration, and rehydration stresses. The RT-qPCR results strongly agreed with RNA-seq analysis, indicating these TFs might play a key role in S. caninervis response to abiotic stress. Our comparative TF characterization and classification provide the foundations for functional investigations of the dominant TF genes involved in S. caninervis stress response, as well as excellent stress tolerance gene resources for plant stress resistance breeding.

Association mapping for broomrape resistance in sunflower

INTRODUCTION

Sunflower breeding for resistance to the parasitic plant sunflower broomrape (Orobanche cumana Wallr.) requires the identification of novel resistance genes. In this research, we conducted a genome-wide association study (GWAS) to identify QTLs associated with broomrape resistance.

METHODS

The marker-trait associations were examined across a germplasm set composed of 104 sunflower accessions. They were genotyped with a 600k AXIOM® genome-wide array and evaluated for resistance to three populations of the parasite with varying levels of virulence (races EFR, FGV, and GTK) in two environments.

RESULTS AND DISCUSSION

The analysis of the genetic structure of the germplasm set revealed the presence of two main groups. The application of optimized treatments based on the general linear model (GLM) and the mixed linear model (MLM) allowed the detection of 14 SNP markers significantly associated with broomrape resistance. The highest number of marker-trait associations were identified on chromosome 3, clustered in two different genomic regions of this chromosome. Other associations were identified on chromosomes 5, 10, 13, and 16. Candidate genes for the main genomic regions associated with broomrape resistance were studied and discussed. Particularly, two significant SNPs on chromosome 3 associated with races EFR and FGV were found at two tightly linked SWEET sugar transporter genes. The results of this study have confirmed the role of some QTL on resistance to sunflower broomrape and have revealed new ones that may play an important role in the development of durable resistance to this parasitic weed in sunflower.

Genome-Wide Characterization and Evolutionary Expansion of Poplar NAC Transcription Factors and Their Tissue-Specific Expression Profiles under Drought

The NAC (NAM, ATAF1/2 and CUC2) is a large gene family of plant-specific transcription factors that play a pivotal role in various physiological processes and abiotic stresses. Due to the lack of genome-wide characterization, intraspecific and interspecific synteny, and drought-responsive expression pattern of NAC genes in poplar, the functional characterization of drought-related NAC genes have been scarcely reported in Populus species. Here, we identified a total of 170 NAC domain-containing genes in the P. trichocarpa genome, 169 of which were unevenly distributed on its nineteen chromosomes. These NAC genes were phylogenetically divided into twenty subgroups, some of which exhibited a similar pattern of exon-intron architecture. The synteny and Ka/Ks analysis indicated that the expansion of NAC genes in poplar was mainly due to gene duplication events occurring before and after the divergence of Populus and Salix. Ten PdNAC (P. deltoids × P. euramericana cv.'Nanlin895') genes were randomly selected and cloned. Their drought-responsive expression profiles showed a tissue-specific pattern. The transcription factor PdNAC013 was verified to be localized in the nucleus. Our research results provide genomic information for the expansion of NAC genes in the poplar genome, and for further characterizing putative poplar NAC genes associated with water-deficit.

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.

Comparative transcriptomic provides novel insights into the soybean response to Colletotrichum truncatum infection

INTRODUCTION

Soybean (Glycine max) is among the most important crops in the world, and its production can be threatened by biotic diseases, such as anthracnose. Soybean anthracnose is a seed-borne disease mainly caused by the hemibiotrophic fungus Colletotrichum truncatum. Typical symptoms are pre- and post-emergence damping off and necrotic lesions on cotyledons, petioles, leaves, and pods. Anthracnose symptoms can appear early in the field, causing major losses to soybean production.

MATERIAL AND METHODS

In preliminary experiments, we observed that the same soybean cultivar can have a range of susceptibility towards different strains of C. truncatum, while the same C. truncatum strain can cause varying levels of disease severity in different soybean cultivars. To gain a better understanding of the molecular mechanisms regulating the early response of different soybean cultivars to different C. truncatum strains, we performed pathogenicity assays to select two soybean cultivars with significantly different susceptibility to two different C. truncatum strains and analyzed their transcriptome profiles at different time points of interaction (0, 12, 48, and 120 h post-inoculation, hpi).

RESULTS AND DISCUSSION

The pathogenicity assays showed that the soybean cultivar Gm1 is more resistant to C. truncatum strain 1080, and it is highly susceptible to strain 1059, while cultivar Gm2 shows the opposite behavior. However, if only trivial anthracnose symptoms appeared in the more resistant phenotype (MRP; Gm1-1080; Gm2-1059) upon 120 hpi, in the more susceptible phenotype (MSP; Gm-1059; Gm2- 1080) plants show mild symptoms already at 72 hpi, after which the disease evolved rapidly to severe necrosis and plant death. Interestingly, several genes related to different cellular responses of the plant immune system (pathogen recognition, signaling events, transcriptional reprogramming, and defense-related genes) were commonly modulated at the same time points only in both MRP. The list of differentially expressed genes (DEGs) specific to the more resistant combinations and related to different cellular responses of the plant immune system may shed light on the important host defense pathways against soybean anthracnose.

Transcriptome Profiling of the Resistance Response of Musa acuminata subsp. burmannicoides, var. Calcutta 4 to Pseudocercospora musae

Banana (Musa spp.), which is one of the world’s most popular and most traded fruits, is highly susceptible to pests and diseases. Pseudocercospora musae, responsible for Sigatoka leaf spot disease, is a principal fungal pathogen of Musa spp., resulting in serious economic damage to cultivars in the Cavendish subgroup. The aim of this study was to characterize genetic components of the early immune response to P. musae in Musa acuminata subsp. burmannicoides, var. Calcutta 4, a resistant wild diploid. Leaf RNA samples were extracted from Calcutta 4 three days after inoculation with fungal conidiospores, with paired-end sequencing conducted in inoculated and non-inoculated controls using lllumina HiSeq 4000 technology. Following mapping to the reference M. acuminata ssp. malaccensis var. Pahang genome, differentially expressed genes (DEGs) were identified and expression representation analyzed on the basis of gene ontology enrichment, Kyoto Encyclopedia of Genes and Genomes orthology and MapMan pathway analysis. Sequence data mapped to 29,757 gene transcript models in the reference Musa genome. A total of 1073 DEGs were identified in pathogen-inoculated cDNA libraries, in comparison to non-inoculated controls, with 32% overexpressed. GO enrichment analysis revealed common assignment to terms that included chitin binding, chitinase activity, pattern binding, oxidoreductase activity and transcription factor (TF) activity. Allocation to KEGG pathways revealed DEGs associated with environmental information processing, signaling, biosynthesis of secondary metabolites, and metabolism of terpenoids and polyketides. With 144 up-regulated DEGs potentially involved in biotic stress response pathways, including genes involved in cell wall reinforcement, PTI responses, TF regulation, phytohormone signaling and secondary metabolism, data demonstrated diverse early-stage defense responses to P. musae. With increased understanding of the defense responses occurring during the incompatible interaction in resistant Calcutta 4, these data are appropriate for the development of effective disease management approaches based on genetic improvement through introgression of candidate genes in superior cultivars.

A chitinase CsChi23 promoter polymorphism underlies cucumber resistance against Fusarium oxysporum f. sp. cucumerinum

Fusarium wilt disease, caused by Fusarium oxysporum f. sp. cucumerinum (Foc), leads to widespread yield loss and quality decline in cucumber. However, the molecular mechanisms underlying Foc resistance remain poorly understood. We report the mapping and functional characterisation of CsChi23, encoding a cucumber class I chitinase with antifungal properties. We assessed sequence variations at CsChi23 and the associated defence response against Foc. We functionally characterised CsChi23 using transgenic assay and expression analysis. The mechanism regulating CsChi23 expression was assessed by genetic and molecular approaches. CsChi23 was induced by Foc infection, which led to rapid upregulation in resistant cucumber lines. Overexpressing CsChi23 enhanced fusarium wilt resistance and reduced fungal biomass accumulation, whereas silencing CsChi23 causes loss of resistance. CsHB15, a homeodomain leucine zipper (HD-Zip) III transcription factor, was found to bind to the CsChi23 promoter region and activate its expression. Furthermore, silencing of CsHB15 reduces CsChi23 expression. A single-nucleotide polymorphism variation -400 bp upstream of CsChi23 abolished the HD-Zip III binding site in a susceptible cucumber line. Collectively, our study indicates that CsChi23 is sufficient to enhance fusarium wilt resistance and reveals a novel function of an HD-Zip III transcription factor in regulating chitinase expression in cucumber defence against fusarium wilt.

Integrated Omic Approaches Reveal Molecular Mechanisms of Tolerance during Soybean and Meloidogyne incognita Interactions

The root-knot nematode (RKN), Meloidogyne incognita, is a devastating soybean pathogen worldwide. The use of resistant cultivars is the most effective method to prevent economic losses caused by RKNs. To elucidate the mechanisms involved in resistance to RKN, we determined the proteome and transcriptome profiles from roots of susceptible (BRS133) and highly tolerant (PI 595099) Glycine max genotypes 4, 12, and 30 days after RKN infestation. After in silico analysis, we described major defense molecules and mechanisms considered constitutive responses to nematode infestation, such as mTOR, PI3K-Akt, relaxin, and thermogenesis. The integrated data allowed us to identify protein families and metabolic pathways exclusively regulated in tolerant soybean genotypes. Among them, we highlighted the phenylpropanoid pathway as an early, robust, and systemic defense process capable of controlling M. incognita reproduction. Associated with this metabolic pathway, 29 differentially expressed genes encoding 11 different enzymes were identified, mainly from the flavonoid and derivative pathways. Based on differential expression in transcriptomic and proteomic data, as well as in the expression profile by RT-qPCR, and previous studies, we selected and overexpressed the GmPR10 gene in transgenic tobacco to assess its protective effect against M. incognita. Transgenic plants of the T₂ generation showed up to 58% reduction in the M. incognita reproduction factor. Finally, data suggest that GmPR10 overexpression can be effective against the plant parasitic nematode M. incognita, but its mechanism of action remains unclear. These findings will help develop new engineered soybean genotypes with higher performance in response to RKN infections.

Silver and copper-oxide nanoparticles prepared with GA3 induced defense in rice plants and caused mortalities to the brown planthopper, Nilaparvata lugens (Stål)

BACKGROUND

Nanoparticles have been employed as nanopesticides for pest control in agriculture. However, the harmful effects of their chemical synthesis on human and environmental health have resulted in increased use of green synthetic approaches, including the use of plant extracts. The brown planthopper, Nilaparvata lugens (Stål) (BPH), is a severe pest of rice plants (Oryza sativa L.), especially in Asia. It is usually controlled chemically but has developed resistance against many insecticides.

RESULTS

In this study, we synthesized metallic silver (Ag-NPs) and copper-oxide (CuO-NPs) nanoparticles using the exogenous phytohormone, gibberellic acid (GA₃), as a reducing agent. We then sprayed them separately on rice plants and BPH together and evaluated their effects on the plants and insects. SEM and TEM images showed that the synthesis was successful, indicated by the sizes (25-60 nm), uniform shape and spherical and cubical structures of Ag-NPs, as well as by the rugby sheet-like of CuO-NPs with lateral sizes of 150-340 nm and thickness of 30-70 nm. Independent applications of the nanoparticles and GA₃ on rice plants induced different volatile profiles, of which the highest number emitted was under Ag-NPs, including the highest emission of linalool. Transcriptome analysis showed that Ag-NPs-treated rice plants showed different transcriptome profiles compared to the control, 24 h after treatment, including the upregulation of the linalool synthase gene, genes of plants transcription factors such as WRKY, bHLH and NAC and other genes involved in plant defense responses. In all treatments, the mortality rate of BPH increased with an increase in NPs concentrations over time but was prominent under Ag-NPs treatment. The LC₅₀ values for Ag-NPs and CuO-NPs decreased with an increase in time. Also, the nanoparticles increased the activities of protective enzymes (POD, SOD and CAT), inhibited that of detoxification enzymes (A-CHE, ACP and AKP), and reduced total protein concentrations in the BPH.

CONCLUSIONS

These results show that synthesizing nanoparticles using phytohormones may be a safer and environmentally friendly option, which also holds promise for controlling the BPH in rice production.

The NIN-Like Protein OsNLP2 Negatively Regulates Ferroptotic Cell Death and Immune Responses to Magnaporthe oryzae in Rice

Nodule inception (NIN)-like proteins (NLPs) have a central role in nitrate signaling to mediate plant growth and development. Here, we report that OsNLP2 negatively regulates ferroptotic cell death and immune responses in rice during Magnaporthe oryzae infection. OsNLP2 was localized to the plant cell nucleus, suggesting that it acts as a transcription factor. OsNLP2 expression was involved in susceptible disease development. ΔOsnlp2 knockout mutants exhibited reactive oxygen species (ROS) and iron-dependent ferroptotic hypersensitive response (HR) cell death in response to M. oryzae. Treatments with the iron chelator deferoxamine, lipid-ROS scavenger ferrostatin-1, actin polymerization inhibitor cytochalasin A, and NADPH oxidase inhibitor diphenyleneiodonium suppressed the accumulation of ROS and ferric ions, lipid peroxidation, and HR cell death, which ultimately led to successful M. oryzae colonization in ΔOsnlp2 mutants. The loss-of-function of OsNLP2 triggered the expression of defense-related genes including OsPBZ1, OsPIP-3A, OsWRKY104, and OsRbohB in ΔOsnlp2 mutants. ΔOsnlp2 mutants exhibited broad-spectrum, nonspecific resistance to diverse M. oryzae strains. These combined results suggest that OsNLP2 acts as a negative regulator of ferroptotic HR cell death and defense responses in rice, and may be a valuable gene source for molecular breeding of rice with broad-spectrum resistance to blast disease.

Gb_ANR-47 Enhances the Resistance of Gossypium barbadense to Fusarium oxysporum f. sp. vasinfectum (FOV) by Regulating the Content of Proanthocyanidins

Anthocyanidin reductase (ANR) is an important regulator of flavonoid metabolism, and proanthocyanidins, the secondary metabolites of flavonoids, play an important role in the response of plants to pathogenic stress. Therefore, in this study, the expression analysis of the ANR gene family of Gossypium barbadense after inoculation with Fusarium oxysporum f. sp. vasinfectum (FOV) was performed at different time points. It was found that Gb_ANR-47 showed significant differences in the disease-resistant cultivar 06-146 and the susceptible cultivar Xinhai 14, as well as in the highest root expression. It was found that the expression of Gb_ANR-47 in the resistant cultivar was significantly higher than that in the susceptible cultivar by MeJA and SA, and different amounts of methyl jasmonate (MeJA) and salicylic acid (SA) response elements were found in the promoter region of Gb_ANR-47. After silencing GbANR-47 in 06-146 material by VIGS technology, its resistance to FOV decreased significantly. The disease severity index (DSI) was significantly increased, and the anthocyanin content was significantly decreased in silenced plants, compared to controls. Our findings suggest that GbANR-47 is a positive regulator of FOV resistance in Gossypium barbadense. The research results provide an important theoretical basis for in-depth analysis of the molecular mechanism of GbANR-47 and improving the anti-FOV of Gossypium barbadense.

The Effector Protein CgNLP1 of Colletotrichum gloeosporioides Affects Invasion and Disrupts Nuclear Localization of Necrosis-Induced Transcription Factor HbMYB8-Like to Suppress Plant Defense Signaling

Fungi secrete numerous effectors to modulate host defense systems. Understanding the molecular mechanisms by which fungal effectors regulate plant defense is of great importance for the development of novel strategies for disease control. In this study, we identified necrosis- and ethylene-inducing protein 1 (Nep1)-like protein (NLP) effector gene, CgNLP1, which contributed to conidial germination, appressorium formation, invasive growth, and virulence of Colletotrichum gloeosporioides to the rubber tree. Transient expression of CgNLP1 in the leaves of Nicotiana benthamiana induced ethylene production in plants. Ectopic expression of CgNLP1 in Arabidopsis significantly enhanced the resistance to Botrytis cinerea and Alternaria brassicicola. An R2R3 type transcription factor HbMYB8-like of rubber tree was identified as the target of CgNLP1.HbMYB8-like, localized on the nucleus, and induced cell death in N. benthamiana. CgNLP1 disrupted nuclear accumulation of HbMYB8-like and suppressed HbMYB8-like induced cell death, which is mediated by the salicylic acid (SA) signal pathway. This study suggested a new strategy whereby C. gloeosporioides exploited the CgNLP1 effector to affect invasion and suppress a host defense regulator HbMYB8-like to facilitate infection.

GWAS Reveals a Novel Candidate Gene CmoAP2/ERF in Pumpkin (Cucurbita moschata) Involved in Resistance to Powdery Mildew

Pumpkin (Cucurbita moschata Duchesne ex Poir.) is a multipurpose cash crop rich in antioxidants, minerals, and vitamins; the seeds are also a good source of quality oils. However, pumpkin is susceptible to the fungus Podosphaera xanthii, an obligate biotrophic pathogen, which usually causes powdery mildew (PM) on both sides of the leaves and reduces photosynthesis. The fruits of infected plants are often smaller than usual and unpalatable. This study identified a novel gene that involves PM resistance in pumpkins through a genome-wide association study (GWAS). The allelic variation identified in the CmoCh3G009850 gene encoding for AP2-like ethylene-responsive transcription factor (CmoAP2/ERF) was proven to be involved in PM resistance. Validation of the GWAS data revealed six single nucleotide polymorphism (SNP) variations in the CmoAP2/ERF coding sequence between the resistant (IT 274039 [PMR]) and the susceptible (IT 278592 [PMS]). A polymorphic marker (dCAPS) was developed based on the allelic diversity to differentiate these two haplotypes. Genetic analysis in the segregating population derived from PMS and PMR parents provided evidence for an incomplete dominant gene-mediated PM resistance. Further, the qRT-PCR assay validated the elevated expression of CmoAP2/ERF during PM infection in the PMR compared with PMS. These results highlighted the pivotal role of CmoAP2/ERF in conferring resistance to PM and identifies it as a valuable molecular entity for breeding resistant pumpkin cultivars.

Low Temperature Plasma Strategies for Xylella fastidiosa Inactivation

The quarantine bacterium Xylella fastidiosa was first detected in Salento (Apulia, Italy) in 2013 and caused severe symptoms in olives, leading to plant death. The disease, named Olive Quick Decline Syndrome (OQDS), is caused by the strain “De Donno” ST53 of the subspecies pauca of this bacterium (XfDD), which is spread by the insect Philaenus spumarius. The epidemic poses a serious threat to the agricultural economy and the landscape, as X. fastidiosa infects several plant species and there is yet no recognized solution. Research on OQDS is focused on finding strategies to control its spread or mitigate its symptoms. As a perspective solution, we investigated the efficacy of the low-temperature plasma and plasma-activated water to kill bacterial cells. Experiments were conducted in vitro to test the biocidal effect of the direct application of a Surface Dielectric Barrier Discharge (SDBD) plasma on bacteria cells and Plasma Activated Water (PAW). PAW activity was tested as a possible biocidal agent that can move freely in the xylem network paving the way to test the strategy on infected plants. The results showed a high decontamination rate even for cells of XfDD embedded in biofilms grown on solid media and complete inactivation in liquid culture medium.

Genome-wide investigation of maize RAD51 binding affinity through phage display

Background

RAD51 proteins, which are conserved in all eukaryotes, repair DNA double-strand breaks. This is critical to homologous chromosome pairing and recombination enabling successful reproduction. Work in Arabidopsis suggests that RAD51 also plays a role in plant defense; the Arabidopsis rad51 mutant is more susceptible to Pseudomonas syringae. However, the defense functions of RAD51 and the proteins interacting with RAD51 have not been thoroughly investigated in maize. Uncovering ligands of RAD51 would help to understand meiotic recombination and possibly the role of RAD51 in defense. This study used phage display, a tool for discovery of protein-protein interactions, to search for proteins interacting with maize RAD51A1.

Results

Maize RAD51A1 was screened against a random phage library. Eleven short peptide sequences were recovered from 15 phages which bound ZmRAD51A1 in vitro; three sequences were found in multiple successfully binding phages. Nine of these phage interactions were verified in vitro through ELISA and/or dot blotting.

BLAST searches did not reveal any maize proteins which contained the exact sequence of any of the selected phage peptides, although one of the selected phages had a strong alignment (E-value = 0.079) to a binding domain of maize BRCA2. Therefore, we designed 32 additional short peptides using amino acid sequences found in the predicted maize proteome. These peptides were not contained within phages. Of these synthesized peptides, 14 bound to ZmRAD51A1 in a dot blot experiment. These 14 sequences are found in known maize proteins including transcription factors putatively involved in defense.

Conclusions

These results reveal several peptides which bind ZmRAD51A1 and support a potential role for ZmRAD51A1 in transcriptional regulation and plant defense. This study also demonstrates the applicability of phage display to basic science questions, such as the search for binding partners of a known protein, and raises the possibility of an iterated approach to test peptide sequences that closely but imperfectly align with the selected phages.

Molecular insights into the jasmonate signaling and associated defense responses against wilt caused by Fusarium oxysporum

Biotic and abiotic stress factors drastically limit plant growth and development as well as alter the physiological, biochemical and cellular processes. This negatively impacts plant productivity, ultimately leading to agricultural and economical loss. Plant defense mechanisms elicited in response to these stressors are crucially regulated by the intricate crosstalk between defense hormones such as jasmonic acid (JA), salicylic acid and ethylene. These hormones orchestrate adaptive responses by modulating the gene regulatory networks leading to sequential changes in the root architecture, cell wall composition, secondary metabolite production and expression of defense-related genes. Fusarium wilt is a widespread vascular disease in plants caused by the soil-borne ascomycete Fusarium oxysporum and is known to attack several economically important plant cultivars. JA along with its conjugated forms methyl jasmonate and jasmonic acid isoleucine critically tunes plant defense mechanisms by regulating the expression of JA-associated genes imparting resistance phenotype. However, it should be noted that some members of F. oxysporum utilize the JA signaling pathway for disease development leading to susceptibility in plants. Therefore, JA signaling pathway becomes one of the important targets amenable for modulation to develop resistance response against Fusarium wilt in plants. In this review, we have emphasized on the physiological and molecular aspects of JA and its significant role in mounting an early defense response against Fusarium wilt disease. Further, utilization of the inherent JA signaling pathway and/or exogenous application of JA in generating Fusarium wilt resistant plants is discussed.

Integrated Metabolomic and Transcriptomic Analyses to Understand the Effects of Hydrogen Water on the Roots of Ficus hirta Vahl

Wuzhimaotao (Ficus hirta Vahl) is an important medicinal and edible plant in China. The extract from the roots of Ficus hirta Vahl contains phenylpropanoid compounds, such as coumarins and flavonoids, which are the main active components of this Chinese herbal medicine. In this study, we analyzed the transcriptomic and metabolomic data of the hydrogen-water-treated roots of Ficus hirta Vahl and a control group. The results showed that many genes and metabolites were regulated in the roots of Ficus hirta Vahl that were treated with hydrogen water. Compared with the control group, 173 genes were downregulated and 138 genes were upregulated in the hydrogen-rich water treatment group. Differential metabolite analysis through LC-MS showed that 168 and 109 metabolites had significant differences in positive and negative ion mode, respectively. In the upregulated metabolites, the main active components of Wuzhimaotao, such as the phenylpropane compounds naringin, bergaptol, hesperidin, and benzofuran, were found. Integrated transcriptomic and metabolomic data analysis showed that four and one of the most relevant pathways were over enriched in positive and negative ion mode, respectively. In the relationship between metabolites and DEGs, phenylpropanoid biosynthesis and metabolism play an important role. This indicates that phenylpropanoid biosynthesis and metabolism may be the main metabolic pathways regulated by hydrogen water. Our transcriptome analysis showed that most of the DEGs with |log2FC| ≥ 1 are transcription factor genes, and most of them are related to plant hormone signal transduction, stress resistance, and secondary metabolism, mainly phenylpropanoid biosynthesis and metabolism. This study provides important evidence and clues for revealing the botanical effect mechanism of hydrogen and a theoretical basis for the application of hydrogen agriculture in the cultivation of Chinese herbal medicine.

Suberin Biosynthesis, Assembly, and Regulation

Suberin is a specialized cell wall modifying polymer comprising both phenolic-derived and fatty acid-derived monomers, which is deposited in below-ground dermal tissues (epidermis, endodermis, periderm) and above-ground periderm (i.e., bark). Suberized cells are largely impermeable to water and provide a critical protective layer preventing water loss and pathogen infection. The deposition of suberin is part of the skin maturation process of important tuber crops such as potato and can affect storage longevity. Historically, the term "suberin" has been used to describe a polyester of largely aliphatic monomers (fatty acids, ω-hydroxy fatty acids, α,ω-dioic acids, 1-alkanols), hydroxycinnamic acids, and glycerol. However, exhaustive alkaline hydrolysis, which removes esterified aliphatics and phenolics from suberized tissue, reveals a core poly(phenolic) macromolecule, the depolymerization of which yields phenolics not found in the aliphatic polyester. Time course analysis of suberin deposition, at both the transcriptional and metabolite levels, supports a temporal regulation of suberin deposition, with phenolics being polymerized into a poly(phenolic) domain in advance of the bulk of the poly(aliphatics) that characterize suberized cells. In the present review, we summarize the literature describing suberin monomer biosynthesis and speculate on aspects of suberin assembly. In addition, we highlight recent advances in our understanding of how suberization may be regulated, including at the phytohormone, transcription factor, and protein scaffold levels.

Defense Strategies: The Role of Transcription Factors in Tomato-Pathogen Interaction

Simple Summary

Tomato is one of the most cultivated and economically important vegetable crops throughout the world. It is affected by a panoply of different pathogens that cause infectious diseases that reduce tomato yield and affect product quality, with the most common symptoms being wilts, leaf spots/blights, fruit spots, and rots. To survive, tomato, as other plants, have developed elaborate defense mechanisms against plant pathogens. Among several genes already identified in tomato response to pathogens, we highlight those encoding the transcription factors (TFs). TFs are regulators of gene expression and are involved in large-scale biological phenomena. Here, we present an overview of recent studies of tomato TFs regarding defense responses to pathogen attack, selected for their abundance, importance, and availability of functionally well-characterized members. Tomato TFs’ roles and the possibilities related to their use for genetic engineering in view of crop breeding are presented.

Abstract

Tomato, one of the most cultivated and economically important vegetable crops throughout the world, is affected by a panoply of different pathogens that reduce yield and affect product quality. The study of tomato–pathogen system arises as an ideal system for better understanding the molecular mechanisms underlying disease resistance, offering an opportunity of improving yield and quality of the products. Among several genes already identified in tomato response to pathogens, we highlight those encoding the transcription factors (TFs). TFs act as transcriptional activators or repressors of gene expression and are involved in large-scale biological phenomena. They are key regulators of central components of plant innate immune system and basal defense in diverse biological processes, including defense responses to pathogens. Here, we present an overview of recent studies of tomato TFs regarding defense responses to biotic stresses. Hence, we focus on different families of TFs, selected for their abundance, importance, and availability of functionally well-characterized members in response to pathogen attack. Tomato TFs’ roles and possibilities related to their use for engineering pathogen resistance in tomato are presented. With this review, we intend to provide new insights into the regulation of tomato defense mechanisms against invading pathogens in view of plant breeding.

Dynamic Transcriptome Profiling of Mungbean Genotypes Unveil the Genes Respond to the Infection of Mungbean Yellow Mosaic Virus

Yellow mosaic disease (YMD), incited by mungbean yellow mosaic virus (MYMV), is a primary viral disease that reduces mungbean production in South Asia, especially in India. There is no detailed knowledge regarding the genes and molecular mechanisms conferring resistance of mungbean to MYMV. Therefore, disclosing the genetic and molecular bases related to MYMV resistance helps to develop the mungbean genotypes with MYMV resistance. In this study, transcriptomes of mungbean genotypes, VGGRU-1 (resistant) and VRM (Gg) 1 (susceptible) infected with MYMV were compared to those of uninfected controls. The number of differentially expressed genes (DEGs) in the resistant and susceptible genotypes was 896 and 506, respectively. Among them, 275 DEGs were common between the resistant and susceptible genotypes. Functional annotation of DEGs revealed that the DEGs belonged to the following categories defense and pathogenesis, receptor-like kinases; serine/threonine protein kinases, hormone signaling, transcription factors, and chaperons, and secondary metabolites. Further, we have confirmed the expression pattern of several DEGs by quantitative real-time PCR (qRT-PCR) analysis. Collectively, the information obtained in this study unveils the new insights into characterizing the MYMV resistance and paved the way for breeding MYMV resistant mungbean in the future.

Deciphering the Crosstalk Mechanisms of Wheat-Stem Rust Pathosystem: Genome-Scale Prediction Unravels Novel Host Targets

Triticum aestivum (wheat), a major staple food grain, is affected by various biotic stresses. Among these, fungal diseases cause about 15-20% of yield loss, worldwide. In this study, we performed a comparative analysis of protein-protein interactions between two Puccinia graminis races (Pgt 21-0 and Pgt Ug99) that cause stem (black) rust in wheat. The available molecular techniques to study the host-pathogen interaction mechanisms are expensive and labor-intensive. We implemented two computational approaches (interolog and domain-based) for the prediction of PPIs and performed various functional analysis to determine the significant differences between the two pathogen races. The analysis revealed that T. aestivum-Pgt 21-0 and T. aestivum-Pgt Ug99 interactomes consisted of ∼90M and ∼56M putative PPIs, respectively. In the predicted PPIs, we identified 115 Pgt 21-0 and 34 Pgt Ug99 potential effectors that were highly involved in pathogen virulence and development. Functional enrichment analysis of the host proteins revealed significant GO terms and KEGG pathways such as O-methyltransferase activity (GO:0008171), regulation of signal transduction (GO:0009966), lignin metabolic process (GO:0009808), plastid envelope (GO:0009526), plant-pathogen interaction pathway (ko04626), and MAPK pathway (ko04016) that are actively involved in plant defense and immune signaling against the biotic stresses. Subcellular localization analysis anticipated the host plastid as a primary target for pathogen attack. The highly connected host hubs in the protein interaction network belonged to protein kinase domain including Ser/Thr protein kinase, MAPK, and cyclin-dependent kinase. We also identified 5,577 transcription factors in the interactions, associated with plant defense during biotic stress conditions. Additionally, novel host targets that are resistant to stem rust disease were also identified. The present study elucidates the functional differences between Pgt 21-0 and Pgt Ug99, thus providing the researchers with strain-specific information for further experimental validation of the interactions, and the development of durable, disease-resistant crop lines.

Interaction between Rag genes results in a unique synergistic transcriptional response that enhances soybean resistance to soybean aphids

BACKGROUND

Pyramiding different resistance genes into one plant genotype confers enhanced resistance at the phenotypic level, but the molecular mechanisms underlying this effect are not well-understood. In soybean, aphid resistance is conferred by Rag genes. We compared the transcriptional response of four soybean genotypes to aphid feeding to assess how the combination of Rag genes enhanced the soybean resistance to aphid infestation.

RESULTS

A strong synergistic interaction between Rag1 and Rag2, defined as genes differentially expressed only in the pyramid genotype, was identified. This synergistic effect in the Rag1/2 phenotype was very evident early (6 h after infestation) and involved unique biological processes. However, the response of susceptible and resistant genotypes had a large overlap 12 h after aphid infestation. Transcription factor (TF) analyses identified a network of interacting TF that potentially integrates signaling from Rag1 and Rag2 to produce the unique Rag1/2 response. Pyramiding resulted in rapid induction of phytochemicals production and deposition of lignin to strengthen the secondary cell wall, while repressing photosynthesis. We also identified Glyma.07G063700 as a novel, strong candidate for the Rag1 gene.

CONCLUSIONS

The synergistic interaction between Rag1 and Rag2 in the Rag1/2 genotype can explain its enhanced resistance phenotype. Understanding molecular mechanisms that support enhanced resistance in pyramid genotypes could facilitate more directed approaches for crop improvement.

Transcriptomics and Metabolomics Changes Triggered by Inflorescence Removal in Panax notoginseng (Burk.)

The root of Panax notoginseng (Burk.), in which saponins are the major active components, is a famous traditional Chinese medicine used to stop bleeding and to decrease inflammation and heart disease. Inflorescence removal increases the yield and quality of P. notoginseng, but the underlying molecular mechanisms are unknown. Here, the differences between inflorescence-removal treatment and control groups of P. notoginseng were compared using transcriptomics and metabolomics analyses. Illumina sequencing of cDNA libraries prepared from the rhizomes, leaves and roots of the two groups independently identified 6,464, 4,584, and 7,220 differentially expressed genes (DEG), respectively. In total, 345 differentially expressed transcription factors (TFs), including MYB and WRKY family members, were induced by the inflorescence-removal treatment. Additionally, 215 DEGs involved in saponin terpenoid backbone biosynthetic pathways were identified. Most genes involved in the mevalonic acid (MVA) and methylerythritol phosphate (MEP) pathways were activated by inflorescence removal. The co-expression analysis showed that the low expression levels of flavonoid biosynthesis-related genes (e.g., C4H and F3H) decreased the biosynthesis and accumulation of some flavonoids after inflorescence removal. The results not only provide new insights into the fundamental mechanisms underlying the poorly studied inflorescence-removal process in P. notoginseng and other rhizome crops, but they also represent an important resource for future research on gene functions during inflorescence-removal treatments and the reproductive stage.

Elucidation of the biochemical and molecular basis of the differential disease expression in two cultivars of chili (Capsicum annuum) in response to Colletotrichum capsici infection

Chili plants are affected by the hemibiotrophic ascomycota fungus Colletotrichum capsici causing Anthracnose. Infection results in yield and marketability loss due to a decrease in the quality of fruits. The study of morphological symptom development in two cultivars, Bullet, and Beldanga, showed very different disease expression pattern. To understand the reasons behind such differential response, we investigated, in a time-dependent manner, biochemical activities of important defense enzymes, PR proteins, like peroxidase, polyphenol-oxidase, phenylalanine ammonia lyase, β-glucanase, chitinase, catalase, as well as phenols, flavonoids, chlorophyll and the key signaling molecule nitric oxide in their leaves. We further performed real-time nitric oxide (NO) detection studies. The results showed striking differences in the activity profile of these defense molecules through the course of the study. We monitored the gene expression levels of 12 important defense-related genes under in vivo condition. The transcription levels were mostly increased in the tolerant cultivar till 7 days post-infection (DPI), while downregulation of some of the genes were observed in the susceptible one. These data indicated that disease manifestation is a simulated response of these defense molecules which can nullify the effect of the pathogen and its products, when resistance occurs. Alternatively, the pathogen suppresses the host defense when the disease develops.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11738-021-03334-x.

Comparative transcriptomic and metabolomic analyses reveal the protective effects of silicon against low phosphorus stress in tomato plants

Phosphorus (P) is an essential nutrient controlling plant growth and development through the regulation of basic metabolic processes. Soil P deficiency is one of the major limiting factors for sustainable crop production worldwide. Previous studies have demonstrated that silicon (Si), as a beneficial element, promotes plant nutrition, growth, development, and responses to low P (LP) stress; however, the molecular mechanisms underlying Si-mediated LP tolerance remain largely unclear. Here, we found that LP + Si treatment increased the net photosynthetic rate and shoot fresh weight by 34.3%, and 121.3%, respectively compared with LP alone. RNA-sequencing and metabolomic analyses were subsequently performed with tomato plants grown under control and P depleted conditions with or without Si amendment. RNA-sequencing showed that Si supply alters not only the expression of genes involved in the metabolism of carbon (C), nitrogen (N), and P but also phosphorylation processes and metabolism of glutathione and reactive active oxygen in tomato roots. Si also affected the expression of genes encoding major transcription factors such as WRKY and MYB under LP stress. Moreover, a set of genes encoding the enzymes or regulators of organic acid (OA) metabolism or secretion were differentially expressed in Si-treated P deficient roots compared with those in LP stress alone. Furthermore, the metabolomic analysis showed that the levels of several OAs were significantly elevated in Si-treated P deficient roots. Taken together, these results indicate that exogenous Si increases the secretion of OAs by modulating C/N metabolism in LP-treated tomato roots and thereby improving plant growth under LP stress.

Proteomic and Transcriptomic Analyses Provide Novel Insights into the Crucial Roles of Host-Induced Carbohydrate Metabolism Enzymes in Xanthomonas oryzae pv. oryzae Virulence and Rice-Xoo Interaction

BACKGROUND

Xanthomonas oryzae pv. oryzae (Xoo) causes bacterial leaf blight, a devastating rice disease. The Xoo-rice interaction, wherein wide ranging host- and pathogen-derived proteins and genes wage molecular arms race, is a research hotspot. Hence, the identification of novel rice-induced Xoo virulence factors and characterization of their roles affecting rice global gene expression profiles will provide an integrated and better understanding of Xoo-rice interactions from the molecular perspective.

RESULTS

Using comparative proteomics and an in vitro interaction system, we revealed that 5 protein spots from Xoo exhibited significantly different expression patterns (|fold change| > 1.5) at 3, 6, 12 h after susceptible rice leaf extract (RLX) treatment. MALDI-TOF MS analysis and pathogenicity tests showed that 4 host-induced proteins, including phosphohexose mutase, inositol monophosphatase, arginase and septum site-determining protein, affected Xoo virulence. Among them, mutants of two host-induced carbohydrate metabolism enzyme-encoding genes, ΔxanA and Δimp, elicited enhanced defense responses and nearly abolished Xoo virulence in rice. To decipher rice differentially expressed genes (DEGs) associated with xanA and imp, transcriptomic responses of ΔxanA-treated and Δimp-treated susceptible rice were compared to those in rice treated with PXO99^A at 1 and 3 dpi. A total of 1521 and 227 DEGs were identified for PXO99^A vs Δimp at 1 and 3 dpi, while for PXO99^A vs ΔxanA, there were 131 and 106 DEGs, respectively. GO, KEGG and MapMan analyses revealed that the DEGs for PXO99^A vs Δimp were mainly involved in photosynthesis, signal transduction, transcription, oxidation-reduction, hydrogen peroxide catabolism, ion transport, phenylpropanoid biosynthesis and metabolism of carbohydrates, lipids, amino acids, secondary metabolites, hormones, and nucleotides, while the DEGs from PXO99^A vs ΔxanA were predominantly associated with photosynthesis, signal transduction, oxidation-reduction, phenylpropanoid biosynthesis, cytochrome P450 and metabolism of carbohydrates, lipids, amino acids, secondary metabolites and hormones. Although most pathways were associated with both the Δimp and ΔxanA treatments, the underlying genes were not the same.

CONCLUSION

Our study identified two novel host-induced virulence factors XanA and Imp in Xoo, and revealed their roles in global gene expression in susceptible rice. These results provide valuable insights into the molecular mechanisms of pathogen infection strategies and plant immunity.

Tolerance mechanisms in maize identified through phenotyping and transcriptome analysis in response to water deficit stress

Water deficit is a key limiting factor for maize (Zea mays L.) productivity. Elucidating the molecular regulatory networks of stress tolerance is crucial for genetic enhancement of drought tolerance. Two genotypes of maize contrasting in their yield response to water deficit were evaluated for tolerance traits of water relations, net CO₂ assimilation rate, antioxidative metabolism and grain yield in relation to the expression levels, based on transcription profiling of genes involved in stress signaling, protein processing and energy metabolism to identify functional tolerance mechanisms. In the genotype SNJ201126 upregulation of calcium mediated signaling, plasma membrane and tonoplast intrinsic proteins and the membrane associated transporters contributed to better maintenance of water relations as evident from the higher relative water content and stomatal conductance at seedling and anthesis stages coupled with robust photosynthetic capacity and antioxidative metabolism. Further the protein folding machinery consisting of calnexin/calreticulin (CNX/CRT) cycle was significantly upregulated only in SNJ201126. While the down regulation of genes involved in photosystems and the enzymes of carbon fixation led to the relative susceptibility of genotype HKI161 in terms of reduced net CO₂ assimilation rate, biomass and grain yield. Our results provide new insight into intrinsic functional mechanisms related to tolerance in maize.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01003-4.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

MicroRNA319a regulates plant resistance to Sclerotinia stem rot

MicroRNA319a (miR319a) controls cell division arrest in plant leaves by inhibiting the expression of TCP (TEOSINTE BRANCHED 1/CYCLOIDEA/PCF) family genes. However, it is unclear whether miR319a influences infection by necrotrophic pathogens and host susceptibility. In this study, we revealed that miR319a affects plant resistance to stem rot disease caused by Sclerotinia sclerotiorum. In Brassica rapa plants infected with S. sclerotiorum, miR319a levels increased while the expression levels of several BraTCP genes significantly decreased compared with those of uninfected plants. Overexpression of BraMIR319a in B. rapa increased the susceptibility of the plants to S. sclerotiorum and aggravated stem rot disease, whereas overexpression of BraTCP4-1 promoted plant resistance. RNA sequencing data revealed a potential relationship between miR319a and pathogen-related WRKY genes. Chromatin immunoprecipitation, electrophoretic mobility shift, and reporter transaction assays showed that BraTCP4-1 could bind to the promoters of WRKY75, WRKY70, and WRKY33 and directly activate these pathogen-related genes. Moreover, the expression levels of WRKY75, WRKY70, and WRKY33 in plants overexpressing BraMIR319a decreased significantly, whereas those of plants overexpressing BraTCP4-1 increased significantly, relative to the wild type. These results suggest that miR319a and its target gene BraTCP4 control stem rot resistance through pathways of WRKY genes.

Genome wide identification and expression analysis of pepper C2H2 zinc finger transcription factors in response to anthracnose pathogen Colletotrichum truncatum

Although, the C₂H₂ zinc finger (ZF) family of plant transcription factors have been implicated in multiple biological processes, they are yet to be characterized in the economically important chilli pepper (Capsicum annuum). In this study, a total of 79 C₂H₂ ZF genes were identified in the pepper genome. Phylogenetic analysis categorized the pepper C₂H₂ ZF (CaZF) members into five subfamilies each with unique conserved domains and functions. Genomic organization revealed that CaZF genes have variable number of introns consistent with the characteristics defined by the evolutionary analysis. Segmental duplication-based purifying selection contributed to the expansion of CaZF genes in pepper. Additionally, 11 CaZF genes were identified as targets for 38 miRNAs indicating their role in post-transcriptional silencing-mediated genetic regulation. Gene expression analysis revealed that 18 CaZF genes were differentially expressed post-infection with the anthrocnose pathogen Colletotrichum truncatum, uncovering their potential function in pepper response to biotic stresses. Moreover, CaZFs were significantly induced post-treatment with methyl jasmonate and ethylene indicating their role in defense signaling. Notably, the MeJA responsive cis-elements were detected in the promoter regions of majority of CaZF genes, suggesting that CaZFs may be implicated in defense-responsive signal cross talking. Additionally, 18 CaZF genes were differentially expressed under drought and heat treatment, indicating their involvement in plant response to abiotic stresses. Overall, a comprehensive analysis of CaZF gene family in pepper provided significant insights into the understanding of C₂H₂ ZF-mediated stress regulation network, which would benefit the genetic improvement of pepper and other allied plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02601-x.

A Novel Role of Salt- and Drought-Induced RING 1 Protein in Modulating Plant Defense Against Hemibiotrophic and Necrotrophic Pathogens

Many plant-encoded E3 ligases are known to be involved in plant defense. Here, we report a novel role of E3 ligase SALT- AND DROUGHT-INDUCED RING FINGER1 (SDIR1) in plant immunity. Even though SDIR1 is reasonably well-characterized, its role in biotic stress response is not known. The silencing of SDIR1 in Nicotiana benthamiana reduced the multiplication of the virulent bacterial pathogen Pseudomonas syringae pv. tabaci. The Arabidopsis sdir1 mutant is resistant to virulent pathogens, whereas SDIR1 overexpression lines are susceptible to both host and nonhost hemibiotrophic bacterial pathogens. However, sdir1 mutant and SDIR1 overexpression lines showed hypersusceptibility and resistance, respectively, against the necrotrophic pathogen Erwinia carotovora. The mutant of SDIR1 target protein, i.e., SDIR-interacting protein 1 (SDIR1P1), also showed resistance to host and nonhost pathogens. In SDIR1 overexpression plants, transcripts of NAC transcription factors were less accumulated and the levels of jasmonic acid (JA) and abscisic acid were increased. In the sdir1 mutant, JA signaling genes JAZ7 and JAZ8 were downregulated. These data suggest that SDIR1 is a susceptibility factor and its activation or overexpression enhances disease caused by P. syringae pv. tomato DC3000 in Arabidopsis. Our results show a novel role of SDIR1 in modulating plant defense gene expression and plant immunity.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Transcription Factors as the "Blitzkrieg" of Plant Defense: A Pragmatic View of Nitric Oxide's Role in Gene Regulation

Plants are in continuous conflict with the environmental constraints and their sessile nature demands a fine-tuned, well-designed defense mechanism that can cope with a multitude of biotic and abiotic assaults. Therefore, plants have developed innate immunity, R-gene-mediated resistance, and systemic acquired resistance to ensure their survival. Transcription factors (TFs) are among the most important genetic components for the regulation of gene expression and several other biological processes. They bind to specific sequences in the DNA called transcription factor binding sites (TFBSs) that are present in the regulatory regions of genes. Depending on the environmental conditions, TFs can either enhance or suppress transcriptional processes. In the last couple of decades, nitric oxide (NO) emerged as a crucial molecule for signaling and regulating biological processes. Here, we have overviewed the plant defense system, the role of TFs in mediating the defense response, and that how NO can manipulate transcriptional changes including direct post-translational modifications of TFs. We also propose that NO might regulate gene expression by regulating the recruitment of RNA polymerase during transcription.

Microarray analysis of Arabidopsis thaliana exposed to single and mixed infections with Cucumber mosaic virus and turnip viruses

Cucumber mosaic virus (CMV), Turnip mosaic virus (TuMV) and Turnip crinkle virus (TCV) are important plant infecting viruses. In the present study, whole transcriptome alteration of Arabidopsis thaliana in response to CMV, TuMV and TCV, individual as well as mixed infections of CMV and TuMV/CMV and TCV were investigated using microarray data. In response to CMV, TuMV and TCV infections, a total of 2517, 3985 and 277 specific differentially expressed genes (DEGs) were up-regulated, while 2615, 3620 and 243 specific DEGs were down-regulated, respectively. The number of 1222 and 30 common DEGs were up-regulated during CMV and TuMV as well as CMV and TCV infections, while 914 and 24 common DEGs were respectively down-regulated. Genes encoding immune response mediators, signal transducer activity, signaling and stress response functions were among the most significantly upregulated genes during CMV and TuMV or CMV and TCV mixed infections. The NAC, C3H, C2H2, WRKY and bZIP were the most commonly presented transcription factor (TF) families in CMV and TuMV infection, while AP2-EREBP and C3H were the TF families involved in CMV and TCV infections. Moreover, analysis of miRNAs during CMV and TuMV and CMV and TCV infections have demonstrated the role of miRNAs in the down regulation of host genes in response to viral infections. These results identified the commonly expressed virus-responsive genes and pathways during plant-virus interaction which might develop novel antiviral strategies for improving plant resistance to mixed viral infections.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00925-3.

Strawberry FaWRKY25 Transcription Factor Negatively Regulated the Resistance of Strawberry Fruits to Botrytis cinerea

WRKY genes and jasmonic acid (JA) play a crucial role in plants’ responses against biotic and abiotic stress. However, the regulating mechanism of WRKY genes on strawberry fruits’ resistance against Botrytis cinerea is largely unknown, and few studies have been performed on their effect on the JA-mediated defense mechanism against B. cinerea. This study explored the effect of FaWRKY25 on the JA-mediated strawberry resistance against B. cinerea. Results showed that the JA content decreased significantly as the fruits matured, whereas the FaWRKY25 expression rose substantially, which led to heightened susceptibility to B. cinerea and in strawberries. External JA treatment significantly increased the JA content in strawberries and reduced the FaWRKY25 expression, thereby enhancing the fruits’ resistance against B. cinerea. FaWRKY25 overexpression significantly lowered the fruits’ resistance against B. cinerea, whereas FaWRKY25 silencing significantly increased resistance. Moreover, FaWRKY25 overexpression significantly lowered the JA content, whereas FaWRKY25 silencing significantly increased it. FaWRKY25 expression level substantially affects the expression levels of genes related to JA biosynthesis and metabolism, other members of the WRKY family, and defense genes. Accordingly, FaWRKY25 plays a crucial role in regulating strawberries’ resistance against B. cinerea and may negatively regulate their JA-mediated resistance mechanism against B. cinerea.

RNA-Seq Analysis of Prickled and Prickle-Free Epidermis Provides Insight into the Genetics of Prickle Development in Red Raspberry (Rubus ideaus L.)

Red raspberry (Rubus idaeus L.) is a globally commercialized specialty crop with growing demand worldwide. The presence of prickles on the stems, petioles and undersides of the leaves complicates both the field management and harvesting of raspberries. An RNA sequencing analysis was used to identify differentially expressed genes in the epidermal tissue of prickled “Caroline” and prickle-free “Joan J.” and their segregating progeny. Expression patterns of differentially expressed genes (DEGs) in prickle-free plants revealed the downregulation of some vital development-related transcription factors (TFs), including a MIXTA-like R2R3-MYB family member; MADS-box; APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) and NAM, ATAF1/2 and CUC2 (NAC) in prickle-free epidermis tissue. The downregulation of these TFs was confirmed by qRT-PCR analysis, indicating a key regulatory role in prickle development. This study adds to the understanding of prickle development mechanisms in red raspberries needed for utilizing genetic engineering strategies for developing prickle-free raspberry cultivars and, possibly, other Rubus species, such as blackberry (Rubus sp.) and black raspberry (R. occidentalis L.).