Novel resistance strategies to soybean cyst nematode (SCN) in wild soybean

Soybean cyst nematode (SCN, Heterodera glycine Ichinohe) is the most damaging soybean pest worldwide and management of SCN remains challenging. The current SCN resistant soybean cultivars, mainly developed from the cultivated soybean gene pool, are losing resistance due to SCN race shifts. The domestication process and modern breeding practices of soybean cultivars often involve strong selection for desired agronomic traits, and thus, decreased genetic variation in modern cultivars, which consequently resulted in limited sources of SCN resistance. Wild soybean (Glycine soja) is the wild ancestor of cultivated soybean (Glycine max) and it’s gene pool is indisputably more diverse than G. max. Our aim is to identify novel resistant genetic resources from wild soybean for the development of new SCN resistant cultivars. In this study, resistance response to HG type 2.5.7 (race 5) of SCN was investigated in a newly identified SCN resistant ecotype, NRS100. To understand the resistance mechanism in this ecotype, we compared RNA seq-based transcriptomes of NRS100 with two SCN-susceptible accessions of G. soja and G. max, as well as an extensively studied SCN resistant cultivar, Peking, under both control and nematode J2-treated conditions. The proposed mechanisms of resistance in NRS100 includes the suppression of the jasmonic acid (JA) signaling pathway in order to allow for salicylic acid (SA) signaling-activated resistance response and polyamine synthesis to promote structural integrity of root cell walls. Our study identifies a set of novel candidate genes and associated pathways involved in SCN resistance and the finding provides insight into the mechanism of SCN resistance in wild soybean, advancing the understanding of resistance and the use of wild soybean-sourced resistance for soybean improvement.

The expression patterns of rhg1 and Rhg4 in NRS100. We analyzed the expression patterns of rhg1 and Rhg4, which are the well-studied loci conferring SCN resistance in cultivated soybean PI88788 and Peking 17,19,48,49 . The rhg1 locus differs between Peking-type (rhg1-a) and PI88788-type (rhg1-b), but shares the AAT (amino acid transporter) gene (Glyma.18G022400, alternative ID: Glyma18g02580) and SNAP gene. Peking-type resistance is conferred by SNAP gene at rhg1 and SHMT at Rhg4 each, and respectively 18,20 . Relative FPKM (Fragments Per Kilobase of transcript per Million mapped reads) of SHMT at Rhg4, AAT at rhg1, and SNAP at rhg1 for all four studied genotypes at control and SCN treatment conditions revealed different expression profiles (Fig. 1). Expression of rhg1-b genes in PI88788-type resistance was not significant in NRS100. A principal component analysis of expression of SHMT and SNAP, confirmed distinct expression profiles (Fig. 1B). Noticeably, NRS100 did not exhibit any significant difference in expression at rhg1 or Rhg4 from the SCN susceptible wild soybean, S-soja, or cultivar, Williams 82.
Analysis of the transcript sequence of SNAP at rhg1 and SHMT at Rhg4 revealed Peking specific variations not found in the other three genotypes. In SNAP, Peking has two non-synonymous SNPs, one at 1,634,660 bp on exon 6 (G/C), resulting in an amino acid change from Aspartic acid to Glutamic acid. The other at 1,645,409 bp on exon 9 (A/C), resulting in an amino acid change from Threonine to Asparagine (Fig. S2). In addition, a G. soja-specific synonymous variation exists on at 1,641,767 bp at exon 2 (G/A). Another synonymous variation exists at 1,641,789 on exon 2, specific to NRS100 (G/A) (Fig. S2). The only non-synonymous variations in SNAP are found in Peking, while NRS100 has a similar susceptible-type SNAP gene similar to Williams 82. A similar trend was seen in the SHMT transcript sequence. A synonymous SNP was found in NRS100 at 8,361,269 bp (C/T) (Fig. S3). A non-synonymous variation was found in Peking at 8,361,924 bp (A/T), resulting in an amino acid change from Asparagine to Tyrosine (Fig. S3), as well as distinct variation in mapping from the reference Scientific Reports | (2021) 11:7967 | https://doi.org/10.1038/s41598-021-86793-z www.nature.com/scientificreports/ (Williams 82) in this region, which was not found in wild soybean NRS100 or S-soja. All sequence variations that resulted in amino acid changes were found only in Peking-type SNAP and SHMT sequences. The resistance conferring gene of rhg1-a in Peking-type resistance, SNAP gene, is located at chromosome 18. Paralogous SNAP genes were analyzed to determine if resistance in NRS100 was conferred by different members of the SNAP gene family. No significant expression or fold change is found on SNAP 2, 9, 11, or 14 by RT-qPCR or RNA-sequencing, as if found in resistant soybean cultivars using SNAP for resistance 50 . The results showed that fold changes for all NRS100 genes < 0.4 between SCN treatment vs control (Fig. S4), which suggested that NRS100 does not have significant differential expression at SHMT, SNAP 18 or a SNAP paralog.
Expression patterns of previously identified candidate genes involved in SCN resistance in wild soybean. Two major QTLs, cqSCN-006 on chromosome 15 and cqSCN-007 on chromosome 18, have been previously identified involved in SCN resistance to race 5 in wild soybean genotype PI468916 [51][52][53] . Interestingly, none of the candidate genes on cqSCN-006 were found significantly upregulated in NRS100, including the y-SNAP gene (Glyma.15g191200), while one candidate gene, Apetala 2 (AP2), on cqSCN-007 was significantly upregulated in NRS100. This region also collocates with SCN resistance loci qSCN18 in soybean landrace PI 567516C 54 . AP2 is a, transcription factor (Glyma.18g244600) 53 , known to regulate multiple developmental pathways and play a role in abiotic and biotic stress response 55 .
A previous RNA-sequencing study of a SCN race 5 resistant wild soybean genotype S54, revealed potential pathways involved in resistance in this genotype 46,47 . Expression patterns of S54 and NRS100 were compared to (C) Significantly up and down regulated genes by genotype. DEGs considered significant if q value < 0.01, log 2 foldchange no less than ± 0.6. (D) Heatmap of all DEGs by log 2 FC (fold change) of induced response to SCN treatment of all genotypes. Cutoff Log 2 FC no less than ± 2 for NRS100, including 2076 genes. Green bar indicates NRS100 specific up regulated genes in response to SCN stress. Red bar indicates NRS100 specific down regulated genes in response to SCN stress. Blue bar indicates Peking specific up regulated genes in response to SCN stress. Orange bar indicates Peking specific down regulated genes in response to SCN stress. www.nature.com/scientificreports/ determine if the mechanism of resistance was the same. Many differences between NRS100 and S54 were found, similarities between these genotypes were limited to the following. Previously identified Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) differentially expressed genes (DEGs) in G. soja s54 SCN response, Brassinosteroid insensitive 1 (BRI1)-associated kinases (BAK1) (Glyma.05G119500, Glyma.05G119600) genes and BAK-1interacting receptor kinase 1 (SOB1R1) genes (Glyma.04G190400, Glyma.06G175100) were also induced in response to SCN in NRS100, but not exclusively to this genotype 46 . Two genes (Glyma.02G270700, Glyma.15G111300) of chitin elicitor receptor kinase 1 (CERK1) were down and upregulated respectively in NRS100, a similar patten was seen in S54 SCN race 5 response. Lectin receptor kinase (Glyma.07G135400) was also induced in NRS100, which was also found in previously studied S54 46 . In addition, significant downregulation of NBS-LRR gene (Glyma.16G209000) and upregulation of NBS-LRR gene (Glyma.17G180000) was found in NRS100, like in S54. Significant upregulated DEGs associated with calmodulin binding (Glyma.05G237200, Glyma.07G093900, Glyma.08G044400) were upregulated in NRS100 and S54. Calcium transport genes: autoinhibited Ca2+-ATPase 9 (ACA9, Glyma.07G004300) and protein calcium exchanger 7 (CAX7, Glyma.19G066700) were upregulated in both NRS100 and S54. MAPKKKs (Glyma.17G245300 and Glyma.05G094400) were significantly up and down regulated in both S54 and NRS100, however, the strongest DEG in the MAPK cascade found in S54 (Glyma.15G048500), was not found to exhibit the same response in NRS100. Jasmonic acid (JA) pathway gene, oxophytodienoate-reductase 3 (DDE1/OPR3, Glyma.13G109800) showed a similar induced effect in NRS100 as in S54, however this was the only one gene of the JA pathway genes to show this similarity between the two genotypes. Members of WRKY40 transcription factor (Glyma.13G370100 and Glyma.17G222500) were found highly induced in NRS100 and S54, although much higher in S54. No similarities were found in the upregulated chitinase genes in S54, however, members of Chitinase-like protein (CLT), Glyma.09G038500 and Glyma.15G143600, were found down regulated in both NRS100 and S54 46 . Overall, similarities between expression in NRS100 and S54 are limited, suggesting a different regulatory mechanism of resistance between NRS100 and S54.

Induced gene expression in response to SCN infection in NRS100.
The entire RNA-seq based transcriptome comparison of NRS100 with the other genotypes, including susceptible wild soybean genotype S-soja, as well as resistant cultivar Peking and susceptible cultivar William82, was conducted to understand the unique resistance mechanism in NRS100. Illumina sequencing generated 17.9-29.7 million raw reads per library. Filtering by quality of reads resulted in 93.8-95% high quality (quality score > 20) reads for alignment. The mapping rate of quality-controlled reads ranged from 84.9 to 92.8% of all 24 libraries to G. max Wm82.a2.v1 (Table S1). On average, 56, 710 expressed genes were found in root tissue from each genotype. The amount of significantly (q < 0.01) differentially expressed genes between control and treatment groups for each genotype were 4641, 2911, 498, and 663 in NRS100, S-soja, Peking, and Williams 82, respectively.
To determine the genes associated with induced resistance response to SCN in NRS100, DEGs of four genotypes between control and SCN treatment were compared. A count of the significant DEGs of all four genotypes, NRS100, S-soja, Peking, Williams 82, in response to SCN infection revealed a greater quantity of expressed genes in NRS100 compared to other genotypes, with 1602 and 3039 genes significantly up and down regulated respectively (Fig. 1C). The SCN susceptible wild soybean (S-soja) has the second most substantial response to SCN infection. Both cultivars, Peking and William 82, exhibit less extensive and more tailored response to SCN infection, with Peking exhibiting a lower number of down regulated genes. This trend is further supported by the expression profiles of highly-significantly expressed genes, with a cutoff Log 2 FC no less than ± 2 for NRS100 ( Fig. 1D).
NRS100 is seen to have a different set of genes induced in response to SCN infection, (green and red bars ( Fig. 1D), when compared to other three genotypes. By contrast, the induced response to SCN infection in Peking (blue and orange bars (Fig. 1D), is controlled by a smaller set of genes, separate from those induced in NRS100. Both shared responses by all genotypes, and species-specific shared responses were seen. GO enrichment analysis of resistance specific, shared response, and stress response genes was done after visualization of all DEG's in a Venn-diagram 56 (Fig. S5).
To elucidate the expression profile associated with resistance and eliminate generalized stress responses, species-specific induced responses were compared ( Fig. 2A,B). Shared induced responses by genotype were eliminated to create a subset of resistance-specific induced DEGs (Fig. 2C). Again, resistance-specific expression in NRS100 is at a much larger scale than that in Peking, consisting of expression profile more than ten-fold to Peking. Comparison of resistance-specific expression profiles of Peking and NRS100 reveal a shared response consisting of 69 genes, contributing to only 1.8% of the overall resistance-specific profile of NRS100. GO functions associated with NRS100 resistance include up regulation of systemic acquired resistance mechanisms, salicylic acid signaling pathway, and regulation of defense response mechanisms (Fig. S5, Tables S2, S3) and down regulation of secondary cell wall biogenesis and light harvesting by photosynthesis (Tables S3, S4).

Constitutive gene expression in NRS100.
To further analyze the potential mechanism of SCN resistance in NRS100, the constitutive expression profile specific to NRS100 was considered due to its important role in plant-pathogen protection [57][58][59][60][61] . To determine the constitutively expressed genes associated with SCN resistance in NRS100, the basal regulation of NRS100 was compared to S-soja and Peking at control conditions. This eliminates specific-specific regulation shared with wild soybeans and resistance-specific regulation shared with cultivars. The significantly induced expression profile from SCN infection in NRS100 was also eliminated from the basally regulated genes. This method ensures that only the genes highly expressed in the control condition in NRS100 specifically, that also did not significantly alter their expression after SCN infection, were captured. Comparison of all expression profiles results in a set of 1656 genes that are constitutively expressed in  (Table S5).
Expression patterns in NRS100-specific SCN resistance pathways. Based on the above findings, putative metabolic pathways conferring resistance were selected based on exclusive constitutive and induced responses, importance in GO enrichment, and previous understanding of plant-pest interactions. The jasmonic acid (JA) signaling pathway is a well-studied pathway, associated with plant-pathogen interactions, and known to have antagonistic effects on the salicylic acid (SA) signaling pathway, and therefore was analyzed further to determine if signal transduction plays a role in SCN resistance in NRS100 [64][65][66][67][68] . Within the JA signaling pathway, NRS100-specific downregulation occurred at JAR1 (jasmonate resistant 1) genes and COI1 (coronatine-insensitive protein 1) gene. JAR1 genes, Glyma.03G256200 and Glyma.16G026900, were exclusively down regulated in NRS100 in response to SCN. Downstream of JAR1, COI1 gene, Glyma.18G030200, was downregulated in NRS100 exclusively in response to nematodes. No other genotypes had a significant induced response at these genes. Regulation at known genes in JAZ (jasmonate ZIM domain-containing protein) and MYC2 transcription factor (Glyma.08G271900) were not induced in NRS100 but were highly expressed in control conditions, or constitutively expressed significantly higher than other genotypes (Fig. S8). Expression at these significant genes in the JA signaling pathway were confirmed in the RT-qPCR results with representatives for JAR1 and JAZ due to modulation of these proteins by multiple genes (Table S6). The polyamine biosynthesis pathway, specifically at spermidine synthesis, was analyzed in further detail to determine if NRS100 uses a polyamine related defense response known to play a role in resistance to other is significantly upregulated in response to nematode infection in both Peking and NRS100. There are two paralogous spermidine synthase genes, SPDS17 (Glyma.17G091123) and SPDS5 (Glyma.05G036300). SPDS5 is significantly upregulated in both NRS100 and Peking. However, SPDS17 is significantly upregulated in response to nematode infection in NRS100 only. These expression results are confirmed by both RNA-sequencing and RT-qPCR results. Downstream of spermidine synthesis, expression of PR-1 (pathogenesis-related protein 1) genes is highly constitutively expressed in NRS100 at Glyma.15G062300, and non-specifically significantly induced at Glyma.15G062400 and Glyma.15G062500. All three PR-1 genes are significantly induced in response to SCN infection in S-soja, but not in either cultivar. Species-specific response at PR-1 is found, with NRS100 specific constitutive expression of Glyma.15G062300 and overall higher level of expression when induced in NRS100. The polyamine related defense response includes direct and indirect resistance mechanisms, where indirect involves the production of an unknown secondary metabolite.

Validation of expression by RT-qPCR.
RNA-sequencing expression of 20 genes (Table S6), including randomly selected and selected from important pathways detailed in PATHVIEW results (Fig. S8), was validated using RT-qPCR for all genotypes and condition replicates. RT-qPCR results of the 20 selected genes significantly correlated with RNA-sequencing results. An expected indirect correlation (R = 0.621) of RT-qPCR derived ΔCT to RNA-sequencing derived FPKM of all biological replicates and averaged technical replicates was found (Fig. 3). In addition, fold changes found in RT-qPCR significantly correlates (R = 0.762) with that from RNA-sequencing data (Fig. S9). RT-qPCR was also used to better understand the expression of genes of interest, including those within rhg1, Rhg4, SNAP paralogs, the JA signaling pathway, and the Polyamine defense response.

Discussion
SCN is the most destructive pest to soybean, and is rapidly evolving enhanced and complex virulence against resistant crops 2,4,9,71 . The majority of efforts to understand SCN resistance mechanisms is focused on PI88788 and Peking type resistance, which use rhg1-b or a combination of rhg1-a and Rhg4 [15][16][17][18][19][20]72 . Little effort has been placed on sources of resistance loci outside rhg1 or Rhg4, but is necessary to combat current and evolving SCN populations. Due to the long history of environmental adaptation in natural populations, as well as the fact that the resistant Rhg4 allele emerged after domestication, it is likely that the wild soybean has different strategies to battle SCN stress 22,23 . However, research to understand the mechanisms of resistance in wild G. soja is lacking, and only a handful of studies describe potential mechanisms [44][45][46][47][51][52][53]73,74 .
The molecular mechanisms found in cultivated soybean cultivars could not explain NRS100 resistance. Our results showed that NRS100 uses a resistance mechanism independent of Peking-type or PI88788-type resistance: first, the genes at rhg1-b, rhg1-a, and Rhg4, were not induced in response to SCN infection in NRS100 (Fig. 1A). Second, analysis of the molecular sequence of exons in SNAP and SHMT revealed that non-synonymous mutations found in Peking were not found in any other genotypes in our study, including www.nature.com/scientificreports/ NRS100. This confirmed that the Peking-type variation of these resistance genes were not found in NRS100. It is known that resistance to SCN is conferred by a complex mechanism in the cultivated soybean, which consists of variation in copy number and genetic composition of genes at rhg1 and Rhg4 and associated promoters [15][16][17][18][19][20] . Our results suggest NRS100 does not make use of rhg1 and Rhg4 based on transcription of these genes constitutively and induced, and molecular sequence of the coding region, but copy number of these genes and composition of their promotors are yet to be determined. It has been suggested that the resistance mechanism involvig SNAP in Peking is likely due to interactions of the SNAP product with proteins within the SCN secretion at the point of infection, such as Hg-SLP-1 72,75,76 . Due to the importance of SNAP in SCN resistance, the SNAP family was analyzed further to determine if another member can help confer resistance by a similar mechanism. Previous characterization of the SNAP family revealed that SNAP 11 contributes to additive resistance to SCN in Peking-type resistance, however our results suggest that none of the SNAP members in the family (SNAP2, SNAP9, SNAP11, SNAP14, and rhg1 SNAP18) are significantly highly expressed during SCN infection in NRS100. In addition, NSF gene, Glyma.07g195900, previously shown to be associated with resistant SNAP at rhg1, is not significantly induced or highly constitutively expressed in NRS100 50,72 . Thus, we suggest that a novel mechanism, different from the previously identified rhg1, Rhg4, and SNAP family, confer resistance in wild soybean NRS100.
Genotype-specific SCN resistance in wild soybean. One of the most well studied wild soybean G. soja genotype for SCN resistance, thus far, is PI468916. The fine mapping studies revealed two candidate loci, cqSCN-006 on chromosome 15 and cqSCN-007 on chromosome 18 51-53 are the two major QTLs conferring SCN resistance to race 5. None of the candidate genes on locus cqSCN-006 were found significantly upregulated in NRS100, including the y-SNAP gene (Glyma.15g191200). However, one candidate gene, Apetala 2 (AP2), on locus cqSCN-007 was found significantly upregulated in NRS100. AP2 is a transcription factor that is related to flower development and stress response in Arabidopsis (Arabidopsis thaliana) 53,55 . AP2 can be activated by either JA or SA pathways, which are established resistance-critical pathways [77][78][79][80] . It is plausible that cqSCN-007 might play a role in SCN resistance in NRS100 at AP2 specifically, however differences in regulation of other candidate genes between PI468916 and NRS100 indicates a genotype-specific SCN resistance mechanism in each ecotype.
Another well studied wild soybean genotype for SCN resistance is S54 44 . Similar responses to SCN race 5 in NRS100 and S54 were limited to RLK-LRRs, Non-LRR domain RLKs, calmodulin binding genes, calcium transport genes, WRKY40, a JA pathway associated gene, MAPKKKs, and downregulated CLTs 44 . Moreover, the expression profile found in NRS100 did not include the important genes involved in the proposed resistance mechanism in S54 46 , which indicates that NRS100 might use a different strategy to confer SCN resistance to race 5. It is not surprising that genotype-specific resistance is found in wild soybean, since the two well-studied soybean cultivars, Peking and PI88788, also have genotype-specific resistance mechanisms to SCN.
Novel mechanism of SCN resistance in NRS100. NRS100 displays a strikingly different expression profile with and without SCN inoculation, compared to SCN resistant soybean cultivar Peking, indicating different mechanisms employed in resistance. The amount of DEGs in response to SCN is much greater in NRS100, with twice as many DEGs than the response of susceptible S-soja, and almost tenfold the amount found in Peking or Williams 82. This strong response in NRS100 suggests that resistance might be more complex, and less tailored than that found in the cultivated Peking-type resistance. NRS100 and Peking have a limited number of shared responsive genes to SCN, consisting of 39 up-and 30 down-regulated genes, which is much smaller amount compared to the 1031 up-and 2604 down-regulated genes specific to NRS100. There might be a regulatory network specific to NRS100 that is non-existent in the SCN resistant cultivar, Peking, or either susceptible genotype, which has not been directly demonstrated in any other G. soja ecotype.
In addition to the clearly induced effect, constitutively expressed genes, or genes that are highly upregulated in control conditions of NRS100 compared to the other genotypes, suggest a resistance mechanism that is already turned on without SCN infection. Constitutively expressed genes are known to play a significant role in plant-pathogen resistance, which can promote toxin production, increase the sensitivity to pathogen stimuli, or support an inducible immunity [57][58][59][60][61] . Intracellular signal transduction is highly upregulated in the control NRS100 roots, suggesting a higher sensitivity to pathogen trigger via signal transducers. Ethylene biosynthesis, another signaling molecule, is also highly expressed, which was reported associated with induced resistance to root knot nematode and SCN, but conversely found to increase susceptibility to SCN in a separate study [81][82][83][84][85][86] . To explain the unique resistance found in NRS100, we hypothesized two modes of resistance in NRS100, (1) induced immunity by suppression of JA signaling pathway, and (2) spermidine synthesis mediated induced defense against SCN growth. Jasmonate acids (JA) suppression in induced immunity. Signaling pathways are widely studied for their role in pest resistance in plants. Signaling pathways, including JA, SA, and ethylene act as detection system to decipher pest and pathogen signals in order to allow a plant to mount a defense 68,87 . High levels of jasmonate acids (JA) are associated with increase wound signaling to insect pests by inducing an effective defense against pests, and induced systemic defense in root knot nematode by PR1 65,85,88,89 . However, JA signaling is antagonistic with SA signaling 87,90 . It has been shown in tomato (Solanum lycopersicum) that JA signaling actually has a negative impact on resistance to aphids and root knot nematode, and JA-deficient mutants exhibit elevated resistance 91,92 . In the interaction between tomato and aphids, susceptibility can be restored by suppressing the SA signaling pathway, demonstrating the direct interplay between the JA and SA signaling pathways 90 . A separate example exists in rice where modification of JA signaling increases SA levels and in turn leads to increased resistance to multiple insect herbivory 93 . It is likely that JA is important for broad resistance and wounding Scientific Reports | (2021) 11:7967 | https://doi.org/10.1038/s41598-021-86793-z www.nature.com/scientificreports/ response, but can be costly due to its negative interaction with SA signaling which is needed for constitutive and inducible resistance mechanisms 67,[94][95][96] . In addition, suppression of the JA signaling pathway can promote plant growth via activation of PIF regulated gibberellin signaling 97 . It is also possible that JA signaling is suppressed by a deactivation of NRS100 defense mechanisms by a decoy produced by the nematode, similar to that found in the tobacco-Helicoverpa zea larvae interaction 67,98 . However, this is unlikely considering this suppression is only found in NRS100 and not the SCN susceptible genotypes. With suppression of JA, SA can accumulate. SA accumulation has been found to be important in specialized defense response in Arabidopsis 99 . In NRS100, JA signaling is suppressed at JAR1 and COI1, important regulators in the pathway 89,100,101 . This suppression reduces the ability of COI1 and JAZ to interact physically, and does not allow a release of JA mediated transcription factors 89 . Upregulation of JAZ and MYC2 act as repressors to the JA pathway 102,103 . Suppression of JA allows SA signaling to induce a defense response, avoiding potential conflicts with SCN-produced decoys to mediate traditional JA signaling in SCN defense response (Fig. 4). Interestingly, the downregulation of JAR1 and COI1 is induced in NRS100 only, but the induced response at JAZ is found in both Peking and NRS100 but regulated by separate genes. This suggests that JA signaling is not only regulated differently at key steps (COI1 and JAR1) in NRS100 but is regulated by different genes at points of similar expression of with Peking (JAZ). It is likely that many different JAZ products exist, especially in the case of regulation by separate genes, which is expected to occur in NRS100 67,96 . Spermidine synthesis mediated induced defense against SCN growth. We hypothesize that spermidine synthesis and secondary metabolite production induces defense against SCN growth and syncytium formation in NRS100. Polyamine metabolism, including spermidine synthesis, is often altered in abiotic and biotic stress responses in plants 104,105 . Spermidine synthesis upregulation can result in accumulation of free spermidine or lead to the production of polyamine derived molecules. This can trigger H 2 O 2 regulated hypersensitive response, cooperative effects with SA signaling to induce PR proteins or other modes of pathogen resistance, lead to production of secondary metabolites like nicotine, and reinforce/maintain structural integrity of cell walls to hinder pathogen invasion 62,63,69,70,[106][107][108][109][110] . In addition, higher expression of spermine synthesis, directly following spermidine synthesis, has been shown to induce pathogen resistance 63,69 . NRS100 exhibits non-exclusive upregulation of spermine synthesis, downstream of exclusively upregulated SPDS. Therefore, accumulation of spermine and spermidine in response to SCN is likely. Conversely, it has been found that upregulation of SPDS can be induced by effector proteins secreted by nematodes, and therefore leads to a compatible plantpathogen interaction 111 . Spermidine synthesis gene, SPDS17, is found to be significantly upregulated in NRS100, www.nature.com/scientificreports/ a response that is specific to the genotype. Interestingly, a separate SPDS gene on chromosome 15 is also highly upregulated in both Peking and NRS100, but not to the same extent as SPDS17 is in NRS100. In addition, upregulation of precursors, at putrescine synthesis (ODC1) is upregulated in both NRS100 and Peking but at a much higher level in NRS100. Overexpression of ODC1 has been shown to lead to elevated levels of putrescine but not spermidine, meaning that the upregulation of spermidine synthesis in NRS100 might activate an independent response and production of spermidine, regardless of the shared ODC1 upregulation with Peking 112,113 .
Although polyamine synthesis can be related to both compatible and incompatible interactions with pathogens, we propose that SPDS17-regulated spermidine synthesis in NRS100 might have a highly incompatible interaction with SCN 104 . The outcome impedes the nematodes' ability to expand a feeding site, leading to decreased nutrition and death.

Conclusion
These results indicated not only species-specific, but also genotype-specific mechanisms involved in SCN resistance. Here, we showed a distinct expression profile unique to a newly identified SCN-resistant wild soybean genotype, NRS100, in response to SCN and proposed candidate mechanisms associated with SCN resistance. The experiment design made use of both wild and cultivated soybeans to identify novel resistance mechanisms. The identification of candidate genes and pathways involved in SCN resistance advances the understanding of molecular mechanisms of plant and pest interactions. Prospectively, novel genes involved in SCN resistance can be incorporated into the cultivated soybean with diverse strategies, such as marker assisted breeding, gene transformation, genome editing technology, toward the long-term goal of developing new and diverse SCN resistant soybean cultivars, which has crucial significance to agriculture and environmental sustainability.

Materials and methods
Validation of NRS100 SCN resistance and experiment design. Our previous study showed high resistance of NRS100 to SCN race 5 44,45 . To verify this finding, we rescreened the NRS100 using race 5, along with indicator lines (Peking, PI88788, PI90763, PI437654, PI209332, PI89772, and PI548316), employed the same protocol described in Zhang et al. 44,45 . We also selected a susceptible wild soybean genotype, named S-soja (FI = 149) 44 , to compare the SCN-infection response differences with NRS100. Of the two main cultivated sources of resistance, Peking and PI88788, the former is the only one that is resistant to SCN race 5. Thus, we used Peking as a resistant soybean cultivar to test species-specific resistance strategies. We also used William82 as a susceptible soybean cultivar for comparison. All the study material were from USDA germplasm and this study comply with relevant institutional, national, and international guidelines. SCN hatching and stress treatment experiment. SCN race 5 were reared on soybean cv. Williams 82 plants in the greenhouse under optimal conditions at 27 °C and 15 h light for over 30 generations. Female SCN cysts were harvested from the roots using a nested sieve collection (850 and 250 µm) method, and females released from cysts by pressure into a 25 µm sieve 46 . Released eggs were purified by sucrose flotation 114 and allowed to hatch on incubation trays over water. second stage juveniles (J2) were collected from the water and diluted to a concentration of 2500 eggs/ml for inoculation.
A controlled SCN stress experiment was performed on SCN-resistant wild soybean accession (NRS100), SCN-susceptible wild soybean accession (S-soja), SCN-resistant cultivated soybean accession (Peking), and SCN-susceptible cultivated soybean accession (Williams 82). Germinated soybeans were planted in cone planters (Stuewe & Sons, Tangent, Oregon, USA) in replicates of 12 per condition, treatment and control for each accession in sterilized sand with nutrients. After 2 days post transplanting, each plant was inoculated with hatched J2 of SCN race 5 (1500/plant) suspended in 0.09% agarose, or a blank (1 ml 0.09% agarose) for control groups. Optimal growing conditions were kept constant at 27 °C and 15 h light at 50% relative humidity in an environmental chamber. Roots were collected 6 days after SCN inoculation and fragmented in liquid nitrogen in preparation for RNA isolation. Each biological replicate contained pooled roots from three individual plants, providing four replicates of each group in total.
Library construction and RNA sequencing. We used the same strategy and methods in root tissue collection, RNA extraction, RNA library preparation, and RNA sequencing as described in our pervious study 46,47 . Technical and biological replicates of total RNA sequences of controlled and treatment group root tissues were obtained using the Illumina 1.9 Hi-seq platform (NEB, Beverly, MA, USA). Each replicate was on two separate lanes and merged, generating 24 total libraries in triplicates per sample group (control 1-3, treatment 1-3 per genotype). Raw fastq reads were checked for quality control with FastQC (version 0.11.6) and filtered for quality, and trimmed (quality score < 30, and adapters) using FASTx toolkit (version 0.0.13). Reads were mapped against reference genome G. max Wm82.a2.v1 115 with Tophat (version 2.1.1) and Bowtie2 (version 2.2.9) for all four genotypes, control and treatment. Assembly of transcripts of each replicate to the reference was done using Cufflinks (Version 2.2.1), a final assembly of each replicate was obtained, and differently expressed genes (DEGs) as the final step in the Cufflinks pipeline 116-119 . Comparative transcriptome analyses. In order to investigate the global gene expression changes, we compare RNA-seq-based transcriptomes of the NRS100, S-soja, Peking, and Williams 82, under both control and nematode-treated conditions. To determine induced response, differentially expressed profiles were found between control and treatment groups for each genotype. To determine constitutive expression, differentially expression profiles were compared between control conditions of each genotype.