The expression of a naturally occurring, truncated allele of an α-SNAP gene suppresses plant parasitic nematode infection

Fortifying nematode resistance through susceptibility gene inactivation

The predominant genetic defense mechanism against soybean cyst nematode (SCN) in 95% of the North America market is under threat by virulent SCN populations. Usovsky et al. identified GmSNAP02 as an SCN susceptibility gene through fine-mapping of unique bi-parental populations. Loss-of-function of GmSNAP02 confers enhanced resistance to more virulent SCN.

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.

Cataloging SCN resistance loci in North American public soybean breeding programs

Soybean cyst nematode (SCN) is a destructive pathogen of soybeans responsible for annual yield loss exceeding $1.5 billion in the United States. Here, we conducted a series of genome-wide association studies (GWASs) to understand the genetic landscape of SCN resistance in the University of Missouri soybean breeding programs (Missouri panel), as well as germplasm and cultivars within the United States Department of Agriculture (USDA) Uniform Soybean Tests-Northern Region (NUST). For the Missouri panel, we evaluated the resistance of breeding lines to SCN populations HG 2.5.7 (Race 1), HG 1.2.5.7 (Race 2), HG 0 (Race 3), HG 2.5.7 (Race 5), and HG 1.3.6.7 (Race 14) and identified seven quantitative trait nucleotides (QTNs) associated with SCN resistance on chromosomes 2, 8, 11, 14, 17, and 18. Additionally, we evaluated breeding lines in the NUST panel for resistance to SCN populations HG 2.5.7 (Race 1) and HG 0 (Race 3), and we found three SCN resistance-associated QTNs on chromosomes 7 and 18. Through these analyses, we were able to decipher the impact of seven major genetic loci, including three novel loci, on resistance to several SCN populations and identified candidate genes within each locus. Further, we identified favorable allelic combinations for resistance to individual SCN HG types and provided a list of available germplasm for integration of these unique alleles into soybean breeding programs. Overall, this study offers valuable insight into the landscape of SCN resistance loci in U.S. public soybean breeding programs and provides a framework to develop new and improved soybean cultivars with diverse plant genetic modes of SCN 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.

The heterologous expression of conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs suppresses Meloidogyne incognita parasitism in Gossypium hirsutum (upland cotton)

Two conserved Glycine max (soybean) mitogen activated protein kinase 3 (MAPK3) paralogs function in defense to the parasitic soybean cyst nematode Heterodera glycines. Gene Ontology analyses of RNA seq data obtained from MAPK3-1-overexpressing (OE) and MAPK3-2-OE roots compared to their control, as well as MAPK3-1-RNA interference (RNAi) and MAPK3-2-RNAi compared to their control, hierarchically orders the induced and suppressed genes, strengthening the hypothesis that their heterologous expression in Gossypium hirsutum (upland cotton) would impair parasitism by the root knot nematode (RKN) Meloidogyne incognita. MAPK3-1 expression (E) in G. hirsutum suppresses the production of M. incognita root galls, egg masses, and second stage juveniles (J2s) by 80.32%, 82.37%, and 88.21%, respectfully. Unexpectedly, egg number increases by 28.99% but J2s are inviable. MAPK3-2-E effects are identical, statistically. MAPK3-1-E and MAPK3-2-E decreases root mass 1.49-fold and 1.55-fold, respectively, as compared to the pRAP15-ccdB-E control. The reproductive factor (RF) of M. incognita for G. hirsutum roots expressing MAPK3-1-E or MAPK3-2-E decreases 60.39% and 50.46%, respectively, compared to controls. The results are consistent with upstream pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) functioning in defense to H. glycines. The experiments showcase the feasibility of employing MAPK3, through heterologous expression, to combat M. incognita parasitism, possibly overcoming impediments otherwise making G. hirsutum's defense platform deficient. MAPK homologs are identified in other important crop species for future functional analyses.

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.

The central circadian regulator CCA1 functions in Glycine max during defense to a root pathogen, regulating the expression of genes acting in effector triggered immunity (ETI) and cell wall metabolism

Expression of the central circadian oscillator components CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), TIMING OF CAB1 (TOC1), GIGANTEA (GI), and CONSTANS (CO) occurs in Glycine max root cells (syncytia) parasitized by the nematode Heterodera glycines while undergoing resistance, indicating a defense role. GmCCA1-1 relative transcript abundance (RTA) in roots experiencing overexpression (OE) or RNA interference (RNAi) is characterized by rhythmic oscillations, compared to a ribosomal protein gene (GmRPS21) control. A GmCCA1-1 RTA change, advancing by 12 h, exists in H. glycines-infected as compared to uninfected controls in wild-type, H. glycines-resistant, G. max[Peking/PI 548402]. The G. max[Peking/PI 548402] transgenic controls exhibit the RTA change by 4 h post infection (hpi), not consistently occurring in the H. glycines-susceptible G. max[Williams 82/PI 518671] until 56 hpi. GmCCA1-1 expression is observed to be reduced in H. glycines-infected GmCCA1-1-OE roots as compared to non-infected transgenic roots with no significant change observed among RNAi roots. The GmCCA1-1 expression in transgenic GmCCA1-1-OE roots remains higher than control and RNAi roots. Decreased GmCCA1-1 mRNA among infected roots shows the altered expression is targeted by H. glycines. Gene expression of proven defense genes including 9 different mitogen activated protein kinases (GmMAPKs), NON-RACE SPECIFIC DISEASE RESISTANCE 1 (GmNDR1-1), RPM1-INTERACTING PROTEIN 4 (GmRIN4-4), and the secreted xyloglucan endotransglycosylase/hydrolase 43 (GmXTH43) in GmCCA1-1-OE and GmCCA1-1-RNAi roots, compared to controls, reveal a significant role of GmCCA1-1 expression in roots undergoing defense to H. glycines parasitism. The observation that GmCCA1-1 regulates GmXTH43 expression links the central circadian oscillator to the functionality of the secretion system.

Glycine max Homologs of DOESN'T MAKE INFECTIONS 1, 2, and 3 Function to Impair Heterodera glycines Parasitism While Also Regulating Mitogen Activated Protein Kinase Expression

Glycine max root cells developing into syncytia through the parasitic activities of the pathogenic nematode Heterodera glycines underwent isolation by laser microdissection (LM). Microarray analyses have identified the expression of a G. max DOESN'T MAKE INFECTIONS3 (DMI3) homolog in syncytia undergoing parasitism but during a defense response. DMI3 encodes part of the common symbiosis pathway (CSP) involving DMI1, DMI2, and other CSP genes. The identified DMI gene expression, and symbiosis role, suggests the possible existence of commonalities between symbiosis and defense. G. max has 3 DMI1, 12 DMI2, and 2 DMI3 paralogs. LM-assisted gene expression experiments of isolated syncytia under further examination here show G. max DMI1-3, DMI2-7, and DMI3-2 expression occurring during the defense response in the H. glycines-resistant genotypes G.max _{[Peking/PI548402]} and G.max _[PI88788] indicating a broad and consistent level of expression of the genes. Transgenic overexpression (OE) of G. max DMI1-3, DMI2-7, and DMI3-2 impairs H. glycines parasitism. RNA interference (RNAi) of G. max DMI1-3, DMI2-7, and DMI3-2 increases H. glycines parasitism. The combined opposite outcomes reveal a defense function for these genes. Prior functional transgenic analyses of the 32-member G. max mitogen activated protein kinase (MAPK) gene family has determined that 9 of them act in the defense response to H. glycines parasitism, referred to as defense MAPKs. RNA-seq analyses of root RNA isolated from the 9 G. max defense MAPKs undergoing OE or RNAi reveal they alter the relative transcript abundances (RTAs) of specific DMI1, DMI2, and DMI3 paralogs. In contrast, transgenically-manipulated DMI1-3, DMI2-7, and DMI3-2 expression influences MAPK3-1 and MAPK3-2 RTAs under certain circumstances. The results show G. max homologs of the CSP, and defense pathway are linked, apparently involving co-regulated gene expression.

Enhanced resistance to soybean cyst nematode in transgenic soybean via host-induced silencing of vital Heterodera glycines genes

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is the most economically damaging pathogen affecting soybean production worldwide. Host-induced gene silencing provides a promising approach to confer resistance to plant parasitic nematodes. In the present study, we produced stable transgenic soybean plants individually harboring the inverted repeats of three essential H. glycines genes, Hg-rps23, Hg-snb1, and Hg-cpn1, and evaluated their resistance to SCN infection. Molecular characterization confirmed the stable integration of the hairpin double stranded (ds) RNA in host plants. Inoculation assays with SCN race 3 showed significant reduction of female index (FI, 11.84 ~ 17.47%) on the roots of T4 transgenic plants, with 73.29 ~ 81.90% reduction for the three RNA interference (RNAi) constructs, compared to non-transformed plants (NT, 65.43%). Enhanced resistance to SCN race 3 was further confirmed in subsequent generations (T5) of transgenic soybean. Moreover, when inoculated with SCN race 4 which was considered highly virulent to most of soybean germplasms and varieties, transgenic soybean plants also exhibited reduced FIs (9.96 ~ 23.67%) and increased resistance, relative to the NT plants (46.46%). Consistently, significant down-regulation in transcript levels of the Hg-rps23, Hg-snb1, Hg-cpn1 genes were observed in the nematodes feeding on the transgenic roots, suggesting a broad-spectrum resistance mediated by the host-mediated silencing of vital H. glycines genes. There were no significant differences in morphological traits between transgenic and NT soybean plants under conditions with negligible SCN infection. In summary, our results demonstrate the effectiveness of host-induced silencing of essential H. glycines genes to enhance broad-spectrum SCN resistance in stable transgenic soybean plants, without negative consequences on the agronomic performance.

Detection of rare nematode resistance Rhg1 haplotypes in Glycine soja and a novel Rhg1 α-SNAP

This study pursued the hypothesis that wild plant germplasm accessions carrying alleles of interest can be identified using available single nucleotide polymorphism (SNP) genotypes for particular alleles of other (unlinked) genes that contribute to the trait of interest. The soybean cyst nematode (SCN, Heterodera glycines [HG]) resistance locus Rhg1 is widely used in farmed soybean [Glycine max (L.) Merr.]. The two known resistance-conferring haplotypes, rhg1-a and rhg1-b, typically contain three or seven to 10 tandemly duplicated Rhg1 segments, respectively. Each Rhg1 repeat carries four genes, including Glyma.18G022500, which encodes unusual isoforms of the vesicle-trafficking chaperone α-SNAP. Using SoySNP50K data for NSFRAN07 allele presence, we discovered a new Rhg1 haplotype, rhg1-ds, in six accessions of wild soybean, Glycine soja Siebold & Zucc. (0.5% of the ∼1,100 G. soja accessions in the USDA collection). The α-SNAP encoded by rhg1-ds is unique at an important site of amino acid variation and shares with the rhg1-a and rhg1-b α-SNAP proteins the traits of cytotoxicity and altered N-ethylmaleimide sensitive factor (NSF) protein interaction. Copy number assays indicate three repeats of rhg1-ds. G. soja PI 507613 and PI 507623 exhibit resistance to HG type 2.5.7 SCN populations, in part because of contributions from other loci. In a segregating F2 population, rhg1-b and rhg1-ds made statistically indistinguishable contributions to resistance to a partially virulent HG type 2.5.7 SCN population. Hence, the unusual multigene copy number variation Rhg1 haplotype was present but rare in ancestral G. soja and was present in accessions that offer multiple traits for SCN resistance breeding. The accessions were initially identified for study based on an unlinked SNP.

A Qa-SNARE complex contributes to soybean cyst nematode resistance via regulation of mitochondria-mediated cell death

The resistance to Heterodera glycines 1 (Rhg1) locus is widely used by soybean breeders to reduce yield loss caused by soybean cyst nematode (SCN). α-SNAP (α-soluble NSF attachment protein) within Rhg1 locus contributes to SCN resistance by modulation of cell status at the SCN feeding site; however, the underlying mechanism is largely unclear. Here, we identified an α-SNAP-interacting protein, GmSYP31A, a Qa-SNARE (soluble NSF attachment protein receptor) protein from soybean. Expression of GmSYP31A significantly induced cell death in Nicotiana benthamiana leaves, and co-expression of α-SNAP and GmSYP31A could accelerate cell death. Overexpression of GmSYP31A increased SCN resistance, while silencing or overexpression of a dominant-negative form of GmSYP31A increased SCN sensitivity. GmSYP31A expression also disrupted endoplasmic reticulum-Golgi trafficking, and the exocytosis pathway. Moreover, α-SNAP was also found to interact with GmVDAC1D (voltage-dependent anion channel). The cytotoxicity induced by the expression of GmSYP31A could be relieved either with the addition of an inhibitor of VDAC protein, or by silencing the VDAC gene. Taken together, our data not only demonstrate that α-SNAP works together with GmSYP31A to increase SCN resistance through triggering cell death, but also highlight the unexplored link between the mitochondrial apoptosis pathway and vesicle trafficking.

Conserved oligomeric Golgi (COG) complex genes functioning in defense are expressed in root cells undergoing a defense response to a pathogenic infection and exhibit regulation my MAPKs

The conserved oligomeric Golgi (COG) complex maintains correct Golgi structure and function during retrograde trafficking. Glycine max has 2 paralogs of each COG gene, with one paralog of each gene family having a defense function to the parasitic nematode Heterodera glycines. Experiments presented here show G. max COG paralogs functioning in defense are expressed specifically in the root cells (syncytia) undergoing the defense response. The expressed defense COG gene COG7-2-b is an alternate splice variant, indicating specific COG variants are important to defense. Transcriptomic experiments examining RNA isolated from COG overexpressing and RNAi roots show some COG genes co-regulate the expression of other COG complex genes. Examining signaling events responsible for COG expression, transcriptomic experiments probing MAPK overexpressing roots show their expression influences the relative transcript abundance of COG genes as compared to controls. COG complex paralogs are shown to be found in plants that are agriculturally relevant on a world-wide scale including Manihot esculenta, Zea mays, Oryza sativa, Triticum aestivum, Hordeum vulgare, Sorghum bicolor, Brassica rapa, Elaes guineensis and Saccharum officinalis and in additional crops significant to U.S. agriculture including Beta vulgaris, Solanum tuberosum, Solanum lycopersicum and Gossypium hirsutum. The analyses provide basic information on COG complex biology, including the coregulation of some COG genes and that MAPKs functioning in defense influence their expression. Furthermore, it appears in G. max and likely other crops that some level of neofunctionalization of the duplicated genes is occurring. The analysis has identified important avenues for future research broadly in plants.

The canonical α-SNAP is essential for gametophytic development in Arabidopsis

The development of male and female gametophytes is a pre-requisite for successful reproduction of angiosperms. Factors mediating vesicular trafficking are among the key regulators controlling gametophytic development. Fusion between vesicles and target membranes requires the assembly of a fusogenic soluble N-ethylmaleimide sensitive factor attachment protein receptors (SNAREs) complex, whose disassembly in turn ensures the recycle of individual SNARE components. The disassembly of post-fusion SNARE complexes is controlled by the AAA+ ATPase N-ethylmaleimide-sensitive factor (Sec18/NSF) and soluble NSF attachment protein (Sec17/α-SNAP) in yeast and metazoans. Although non-canonical α-SNAPs have been functionally characterized in soybeans, the biological function of canonical α-SNAPs has yet to be demonstrated in plants. We report here that the canonical α-SNAP in Arabidopsis is essential for male and female gametophytic development. Functional loss of the canonical α-SNAP in Arabidopsis results in gametophytic lethality by arresting the first mitosis during gametogenesis. We further show that Arabidopsis α-SNAP encodes two isoforms due to alternative splicing. Both isoforms interact with the Arabidopsis homolog of NSF whereas have distinct subcellular localizations. The presence of similar alternative splicing of human α-SNAP indicates that functional distinction of two α-SNAP isoforms is evolutionarily conserved.

Xyloglucan endotransglycosylase/hydrolase increases tightly-bound xyloglucan and chain number but decreases chain length contributing to the defense response that Glycine max has to Heterodera glycines

The Glycine max xyloglucan endotransglycosylase/hydrolase (EC 2.4.1.207), GmXTH43, has been identified through RNA sequencing of RNA isolated through laser microdissection of Heterodera glycines-parasitized root cells (syncytia) undergoing the process of defense. Experiments reveal that genetically increasing XTH43 transcript abundance in the H. glycines-susceptible genotype G. max[Williams 82/PI 518671] decreases parasitism. Experiments presented here show decreasing XTH43 transcript abundance through RNA interference (RNAi) in the H. glycines-resistant G. max[Peking/PI 548402] increases susceptibility, but it is unclear what role XTH43 performs. The experiments presented here show XTH43 overexpression decreases the relative length of xyloglucan (XyG) chains, however, there is an increase in the amount of those shorter chains. In contrast, XTH43 RNAi increases XyG chain length. The experiments show that XTH43 has the capability to function, when increased in its expression, to limit XyG chain extension. This outcome would likely impair the ability of the cell wall to expand. Consequently, XTH43 could provide an enzymatically-driven capability to the cell that would allow it to limit the ability of parasitic nematodes like H. glycines to develop a feeding structure that, otherwise, would facilitate parasitism. The experiments presented here provide experimentally-based proof that XTHs can function in ways that could be viewed as being able to limit the expansion of the cell wall.

Mitogen activated protein kinase (MAPK)-regulated genes with predicted signal peptides function in the Glycine max defense response to the root pathogenic nematode Heterodera glycines

Glycine max has 32 mitogen activated protein kinases (MAPKs), nine of them exhibiting defense functions (defense MAPKs) to the plant parasitic nematode Heterodera glycines. RNA seq analyses of transgenic G. max lines overexpressing (OE) each defense MAPK has led to the identification of 309 genes that are increased in their relative transcript abundance by all 9 defense MAPKs. Here, 71 of those genes are shown to also have measurable amounts of transcript in H. glycines-induced nurse cells (syncytia) produced in the root that are undergoing a defense response. The 71 genes have been grouped into 7 types, based on their expression profile. Among the 71 genes are 8 putatively-secreted proteins that include a galactose mutarotase-like protein, pollen Ole e 1 allergen and extensin protein, endomembrane protein 70 protein, O-glycosyl hydrolase 17 protein, glycosyl hydrolase 32 protein, FASCICLIN-like arabinogalactan protein 17 precursor, secreted peroxidase and a pathogenesis-related thaumatin protein. Functional transgenic analyses of all 8 of these candidate defense genes that employ their overexpression and RNA interference (RNAi) demonstrate they have a role in defense. Overexpression experiments that increase the relative transcript abundance of the candidate defense gene reduces the ability that the plant parasitic nematode Heterodera glycines has in completing its life cycle while, in contrast, RNAi of these genes leads to an increase in parasitism. The results provide a genomic analysis of the importance of MAPK signaling in relation to the secretion apparatus during the defense process defense in the G. max-H. glycines pathosystem and identify additional targets for future studies.

Exocyst components promote an incompatible interaction between Glycine max (soybean) and Heterodera glycines (the soybean cyst nematode)

Vesicle and target membrane fusion involves tethering, docking and fusion. The GTPase SECRETORY4 (SEC4) positions the exocyst complex during vesicle membrane tethering, facilitating docking and fusion. Glycine max (soybean) Sec4 functions in the root during its defense against the parasitic nematode Heterodera glycines as it attempts to develop a multinucleate nurse cell (syncytium) serving to nourish the nematode over its 30-day life cycle. Results indicate that other tethering proteins are also important for defense. The G. max exocyst is encoded by 61 genes: 5 EXOC1 (Sec3), 2 EXOC2 (Sec5), 5 EXOC3 (Sec6), 2 EXOC4 (Sec8), 2 EXOC5 (Sec10) 6 EXOC6 (Sec15), 31 EXOC7 (Exo70) and 8 EXOC8 (Exo84) genes. At least one member of each gene family is expressed within the syncytium during the defense response. Syncytium-expressed exocyst genes function in defense while some are under transcriptional regulation by mitogen-activated protein kinases (MAPKs). The exocyst component EXOC7-H4-1 is not expressed within the syncytium but functions in defense and is under MAPK regulation. The tethering stage of vesicle transport has been demonstrated to play an important role in defense in the G. max-H. glycines pathosystem, with some of the spatially and temporally regulated exocyst components under transcriptional control by MAPKs.

The heterologous expression of a soybean (Glycine max) xyloglucan endotransglycosylase/hydrolase (XTH) in cotton (Gossypium hirsutum) suppresses parasitism by the root knot nematode Meloidogyne incognita

A Glycine max (soybean) hemicellulose modifying gene, xyloglucan endotransglycoslase/hydrolase (XTH43), has been identified as being expressed within a nurse cell known as a syncytium developing within the soybean root undergoing the process of defense to infection by the parasitic nematode, Heterodera glycines. The highly effective nature of XTH43 overexpression in suppressing H. glycines parasitism in soybean has led to experiments examining whether the heterologous expression of XTH43 in Gossypium hirsutum (upland cotton) could impair the parasitism of Meloidogyne incognita, that form a different type of nurse cell called a giant cell that is enclosed within a swollen root structure called a gall. The heterologous transgenic expression of XTH43 in cotton resulted in an 18% decrease in the number of galls, 70% decrease in egg masses, 64% decrease in egg production and a 97% decrease in second stage juvenile (J2) production as compared to transgenic controls. The heterologous XTH43 expression does not significantly affect root mass. The results demonstrate XTH43 expression functions effectively in impairing the development of M. incognita at numerous life cycle stages occurring within the cotton root. The experiments reveal that there are highly conserved aspects of the defense response of G. max that can function effectively in G. hirsutum to impair M. incognita having a different method of parasitism.

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.

The Glycine max Conserved Oligomeric Golgi (COG) Complex Functions During a Defense Response to Heterodera glycines

The conserved oligomeric Golgi (COG) complex, functioning in retrograde trafficking, is a universal structure present among eukaryotes that maintains the correct Golgi structure and function. The COG complex is composed of eight subunits coalescing into two sub-complexes. COGs1-4 compose Sub-complex A. COGs5-8 compose Sub-complex B. The observation that COG interacts with the syntaxins, suppressors of the erd2-deletion 5 (Sed5p), is noteworthy because Sed5p also interacts with Sec17p [alpha soluble NSF attachment protein (α-SNAP)]. The α-SNAP gene is located within the major Heterodera glycines [soybean cyst nematode (SCN)] resistance locus (rhg1) and functions in resistance. The study presented here provides a functional analysis of the Glycine max COG complex. The analysis has identified two paralogs of each COG gene. Functional transgenic studies demonstrate at least one paralog of each COG gene family functions in G. max during H. glycines resistance. Furthermore, treatment of G. max with the bacterial effector harpin, known to function in effector triggered immunity (ETI), leads to the induced transcription of at least one member of each COG gene family that has a role in H. glycines resistance. In some instances, altered COG gene expression changes the relative transcript abundance of syntaxin 31. These results indicate that the G. max COG complex functions through processes involving ETI leading to H. glycines resistance.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid–nematode interactions, particularly in the soybean–SBA–SCN system.

The rhg1-a (Rhg1 low-copy) nematode resistance source harbors a copia-family retrotransposon within the Rhg1-encoded α-SNAP gene

Soybean growers widely use the Resistance to Heterodera glycines 1 (Rhg1) locus to reduce yield losses caused by soybean cyst nematode (SCN). Rhg1 is a tandemly repeated four gene block. Two classes of SCN resistance-conferring Rhg1 haplotypes are recognized: rhg1-a ("Peking-type," low-copy number, three or fewer Rhg1 repeats) and rhg1-b ("PI 88788-type," high-copy number, four or more Rhg1 repeats). The rhg1-a and rhg1-b haplotypes encode α-SNAP (alpha-Soluble NSF Attachment Protein) variants α-SNAP Rhg1 LC and α-SNAP Rhg1 HC, respectively, with differing atypical C-terminal domains, that contribute to SCN resistance. Here we report that rhg1-a soybean accessions harbor a copia retrotransposon within their Rhg1 Glyma.18G022500 (α-SNAP-encoding) gene. We termed this retrotransposon "RAC," for Rhg1 alpha-SNAP copia. Soybean carries multiple RAC-like retrotransposon sequences. The Rhg1 RAC insertion is in the Glyma.18G022500 genes of all true rhg1-a haplotypes we tested and was not detected in any examined rhg1-b or Rhg1WT (single-copy) soybeans. RAC is an intact element residing within intron 1, anti-sense to the rhg1-a α-SNAP open reading frame. RAC has intrinsic promoter activities, but overt impacts of RAC on transgenic α-SNAP Rhg1 LC mRNA and protein abundance were not detected. From the native rhg1-a RAC+ genomic context, elevated α-SNAP Rhg1 LC protein abundance was observed in syncytium cells, as was previously observed for α-SNAP Rhg1 HC (whose rhg1-b does not carry RAC). Using a SoySNP50K SNP corresponding with RAC presence, just ~42% of USDA accessions bearing previously identified rhg1-a SoySNP50K SNP signatures harbor the RAC insertion. Subsequent analysis of several of these putative rhg1-a accessions lacking RAC revealed that none encoded α-SNAPRhg1LC, and thus, they are not rhg1-a. rhg1-a haplotypes are of rising interest, with Rhg4, for combating SCN populations that exhibit increased virulence against the widely used rhg1-b resistance. The present study reveals another unexpected structural feature of many Rhg1 loci, and a selectable feature that is predictive of rhg1-a haplotypes.

The mitogen activated protein kinase (MAPK) gene family functions as a cohort during the Glycine max defense response to Heterodera glycines

Mitogen activated protein kinases (MAPKs) play important signal transduction roles. However, little is known regarding how they influence the gene expression of other family members and the relationship to a biological process, including the Glycine max defense response to Heterodera glycines. Transcriptomics have identified MAPK gene expression occurring within root cells undergoing a defense response to a pathogenic event initiated by H. glycines in the allotetraploid Glycine max. Functional analyses are presented for its 32 MAPKs revealing 9 have a defense role, including homologs of Arabidopsis thaliana MAPK (MPK) MPK2, MPK3, MPK4, MPK5, MPK6, MPK13, MPK16 and MPK20. Defense signaling occurring through pathogen activated molecular pattern (PAMP) triggered immunity (PTI) and effector triggered immunity (ETI) have been determined in relation to these MAPKs. Five different types of gene expression relate to MAPK expression, influencing PTI and ETI gene expression and proven defense genes including an ABC-G transporter, 20S membrane fusion particle components, glycoside biosynthesis, carbon metabolism, hemicellulose modification, transcription and secretion. The experiments show MAPKs broadly influence defense MAPK gene expression, including the co-regulation of parologous MAPKs and reveal its relationship to proven defense genes. The experiments reveal each defense MAPK induces the expression of a G. max homolog of a PATHOGENESIS RELATED1 (PR1), itself shown to function in defense in the studied pathosystem.

Comparative Transcriptome Analysis of the Pinewood Nematode Bursaphelenchus xylophilus Reveals the Molecular Mechanism Underlying Its Defense Response to Host-Derived α-pinene

Bursaphelenchus xylophilus is fatal to the pine trees around the world. The production of the pine tree secondary metabolite gradually increases in response to a B. xylophilus infestation, via a stress reaction mechanism(s). α-pinene is needed to combat the early stages of B. xylophilus infection and colonization, and to counter its pathogenesis. Therefore, research is needed to characterize the underlying molecular response(s) of B. xylophilus to resist α-pinene. We examined the effects of different concentrations of α-pinene on the mortality and reproduction rate of B. xylophilus in vitro. The molecular response by which B. xylophilus resists α-pinene was examined via comparative transcriptomics of the nematode. Notably, B. xylophilus genes involved in detoxification, transport, and receptor activities were differentially expressed in response to two different concentrations of α-pinene compared with control. Our results contribute to our understanding of the molecular mechanisms by which B. xylophilus responds to monoterpenes in general, and the pathogenesis of B. xylophilus.

Harpin-inducible defense signaling components impair infection by the ascomycete Macrophomina phaseolina

Soybean (Glycine max) infection by the charcoal rot (CR) ascomycete Macrophomina phaseolina is enhanced by the soybean cyst nematode (SCN) Heterodera glycines. We hypothesized that G. max genetic lines impairing infection by M. phaseolina would also limit H. glycines parasitism, leading to resistance. As a part of this M. phaseolina resistance process, the genetic line would express defense genes already proven to impair nematode parasitism. Using G. max[DT97-4290/PI 642055], exhibiting partial resistance to M. phaseolina, experiments show the genetic line also impairs H. glycines parasitism. Furthermore, comparative studies show G. max[DT97-4290/PI 642055] exhibits induced expression of the effector triggered immunity (ETI) gene NON-RACE SPECIFIC DISEASE RESISTANCE 1/HARPIN INDUCED1 (NDR1/HIN1) that functions in defense to H. glycines as compared to the H. glycines and M. phaseolina susceptible line G. max[Williams 82/PI 518671]. Other defense genes that are induced in G. max[DT97-4290/PI 642055] include the pathogen associated molecular pattern (PAMP) triggered immunity (PTI) genes ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1), NONEXPRESSOR OF PR1 (NPR1) and TGA2. These observations link G. max defense processes that impede H. glycines parasitism to also potentially function toward impairing M. phaseolina pathogenicity. Testing this hypothesis, G. max[Williams 82/PI 518671] genetically engineered to experimentally induce GmNDR1-1, EDS1-2, NPR1-2 and TGA2-1 expression leads to impaired M. phaseolina pathogenicity. In contrast, G. max[DT97-4290/PI 642055] engineered to experimentally suppress the expression of GmNDR1-1, EDS1-2, NPR1-2 and TGA2-1 by RNA interference (RNAi) enhances M. phaseolina pathogenicity. The results show components of PTI and ETI impair both nematode and M. phaseolina pathogenicity.

An atypical N-ethylmaleimide sensitive factor enables the viability of nematode-resistant Rhg1 soybeans

N-ethylmaleimide sensitive factor (NSF) and α-soluble NSF attachment protein (α-SNAP) are essential eukaryotic housekeeping proteins that cooperatively function to sustain vesicular trafficking. The "resistance to Heterodera glycines 1" (Rhg1) locus of soybean (Glycine max) confers resistance to soybean cyst nematode, a highly damaging soybean pest. Rhg1 loci encode repeat copies of atypical α-SNAP proteins that are defective in promoting NSF function and are cytotoxic in certain contexts. Here, we discovered an unusual NSF allele (Rhg1-associated NSF on chromosome 07; NSFRAN07 ) in Rhg1+ germplasm. NSFRAN07 protein modeling to mammalian NSF/α-SNAP complex structures indicated that at least three of the five NSFRAN07 polymorphisms reside adjacent to the α-SNAP binding interface. NSFRAN07 exhibited stronger in vitro binding with Rhg1 resistance-type α-SNAPs. NSFRAN07 coexpression in planta was more protective against Rhg1 α-SNAP cytotoxicity, relative to WT NSFCh07 Investigation of a previously reported segregation distortion between chromosome 18 Rhg1 and a chromosome 07 interval now known to contain the Glyma.07G195900 NSF gene revealed 100% coinheritance of the NSFRAN07 allele with disease resistance Rhg1 alleles, across 855 soybean accessions and in all examined Rhg1+ progeny from biparental crosses. Additionally, we show that some Rhg1-mediated resistance is associated with depletion of WT α-SNAP abundance via selective loss of WT α-SNAP loci. Hence atypical coevolution of the soybean SNARE-recycling machinery has balanced the acquisition of an otherwise disruptive housekeeping protein, enabling a valuable disease resistance trait. Our findings further indicate that successful engineering of Rhg1-related resistance in plants will require a compatible NSF partner for the resistance-conferring α-SNAP.

A harpin elicitor induces the expression of a coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene and others functioning during defense to parasitic nematodes

The bacterial effector harpin induces the transcription of the Arabidopsis thaliana NON-RACE SPECIFIC DISEASE RESISTANCE 1/HARPIN INDUCED1 (NDR1/HIN1) coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene. In Glycine max, Gm-NDR1-1 transcripts have been detected within root cells undergoing a natural resistant reaction to parasitism by the syncytium-forming nematode Heterodera glycines, functioning in the defense response. Expressing Gm-NDR1-1 in Gossypium hirsutum leads to resistance to Meloidogyne incognita parasitism. In experiments presented here, the heterologous expression of Gm-NDR1-1 in G. hirsutum impairs Rotylenchulus reniformis parasitism. These results are consistent with the hypothesis that Gm-NDR1-1 expression functions broadly in generating a defense response. To examine a possible relationship with harpin, G. max plants topically treated with harpin result in induction of the transcription of Gm-NDR1-1. The result indicates the topical treatment of plants with harpin, itself, may lead to impaired nematode parasitism. Topical harpin treatments are shown to impair G. max parasitism by H. glycines, M. incognita and R. reniformis and G. hirsutum parasitism by M. incognita and R. reniformis. How harpin could function in defense has been examined in experiments showing it also induces transcription of G. max homologs of the proven defense genes ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1), TGA2, galactinol synthase, reticuline oxidase, xyloglucan endotransglycosylase/hydrolase, alpha soluble N-ethylmaleimide-sensitive fusion protein (α-SNAP) and serine hydroxymethyltransferase (SHMT). In contrast, other defense genes are not directly transcriptionally activated by harpin. The results indicate harpin induces pathogen associated molecular pattern (PAMP) triggered immunity (PTI) and effector-triggered immunity (ETI) defense processes in the root, activating defense to parasitic nematodes.

A Glycine max homolog of NON-RACE SPECIFIC DISEASE RESISTANCE 1 (NDR1) alters defense gene expression while functioning during a resistance response to different root pathogens in different genetic backgrounds

A Glycine max homolog of the Arabidopsis thaliana NON-RACE SPECIFIC DISEASE RESISTANCE 1 (NDR1) coiled-coil nucleotide binding leucine rich repeat (CC-NB-LRR) defense signaling gene (Gm-NDR1-1) is expressed in root cells undergoing a defense response to the root pathogenic nematode, Heterodera glycines. Gm-NDR1-1 overexpression in the H. glycines-susceptible genotype G. max[Williams 82/PI 518671] impairs parasitism. In contrast, Gm-NDR1-1 RNA interference (RNAi) in the H. glycines-resistant genotype G. max[Peking/PI 548402] facilitates parasitism. The broad effectiveness of Gm-NDR1-1 in impairing parasitism has then been examined by engineering its heterologous expression in Gossypium hirsutum which is susceptible to the root pathogenic nematode Meloidogyne incognita. The heterologous expression of Gm-NDR1-1 in G. hirsutum effectively impairs M. incognita parasitism, reducing gall, egg mass, egg and juvenile numbers. In contrast to our prior experiments examining the effectiveness of the heterologous expression of a G. max homolog of the A. thaliana salicyclic acid signaling (SA) gene NONEXPRESSOR OF PR1 (Gm-NPR1-2), no cumulative negative effect on M. incognita parasitism has been observed in G. hirsutum expressing Gm-NDR1-1. The results indicate a common genetic basis exists for plant resistance to parasitic nematodes that involves Gm-NDR1. However, the Gm-NDR1-1 functions in ways that are measurably dissimilar to Gm-NPR1-2. Notably, Gm-NDR1-1 overexpression leads to increased relative transcript levels of its homologs of A. thaliana genes functioning in SA signaling, including NPR1-2, TGA2-1 and LESION SIMULATING DISEASE1 (LSD1-2) that is lost in Gm-NDR1-1 RNAi lines. Similar observations have been made regarding the expression of other defense genes.

Characterization of the Soluble NSF Attachment Protein gene family identifies two members involved in additive resistance to a plant pathogen

Proteins with Tetratricopeptide-repeat (TPR) domains are encoded by large gene families and distributed in all plant lineages. In this study, the Soluble NSF-Attachment Protein (SNAP) subfamily of TPR containing proteins is characterized. In soybean, five members constitute the SNAP gene family: GmSNAP18, GmSNAP11, GmSNAP14, GmSNAP02, and GmSNAP09. Recently, GmSNAP18 has been reported to mediate resistance to soybean cyst nematode (SCN). Using a population of recombinant inbred lines from resistant and susceptible parents, the divergence of the SNAP gene family is analysed over time. Phylogenetic analysis of SNAP genes from 22 diverse plant species showed that SNAPs were distributed in six monophyletic clades corresponding to the major plant lineages. Conservation of the four TPR motifs in all species, including ancestral lineages, supports the hypothesis that SNAPs were duplicated and derived from a common ancestor and unique gene still present in chlorophytic algae. Syntenic analysis of regions harbouring GmSNAP genes in soybean reveals that this family expanded from segmental and tandem duplications following a tetraploidization event. qRT-PCR analysis of GmSNAPs indicates a co-regulation following SCN infection. Finally, genetic analysis demonstrates that GmSNAP11 contributes to an additive resistance to SCN. Thus, GmSNAP11 is identified as a novel minor gene conferring resistance to SCN.

Components of the SNARE-containing regulon are co-regulated in root cells undergoing defense

The term regulon has been coined in the genetic model plant Arabidopsis thaliana, denoting a structural and physiological defense apparatus defined genetically through the identification of the penetration (pen) mutants. The regulon is composed partially by the soluble N-ethylmaleimide-sensitive fusion protein attachment protein receptor (SNARE) syntaxin PEN1. PEN1 has homology to a Saccharomyces cerevisae gene that regulates a Secretion (Sec) protein, Suppressor of Sec 1 (Sso1p). The regulon is also composed of the β-glucosidase (PEN2) and an ATP binding cassette (ABC) transporter (PEN3). While important in inhibiting pathogen infection, limited observations have been made regarding the transcriptional regulation of regulon genes until now. Experiments made using the model agricultural Glycine max (soybean) have identified co-regulated gene expression of regulon components. The results explain the observation of hundreds of genes expressed specifically in the root cells undergoing the natural process of defense. Data regarding additional G. max genes functioning within the context of the regulon are presented here, including Sec 14, Sec 4 and Sec 23. Other examined G. max homologs of membrane fusion genes include an endosomal bromo domain-containing protein1 (Bro1), syntaxin6 (SYP6), SYP131, SYP71, SYP8, Bet1, coatomer epsilon (ϵ-COP), a coatomer zeta (ζ-COP) paralog and an ER to Golgi component (ERGIC) protein. Furthermore, the effectiveness of biochemical pathways that would function within the context of the regulon ave been examined, including xyloglucan xylosyltransferase (XXT), reticuline oxidase (RO) and galactinol synthase (GS). The experiments have unveiled the importance of the regulon during defense in the root and show how the deposition of callose relates to the process.

Combined transcriptome and metabolome analyses to understand the dynamic responses of rice plants to attack by the rice stem borer Chilo suppressalis (Lepidoptera: Crambidae)

BACKGROUND

Rice (Oryza sativa L.), which is a staple food for more than half of the world's population, is frequently attacked by herbivorous insects, including the rice stem borer, Chilo suppressalis. C. suppressalis substantially reduces rice yields in temperate regions of Asia, but little is known about how rice plants defend themselves against this herbivore at molecular and biochemical level.

RESULTS

In the current study, we combined next-generation RNA sequencing and metabolomics techniques to investigate the changes in gene expression and in metabolic processes in rice plants that had been continuously fed by C. suppressalis larvae for different durations (0, 24, 48, 72, and 96 h). Furthermore, the data were validated using quantitative real-time PCR. There were 4,729 genes and 151 metabolites differently regulated when rice plants were damaged by C. suppressalis larvae. Further analyses showed that defense-related phytohormones, transcript factors, shikimate-mediated and terpenoid-related secondary metabolism were activated, whereas the growth-related counterparts were suppressed by C. suppressalis feeding. The activated defense was fueled by catabolism of energy storage compounds such as monosaccharides, which meanwhile resulted in the increased levels of metabolites that were involved in rice plant defense response. Comparable analyses showed a correspondence between transcript patterns and metabolite profiles.

CONCLUSION

The current findings greatly enhance our understanding of the mechanisms of induced defense response in rice plants against C. suppressalis infestation at molecular and biochemical levels, and will provide clues for development of insect-resistant rice varieties.

Disease resistance through impairment of α-SNAP-NSF interaction and vesicular trafficking by soybean Rhg1

α-SNAP [soluble NSF (N-ethylmaleimide-sensitive factor) attachment protein] and NSF proteins are conserved across eukaryotes and sustain cellular vesicle trafficking by mediating disassembly and reuse of SNARE protein complexes, which facilitate fusion of vesicles to target membranes. However, certain haplotypes of the Rhg1 (resistance to Heterodera glycines 1) locus of soybean possess multiple repeat copies of an α-SNAP gene (Glyma.18G022500) that encodes atypical amino acids at a highly conserved functional site. These Rhg1 loci mediate resistance to soybean cyst nematode (SCN; H. glycines), the most economically damaging pathogen of soybeans worldwide. Rhg1 is widely used in agriculture, but the mechanisms of Rhg1 disease resistance have remained unclear. In the present study, we found that the resistance-type Rhg1 α-SNAP is defective in interaction with NSF. Elevated in planta expression of resistance-type Rhg1 α-SNAPs depleted the abundance of SNARE-recycling 20S complexes, disrupted vesicle trafficking, induced elevated abundance of NSF, and caused cytotoxicity. Soybean, due to ancient genome duplication events, carries other loci that encode canonical (wild-type) α-SNAPs. Expression of these α-SNAPs counteracted the cytotoxicity of resistance-type Rhg1 α-SNAPs. For successful growth and reproduction, SCN dramatically reprograms a set of plant root cells and must sustain this sedentary feeding site for 2-4 weeks. Immunoblots and electron microscopy immunolocalization revealed that resistance-type α-SNAPs specifically hyperaccumulate relative to wild-type α-SNAPs at the nematode feeding site, promoting the demise of this biotrophic interface. The paradigm of disease resistance through a dysfunctional variant of an essential gene may be applicable to other plant-pathogen interactions.

A SNARE-Like Protein and Biotin Are Implicated in Soybean Cyst Nematode Virulence

Phytoparasitic nematodes that are able to infect and reproduce on plants that are considered resistant are referred to as virulent. The mechanism(s) that virulent nematodes employ to evade or suppress host plant defenses are not well understood. Here we report the use of a genetic strategy (allelic imbalance analysis) to associate single nucleotide polymorphisms (SNPs) with nematode virulence genes in Heterodera glycines, the soybean cyst nematode (SCN). To accomplish this analysis, a custom SCN SNP array was developed and used to genotype SCN F3-derived populations grown on resistant and susceptible soybean plants. Three SNPs reproducibly showed allele imbalances between nematodes grown on resistant and susceptible plants. Two candidate SCN virulence genes that were tightly linked to the SNPs were identified. One SCN gene encoded biotin synthase (HgBioB), and the other encoded a bacterial-like protein containing a putative SNARE domain (HgSLP-1). The two genes mapped to two different linkage groups. HgBioB contained sequence polymorphisms between avirulent and virulent nematodes. However, the gene encoding HgSLP-1 had reduced copy number in virulent nematode populations and appears to produce multiple forms of the protein via intron retention and alternative splicing. We show that HgSLP-1 is an esophageal-gland protein that is secreted by the nematode during plant parasitism. Furthermore, in bacterial co-expression experiments, HgSLP-1 co-purified with the SCN resistance protein Rhg1 α-SNAP, suggesting that these two proteins physically interact. Collectively our data suggest that multiple SCN genes are involved in SCN virulence, and that HgSLP-1 may function as an avirulence protein and when absent it helps SCN evade host defenses.

SNP identification and marker assay development for high-throughput selection of soybean cyst nematode resistance

BACKGROUND

Soybean cyst nematode (SCN) is the most economically devastating pathogen of soybean. Two resistance loci, Rhg1 and Rhg4 primarily contribute resistance to SCN race 3 in soybean. Peking and PI 88788 are the two major sources of SCN resistance with Peking requiring both Rhg1 and Rhg4 alleles and PI 88788 only the Rhg1 allele. Although simple sequence repeat (SSR) markers have been reported for both loci, they are linked markers and limited to be applied in breeding programs due to accuracy, throughput and cost of detection methods. The objectives of this study were to develop robust functional marker assays for high-throughput selection of SCN resistance and to differentiate the sources of resistance.

RESULTS

Based on the genomic DNA sequences of 27 soybean lines with known SCN phenotypes, we have developed Kompetitive Allele Specific PCR (KASP) assays for two Single nucleotide polymorphisms (SNPs) from Glyma08g11490 for the selection of the Rhg4 resistance allele. Moreover, the genomic DNA of Glyma18g02590 at the Rhg1 locus from 11 soybean lines and cDNA of Forrest, Essex, Williams 82 and PI 88788 were fully sequenced. Pairwise sequence alignment revealed seven SNPs/insertion/deletions (InDels), five in the 6th exon and two in the last exon. Using the same 27 soybean lines, we identified one SNP that can be used to select the Rhg1 resistance allele and another SNP that can be employed to differentiate Peking and PI 88788-type resistance. These SNP markers have been validated and a strong correlation was observed between the SNP genotypes and reactions to SCN race 3 using a panel of 153 soybean lines, as well as a bi-parental population, F5-derived recombinant inbred lines (RILs) from G00-3213xLG04-6000.

CONCLUSIONS

Three functional SNP markers (two for Rhg1 locus and one for Rhg4 locus) were identified that could provide genotype information for the selection of SCN resistance and differentiate Peking from PI 88788 source for most germplasm lines. The robust KASP SNP marker assays were developed. In most contexts, use of one or two of these markers is sufficient for high-throughput marker-assisted selection of plants that will exhibit SCN resistance.

The syntaxin 31-induced gene, LESION SIMULATING DISEASE1 (LSD1), functions in Glycine max defense to the root parasite Heterodera glycines

Experiments show the membrane fusion genes α soluble NSF attachment protein (α-SNAP) and syntaxin 31 (Gm-SYP38) contribute to the ability of Glycine max to defend itself from infection by the plant parasitic nematode Heterodera glycines. Accompanying their expression is the transcriptional activation of the defense genes ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) and NONEXPRESSOR OF PR1 (NPR1) that function in salicylic acid (SA) signaling. These results implicate the added involvement of the antiapoptotic, environmental response gene LESION SIMULATING DISEASE1 (LSD1) in defense. Roots engineered to overexpress the G. max defense genes Gm-α-SNAP, SYP38, EDS1, NPR1, BOTRYTIS INDUCED KINASE1 (BIK1) and xyloglucan endotransglycosylase/hydrolase (XTH) in the susceptible genotype G. max[Williams 82/PI 518671] have induced Gm-LSD1 (Gm-LSD1-2) transcriptional activity. In reciprocal experiments, roots engineered to overexpress Gm-LSD1-2 in the susceptible genotype G. max[Williams 82/PI 518671] have induced levels of SYP38, EDS1, NPR1, BIK1 and XTH, but not α-SNAP prior to infection. In tests examining the role of Gm-LSD1-2 in defense, its overexpression results in ∼52 to 68% reduction in nematode parasitism. In contrast, RNA interference (RNAi) of Gm-LSD1-2 in the resistant genotype G. max[Peking/PI 548402] results in an 3.24-10.42 fold increased ability of H. glycines to parasitize. The results identify that Gm-LSD1-2 functions in the defense response of G. max to H. glycines parasitism. It is proposed that LSD1, as an antiapoptotic protein, may establish an environment whereby the protected, living plant cell could secrete materials in the vicinity of the parasitizing nematode to disarm it. After the targeted incapacitation of the nematode the parasitized cell succumbs to its targeted demise as the infected root region is becoming fortified.

Manipulation of two α-endo-β-1,4-glucanase genes, AtCel6 and GmCel7, reduces susceptibility to Heterodera glycines in soybean roots

Plant endo-β-1,4-glucanases (EGases) include cell wall-modifying enzymes that are involved in nematode-induced growth of syncytia (feeding structures) in nematode-infected roots. EGases in the α- and β-subfamilies contain signal peptides and are secreted, whereas those in the γ-subfamily have a membrane-anchoring domain and are not secreted. The Arabidopsis α-EGase At1g48930, designated as AtCel6, is known to be down-regulated by beet cyst nematode (Heterodera schachtii) in Arabidopsis roots, whereas another α-EGase, AtCel2, is up-regulated. Here, we report that the ectopic expression of AtCel6 in soybean roots reduces susceptibility to both soybean cyst nematode (SCN; Heterodera glycines) and root knot nematode (Meloidogyne incognita). Suppression of GmCel7, the soybean homologue of AtCel2, in soybean roots also reduces the susceptibility to SCN. In contrast, in studies on two γ-EGases, both ectopic expression of AtKOR2 in soybean roots and suppression of the soybean homologue of AtKOR3 had no significant effect on SCN parasitism. Our results suggest that secreted α-EGases are likely to be more useful than membrane-bound γ-EGases in the development of an SCN-resistant soybean through gene manipulation. Furthermore, this study provides evidence that Arabidopsis shares molecular events of cyst nematode parasitism with soybean, and confirms the suitability of the Arabidopsis-H. schachtii interaction as a model for the soybean-H. glycines pathosystem.

Regulatory interplay between soybean root and soybean cyst nematode during a resistant and susceptible reaction

BACKGROUND

Plant-parasitic nematodes (PPNs) are obligate parasites that feed on the roots of living host plants. Often, these nematodes can lay hundreds of eggs, each capable of surviving without a host for as long as 12 years. When it comes to wreaking havoc on agricultural yield, few nematodes can compare to the soybean cyst nematode (SCN). Quantifying soybean (Glycine max) transcription factor binding sites (TFBSs) during a late-stage SCN resistant and susceptible reaction can shed light onto the systematic interplay between host and pathogen, thereby elucidating underlying cis-regulatory mechanisms.

RESULTS

We sequenced the soybean root transcriptome at 6 and 8 days upon independent inoculation with a virulent and avirulent SCN population. Genes such as β-1,4 glucanase, chalcone synthase, superoxide dismutase and various heat shock proteins (HSPs) exhibited reaction-specific expression profiles. Several likely defense-response genes candidates were also identified which are believed to confer SCN resistance. To explore magnitude of TFBS representation during SCN pathogenesis, a multivariate statistical software identified 46 over-represented TFBSs which capture soybean regulatory dynamics across both reactions.

CONCLUSIONS

Our results reveal a set of soybean TFBSs which are over-represented solely throughout a resistant and susceptible SCN reaction. This set furthers our understanding of soybean cis-regulatory dynamics by providing reaction-specific levels of over-representation at 6 and 8 days after inoculation (dai) with SCN.

Modification of the expression of two NPR1 suppressors, SNC1 and SNI1, in soybean confers partial resistance to the soybean cyst nematode, Heterodera glycines

Systemic acquired resistance (SAR) is an enhanced defence response triggered when plants detect a pathogen. The response is extended to uninfected organs to protect against future attack. NPR1 is a nuclear leucine-rich repeat protein with a key role in SAR. It binds specifically to salicylic acid, and acts as a transcriptional coregulator of SAR activators and an inhibitor of transcriptional repressors. The proteins encoded by Suppressor of NPR1, Constitutive (SNC1) and Suppressor of NPR1, Inducible (SNI1) interact with NPR1 to regulate the expression of pathogenesis-related genes. The Arabidopsis thaliana (L.) Heynh. snc1 mutant exhibits a constitutive resistance response, but in the sni1 mutant, the SNI1 protein is rendered incapable of suppressing pathogen resistance genes. To study the influence of SNC1 and SNI1 on resistance to the soybean cyst nematode (Heterodera glycines), soybean (Glycine max (L.) Merr.) roots were separately transformed with four constructs designed to: (i) overexpress GmSNC1, the soybean orthologue of AtSNC1; (ii) overexpress AtSNI1; (iii) silence GmSNC1 and (iv) silence GmSNI1. A significant reduction of the female nematode population was observed in Treatments (i) and (iv). The expression of SAR marker genes was analysed in these treatments. The unusual pattern of expression of pathogen resistance genes shows there are differences in the effect resistance genes have on soybean and A. thaliana. Although NPR1 is involved in the cross-talk between the salicylic acid, jasmonic acid and ethylene pathways, understanding the nematode resistance mechanism in plants is still imprecise. These results provide further insights into the soybean defence response.

Distinct Copy Number, Coding Sequence, and Locus Methylation Patterns Underlie Rhg1-Mediated Soybean Resistance to Soybean Cyst Nematode

Copy number variation of kilobase-scale genomic DNA segments, beyond presence/absence polymorphisms, can be an important driver of adaptive traits. Resistance to Heterodera glycines (Rhg1) is a widely utilized quantitative trait locus that makes the strongest known contribution to resistance against soybean cyst nematode (SCN), Heterodera glycines, the most damaging pathogen of soybean (Glycine max). Rhg1 was recently discovered to be a complex locus at which resistance-conferring haplotypes carry up to 10 tandem repeat copies of a 31-kb DNA segment, and three disparate genes present on each repeat contribute to SCN resistance. Here, we use whole-genome sequencing, fiber-FISH (fluorescence in situ hybridization), and other methods to discover the genetic variation at Rhg1 across 41 diverse soybean accessions. Based on copy number variation, transcript abundance, nucleic acid polymorphisms, and differentially methylated DNA regions, we find that SCN resistance is associated with multicopy Rhg1 haplotypes that form two distinct groups. The tested high-copy-number Rhg1 accessions, including plant introduction (PI) 88788, contain a flexible number of copies (seven to 10) of the 31-kb Rhg1 repeat. The identified low-copy-number Rhg1 group, including PI 548402 (Peking) and PI 437654, contains three copies of the Rhg1 repeat and a newly identified allele of Glyma18g02590 (a predicted α-SNAP [α-soluble N-ethylmaleimide-sensitive factor attachment protein]). There is strong evidence for a shared origin of the two resistance-conferring multicopy Rhg1 groups and subsequent independent evolution. Differentially methylated DNA regions also were identified within Rhg1 that correlate with SCN resistance. These data provide insights into copy number variation of multigene segments, using as the example a disease resistance trait of high economic importance.

Syntaxin 31 functions in Glycine max resistance to the plant parasitic nematode Heterodera glycines

A Glycine max syntaxin 31 homolog (Gm-SYP38) was identified as being expressed in nematode-induced feeding structures known as syncytia undergoing an incompatible interaction with the plant parasitic nematode Heterodera glycines. The observed Gm-SYP38 expression was consistent with prior gene expression analyses that identified the alpha soluble NSF attachment protein (Gm-α-SNAP) resistance gene because homologs of these genes physically interact and function together in other genetic systems. Syntaxin 31 is a protein that resides on the cis face of the Golgi apparatus and binds α-SNAP-like proteins, but has no known role in resistance. Experiments presented here show Gm-α-SNAP overexpression induces Gm-SYP38 transcription. Overexpression of Gm-SYP38 rescues G. max [Williams 82/PI 518671], genetically rhg1 (-/-), by suppressing H. glycines parasitism. In contrast, Gm-SYP38 RNAi in the rhg1 (+/+) genotype G. max [Peking/PI 548402] increases susceptibility. Gm-α-SNAP and Gm-SYP38 overexpression induce the transcriptional activity of the cytoplasmic receptor-like kinase BOTRYTIS INDUCED KINASE 1 (Gm-BIK1-6) which is a family of defense proteins known to anchor to membranes through a 5' MGXXXS/T(R) N-myristoylation sequence. Gm-BIK1-6 had been identified previously by RNA-seq experiments as expressed in syncytia undergoing an incompatible reaction. Gm-BIK1-6 overexpression rescues the resistant phenotype. In contrast, Gm-BIK1-6 RNAi increases parasitism. The analysis demonstrates a role for syntaxin 31-like genes in resistance that until now was not known.

Quantitative field testing Heterodera glycines from metagenomic DNA samples isolated directly from soil under agronomic production

A quantitative PCR procedure targeting the Heterodera glycines ortholog of the Caenorhabditis elegans uncoordinated-78 gene was developed. The procedure estimated the quantity of H. glycines from metagenomic DNA samples isolated directly from field soil under agronomic production. The estimation of H. glycines quantity was determined in soil samples having other soil dwelling plant parasitic nematodes including Hoplolaimus, predatory nematodes including Mononchus, free-living nematodes and biomass. The methodology provides a framework for molecular diagnostics of nematodes from metagenomic DNA isolated directly from field soil.

The activation and suppression of plant innate immunity by parasitic nematodes

Plant-parasitic nematodes engage in prolonged and intimate relationships with their host plants, often involving complex alterations in host cell morphology and function. It is puzzling how nematodes can achieve this, seemingly without activating the innate immune system of their hosts. Secretions released by infective juvenile nematodes are thought to be crucial for host invasion, for nematode migration inside plants, and for feeding on host cells. In the past, much of the research focused on the manipulation of developmental pathways in host plants by plant-parasitic nematodes. However, recent findings demonstrate that plant-parasitic nematodes also deliver effectors into the apoplast and cytoplasm of host cells to suppress plant defense responses. In this review, we describe the current insights in the molecular and cellular mechanisms underlying the activation and suppression of host innate immunity by plant-parasitic nematodes along seven critical evolutionary and developmental transitions in plant parasitism.

Characterization of the Pinus massoniana transcriptional response to Bursaphelenchus xylophilus infection using suppression subtractive hybridization

Pine wilt disease (PWD) caused by pine wood nematode (PWN), Bursaphelenchus xylophilus, is the most destructive diseases of pine and poses a threat of serious economic losses worldwide. Although several of the mechanisms involved in disease progression have been discovered, the molecular response of Pinus massoniana to PWN infection has not been explored. We constructed four subtractive suppression hybridization cDNA libraries by taking time-course samples from PWN-inoculated Masson pine trees. One-hundred forty-four significantly differentially expressed sequence tags (ESTs) were identified, and 124 high-quality sequences with transcriptional features were selected for gene ontology (GO) and individual gene analyses. There were marked differences in the types of transcripts, as well as in the timing and levels of transcript expression in the pine trees following PWN inoculation. Genes involved in signal transduction, transcription and translation and secondary metabolism were highly expressed after 24 h and 72 h, while stress response genes were highly expressed only after 72 h. Certain transcripts responding to PWN infection were discriminative; pathogenesis and cell wall-related genes were more abundant, while detoxification or redox process-related genes were less abundant. This study provides new insights into the molecular mechanisms that control the biochemical and physiological responses of pine trees to PWN infection, particularly during the initial stage of infection.

Engineered resistance and hypersusceptibility through functional metabolic studies of 100 genes in soybean to its major pathogen, the soybean cyst nematode

During pathogen attack, the host plant induces genes to ward off the pathogen while the pathogen often produces effector proteins to increase susceptibility of the host. Gene expression studies of syncytia formed in soybean root by soybean cyst nematode (Heterodera glycines) identified many genes altered in expression in resistant and susceptible roots. However, it is difficult to assess the role and impact of these genes on resistance using gene expression patterns alone. We selected 100 soybean genes from published microarray studies and individually overexpressed them in soybean roots to determine their impact on cyst nematode development. Nine genes reduced the number of mature females by more than 50 % when overexpressed, including genes encoding ascorbate peroxidase, β-1,4-endoglucanase, short chain dehydrogenase, lipase, DREPP membrane protein, calmodulin, and three proteins of unknown function. One gene encoding a serine hydroxymethyltransferase decreased the number of mature cyst nematode females by 45 % and is located at the Rhg4 locus. Four genes increased the number of mature cyst nematode females by more than 200 %, while thirteen others increased the number of mature cyst nematode females by more than 150 %. Our data support a role for auxin and ethylene in susceptibility of soybean to cyst nematodes. These studies highlight the contrasting gene sets induced by host and nematode during infection and provide new insights into the interactions between host and pathogen at the molecular level. Overexpression of some of these genes result in a greater decrease in the number of cysts formed than recognized soybean cyst nematode resistance loci.

Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean

The rhg1-b allele of soybean is widely used for resistance against soybean cyst nematode (SCN), the most economically damaging pathogen of soybeans in the United States. Gene silencing showed that genes in a 31-kilobase segment at rhg1-b, encoding an amino acid transporter, an α-SNAP protein, and a WI12 (wound-inducible domain) protein, each contribute to resistance. There is one copy of the 31-kilobase segment per haploid genome in susceptible varieties, but 10 tandem copies are present in an rhg1-b haplotype. Overexpression of the individual genes in roots was ineffective, but overexpression of the genes together conferred enhanced SCN resistance. Hence, SCN resistance mediated by the soybean quantitative trait locus Rhg1 is conferred by copy number variation that increases the expression of a set of dissimilar genes in a repeated multigene segment.