The Soybean Rhg1 locus for resistance to the soybean cyst nematode Heterodera glycines regulates the expression of a large number of stress- and defense-related genes in degenerating feeding cells

A novel toolbox of GATEWAY-compatible vectors for rapid functional gene analysis in soybean composite plants

KEY MESSAGE

We developed a set of GATEWAY vectors to accelerate gene function analysis in soybean composite plants to rapidly screen transgenic roots and investigate subcellular localization, protein-protein interactions, and root-pathogen interactions. The generation of transgenic plants is essential for plant biology research to investigate plant physiology, pathogen interactions, and gene function. However, producing stable transgenic plants for plants such as soybean is a laborious and time-consuming process, which can impede research progress. Composite plants consisting of wild-type shoots and transgenic roots are an alternative method for generating transgenic plant tissues that can facilitate functional analysis of genes-of-interest involved in root development or root-microbe interactions. In this report, we introduce a novel set of GATEWAY-compatible vectors that enable a wide range of molecular biology uses in roots of soybean composite plants. These vectors incorporate in-frame epitope fusions of green fluorescent protein, 3x-HA, or miniTurbo-ID, which can be easily fused to a gene-of-interest using the GATEWAY cloning system. Moreover, these vectors allow for the identification of transgenic roots using either mCherry fluorescence or the RUBY marker. We demonstrate the functionality of these vectors by expressing subcellular markers in soybean, providing evidence of their effectiveness in generating protein fusions in composite soybean plants. Furthermore, we show how these vectors can be used for gene function analysis by expressing the bacterial effector, AvrPphB in composite roots, enabling the identification of soybean targets via immunoprecipitation followed by mass spectrometry. Additionally, we demonstrate the successful expression of stable miniTurbo-ID fusion proteins in composite roots. Overall, this new set of vectors is a powerful tool that can be used to assess subcellular localization and perform gene function analyses in soybean roots without the need to generate stable transgenic plants.

Discovery of two tightly linked soybean genes at the qSCN10 (O) locus conferring broad-spectrum resistance to soybean cyst nematode

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is a major threat to global soybean yield. Resistance genes at the rhg1 locus from PI 88788 are majorly utilized in 95% of the U.S. breeding programs. Continuous use of this resistance source leads to a shift in the virulence of SCN populations and overcomes host resistance. Therefore, it is necessary to identify alternative SCN resistance sources to combat this ever-changing pest. Previously, we identified an exotic soybean line, PI 567516C, which carries a novel qSCN10 (O) locus for SCN resistance demonstrating different resistance responses compared to the known rhg1 and Rhg4 loci. Here, we narrowed the qSCN10 QTL region to 142-kb (containing 20 genes). Based on gene expression, gene ontology, in-silico analysis, and QTL-based haplotyping, two genes were identified for functional characterization. Overexpression of the transcription factor TGA1-related and Shugoshin C-terminus in the SCN-susceptible Williams 82 reduced the cyst number by 6.4-fold (84.6%) and 5.3-fold (81.2%), respectively. GmTGA1-10 and GmSCT-10 Tilling mutants showed high cyst numbers. The two genes associated with the qSCN10 QTL have significant potential to reduce the SCN population. They also offer an alternative source of durable SCN resistance that is independent of rhg1 and Rhg4.

Protein Dynamics in Plant Immunity: Insights into Plant-Pest Interactions

All living organisms regulate biological activities by proteins. When plants encounter pest invasions, the delicate balance between protein synthesis and degradation becomes even more pivotal for mounting an effective defense response. In this review, we summarize the mechanisms by which plants regulate their proteins to effectively coordinate immune responses during plant-pest interactions. Additionally, we discuss the main pathway proteins through which pest effectors manipulate host protein homeostasis in plants to facilitate their infestation. Understanding these processes at the molecular level not only deepens our knowledge of plant immunity but also holds the potential to inform strategies for developing pest-resistant crops, contributing to sustainable and resilient agriculture.

Differential Transcriptome Reprogramming Induced by the Soybean Cyst Nematode Type 0 and Type 1.2.5.7 During Resistant and Susceptible Interactions

Soybean cyst nematode (SCN, Heterodera glycines [Hg]) is a serious root parasite of soybean (Glycine max) that induces extensive gene expression changes associated with pleiotropic biological activities in infected cells. However, the impacts of various SCN Hg types on host transcriptome reprogramming remain largely unknown. Here, we developed and used two recombinant inbred lines (RIL; RIL-72 and RIL-137) to profile transcriptome reprogramming in the infection sites during the resistant and susceptible interactions with SCN Hg Type 1.2.5.7 and Type 0. SCN bioassays indicated that RIL-72 was susceptible to Type 1.2.5.7 but resistant to Type 0, whereas RIL-137 was resistant to both types. Comparative analysis of gene expression changes induced by Type 1.2.5.7 in the resistant and susceptible lines revealed distinct transcriptome regulation with a number of similarly and oppositely regulated genes. The expression levels of similarly regulated genes in the susceptible line appeared to be insufficient to mount an effective defense against SCN. The functional importance of oppositely regulated genes was confirmed using virus-induced gene silencing (VIGS) and overexpression approaches. Further transcriptome comparisons revealed shared as well as Hg type- and genotype-specific transcriptome reprogramming. Shared transcriptome responses were mediated through common SCN-responsive genes and conserved immune signaling, whereas genotype-specific responses were derived from genetic variability, metabolic and hormonal differences, and varied regulation of protein phosphorylation and ubiquitination. The conserved defense mechanisms together with genotype-specific responses would enable plants to trigger effective and tailored immune responses to various Hg types and adapt the defense response to their genetic backgrounds. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Soybean Cyst Nematode Effector Cysteine Protease 1 (CPR1) Targets a Mitochondrial Soybean Branched-Chain Amino Acid Aminotransferase (GmBCAT1)

The soybean cyst nematode (SCN; Heterodera glycines) facilitates infection by secreting a repertoire of effector proteins into host cells to establish a permanent feeding site composed of a syncytium of root cells. Among the diverse proteins secreted by the nematode, we were specifically interested in identifying proteases to pursue our goal of engineering decoy substrates that elicit an immune response when cleaved by an SCN protease. We identified a cysteine protease that we named Cysteine Protease 1 (CPR1), which was predicted to be a secreted effector based on transcriptomic data obtained from SCN esophageal gland cells, the presence of a signal peptide, and the lack of transmembrane domains. CPR1 is conserved in all isolates of SCN sequenced to date, suggesting it is critical for virulence. Transient expression of CPR1 in Nicotiana benthamiana leaves suppressed cell death induced by a constitutively active nucleotide binding leucine-rich repeat protein, RPS5, indicating that CPR1 inhibits effector-triggered immunity. CPR1 localizes in part to the mitochondria when expressed in planta. Proximity-based labeling in transgenic soybean roots, co-immunoprecipitation, and cleavage assays identified a branched-chain amino acid aminotransferase from soybean (GmBCAT1) as a substrate of CPR1. Consistent with this, GmBCAT1 also localizes to mitochondria. Silencing of the CPR1 transcript in the nematode reduced penetration frequency in soybean roots, while the expression of CPR1 in soybean roots enhanced susceptibility. Our data demonstrates that CPR1 is a conserved effector protease with a direct target in soybean roots, highlighting it as a promising candidate for decoy engineering. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The molecular dynamics between reactive oxygen species (ROS), reactive nitrogen species (RNS) and phytohormones in plant's response to biotic stress

Reactive oxygen species (ROS) and reactive nitrogen species (RNS) are critical for plant development as well as for its stress response. They can function as signaling molecules to orchestrate a well-defined response of plants to biotic stress. These responses are further fine-tuned by phytohormones, such as salicylic acid, jasmonic acid, and ethylene, to modulate immune response. In the past decades, the intricacies of redox and phytohormonal signaling have been uncovered during plant-pathogen interactions. This review explores the dynamic interplay of these components, elucidating their roles in perceiving biotic threats and shaping the plant's defense strategy. Molecular regulators and sites of oxidative burst have been explored during pathogen perception. Further, the interplay between various components of redox and phytohormonal signaling has been explored during bacterial, fungal, viral, and nematode infections as well as during insect pest infestation. Understanding these interactions highlights gaps in the current knowledge and provides insights into engineering crop varieties with enhanced resistance to pathogens and pests. This review also highlights potential applications of manipulating regulators of redox signaling to bolster plant immunity and ensure global food security. Future research should explore regulators of these signaling pathways as potential target to develop biotic stress-tolerant crops. Further insights are also needed into roles of endophytes and host microbiome modulating host ROS and RNS pool for exploiting them as biocontrol agents imparting resistance against pathogens in plants.

Impact of Rhg1 copy number variation on a soybean cyst nematode resistance transcriptional network

Soybean yield loss due to soybean cyst nematode (SCN) infestation has a negative impact on the U.S. economy. Most SCN-resistant soybeans carry a common resistance locus (Rhg1), conferred by copy number variation of a 31.2-kb segment at the Rhg1 locus. To identify the effects of Rhg1 copy number on the plant prior to SCN infection, we investigated genome-wide expression profiles in isogenic Fayette plants carrying different copy numbers at the Rhg1 locus (9-11 copies), that confer different levels of resistance to SCN. We found that even small differences in copy number lead to large changes in expression of downstream defense genes. The co-expression network constructed from differentially expressed genes (DEGs) outside the Rhg1 locus revealed complex effects of Rhg1 copy number on transcriptional regulation involving signal transduction and ethylene-mediated signaling pathways. Moreover, we report a variation in expression levels of phytoalexin biosynthesis-related genes that is correlated with copy number, and the activation of different NBS-LRR gene sets, indicating a broad effect of copy number on defense responses. Using qRT-PCR time series during SCN infection, we validated the SCN responses of DEGs detected in the copy number comparison and showed a stable upregulation of genes related to phytoalexin biosynthesis in resistant Fayette lines during the early stages of the incompatible interaction between soybeans and SCN, before syncytium formation. These results suggest additional genes that could enhance Rhg1-mediated SCN resistance.

CRISPR-Cas9-mediated editing of GmARM improves resistance to multiple stresses in soybean

The growth and development of soybean plants can be affected by both abiotic and biotic stressors, such as saline-alkali stress and Phytophthora root rot. In this study, we identified a stress-related gene-GmARM-whose promoter contained several hormone-response and stress-regulatory elements, including ABRE, TCA element, STRE, and MBS. qRT-PCR analysis showed that the expression of GmARM was the highest in seeds at 55 days after flowering. Furthermore, this gene was upregulated after exposure to saline-alkali stress and Phytophthora root rot infection at the seedling stage. Thus, we generated GmARM mutants using the CRISPR-Cas9 system to understand the role of this gene in stress response. T3 plants showed significantly improved salt tolerance, alkali resistance, and disease resistance, with a significantly higher survival rate than the wildtype plants. Moreover, mutations in GmARM affected the expression of related stress-resistance genes, indicating that GmARM mutants achieved multiple stress tolerance. Therefore, this study provides a foundation for further exploration of the genes involved in resistance to multiple stresses in soybean that can be used for breeding multiple stress-resistance soybean varieties.

Genome scans for selection signatures identify candidate virulence genes for adaptation of the soybean cyst nematode to host resistance

Plant pathogens are constantly under selection pressure for host resistance adaptation. Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean primarily managed through resistant cultivars; however, SCN populations have evolved virulence in response to selection pressures driven by repeated monoculture of the same genetic resistance. Resistance to SCN is mediated by multiple epistatic interactions between Rhg (for resistance to H. glycines) genes. However, the identity of SCN virulence genes that confer the ability to overcome resistance remains unknown. To identify candidate genomic regions showing signatures of selection for increased virulence, we conducted whole genome resequencing of pooled individuals (Pool-Seq) from two pairs of SCN populations adapted on soybeans with Peking-type (rhg1-a, rhg2, and Rhg4) resistance. Population differentiation and principal component analysis-based approaches identified approximately 0.72-0.79 million SNPs, the frequency of which showed potential selection signatures across multiple genomic regions. Chromosomes 3 and 6 between population pairs showed the greatest density of outlier SNPs with high population differentiation. Conducting multiple outlier detection tests to identify overlapping SNPs resulted in a total of 966 significantly differentiated SNPs, of which 285 exon SNPs were mapped to 97 genes. Of these, six genes encoded members of known stylet-secreted effector protein families potentially involved in host defence modulation including venom-allergen-like, annexin, glutathione synthetase, SPRYSEC, chitinase, and CLE effector proteins. Further functional analysis of identified candidate genes will provide new insights into the genetic mechanisms by which SCN overcomes soybean resistance and inform the development of molecular markers for rapidly screening the virulence profile of an SCN-infested field.

Glycine max polygalacturonase inhibiting protein 11 (GmPGIP11) functions in the root to suppress Heterodera glycines parasitism

Pathogen-secreted polygalacturonases (PGs) alter plant cell wall structure by cleaving the α-(1 → 4) linkages between D-galacturonic acid residues in homogalacturonan (HG), macerating the cell wall, facilitating infection. Plant PG inhibiting proteins (PGIPs) disengage pathogen PGs, impairing infection. The soybean cyst nematode, Heterodera glycines, obligate root parasite produces secretions, generating a multinucleate nurse cell called a syncytium, a byproduct of the merged cytoplasm of 200-250 root cells, occurring through cell wall maceration. The common cytoplasmic pool, surrounded by an intact plasma membrane, provides a source from which H. glycines derives nourishment but without killing the parasitized cell during a susceptible reaction. The syncytium is also the site of a naturally-occurring defense response that happens in specific G. max genotypes. Transcriptomic analyses of RNA isolated from the syncytium undergoing the process of defense have identified that one of the 11 G. max PGIPs, GmPGIP11, is expressed during defense. Functional transgenic analyses show roots undergoing GmPGIP11 overexpression (OE) experience an increase in its relative transcript abundance (RTA) as compared to the ribosomal protein 21 (GmRPS21) control, leading to a decrease in H. glycines parasitism as compared to the overexpression control. The GmPGIP11 undergoing RNAi experiences a decrease in its RTA as compared to the GmRPS21 control with transgenic roots experiencing an increase in H. glycines parasitism as compared to the RNAi control. Pathogen associated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) components are shown to influence GmPGIP11 expression while numerous agricultural crops are shown to have homologs.

GmSABP2-1 encodes methyl salicylate esterase and functions in soybean defense against soybean cyst nematode

KEY MESSAGE

The soybean gene GmSABP2-1 encodes methyl salicylate esterase and its overexpression led to significant reduction in development of pathogenic soybean cyst nematode. Soybean cyst nematode (SCN, Heterodera glycines) is one of the most devastating pests of soybean (Glycine max L. Merr.). In searching for SCN-defense genes, a soybean gene of the methylesterase (MES) family was found to be upregulated in an SCN-resistant soybean line and downregulated in an SCN-susceptible line upon SCN infection. This gene was designated as GmSABP2-1. Here, we report on biochemical and overexpression studies of GmSABP2-1 to examine its possible function in SCN resistance. The protein encoded by GmSABP2-1 is closely related to known methyl salicylate esterases. To determine the biochemical function of GmSABP2-1, a full-length cDNA of GmSABP2-1 was cloned into a protein expression vector and expressed in Escherichia coli. The resulting recombinant GmSABP2-1 was demonstrated to catalyze the demethylation of methyl salicylate. The biochemical properties of GmSABP2-1 were determined. Its apparent Km value was 46.2 ± 2.2 μM for methyl salicylate, comparable to those of the known methyl salicylate esterases. To explore the biological significance of GmSABP2-1 in soybean defense against SCN, we first overexpressed GmSABP2-1 in transgenic hairy roots of an SCN-susceptible soybean line. When infected with SCN, GmSABP2-1-overexpressing hairy roots showed 84.5% reduction in the development of SCN beyond J2 stage. To provide further genetic evidence for the role of GmSABP2-1 in SCN resistance, stable transgenic soybean plants overexpressing GmSABP2-1 were produced. Analysis of the GmSABP2-1-overexpressing lines showed a significant reduction in SCN development compared to non-transgenic plants. In conclusion, we demonstrated that GmSABP2-1 encodes methyl salicylate esterase and functions as a resistance-related gene against SCN.

Syncytium Induced by Plant-Parasitic Nematodes

Plant-parasitic nematodes from the genera Globodera, Heterodera (cyst-forming nematodes), and Meloidogyne (root-knot nematodes) are notorious and serious pests of crops. They cause tremendous economic losses between US $80 and 358 billion a year. Nematodes infect the roots of plants and induce the formation of specialised feeding structures (syncytium and giant cells, respectively) that nourish juveniles and adults of the nematodes. The specialised secretory glands enable nematodes to synthesise and secrete effectors that facilitate migration through root tissues and alter the morphogenetic programme of host cells. The formation of feeding sites is associated with the suppression of plant defence responses and deep reprogramming of the development and metabolism of plant cells.In this chapter, we focus on syncytia induced by the sedentary cyst-forming nematodes and provide an overview of ultrastructural changes that occur in the host roots during syncytium formation in conjunction with the most important molecular changes during compatible and incompatible plant responses to infection with nematodes.

Loss-of-function of an α-SNAP gene confers resistance to soybean cyst nematode

Plant-parasitic nematodes are one of the most economically impactful pests in agriculture resulting in billions of dollars in realized annual losses worldwide. Soybean cyst nematode (SCN) is the number one biotic constraint on soybean production making it a priority for the discovery, validation and functional characterization of native plant resistance genes and genetic modes of action that can be deployed to improve soybean yield across the globe. Here, we present the discovery and functional characterization of a soybean resistance gene, GmSNAP02. We use unique bi-parental populations to fine-map the precise genomic location, and a combination of whole genome resequencing and gene fragment PCR amplifications to identify and confirm causal haplotypes. Lastly, we validate our candidate gene using CRISPR-Cas9 genome editing and observe a gain of resistance in edited plants. This demonstrates that the GmSNAP02 gene confers a unique mode of resistance to SCN through loss-of-function mutations that implicate GmSNAP02 as a nematode virulence target. We highlight the immediate impact of utilizing GmSNAP02 as a genome-editing-amenable target to diversify nematode resistance in commercially available cultivars.

Transcriptome Analysis of GmPUB20A Overexpressing and RNA-Interferencing Transgenic Hairy Roots Reveals Underlying Negative Role in Soybean Resistance to Cyst Nematode

Ubiquitination genes are key components of plant responses to biotic stress. GmPUB20A, a ubiquitination gene, plays a negative role in soybean resistance to soybean cyst nematode (SCN). In this study, we employed high-throughput sequencing to investigate transcriptional changes in GmPUB20A overexpressing and RNA-interfering transgenic hairy roots. Totally, 7661 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that DEGs were significantly enriched in disease resistance and signal transduction pathways. In addition, silencing Glyma.15G021600 and Glyma.09G284700 by siRNA, the total number of nematodes was decreased by 33.48% and 27.47% than control plants, respectively. Further, GUS activity and reactive oxygen species (ROS) assays revealed that GmPUB20A, Glyma.15G021600, and Glyma.09G284700 respond to SCN parasitism and interfere with the accumulation of ROS in plant roots, respectively. Collectively, our study provides insights into the molecular mechanism of GmPUB20A in soybean resistance to SCN.

Identification of genes governing resistance to PCN (Globodera rostochiensis) through transcriptome analysis in Solanum tuberosum

Potato cyst nematodes (PCNs) are major pests worldwide that affect potato production. The molecular changes happening in the roots upon PCN infection are still unknown. Identification of transcripts and genes governing PCN resistance will help in the development of resistant varieties. Hence, differential gene expression of compatible (Kufri Jyoti) and incompatible (JEX/A-267) potato genotypes was studied before (0 DAI) and after (10 DAI) inoculation of Globodera rostochiensis J2s through RNA sequencing (RNA-Seq). Total sequencing reads generated ranged between 33 and 37 million per sample, with a read mapping of 48-84% to the potato reference genome. In the infected roots of the resistant genotype JEX/A-267, 516 genes were downregulated, and 566 were upregulated. In comparison, in the susceptible genotype Kufri Jyoti, 316 and 554 genes were downregulated and upregulated, respectively. Genes encoding cell wall proteins, zinc finger protein, WRKY transcription factors, MYB transcription factors, disease resistance proteins, and pathogenesis-related proteins were found to be majorly involved in the incompatible reaction after PCN infection in the resistant genotype, JEX/A-267. Furthermore, RNA-Seq results were validated through quantitative real-time PCR (qRT-PCR), and it was observed that ATP, FLAVO, CYTO, and GP genes were upregulated at 5 DAI, which was subsequently downregulated at 10 DAI. The genes encoding ATP, FLAVO, LBR, and GP were present in > 1.5 fold before infection in JEX-A/267 and upregulated 7.9- to 27.6-fold after 5 DAI; subsequently, most of these genes were downregulated to 0.9- to 2.8-fold, except LBR, which was again upregulated to 44.4-fold at 10 DAI.

Single Nematode Transcriptomic Analysis, Using Long-Read Technology, Reveals Two Novel Virulence Gene Candidates in the Soybean Cyst Nematode, Heterodera glycines

The soybean cyst nematode (Heterodera glycines, SCN), is the most damaging disease of soybean in North America. While management of this pest using resistant soybean is generally still effective, prolonged exposure to cultivars derived from the same source of resistance (PI 88788) has led to the emergence of virulence. Currently, the underlying mechanisms responsible for resistance breakdown remain unknown. In this study, we combined a single nematode transcriptomic profiling approach with long-read sequencing to reannotate the SCN genome. This resulted in the annotation of 1932 novel transcripts and 281 novel gene features. Using a transcript-level quantification approach, we identified eight novel effector candidates overexpressed in PI 88788 virulent nematodes in the late infection stage. Among these were the novel gene Hg-CPZ-1 and a pioneer effector transcript generated through the alternative splicing of the non-effector gene Hetgly21698. While our results demonstrate that alternative splicing in effectors does occur, we found limited evidence of direct involvement in the breakdown of resistance. However, our analysis highlighted a distinct pattern of effector upregulation in response to PI 88788 resistance indicative of a possible adaptation process by SCN to host resistance.

Soybean transporter AAT Rhg1 abundance increases along the nematode migration path and impacts vesiculation and ROS

Rhg1 (Resistance to Heterodera glycines 1) mediates soybean (Glycine max) resistance to soybean cyst nematode (SCN; H. glycines). Rhg1 is a 4-gene, ∼30-kb block that exhibits copy number variation, and the common PI 88788-type rhg1-b haplotype carries 9 to 10 tandem Rhg1 repeats. Glyma.18G022400 (Rhg1-GmAAT), 1 of 3 resistance-conferring genes at the complex Rhg1 locus, encodes the putative amino acid transporter AATRhg1 whose mode of action is largely unknown. We discovered that AATRhg1 protein abundance increases 7- to 15-fold throughout root cells along the migration path of SCN. These root cells develop an increased abundance of vesicles and large vesicle-like bodies (VLB) as well as multivesicular and paramural bodies. AATRhg1 protein is often present in these structures. AATRhg1 abundance remained low in syncytia (plant cells reprogrammed by SCN for feeding), unlike the Rhg1 α-SNAP protein, whose abundance has previously been shown to increase in syncytia. In Nicotiana benthamiana, if soybean AATRhg1 was present, oxidative stress promoted the formation of large VLB, many of which contained AATRhg1. AATRhg1 interacted with the soybean NADPH oxidase GmRBOHG, the ortholog of Arabidopsis thaliana RBOHD previously found to exhibit upregulated expression upon SCN infection. AATRhg1 stimulated reactive oxygen species (ROS) generation when AATRhg1 and GmRBOHG were co-expressed. These findings suggest that AATRhg1 contributes to SCN resistance along the migration path as SCN invades the plant and does so, at least in part, by increasing ROS production. In light of previous findings about α-SNAPRhg1, this study also shows that different Rhg1 resistance proteins function via at least 2 spatially and temporally separate modes of action.

Functional insights of plant bcl-2-associated ahanogene (BAG) proteins: Multi-taskers in diverse cellular signal transduction pathways

Bcl-2-associated athanogene (BAG) gene family is a highly conserved molecular chaperone cofactor in evolution from yeast to humans and plants playing important roles in a variety of signal pathways. Plant BAG proteins have special structures, especially those containing CaM-binding IQ motifs which are unique to plants. While early studies focused more on the structure and physiological function of plant BAGs, recent studies have revealed many novel functional mechanisms involved in multiple cellular processes. How to achieve signal specificity has become an interesting topic of plant BAG research. In this review, we have provided a historic view of plant BAG research and summarized recent advances in the establishment of BAG as essential components in normal plant growth, environmental stress response, and plant immunity. Based on the relationship between BAG proteins and their newly interacting proteins, this review highlights the functional mechanisms of various cellular signals mediated by plant BAGs. Future work needs to focus on the post-translational modification of BAG proteins, and on understanding how specificity is achieved among BAG signaling pathways.

Understandings and future challenges in soybean functional genomics and molecular breeding

Soybean (Glycine max) is a major source of plant protein and oil. Soybean breeding has benefited from advances in functional genomics. In particular, the release of soybean reference genomes has advanced our understanding of soybean adaptation to soil nutrient deficiencies, the molecular mechanism of symbiotic nitrogen (N) fixation, biotic and abiotic stress tolerance, and the roles of flowering time in regional adaptation, plant architecture, and seed yield and quality. Nevertheless, many challenges remain for soybean functional genomics and molecular breeding, mainly related to improving grain yield through high-density planting, maize-soybean intercropping, taking advantage of wild resources, utilization of heterosis, genomic prediction and selection breeding, and precise breeding through genome editing. This review summarizes the current progress in soybean functional genomics and directs future challenges for molecular breeding of soybean.

Weighted gene co-expression network analysis identifies genes related to HG Type 0 resistance and verification of hub gene GmHg1

INTRODUCTION

The soybean cyst nematode (SCN) is a major disease in soybean production thatseriously affects soybean yield. At present, there are no studies on weighted geneco-expression network analysis (WGCNA) related to SCN resistance.

METHODS

Here, transcriptome data from 36 soybean roots under SCN HG Type 0 (race 3) stresswere used in WGCNA to identify significant modules.

RESULTS AND DISCUSSION

A total of 10,000 differentially expressed genes and 21 modules were identified, of which the module most related to SCN was turquoise. In addition, the hub gene GmHg1 with high connectivity was selected, and its function was verified. GmHg1 encodes serine/threonine protein kinase (PK), and the expression of GmHg1 in SCN-resistant cultivars ('Dongnong L-204') and SCN-susceptible cultivars ('Heinong 37') increased significantly after HG Type 0 stress. Soybean plants transformed with GmHg1-OX had significantly increased SCN resistance. In contrast, the GmHg1-RNAi transgenic soybean plants significantly reduced SCN resistance. In transgenic materials, the expression patterns of 11 genes with the same expression trend as the GmHg1 gene in the 'turquoise module' were analyzed. Analysis showed that 11genes were co-expressed with GmHg1, which may be involved in the process of soybean resistance to SCN. Our work provides a new direction for studying the Molecular mechanism of soybean resistance to SCN.

Establishment and maintenance of DNA methylation in nematode feeding sites

A growing body of evidence indicates that epigenetic mechanisms, particularly DNA methylation, play key regulatory roles in plant-nematode interactions. Nevertheless, the transcriptional activity of key genes mediating DNA methylation and active demethylation in the nematode feeding sites remains largely unknown. Here, we profiled the promoter activity of 12 genes involved in maintenance and de novo establishment of DNA methylation and active demethylation in the syncytia and galls induced respectively by the cyst nematode Heterodera schachtii and the root-knot nematode Meloidogyne incognita in Arabidopsis roots. The promoter activity assays revealed that expression of the CG-context methyltransferases is restricted to feeding site formation and development stages. Chromomethylase1 (CMT1), CMT2, and CMT3 and Domains Rearranged Methyltransferase2 (DRM2) and DRM3, which mediate non-CG methylation, showed similar and distinct expression patterns in the syncytia and galls at various time points. Notably, the promoters of various DNA demethylases were more active in galls as compared with the syncytia, particularly during the early stage of infection. Mutants impaired in CG or CHH methylation similarly enhanced plant susceptibility to H. schachtii and M. incognita, whereas mutants impaired in CHG methylation reduced plant susceptibility only to M. incognita. Interestingly, hypermethylated mutants defective in active DNA demethylation exhibited contrasting responses to infection by H. schachtii and M. incognita, a finding most likely associated with differential regulation of defense-related genes in these mutants upon nematode infection. Our results point to methylation-dependent mechanisms regulating plant responses to infection by cyst and root-knot nematodes.

Identification of QTL, QTL-by-environment interactions, and their candidate genes for resistance HG Type 0 and HG Type 1.2.3.5.7 in soybean using 3VmrMLM

Introduction

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is an important disease affecting soybean yield in the world. Potential SCN-related QTLs and QTL-by-environment interactions (QEIs) have been used in SCN-resistant breeding.

Methods

In this study, a compressed variance component mixed model, 3VmrMLM, in genome-wide association studies was used to detect QTLs and QEIs for resistance to SCN HG Type 0 and HG Type 1.2.3.5.7 in 156 different soybean cultivars materials.

Results and discussion

The results showed that 53 QTLs were detected in single environment analysis; 36 QTLs and 9 QEIs were detected in multi-environment analysis. Based on the statistical screening of the obtained QTLs, we obtained 10 novel QTLs and one QEI which were different from the previous studies. Based on previous studies, we identified 101 known genes around the significant/suggested QTLs and QEIs. Furthermore, used the transcriptome data of SCN-resistant (Dongnong L-10) and SCN-susceptible (Suinong 14) cultivars, 10 candidate genes related to SCN resistance were identified and verified by Quantitative real time polymerase chain reaction (qRT-PCR) analysis. Haplotype difference analysis showed that Glyma.03G005600 was associated with SCN HG Type 0 and HG Type 1.2.3.5.7 resistance and had a haplotype beneficial to multi-SCN-race resistance. These results provide a new idea for accelerating SCN disease resistance breeding.

Time Course RNA-seq Reveals Soybean Responses against Root-Lesion Nematode and Resistance Players

Pratylenchus brachyurus causes serious damage to soybean production and other crops worldwide. Plant molecular responses to RLN infection remain largely unknown and no resistance genes have been identified in soybean. In this study, we analyzed molecular responses to RLN infection in moderately resistant BRSGO (Chapadões-BRS) and susceptible TMG115 RR (TMG) Glycine max genotypes. Differential expression analysis revealed two stages of response to RLN infection and a set of differentially expressed genes (DEGs) in the first stage suggested a pattern-triggered immunity (PTI) in both genotypes. The divergent time-point of DEGs between genotypes was observed four days post-infection, which included the activation of mitogen-activated protein kinase (MAPK) and plant-pathogen interaction genes in the BRS, suggesting the occurrence of an effector-triggered immunity response (ETI) in BRS. The co-expression analyses combined with single nucleotide polymorphism (SNP) uncovered a key element, a transcription factor phytochrome-interacting factor (PIF7) that is a potential regulator of moderate resistance to RLN infection. Two genes for resistance-related leucine-rich repeat (LRR) proteins were found as BRS-specific expressed genes. In addition, alternative splicing analysis revealed an intron retention in a myo-inositol oxygenase (MIOX) transcript, a gene related to susceptibility, may cause a loss of function in BRS.

Genome-wide identification, characterization and expression of HSP 20 gene family in dove

Davidia involucrata is a significant living fossil with high abiotic stress tolerance. Although heat shock protein 20 (HSP20) has already been linked to heat stress, nothing is known about HSP20 family protein activities in D. involucrata. The functional dynamics of the D. involucrata HSP20 (DiHSP20) gene family were identified and characterized using a thorough genome-wide investigation. From the genome of D. involucrata, a total of 42 HSP20 genes were identified, which are distributed across 16 chromosomes. The DiHSP20 proteins were grouped into seven separate subfamilies by our phylogenetic analysis, which was validated by the conserved motif composition and gene structure studies. Segmental duplication events were shown to play a crucial role in the expansion of the DiHSP20 gene family. Synteny analysis revealed that 19 DiHSP20 genes of D. involucrata shared a syntenic connection with Arabidopsis genes, 39 with C. acuminata genes, and just 6 with O. sativa genes. Additionally, heat stress differently enhanced the expression levels of D. involucrata HSP20 genes. After 1 hour of heat treatment, the expression levels of most DiHSP20 genes, particularly DiHSP20-7, DiHSP20-29, DiHSP20-30, DiHSP20-32, and DiHSP20-34, were dramatically increased, suggestted that they might be employed as heat tolerance candidate genes. Overall, these findings add to our knowledge of the HSP20 family genes and provide helpful information for breeding heat stress resistance in D. involucrata.

Proteomic, Transcriptomic, Mutational, and Functional Assays Reveal the Involvement of Both THF and PLP Sites at the GmSHMT08 in Resistance to Soybean Cyst Nematode

The serine hydroxymethyltransferase (SHMT; E.C. 2.1.2.1) is involved in the interconversion of serine/glycine and tetrahydrofolate (THF)/5,10-methylene THF, playing a key role in one-carbon metabolism, the de novo purine pathway, cellular methylation reactions, redox homeostasis maintenance, and methionine and thymidylate synthesis. GmSHMT08 is the soybean gene underlying soybean cyst nematode (SCN) resistance at the Rhg4 locus. GmSHMT08 protein contains four tetrahydrofolate (THF) cofactor binding sites (L129, L135, F284, N374) and six pyridoxal phosphate (PLP) cofactor binding/catalysis sites (Y59, G106, G107, H134, S190A, H218). In the current study, proteomic analysis of a data set of protein complex immunoprecipitated using GmSHMT08 antibodies under SCN infected soybean roots reveals the presence of enriched pathways that mainly use glycine/serine as a substrate (glyoxylate cycle, redox homeostasis, glycolysis, and heme biosynthesis). Root and leaf transcriptomic analysis of differentially expressed genes under SCN infection supported the proteomic data, pointing directly to the involvement of the interconversion reaction carried out by the serine hydroxymethyltransferase enzyme. Direct site mutagenesis revealed that all mutated THF and PLP sites at the GmSHMT08 resulted in increased SCN resistance. We have shown the involvement of PLP sites in SCN resistance. Specially, the effect of the two Y59 and S190 PLP sites was more drastic than the tested THF sites. This unprecedented finding will help us to identify the biological outcomes of THF and PLP residues at the GmSHMT08 and to understand SCN resistance mechanisms.

Functional analysis of soybean cyst nematode-inducible synthetic promoters and their regulation by biotic and abiotic stimuli in transgenic soybean (Glycine max)

We previously identified cis-regulatory motifs in the soybean (Glycine max) genome during interaction between soybean and soybean cyst nematode (SCN), Heterodera glycines. The regulatory motifs were used to develop synthetic promoters, and their inducibility in response to SCN infection was shown in transgenic soybean hairy roots. Here, we studied the functionality of two SCN-inducible synthetic promoters; 4 × M1.1 (TAAAATAAAGTTCTTTAATT) and 4 × M2.3 (ATATAATTAAGT) each fused to the -46 CaMV35S core sequence in transgenic soybean. Histochemical GUS analyses of transgenic soybean plants containing the individual synthetic promoter::GUS construct revealed that under unstressed condition, no GUS activity is present in leaves and roots. While upon nematode infection, the synthetic promoters direct GUS expression to roots predominantly in the nematode feeding structures induced by the SCN and by the root-knot nematode (RKN), Meloidogyne incognita. There were no differences in GUS activity in leaves between nematode-infected and non-infected plants. Furthermore, we examined the specificity of the synthetic promoters in response to various biotic (insect: fall armyworm, Spodoptera frugiperda; and bacteria: Pseudomonas syringe pv. glycinea, P. syringe pv. tomato, and P. marginalis) stresses. Additionally, we examined the specificity to various abiotic (dehydration, salt, cold, wounding) as well as to the signal molecules salicylic acid (SA), methyl jasmonate (MeJA), and abscisic acid (ABA) in the transgenic plants. Our wide-range analyses provide insights into the potential applications of synthetic promoter engineering for conditional expression of transgenes leading to transgenic crop development for resistance improvement in plant.

Exploring Soybean Resistance to Soybean Cyst Nematode

Resistance to the soybean cyst nematode (SCN) is a topic incorporating multiple mechanisms and multiple types of science. It is also a topic of substantial agricultural importance, as SCN is estimated to cause more yield damage than any other pathogen of soybean, one of the world's main food crops. Both soybean and SCN have experienced jumps in experimental tractability in the past decade, and significant advances have been made. The rhg1-b locus, deployed on millions of farm acres, has been durable and will remain important, but local SCN populations are gradually evolving to overcome rhg1-b. Multiple other SCN resistance quantitative trait loci (QTL) of proven value are now in play with soybean breeders. QTL causal gene discovery and mechanistic insights into SCN resistance are contributing to both basic and applied disciplines. Additional understanding of SCN and other cyst nematodes will also grow in importance and lead to novel disease control strategies.

Essential Roles of Cupredoxin Family Proteins in Soybean Cyst Nematode Resistance

Soybean cyst nematode (SCN, Heterodera glycines), one of the most devastating soybean pathogens, causes a significant yield loss in soybean production. One of the most effective ways to manage SCN is to grow resistant cultivars. Therefore, comparative study using resistant and susceptible soybean cultivars provides a powerful tool to identify new genes involved in soybean SCN resistance. In the present study, a transcriptome analysis was carried out using both the resistant (PI88788) and susceptible (Williams 82) soybean cultivars to characterize the responses to nematode infection. Various defense-related genes and different pathways involved in nematode resistance were recognized as being highly expressed in resistant cultivar. Promoter-GUS analysis was conducted to monitor the spatial expression pattern of the genes highly induced by nematode infection. Two nematode-inducible promoters for Glyma.05g147000 (encoding caffeoyl-CoA O-methyltransferase) and Glyma.06g036700 (encoding cupredoxin superfamily protein) were characterized, and the promoters could efficiently drive the expression of known nematode resistance genes (α-SNAPRhg1HC or GmSHMT) to affect soybean SCN resistance. Interestingly, expression of the cupredoxin family genes was upregulated not only by SCN, but also by jasmonic acid treatment. DNA sequence analysis identified that a conserved motif (GGTGCATG) with high similarity to SCNbox1 and GC-rich element is enriched in their promoter regions, suggesting its potential to serve as a nematode-responsive regulatory element. Overexpression of Glyma.06g036700 significantly enhanced soybean resistance to cyst nematode. Overall, our findings not only highlight the essential role of cupredoxin family genes in SCN resistance, but also offer potential functional tools to develop nematode resistance in crops.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

RNA-Seq of Cyst Nematode Infestation of Potato (Solanum tuberosum L.): A Comparative Transcriptome Analysis of Resistant and Susceptible Cultivars

Potato (Solanum tuberosum L.) is an important food crop worldwide, and potato cyst nematodes (PCNs) are among the most serious pests. The identification of disease resistance genes and molecular markers for PCN infestation can aid in crop improvement research programs against PCN infestation. In the present study, we used high-throughput RNA sequencing to investigate the comprehensive resistance mechanisms induced by PCN infestation in the resistant cultivar Kufri Swarna and the susceptible cultivar Kufri Jyoti. PCN infestation induced 791 differentially expressed genes in resistant cultivar Kufri Swarna, comprising 438 upregulated and 353 downregulated genes. In susceptible cultivar Kufri Jyoti, 2225 differentially expressed genes were induced, comprising 1247 upregulated and 978 downregulated genes. We identified several disease resistance genes (KIN) and transcription factors (WRKY, HMG, and MYB) that were upregulated in resistant Kufri Swarna. The differentially expressed genes from several enriched KEGG pathways, including MAPK signaling, contributed to the disease resistance in Kufri Swarna. Functional network analysis showed that several cell wall biogenesis genes were induced in Kufri Swarna in response to infestation. This is the first study to identify underlying resistance mechanisms against PCN and host interaction in Indian potato varieties.

Kinase-dead mutation: A novel strategy for improving soybean resistance to soybean cyst nematode Heterodera glycines

Protein kinases phosphorylate proteins for functional changes and are involved in nearly all cellular processes, thereby regulating almost all aspects of plant growth and development, and responses to biotic and abiotic stresses. We generated two independent co-expression networks of soybean genes using control and stress response gene expression data and identified 392 differentially highly interconnected kinase hub genes among the two networks. Of these 392 kinases, 90 genes were identified as "syncytium highly connected hubs", potentially essential for activating kinase signalling pathways in the nematode feeding site. Overexpression of wild-type coding sequences of five syncytium highly connected kinase hub genes using transgenic soybean hairy roots enhanced plant susceptibility to soybean cyst nematode (SCN; Heterodera glycines) Hg Type 0 (race 3). In contrast, overexpression of kinase-dead variants of these five syncytium kinase hub genes significantly enhanced soybean resistance to SCN. Additionally, three of the five tested kinase hub genes enhanced soybean resistance to SCN Hg Type 1.2.5.7 (race 2), highlighting the potential of the kinase-dead approach to generate effective and durable resistance against a wide range of SCN Hg types. Subcellular localization analysis revealed that kinase-dead mutations do not alter protein cellular localization, confirming the structure-function of the kinase-inactive variants in producing loss-of-function phenotypes causing significant decrease in nematode susceptibility. Because many protein kinases are highly conserved and are involved in plant responses to various biotic and abiotic stresses, our approach of identifying kinase hub genes and their inactivation using kinase-dead mutation could be translated for biotic and abiotic stress tolerance.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Pumpkin (Cucurbita moschata) HSP20 Gene Family Identification and Expression Under Heat Stress

Pumpkin (Cucurbita moschata) is an important cucurbit vegetable crop that has strong resistance to abiotic stress. While heat shock protein 20 (HSP20) has been implicated in vegetable response to heat stress, little is known regarding activity of HSP20 family proteins in C. moschata. Here, we performed a comprehensive genome-wide analysis to identify and characterize the functional dynamics of the Cucurbita moschata HSP20 (CmoHSP20) gene family. A total of 33 HSP20 genes distributed across 13 chromosomes were identified from the pumpkin genome. Our phylogenetic analysis determined that the CmoHSP20 proteins fell into nine distinct subfamilies, a division supported by the conserved motif composition and gene structure analyses. Segmental duplication events were shown to play a key role in expansion of the CmoHSP20 gene family. Synteny analysis revealed that 19 and 18 CmoHSP20 genes were collinear with those in the cucumber and melon genomes, respectively. Furthermore, the expression levels of pumpkin HSP20 genes were differentially induced by heat stress. The transcript level of CmoHSP20-16, 24 and 25 were down-regulated by heat stress, while CmoHSP20-7, 13, 18, 22, 26 and 32 were up-regulated by heat stress, which could be used as heat tolerance candidate genes. Overall, these findings contribute to our understanding of vegetable HSP20 family genes and provide valuable information that can be used to breed heat stress resistance in cucurbit vegetable crops.

Combining targeted metabolite analyses and transcriptomics to reveal the specific chemical composition and associated genes in the incompatible soybean variety PI437654 infected with soybean cyst nematode HG1.2.3.5.7

BACKGROUND

Soybean cyst nematode, Heterodera glycines, is one of the most devastating pathogens of soybean and causes severe annual yield losses worldwide. Different soybean varieties exhibit different responses to H. glycines infection at various levels, such as the genomic, transcriptional, proteomic and metabolomic levels. However, there have not yet been any reports of the differential responses of incompatible and compatible soybean varieties infected with H. glycines based on combined metabolomic and transcriptomic analyses.

RESULTS

In this study, the incompatible soybean variety PI437654 and three compatible soybean varieties, Williams 82, Zhonghuang 13 and Hefeng 47, were used to clarify the differences in metabolites and transcriptomics before and after the infection with HG1.2.3.5.7. A local metabolite-calibrated database was used to identify potentially differential metabolites, and the differences in metabolites and metabolic pathways were compared between the incompatible and compatible soybean varieties after inoculation with HG1.2.3.5.7. In total, 37 differential metabolites and 20 KEGG metabolic pathways were identified, which were divided into three categories: metabolites/pathways overlapped in the incompatible and compatible soybeans, and metabolites/pathways specific to either the incompatible or compatible soybean varieties. Twelve differential metabolites were found to be involved in predicted KEGG metabolite pathways. Moreover, 14 specific differential metabolites (such as significantly up-regulated nicotine and down-regulated D-aspartic acid) and their associated KEGG pathways (such as the tropane, piperidine and pyridine alkaloid biosynthesis, alanine, aspartate and glutamate metabolism, sphingolipid metabolism and arginine biosynthesis) were significantly altered and abundantly enriched in the incompatible soybean variety PI437654, and likely played pivotal roles in defending against HG1.2.3.5.7 infection. Three key metabolites (N-acetyltranexamic acid, nicotine and D,L-tryptophan) found to be significantly up-regulated in the incompatible soybean variety PI437654 infected by HG1.2.3.5.7 were classified into two types and used for combined analyses with the transcriptomic expression profiling. Associated genes were predicted, along with the likely corresponding biological processes, cellular components, molecular functions and pathways.

CONCLUSIONS

Our results not only identified potential novel metabolites and associated genes involved in the incompatible response of PI437654 to soybean cyst nematode HG1.2.3.5.7, but also provided new insights into the interactions between soybeans and soybean cyst nematodes.

Chromosome-scale genome assembly of Cucumis hystrix-a wild species interspecifically cross-compatible with cultivated cucumber

Cucumis hystrix Chakr. (2n = 2x = 24) is a wild species that can hybridize with cultivated cucumber (C. sativus L., 2n = 2x = 14), a globally important vegetable crop. However, cucumber breeding is hindered by its narrow genetic base. Therefore, introgression from C. hystrix has been anticipated to bring a breakthrough in cucumber improvement. Here, we report the chromosome-scale assembly of C. hystrix genome (289 Mb). Scaffold N50 reached 14.1 Mb. Over 90% of the sequences were anchored onto 12 chromosomes. A total of 23,864 genes were annotated using a hybrid method. Further, we conducted a comprehensive comparative genomic analysis of cucumber, C. hystrix, and melon (C. melo L., 2n = 2x = 24). Whole-genome comparisons revealed that C. hystrix is phylogenetically closer to cucumber than to melon, providing a molecular basis for the success of its hybridization with cucumber. Moreover, expanded gene families of C. hystrix were significantly enriched in "defense response," and C. hystrix harbored 104 nucleotide-binding site-encoding disease resistance gene analogs. Furthermore, 121 genes were positively selected, and 12 (9.9%) of these were involved in responses to biotic stimuli, which might explain the high disease resistance of C. hystrix. The alignment of whole C. hystrix genome with cucumber genome and self-alignment revealed 45,417 chromosome-specific sequences evenly distributed on C. hystrix chromosomes. Finally, we developed four cucumber-C. hystrix alien addition lines and identified the exact introgressed chromosome using molecular and cytological methods. The assembled C. hystrix genome can serve as a valuable resource for studies on Cucumis evolution and interspecific introgression breeding of cucumber.

Screening soybean cyst nematode effectors for their ability to suppress plant immunity

The soybean cyst nematode (SCN), Heterodera glycines, is one of the most destructive pathogens of soybeans. SCN is an obligate and sedentary parasite that transforms host plant root cells into an elaborate permanent feeding site, a syncytium. Formation and maintenance of a viable syncytium is an absolute requirement for nematode growth and reproduction. In turn, sensing pathogen attack, plants activate defence responses and may trigger programmed cell death at the sites of infection. For successful parasitism, H. glycines must suppress these host defence responses to establish and maintain viable syncytia. Similar to other pathogens, H. glycines engages in these molecular interactions with its host via effector proteins. The goal of this study was to conduct a comprehensive screen to identify H. glycines effectors that interfere with plant immune responses. We used Nicotiana benthamiana plants infected by Pseudomonas syringae and Pseudomonas fluorescens strains. Using these pathosystems, we screened 51 H. glycines effectors to identify candidates that could inhibit effector-triggered immunity (ETI) and/or pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI). We identified three effectors as ETI suppressors and seven effectors as PTI suppressors. We also assessed expression modulation of plant immune marker genes as a function of these suppressors.

Targeted suppression of soybean BAG6-induced cell death in yeast by soybean cyst nematode effectors

While numerous effectors that suppress plant immunity have been identified from bacteria, fungi, and oomycete pathogens, relatively little is known for nematode effectors. Several dozen effectors have been reported from the soybean cyst nematode (SCN). Previous studies suggest that a hypersensitive response-like programmed cell death is triggered at nematode feeding sites in soybean during an incompatible interaction. However, virulent SCN populations overcome this incompatibility using unknown mechanisms. A soybean BAG6 (Bcl-2 associated anthanogene 6) gene previously reported by us to be highly up-regulated in degenerating feeding sites induced by SCN in a resistant soybean line was attenuated in response to a virulent SCN population. We show that GmBAG6-1 induces cell death in yeast like its Arabidopsis homolog AtBAG6 and also in soybean. This led us to hypothesize that virulent SCN may target GmBAG6-1 as part of their strategy to overcome soybean defence responses during infection. Thus, we used a yeast viability assay to screen SCN effector candidates for their ability to specifically suppress GmBAG6-1-induced cell death. We identified several effectors that strongly suppressed cell death mediated by GmBAG6-1. Two effectors identified as suppressors showed direct interaction with GmBAG6-1 in yeast, suggesting that one mechanism of cell death suppression may occur through an interaction with this host protein.

A pathogenesis-related protein GmPR08-Bet VI promotes a molecular interaction between the GmSHMT08 and GmSNAP18 in resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is the most devastating pest affecting soybean production worldwide. SCN resistance requires both the GmSHMT08 and the GmSNAP18 in 'Peking'-type resistance. Here, we describe the molecular interaction between GmSHMT08 and GmSNAP18, which is potentiated by a pathogenesis-related protein GmPR08-Bet VI. Like GmSNAP18 and GmSHMT08, GmPR08-Bet VI expression was induced in response to SCN and its overexpression decreased SCN cysts by 65% in infected transgenic soybean roots. Overexpression of GmPR08-Bet VI did not have an effect on SCN resistance when the two cytokinin-binding sites in GmPR08-Bet VI were mutated, indicating a new role of GmPR08-Bet VI in SCN resistance. GmPR08-Bet VI was mapped to a QTL for resistance to SCN using different mapping populations. GmSHMT08, GmSNAP18 and GmPR08-Bet VI localize to the cytosol and plasma membrane. GmSNAP18 expression and localization hyper-accumulated at the plasma membrane and was specific to the root cells surrounding the nematode in SCN-resistant soybeans. Genes encoding key components of the salicylic acid signalling pathway were induced under SCN infection. GmSNAP18 and GmPR08-Bet VI were also induced under salicylic acid and cytokinin exogenous treatments, while GmSHMT08 was induced only when the resistant GmSNAP18 was present, pointing to the presence of a molecular crosstalk between SCN-resistant genes and defence genes. Expression analysis of GmSHMT08 and GmSNAP18 identified the need of a minimum expression requirement to trigger the SCN resistance reaction. These results provide insight into a new response mechanism towards plant nematode resistance involving haplotype compatibility, gene dosage and hormone signalling.

Identification of introduced and stably inherited DNA methylation variants in soybean associated with soybean cyst nematode parasitism

DNA methylation is a widespread epigenetic mark that contributes to transcriptome reprogramming during plant-pathogen interactions. However, the distinct role of DNA methylation in establishing resistant and susceptible responses remains largely unexplored. Here, we developed and used a pair of near-isogenic lines (NILs) to characterize DNA methylome landscapes of soybean roots during the susceptible and resistant interactions with soybean cyst nematode (SCN; Heterodera glycines). We also compared the methylomes of the NILs and their parents to identify introduced and stably inherited methylation variants. The genomes of the NILs were substantially differentially methylated under uninfected conditions. This difference was associated with differential gene expression that may prime the NIL responses to SCN infection. In response to SCN infection, the susceptible line exhibited reduced global methylation levels in both protein-coding genes and transposable elements, whereas the resistant line showed the opposite response, increased global methylation levels. Heritable and novel nonparental differentially methylated regions overlapping with genes associated with soybean response to SCN infection were identified and validated using transgenic hairy root system. Our analyses indicate that DNA methylation patterns associated with the susceptible and resistant interactions are highly specific and that novel and stably inherited methylation variants are of biological significance.

Mutations at the Serine Hydroxymethyltransferase Impact its Interaction with a Soluble NSF Attachment Protein and a Pathogenesis-Related Protein in Soybean

Resistance to soybean cyst nematodes (SCN) in “Peking-type” resistance is bigenic, requiring Rhg4-a and rhg1-a. Rhg4-a encodes a serine hydroxymethyltransferase (GmSHMT08) and rhg1-a encodes a soluble NSF attachment protein (GmSNAP18). Recently, it has been shown that a pathogenesis-related protein, GmPR08-Bet VI, potentiates the interaction between GmSHMT08 and GmSNAP18. Mutational analysis using spontaneously occurring and ethyl methanesulfonate (EMS)-induced mutations was carried out to increase our knowledge of the interacting GmSHMT08/GmSNAP18/GmPR08-Bet VI multi-protein complex. Mutations affecting the GmSHMT08 protein structure (dimerization and tetramerization) and interaction sites with GmSNAP18 and GmPR08-Bet VI proteins were found to impact the multi-protein complex. Interestingly, mutations affecting the PLP/THF substrate binding and catalysis did not affect the multi-protein complex, although they resulted in increased susceptibility to SCN. Most importantly, GmSHMT08 and GmSNAP18 from PI88788 were shown to interact within the cell, being potentiated in the presence of GmPR08-Bet VI. In addition, we have shown the presence of incompatibility between the GmSNAP18 (rhg1-b) of PI88788 and GmSHMT08 (Rhg4-a) from Peking. Components of the reactive oxygen species (ROS) pathway were shown to be induced in the SCN incompatible reaction and were mapped to QTLs for resistance to SCN using different mapping populations.

Transcriptomic profiling of Solanum peruvianum LA3858 revealed a Mi-3-mediated hypersensitive response to Meloidogyne incognita

Background

The Mi-1 gene was the first identified and cloned gene that provides resistance to root-knot nematodes (RKNs) in cultivated tomato. However, owing to its temperature sensitivity, this gene does not meet the need for breeding disease-resistant plants that grow under high temperature. In this study, Mi-3 was isolated from the wild species PI 126443 (LA3858) and was shown to display heat-stable resistance to RKNs. However, the mechanism that regulates this resistance remains unknown.

Results

In this study, 4760, 1024 and 137 differentially expressed genes (DEGs) were enriched on the basis of pairwise comparisons (34 °C vs. 25 °C) at 0 (before inoculation), 3 and 6 days post-inoculation (dpi), respectively. A total of 7035 DEGs were identified from line LA3858 in the respective groups under the different soil temperature treatments. At 3 dpi, most DEGs were enriched in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to plant biotic responses, such as “plant-pathogen interaction” and “plant hormone signal transduction”. Significantly enriched DEGs were found to encode key proteins such as R proteins and heat-shock proteins (HSPs). Moreover, other DEGs were found to participate in Ca²⁺ signal transduction; the production of ROS; DEGs encoding transcription factors (TFs) from the bHLH, TGA, ERF, heat-shock transcription factor (HSF) and WRKY families were highly expressed, which contribute to be involved into the formation of phytohormones, such as salicylic acid (SA), jasmonic acid (JA) and ethylene (ET), the expression of most was upregulated at 3 dpi at the 25 °C soil temperature compared with the 34 °C soil temperature.

Conclusion

Taken together, the results of our study revealed reliable candidate genes from wild materials LA3858, that are related to Mi-3-mediate resistance to Meloidogyne incognita. A large number of vital pathways and DEGs were expressed specifically in accession LA3858 grown at 34 °C and 25 °C soil temperatures at 3 dpi. Upon infection by RKNs, pattern-recognition receptors (PRRs) specifically recognized conserved pathogen-associated molecular patterns (PAMPs) as a result of pathogen-triggered immunity (PTI), and the downstream defensive signal transduction pathway was likely activated through Ca²⁺ signal channels. The expression of various TFs was induced to synthesize phytohormones and activate R proteins related to resistance, resulting in the development of effector-triggered immunity (ETI). Last, a hypersensitive response in the roots occurred, which was probably induced by the accumulation of ROS.

Profiling of defense responsive pathway regulatory genes in Asian rice (Oryza sativa) against infection of Meloidogyne graminicola (Nematoda:Meloidogynidae)

In the present study, 36 Asian rice cultivars/landraces were evaluated against M. graminicola under in vitro conditions using soilless Pluronic gel medium. The cultivars/genotypes Phule Radha, EK 70, LK 248 and Khalibagh showed significantly reduced nematode infection, endoparasitic development, and derived multiplication factor indicating the presence of resistance, while Halvi Sal 17 was found to be most susceptible. Performance of selected genotypes showing resistance/susceptibility under in vitro conditions was further confirmed in soil which also revealed Phule Radha to be highly resistant and Halvi Sal 17 as the most susceptible genotype. Further, expression profile of plant defense responsive genes related to MAPK pathway, phytohormones, PR-proteins and callose and lignin synthesis were quantified in Phule Radha (the most resistant) and Halvi Sal 17 (the most susceptible) at 2 and 6 days post nematode inoculation. Significant upregulated expression of several defensive genes was observed in the resistant cultivar Phule Radha in contrast to insignificant expression in the susceptible varieties. The resistant genotype identified in the present study will be highly promising for resistance breeding in rice against M. graminicola.

Metabolomics as an Emerging Tool for the Study of Plant-Pathogen Interactions

Plants defend themselves from most microbial attacks via mechanisms including cell wall fortification, production of antimicrobial compounds, and generation of reactive oxygen species. Successful pathogens overcome these host defenses, as well as obtain nutrients from the host. Perturbations of plant metabolism play a central role in determining the outcome of attempted infections. Metabolomic analyses, for example between healthy, newly infected and diseased or resistant plants, have the potential to reveal perturbations to signaling or output pathways with key roles in determining the outcome of a plant-microbe interaction. However, application of this -omic and its tools in plant pathology studies is lagging relative to genomic and transcriptomic methods. Thus, it is imperative to bring the power of metabolomics to bear on the study of plant resistance/susceptibility. This review discusses metabolomics studies that link changes in primary or specialized metabolism to the defense responses of plants against bacterial, fungal, nematode, and viral pathogens. Also examined are cases where metabolomics unveils virulence mechanisms used by pathogens. Finally, how integrating metabolomics with other -omics can advance plant pathology research is discussed.

Early transcriptional responses to soybean cyst nematode HG Type 0 show genetic differences among resistant and susceptible soybeans

KEY MESSAGE

Root transcriptome profiling of three soybean cultivars and a wild relative infected with soybean cyst nematode at migratory phase revealed differential resistance pathway responses between resistant and susceptible genotypes. The soybean cyst nematode (SCN), Heterodera glycines, is the most serious pathogen of soybean production throughout the world. Using resistant cultivars is the primary management strategy against SCN infestation. To gain insight into the still obscure mechanisms of genetic resistance to nematodes in different soybean genotypes, RNA-Seq profiling of the roots of Glycine max cv. Peking, Fayette, Williams 82, and a wild relative (Glycine soja PI 468916) was performed during SCN infection at the migratory phase. The analysis showed statistically significant changes of expression beginning at eight hours after inoculation in genes associated with defense mechanisms and pathways, such as the phenylpropanoid biosynthesis pathway, plant innate immunity and hormone signaling. Our results indicate the importance of the early plant response to migratory phase nematodes in pathogenicity determination. The transcriptome changes occurring during early SCN infection included a number of genes and pathways specific to the different resistant genotypes. We observed the most extensive resistant transcriptome reaction in PI 468916, where the resistant response was qualitatively different from that of commonly used G. max varieties.

Genomics of Plant Disease Resistance in Legumes

The constant interactions between plants and pathogens in the environment and the resulting outcomes are of significant importance for agriculture and agricultural scientists. Disease resistance genes in plant cultivars can break down in the field due to the evolution of pathogens under high selection pressure. Thus, the protection of crop plants against pathogens is a continuous arms race. Like any other type of crop plant, legumes are susceptible to many pathogens. The dawn of the genomic era, in which high-throughput and cost-effective genomic tools have become available, has revolutionized our understanding of the complex interactions between legumes and pathogens. Genomic tools have enabled a global view of transcriptome changes during these interactions, from which several key players in both the resistant and susceptible interactions have been identified. This review summarizes some of the large-scale genomic studies that have clarified the host transcriptional changes during interactions between legumes and their plant pathogens while highlighting some of the molecular breeding tools that are available to introgress the traits into breeding programs. These studies provide valuable insights into the molecular basis of different levels of host defenses in resistant and susceptible interactions.

Heat Shock Proteins: Dynamic Biomolecules to Counter Plant Biotic and Abiotic Stresses

Due to the present scenario of climate change, plants have to evolve strategies to survive and perform under a plethora of biotic and abiotic stresses, which restrict plant productivity. Maintenance of plant protein functional conformation and preventing non-native proteins from aggregation, which leads to metabolic disruption, are of prime importance. Plant heat shock proteins (HSPs), as chaperones, play a pivotal role in conferring biotic and abiotic stress tolerance. Moreover, HSP also enhances membrane stability and detoxifies the reactive oxygen species (ROS) by positively regulating the antioxidant enzymes system. Additionally, it uses ROS as a signal to molecules to induce HSP production. HSP also enhances plant immunity by the accumulation and stability of pathogenesis-related (PR) proteins under various biotic stresses. Thus, to unravel the entire plant defense system, the role of HSPs are discussed with a special focus on plant response to biotic and abiotic stresses, which will be helpful in the development of stress tolerance in plant crops.

Re-targeting of a plant defense protease by a cyst nematode effector

Plants mount defense responses during pathogen attacks, and robust host defense suppression by pathogen effector proteins is essential for infection success. 4E02 is an effector of the sugar beet cyst nematode Heterodera schachtii. Arabidopsis thaliana lines expressing the effector-coding sequence showed altered expression levels of defense response genes, as well as higher susceptibility to both the biotroph H. schachtii and the necrotroph Botrytis cinerea, indicating a potential suppression of defenses by 4E02. Yeast two-hybrid analyses showed that 4E02 targets A. thaliana vacuolar papain-like cysteine protease (PLCP) 'Responsive to Dehydration 21A' (RD21A), which has been shown to function in the plant defense response. Activity-based protein profiling analyses documented that the in planta presence of 4E02 does not impede enzymatic activity of RD21A. Instead, 4E02 mediates a re-localization of this protease from the vacuole to the nucleus and cytoplasm, which is likely to prevent the protease from performing its defense function and at the same time, brings it in contact with novel substrates. Yeast two-hybrid analyses showed that RD21A interacts with multiple host proteins including enzymes involved in defense responses as well as carbohydrate metabolism. In support of a role in carbohydrate metabolism of RD21A after its effector-mediated re-localization, we observed cell wall compositional changes in 4E02 expressing A. thaliana lines. Collectively, our study shows that 4E02 removes RD21A from its defense-inducing pathway and repurposes this enzyme by targeting the active protease to different cell compartments.

Soybean root transcriptome profiling reveals a nonhost resistant response during Heterodera glycines infection

Heterodera glycines (soybean cyst nematode, SCN) is one of the most devastating pathogens of soybean worldwide. The compatible and in compatible interactions between soybean and SCN have well documented. Nevertheless, the molecular mechanism of a nonhost resistant response in soybean against SCN infection remains obscure. Toward this end, a global transcriptional comparison was conducted between susceptible and resistant reactions of soybean roots infected by taking advantage of finding a new pathotype of SCN (SCNT). The soybean cultivar Lee, which exhibits resistant to SCNT and susceptible to HG 1.2.3.4.7 (SCNs) was utilized in the expriments. The results highlighted a nonhost resistant response of soybean. Transcriptome analysis indicated that the number of differentially expressed genes (DEGs) in the resistant interaction (3746) was much larger than that in the susceptible interaction (602). A great number of genes acting as intrinsic component of membrane, integral component of membrane, cell periphery and plasma membrance were remarkably enriched only in the resistant interaction, while the taurine and hypotaurine, phenylpropanoid pathway, plant-pathogen interaction and transcript factors were modulated in both interactions. This is the first study to examine genes expression patterns in a soybean genotype in response to invasion by a virulent and avirulent SCN population at the transcriptional level, which will provide insights into the complicate molecular mechanism of the nonhost resistant interaction.