The receptor like kinase at Rhg1-a/Rfs2 caused pleiotropic resistance to sudden death syndrome and soybean cyst nematode as a transgene by altering signaling responses
© Srour et al.; licensee BioMed Central Ltd. 2012
Received: 17 November 2011
Accepted: 12 June 2012
Published: 2 August 2012
Soybean (Glycine max (L. Merr.)) resistance to any population of Heterodera glycines (I.), or Fusarium virguliforme (Akoi, O’Donnell, Homma & Lattanzi) required a functional allele at Rhg1/Rfs2. H. glycines, the soybean cyst nematode (SCN) was an ancient, endemic, pest of soybean whereas F. virguliforme causal agent of sudden death syndrome (SDS), was a recent, regional, pest. This study examined the role of a receptor like kinase (RLK) GmRLK18-1 (gene model Glyma_18_02680 at 1,071 kbp on chromosome 18 of the genome sequence) within the Rhg1/Rfs2 locus in causing resistance to SCN and SDS.
A BAC (B73p06) encompassing the Rhg1/Rfs2 locus was sequenced from a resistant cultivar and compared to the sequences of two susceptible cultivars from which 800 SNPs were found. Sequence alignments inferred that the resistance allele was an introgressed region of about 59 kbp at the center of which the GmRLK18-1 was the most polymorphic gene and encoded protein. Analyses were made of plants that were either heterozygous at, or transgenic (and so hemizygous at a new location) with, the resistance allele of GmRLK18-1. Those plants infested with either H. glycines or F. virguliforme showed that the allele for resistance was dominant. In the absence of Rhg4 the GmRLK18-1 was sufficient to confer nearly complete resistance to both root and leaf symptoms of SDS caused by F. virguliforme and provided partial resistance to three different populations of nematodes (mature female cysts were reduced by 30–50%). In the presence of Rhg4 the plants with the transgene were nearly classed as fully resistant to SCN (females reduced to 11% of the susceptible control) as well as SDS. A reduction in the rate of early seedling root development was also shown to be caused by the resistance allele of the GmRLK18-1. Field trials of transgenic plants showed an increase in foliar susceptibility to insect herbivory.
The inference that soybean has adapted part of an existing pathogen recognition and defense cascade (H.glycines; SCN and insect herbivory) to a new pathogen (F. virguliforme; SDS) has broad implications for crop improvement. Stable resistance to many pathogens might be achieved by manipulation the genes encoding a small number of pathogen recognition proteins.
KeywordsSegregation Pleiotropy Rhg1/Rfs2 Soybean Resistance Soybean cyst nematode (SCN) Sudden death syndrome (SDS) Insect herbivory
Fusarium virguliforme (Akoi, O’Donnell, Homma & Lattanzi), causal agent of soybean (Glycine max L. Merr.) sudden death syndrome (SDS), first caused a significant disease loss in 1987 . F. virguliforme was not prevalent in Asia by 2011 but had spread quickly across the Americas from about 1980–2011. SDS has become a major pest problem for soybean growers and breeders in the Americas . The origins of the disease remain unclear but F. virguliforme may be a new pathogen of soybean since no complete resistance has been reported.
F. virguliforme, like many plant pathogenic Fusaria, were facultative hemi-biotrophic pathogens of plant roots with many host species [1, 3, 4]. However, only soybean among known hosts showed the leaf scorch when infected by F. virguliforme. Soybean cultivars showed a wide range of susceptibility to both leaf scorch and root rot suggesting cultivar-specific partial resistance existed F. virguliforme appeared to be a clonal pathogen [3, 6]. There were some variations in aggressiveness among field isolates and maintained strains but there were no races reported, by 2011.
Soybean resistance to SDS was multi-geneic and had two components; a partial resistance to root infection and rot caused directly at the site of infection by the fungus; and a partial resistance to leaf scorch caused indirectly by translocated fungal toxins [5, 7, 8]. The bases of resistance might include resistances to one or more toxins ; and both local and systemic resistances following pathogen recognition [10–12].
Heterodera glycines I., the soybean cyst nematode (SCN) was probably an ancient pest of soybean since complete resistances to some Hg Types of SCN was found in about 1% of pre-domesticated and early domesticated Plant Introductions (PIs) . Interestingly most of these PIs were also partially resistant to SDS . One locus, the Rfs2/Rhg1 locus on chromosome 18, was shown to underlie coinheritance of resistance to SDS in the roots and also reduce root infestation by SCN [10, 11, 14–17]. Fine map development did not resolve Rfs2 from Rhg1 suggesting the underlying gene(s) were either very closely linked or pleiotropic [11, 16].
Rfs2 underlay partial resistance to the spread of root infections by F. virguliforme [5, 7, 11]. The site of infection did not rot as rapidly when this allele was present and rates of root growth nearly equal to non-infested plants were maintained. Toxin translocation to leaves appeared reduced because leaf scorches did not develop or were less severe. The sudden plant death characteristic of SDS was manifested as both early senescence and an unusual abscission, basal to the leaflets instead of the petiole. Neither occurred if the Rfs2 allele was present.
Equally the Rhg1 locus underlay partial resistance to SCN [13, 18]. H.glycines, like many plant parasitic nematodes, were obligate endoparasites of plant roots. Like F. virguliforme H. glycines has many alternate hosts. Over the past 50 years, the number of Hg Types (ex. races) has expanded from 4 in the 1960’s to 16–20 [19, 20]. However, the Rhg1 locus was constant, being required for partial resistance to all Hg Types in most PIs and cultivars. Full resistance to SCN required 1–4 loci in addition to Rhg1, the number depending on the nature of the cyst population parasitizing the roots [17, 21–25]. Genetic diversity was found among SCN isolates, even inbred cyst populations like PA3, Hg type 0 [13, 26]. Further, variation among the host plant roots response to SCN has been associated with temperature  such that environmental conditions must be rigorously controlled during assays [13, 16].
The resistance or susceptible interaction(s) between the nematode and soybean affected by Rhg1 was not induced until females stopped moving through the roots and established a feeding site comprising several giant cells [28–30]. Full resistance to SCN, based on the combined action of the genes at Rhg1 and one or more additional Rhg loci, was manifest as; cell wall appositions to surround the feeding site; failure to supply the feeding site a tracheary element; and a necrosis as the feeding site develops. If the resistance allele at Rhg1 was present normal rates of root growth were slightly depressed but the above ground stunting, yellowing and early senescence did not occur.
Inheritance of resistance to SCN was first reported in the PI ‘Peking’ . Three recessive loci (rhg1 rhg3) and a dominant locus (Rhg4) were assigned gene names by parsimony though other dominance models were equally likely. The Peking derived resistance alleles of rhg1 and Rhg4 were introgressed into the cv. ‘Forrest’ [31–33]. In crosses based on Forrest and with SCN isolate PA3 the rhg1-a was shown to be codominant with the susceptibility allele of ‘Essex’ (rhg1-e) and alone capable of providing partial resistance. Consequently, the Forrest, ‘Hartwig’, Peking and ‘PI 437654’ allele was renamed to Rhg1-a ( and hereafter) by the Soybean Genetics Committee.
The Rhg1 locus was located to a sub-telomeric region of the soybean chromosome 18 (molecular linkage group G; Lg G) by many studies [17, 21, 23–25]. All of these segregating populations that were later tested with F. virguliforme also had an Rfs2- like activity against SDS [8, 15, 34]. However, an Rhg1-like locus was found at other locations in a few SCN resistant PIs, including Lg B1 (chromosome 11) , mid LgG [35, 36] and Lg B2 (chromosome 14) [37, 38]. The effects of the Rhg1-like loci found at other locations than chromosome 18 on resistance to SDS were not reported by 2012.
It has been shown the resistance allele of Rhg1-a/Rfs2 (from Peking) was associated with reduced seed yield when SCN is not present in the fields [39–41]. That phenomenon might be related to delayed seedling development and stand formation . How the root reduction contributes to resistance may involve a locus on chromosome 7 (LG M) where Rzd, an interacting allele needed for r esistance to z ygote d eath, was strictly co-inherited in phase with Rhg1/Rfs2-a . Therefore, effects on development were predicted for the gene(s) underlying Rhg1-a/Rfs2-a [16, 43–47].
Fine scale genetic maps and BAC based genomic analysis identified a 42 kbp region encompassed in BAC B73p06 as the Rhg1/Rfs2 locus [16, 45, 48]. The region encoded an RLK (Glyma18g02680 named GmRLK18-1 hereafter), a variant laccase (Glyma18g02690) and a predicted Na/H ion antiporter (Glyma18g0270). The GmRLK18-1 and Gmlaccase18-1 were expressed in roots, shoots and flowers but not nodules or seeds. The antiporter was not expressed in any organ tested in Essex or Forrest. However, it might be expressed in flowers and seed at very low abundance . Many nucleotide differences were found in the region encompassing the GmRLK18-1 and Gmlaccase18-1 genes from fragments of sequences from the resistant allele in Forrest compared to susceptible genotypes ‘Asgrow 3244’ and ‘Williams 82’ [16, 29, 48, 50, 51]. However, the GmRLK18-1 was considered the most likely candidate gene underlying Rhg1-a based on fine maps and association analyses.
RLKs are part of the eukaryotic tumor necrosis factor beta receptor super-family [44, 52]. In plants they represent a major class of resistance genes and also a major class of developmental regulators. The GmRLK18-1 gene at Rhg1/Rfs2 had nine alleles recognized among PIs and cultivars [16, 44–48, 50]. Of those four alleles were associated with partial resistance to SCN and five with susceptibility. The second most important resistance allele was from PI88788 and named rhg1-b because it was recessive and discovered second . However, a small scale experiment with RNAi to the RLK (named GmRLK18-1-b) at the rhg1-b in transgenic hairy roots did not show a large effect on cyst numbers . Either dominance or incomplete inhibition of the protein might have occurred. Here, additional experiments were undertaken with the GmRLK18-1 at Rhg1-a in stable transformed soybean lines. Here, a molecular basis for resistance to SDS and SCN was inferred from functional analyses of the Rhg1/Rfs2-a locus in near isogenic lines (NILs) and the Forrest allele of the GmRLK18-1 gene in transgenic plants.
Allelic variations at the Rhg1/Rfs2 locus
Genes predicted from DNA sequence within the Rhg1 -a/ Rfs2 locus encoded by BAC 73p06 (81, 157 bp) from Forrest (resistant) and sequence contig from A3244 (susceptible; AX196295)
SIUC_Sat_−35 (35 bp)
NADP redox coenzyme-like
SIUC_Scaa-23 (19 bp)
SIUC_Sat_001 (2 bp)
SIUC_InD + 3.5 (53 bp)
SIUC_Sca_005/Sac_013 (4 bp)
Receptor like kinase ( Rhg1 -a/ Rfs2 )
SIUC_Indel-3.5 (22 bp)
STS 10893 (12 bp)
SIUC-Satt9.0 (17 bp)
SIUC-Satt9.5 (21 bp)
Variant diphenol oxidase
Na/H antiporter-like 1
SIUC-Satt22/Indel22 (20 bp)
BARC-Satt309 (17 bp)
Na/H antiporter-like 2
DNA helicase-like 1
DNA helicase-like 2
Minisat 2/SIUC-Sctt39 (45 bp)
SIUC-Sat_40 (38 bp)
There was evidence for a large and highly polymorphic region within the BAC (743 SNPs in 59 Kbp; from 1,500–60,500 bp). A highly polymorphic region was expected to be a characteristic of the region introgressed into Forrest from Peking. Equally, relatively monomorphic regions were common when comparing sequences of US cultivars [54, 55]. In fact, the three alleles (Forrest, Williams 82 and A3244) were nearly identical outside the 59 kbp central region. For example the region that encompassed gene 10 had just 5 SNPs across 18.25 kbp among the 3 alleles and no alloproteins.
The RLK contained 2 of the 6 SNPs in the BAC predicted to cause amino acid changes (Ruben et al. 2006; Afzal et al. 2012) when comparing Forrest allele to Williams 82 and Asgrow 3244. However, it contained 7 synonymous SNPs. In comparison the laccase contained 3 SNPs predicted to cause amino acid changes and only one that was synonymous (Figure 1). The antiporter alloproteins contained one non-synonymous and two synonymous SNPs.
The greater number of differences between alleles of the GmRLK18-1 and the other candidate genes was taken as evidence this should be the first gene tested in transgenic plants. The RLK was subcloned on a 9,772 bp insert for plant transformation (Figure 1D) that contained 6 Kbp of 5′ sequence, 2.7 kbp of geneic sequence and 1 kbp of 3′sequence. The 5′ sequence extended to the TATA box of the neighboring gene (#4). The region contained 7 potential cis-regulatory elements (CREs) of 8 bp that were identical to motifs found in the 5′ regions hundreds of plant genes. The four closest to the RLK are shown in Figure 1D. Only the most distant was polymorphic between cultivars. The Forrest allele of this CRE was the identical to sequences found 5′ to 1,091 rice genes and 108 Arabidopsis genes. There was a 53 bp deletion from Forrest (at 33,942 bp) compared to Essex, Williams 82 and Asgrow 3244 in the enhancer region of the gene that was named SIUC-indel-(plus) + 3.5. There was also a 22 bp deletion 3′ to the gene at (at 40,215 bp) named SIUC-indel (minus)-3.5. This deletion was within a complex microsatellite repeat mainly composed of AAAG motifs. Neither deletion appeared to encompass any previously characterized CREs. In the intergeneic region between RLK and laccase were four polymorphic small indels or microsatellite-like polymorphisms.
Syntenic paralogs of the Rhg1/Rfs2 locus
The soybean genome was hypothesized to be the product of a diploidized tetraploid [54, 55]. Therefore, a detailed molecular analysis of the Rhg1/Rfs2 locus required that paralogs and syntenic gene clusters be identified. The most identical paralog (the homeolog) was found on BAC H38f23 (Additional file 1: Figure S1) that was from Lg B1 (chromosome 11) in a region where loci with functions similar to Rhg1/Rfs2 were mapped [22, 36, 37]. Significantly, the sequence of H38f23 contained a complete set of syntenic genes for a second Rhg1/Rfs2 locus and surrounding genes (an RLK, laccase, both antiporters, the kinase and the helicase; Glyma11g35710- Glyma11g35660; Figure 1; HQ008940). Interrogation of Soybase showed the transcript abundance patterns for the syntenic RLK pair and laccase pair were both abundant in roots as reported in the RNA sequence atlas generated from a Williams 82 sister line derived NIL (cultivar ‘P-C609-45-2-2’ a BC5F5 plant derived from G. soja (‘PI 468916’) backcrossed into G. max (‘A81-356022’) . The laccase transcripts were 15 fold more than either of the RLK transcripts in roots. Though it must be noted this genotype was either the rhg1-d or Rhg1-e allele.
Comparing the syntenic BACs the DNA sequence identity was high (~97% in geneic regions; Additional file 1: Figure S1E; Additional file 2: Table S1). At the RLKs, GmRLK18-1 and GmRLK11-1, with intron boundaries as described by  and , the amino acid identity was 93% overall with 94% in the LRR, 93% in the trans-membrane domain and 97% in the kinase domain (Additional file 1: Table S1). The N terminal signal peptide was most divergent (11 changes in 50 amino acids). The laccase and the antiporter also showed 85–96% amino acid identity with their syntenic paralogs.
GmRLK11-1’s amino acid sequence was less identical to GmRLK18-1 than the alleles (98–99%) of the gene. However, four of the six residues that differ among alloproteins of GmRLK18-1 were identical to the Forrest allele in the homeologous protein GmRLK11-1-a. At the other 2 residues the changes were identical to the susceptible alleles of Essex and Williams 82. From sequence comparison and transcript abundances it appears likely the paralogs GmRLK18-1 and GmRLK11-1 would share functions and should hetero-dimerize .
Allele and paralog discrimination in NILs and transgenic plants
Dominant, recessive or co-dominant nature of Rhg1/Rfs2 in NILs
In the NIL line 34–33 that was heterogenous at the Rhg1/Rfs2 locus, about 12% of plants (4/34) were still heterozygous by the F5:7:13 (thirteen generations after the original cross and single seed selection at F5; NIL selection at F7; Figure 2). Those plants were heterozygous at the TMD1 marker because they had not yet been fixed by recombination. The frequency of heterozygous plants among F5:7:13 generation seed was surprisingly high (4/34 or 12%) compared to the expected frequencies at the F5 (6.25%), F7 (1.56%) and F13 (0.025%) generations. Even calculating the expected heterozygous plant frequency starting from the F5 the occurrence was higher than expected (1/256 or 0.39%). The existence of these plants suggested that fixation by re-assortment was selected against in viable zygotes (or gametes) or that the heterozygous state was under positive selection. The remaining 30 plants (88%) were fixed equally to Forrest or Essex alleles (15/34 or 44% each).
The cyst scores for all plants in the NIL population corresponded with the respective alleles at the Rhg1/Rfs2 locus so plants with resistance alleles had female indices (FIs) less than 8% and plants with susceptibility alleles had FIs greater than 84% (Figure 2). For the four heterozygous plants, polymorphic at TMD1, the cyst scores ranged from 5–12% which corresponded to those for resistant or moderately resistant plants. Therefore, the Rhg1/Rfs2 locus was dominant in this set of NILs infested with HgType 0 population JB3. However, heterozygous NILs were reported co-dominant with HgType 0 population PA3 in earlier tests where FIs ranged from 5–40% . Co-dominant and recessive roles of plant disease resistance loci were previously associated with factors needed by the pest for successful parasitism [56, 57]. The discrepancies in dominance among different populations may be associated with the genetic background in which the gene resides or may result from interactions among genes at the Rhg1/Rfs2 locus and/or modifier genes at other loci .
Inhibition of root growth by alleles of Rhg1/Rfs2 in the NILs
When counting the cysts with prior knowledge of the allele at Rhg1/Rfs2 in the standard assay of SCN it was noted that root mass and vigor appeared to differ among genotypes. Measurement of root masses showed a significant difference (P < 0.0015) among NILs that were associated with the allele at Rhg1/Rfs2 (Figure 2). Across several experiments, both NILs that were susceptible and NILs that segregated most susceptible lines had higher root masses (mean 2.24 ± 0.19 g) than their SCN resistant (1.22 ± 0.22 g) or heterozygous (1.79 ± 0.18 g) counterparts. This phenomenon might underlie the global association of resistance to SCN with low seed germination, seedling vigor and ultimately seed yield noted previously [40–42]. In addition, the lower root growth might provide an avoidance mechanism that is part of a broad resistance effective against all Hg Types.
The distribution of recombination events found previously among the six Hg Type 0 susceptible PIs  suggested that the action of Rhg1/Rfs2 required elements to the 3′ side of the RLK intron or in the C terminal portion of the protein. One of the 2 amino acid substitutions found in the intracellular kinase of the complete RLK (Gly539Ala and Ser770Pro) may be key to transducing intracellular signaling leading to the resistance response . Some mutations in the kinases of other plant RLKs are known to be lethal [58–60]. Therefore, it may be the kinase at the Rhg1/Rfs2 locus that underlies restricted root growth in resistant genotypes. The effect may be direct or occur after protein-protein interaction(s). The signal transduced is likely to result in a negative reponse since the kinase domain lacks the RD amino acid motif needed for ATP binding  and so protein phosphorylation.
Transcripts of the GmRLK18-1 alleles at Rhg1/Rfs2 were found in soybean roots
Analyses of disease responses in plants transgenic with the resistance allele of the GmRLK18-1 in the greenhouse
In order to test the hypothesis that GmRLK18-1-a underlay part of the activity of the Rhg1/Rfs2-a locus susceptible plants were made transgenic with this allele. Such plants had a resistance allele at a new location and a susceptibility allele (GmRLK18-1-a; rhg1/rfs2-e) on chromosome 18. Several primary transgenic lines were created by biolistics and fertile lines selected in two different cultivars (‘X5’ and ‘Westag 97’) for analyses of SDS and SCN responses.
Association of mean root growth in transgenic lines with pleiotropic resistance to two pests in two different greenhouse assays and insect herbivory in field tests
Root mass (g)
SCN FI (%)
0 ± 0.0
0 ± 0.0
100 ± 13
0 ± 0.0
60 ± 11
11 ± 3
38 ± 6
130 ± 11
12 ± 3
Root mass (g)
1.0 ± 0.1
3.2 ± 0.7
4.5 ± 0.9
4.3 ± 0.8
1.0 ± 0.1
1.0 ± 0.3
1.5 ± 0.5
1.5 ± 0.4
6.5 ± 1.9
6.3 ± 1.3
0.5 ± 0.3
1.8 ± 0.5
Insect PI (%)
Mean shoot dry weight (g)
Leaf defoliation (%)
22 + 13
33 ± 4.0
40 ± 9
11.6 ± 0.3
59 + 15
90 ± 5.0
36 ± 8
11.2 ± 0.3
The results from X5 transgenics were comaperd to results from Westag97 transgenics in the SDS assay with infested plants (Table 2). Westag97 plant were bigger and more stress resistant than X5. However, the Westag plants transgenic with the RLK again proved more significantly more resistant to SDS caused by F. virguliforme. Therefore, the resistance was not transgenic event or cultivar specific.
In separate assays with plants from the same lines the resistance to SCN was partial in the GmRLK18-1-a transgenic plants (Table 2A) judged by female index (FI). Partial resistance was expected as rhg4-e was present and Rhg4-a absent. SCN FI was reduced by 30–50% across four experiments using three Hg Types (P < 0.01). Because the X5 selfed transgenic plants reported here had a susceptible allele at Rhg4, partial resistance was the expected outcome. The partial resistance was confirmed with another isolate of SCN at Harbin University (China). A cross was made to the RIL EF2 via a susceptible NIL WE1. EF2 was rhg1-e/rfs2, rhg1-e/rfs2, Rhg4, Rhg4. The progeny segregated for the RLK and so were rhg1-e/rfs2, rhg1-e/rfs2, Rhg4, rhg4:: RLK18-1 or rhg1-e/rfs2, rhg1-e/rfs2, Rhg4, rhg4 with no transgene. Plants with the RLK had a 20% lower FI and were almost in the resistant class (FI < 10; Table 2). Again the Westag 97 plants were tested with the same assay. Westag 97 appeared to encode a functional Rhg4 resistance allele. Consequently the RLK transgenic Westag 97 showed nearly complet resistance. Therefore, the GmRLK18-1-a was alone sufficient to provide for an Rhg1–a like activity at different locations in different cultivars. The linked genes may reduce SCN numbers but did require the action of a second locus (Rhg4) to provide full resistance to Hg Type 0 as expected from [23, 61].
Analyses of plants transgenic with the resistance allele of the GmRLK18-1 in the field
Field grown plants showed that the RLK was associated with increased insect herbivory (Table 2C; Figure 5). The primary damaging pest was the Japanese beetle (Popillia japonica, Newman). Insects were attracted to the RLK transgenic plants judged by pest incidence (PI) measured from R1-R7 growth stages. Herbivory resulted in plants with less leaf area that produced less shoot biomass at harvest compared to isogenic plants lacking the GmRLK18-a allele. There was no incidence of SDS or SCN in the field during the 2010 or 2011 seasons.
Effect of the Forrest allele of GmRLK18-1 on transgenic plant development
In both SCN greenhouse assay and field the Forrest allele of GmRLK18-1 caused a reduction in root and shoot mass (Table 2). The reduction was significant even though the X5 cultivar was innately smaller than the NILs due to a very much earlier maturity date. Reductions occurred in seedlings in the SCN assay (Table 2A), the SDS assay (Table 2B) and whole plants in the field (Table 2C) at harvest maturity. The reduction in growth in the field lead to a reduction in biomass but not seed yield and seed number per plant. Mid- to late-season defoliation by insect pests rarely reduces soybean yield. In both SCN and field assays plant growth was limited by available water. A reduction in root growth might be related to reduced water uptake.
However, in the water saturated assays of SDS resistance the Forrest allele of GmRLK18-1 caused an increase in root and shoot mass (Table 2) in the presence of F. virguliforme infestations. The increase in growth lead to a increase in biomass and seed yield (Figure 5) and seed number per plant in plants grown to maturity in the greenhouse (1.7 ± 0.3 compared to 5.7 ± 1.1). Water sufficiency appeared to make the Forrest allele of GmRLK18-1 beneficial to growth and yield when F. virguliforme was present. In contrast, plants that were not infested in the SDS assays showed the Forrest allele of GmRLK18-1 caused a decrease in root and shoot mass.
In field trials of plants transgenic with GmRLK18-1 further evidence of pleiotropy was discovered when leaf herbivory was shown to be significantly worsened. Insect pests are separated into the chewing and sucking guilds. Resistance mechanisms to the guilds were known to operate by different pathways but some involved the Rhg4 locus or region . Arabidopsis resistance to F. graminearum was recently shown to require operation of both SA and JA pathways . The same Nils used here when infested with H. glycines decreased proteins involved in salicylic acid responses and increased those in the jasmonic acid signaling pathways . Therefore, plants transgenic with Gm18RLK18-1 might be altered in their responses to a broad range of pests and pathogens. Additional pathogens will be tested in future studies.
The Rhg1/Rfs2 locus was shown to center on three genes where recombination was suppressed within a wider region that was high in recombination. Perhaps the failure to recombine within the central 3 gene cluster caused the unusually high frequency of heterozygous plants in the neighboring region. A potential suppressor locus acting on the resistance allele of Rhg1 was identified earlier [21, 45, 47]. The allele of the gene on chromosome 7 (Lg M) had to be inherited from the resistant parent to prevent zygote or embryo lethality in SCN resistant plants. Consequently, all resistant plants co-inherited Rhg1/Rfs2 on chromosome 18 (Lg G) and the resistance allele on chromosome 7 (Lg M). That gene may have had a functional homeolog in this study. However, the locus was fixed to the R haplotype in all the NILs, so fine mapping was not possible. Plant transgenic with the RLK were fertile and produced seed and they would lack the resistance allele of the gene on chromosome 7 (Lg M). Therefore, GmRLK18-1 does not cause zygote or embryo lethality. Additional genetic elements must be involved in that phenomenon.
An absence of recombination events in a region can be caused by; deletions; insertions; inversions; the condensed heterochromatin found near centromeres; and recombinant allele lethality. The first four phenomena were not occurring at the Rhg1/Rfs2 locus since DNA sequences in resistant and susceptible cultivars were co-linear over at least 87 kbp. Further, the Rhg1/Rfs2 locus mapped into a region within 1 Mbp of the predicted telomere region of chromosome 18. Therefore, recombinant allele lethality may be the cause and a hypothetical model for its action was developed. The kinase domain of the RLK at Rhg1/Rfs2 was proposed as the element causing lethality to gametes or the zygote carrying a recombination event. The hypothetical anti-porter protein was proposed as the target locus kept in the resistant state when the RLK is conferring SCN resistance (Figure 6). The laccase that was located between these two genes, along with the intergenic regions, are held in phase by the locus. It may be possible that the lethal nature of the three genes linked in association at the resistance locus are not fully suppressed by the modifier locus (they are leaky) and results in the inhibition of root growth observed in SCN resistant cultivars (Figure 2).
The linkage disequilibrium studies reported with probes in, and around, the RLK inferred a larger more complex structure for Rhg1 than simply the RLK polymorphisms [15, 16, 25, 64]. Possible roles included contributions to additive resistance; contributions to resistance in other resistance types (eg PI88788 and Toyosuzu; R types 2 and 3) or contributions to the resistance to other Hg types. The hypothesis that the linked genes were factors necessary in susceptible genotypes for SCN parasitism was unlikely [16, 25, 53]. The dominance of Rhg1/Rfs2 to SCN in NIL segregation also suggested the genes were active in resistant types and inactive in susceptible genotypes.
The Rhg1/Rfs2 locus was shown to include the GmRLK18-1. Previously  also reported inhibition of the resistance allele increased FI in composite roots, but not to a significant degree in small experiments with few replications. Here Rhg1/Rfs2 was inferred to be a complex of three genes assembled and co-inherited over long periods of selection for resistance to a endemic pest, root parasitic nematodes. However, the GmRLK18-1 alone at a new genomic locations in transgenic plants was sufficient for resistance to a relatively new pathogen of soybean F. virguliforme. Revisions to the sequence of the RLK made by re-sequencing PCR products  were confirmed by BAC fragment sequencing here. Fortunately the amino acid changes made to  lay outside of regions under detailed analysis for ligand and antibody binding [65, 66]. Previously CLE peptide in nematode secretions were shown to be perceived by RLKs in the CLAVATA1, CLAVATA2 and CORYNE families . In vitro, the purified GmRLK18-1 LRR domain has been shown to bind strongly to the CLE peptide found in nematode secretions as well as two proteins and two CLE peptide of plant origin involved in tracheary element inhibition [43, 66]. It will be of interest to discover whether the ligand profiles are similar for the two resistances. The biochemical analysis of the RLK protein will lead to refinements of models for modes of action for the pleiotropic resistance.
The allelic discrimination probes developed provide a high throughput alternative to satellite markers for marker assisted selection. These tools will facilitate molecular breeding for resistance to two or more important pests and diseases of soybean . The discovery of a syntenic set of paralogs to the RLK, laccase and transporter gene cluster at Rhg1/Rfs2 may also assist marker assisted breeding. The discovery of the syntenic homeologs raised significant barriers to reverse genetic approaches to analyses of the RLK at Rhg1/Rfs2. Another barrier to reverse genetic approaches was the role of Rhg1/Rfs2 in normal plant development that can be inferred from the restricted root growth of NILs (Figure 2) and transgenic plants (Table 2). A third barrier to reverse genetics was the uncoupling in eukaryotes of transcript and protein abundances. Both effects may have contributed to the recent report of no significant difference between the SCN counts on composite plant roots with the GmRLK18-1-b alone targeted for inhibition . It must be noted that the GmRLK18-1-b protein abundance was not measured in that study in spite of the tools available to do so [43, 65, 66]. Further, note that the numbers of cyst did differ between treatments  and that might be significant in larger, replicated studies or under different conditions. Further note that allele was recessive. However, further proofs of the nature of Rhg1/Rfs2 locus function may require knock-outs of each of the paralogs, or stable transformation to several new locations, followed by measurements of genetic segregation. In each case, the analysis will be complicated by the co-dominant nature of the resistance gene in certain experiments. In fact, the possibility that the susceptible Essex allele of Rhg1/Rfs2 was functional by promoting root growth and susceptibility, and so actively promoted the establishment of parasitism by SCN and F. virguliforme, should be explored with transgenic plants. The model inferring that soybean has adapted part of an existing pathogen recognition and defense cascade (SCN parasitism and insect herbivory), to a new pathogen (F. virguliforme) disease (root rot and leaf scorch and syndrome (SDS) has broad implications for crop improvement. Stable resistance to many pathogens might be achieved by manipulation of the genes encoding a small number of pathogen recognition proteins in the RLK family.
Many of the NILs were described previously [7, 11, 16]. All lines were released and consequently are available on request as seed. Seeds of NIL 34–23 (resistant haplotype between markers Satt 214 to Satt 570) and NIL 34–3 (susceptible haplotype from the marker Satt 214 to the Sat122-Satt 570 interval) were obtained at the F5:13 generation. Genotypes were rhg1-e/rfs2, rhg1-e/rfs2, Rhg4, Rhg4 for NIL 34–3 and Rhg1-a/Rfs2, Rhg1-a/Rfs2, Rhg4, Rhg4 for NIL 34–23 whereas NIL34-33 contained both those and rhg1-e/rfs2, Rhg1- a/Rfs2, Rhg4, Rhg4 in different plants (NIL34-33-1 to −34).
Soybean cv. X5 and Westag 97 were used for transformation because they could be regenerated to plant efficiently from embryo cultures . They were susceptible to both SCN and SDS. X5 was judged to be rhg1-e/rfs2, rhg1-e/rfs2, rhg4, rhg4 and Westag 97 was judged rhg1-e/rfs2rhg1-e/rfs2, Rhg4, Rhg4 and based on the assays reported here. Crosses were made to RIL EF2 for SCN tests which was rhg1-e/rfs2, rhg1-e/rfs2, Rhg4, Rhg4.
Soybean plants were grown in cones containing a 1:1 ratio of sand soil mix, placed in a water bath to maintain the root zone at 26°C or in a growth chamber at 26°C. Infection with the three different Hg Type 0 SCN populations (PA3, JB3 and YC3) in separate assays consisted of inoculating 2,000 eggs to each 14 day old seedling. Plants were watered daily with 100 ml per pot. Infested soybean plants were removed from the cones at 28 days post infestation (dpi) and cyst numbers counted and compared to ‘Essex’ susceptible controls . Some plant roots were harvested at 10 and 20 days for protein analyses. Root masses of the NILs and transgenic plants were noted. Heterozygous NILs were replanted in non-infested soil and grown to maturity. Seed were harvested. The indicator lines female indices (FI) for nematode population JB3 were ‘PI54840’ (FI 7%), PI 88788 (FI 2%), PI90763 (FI 1%), PI437654 (FI 0%), ‘PI 209332’ (FI 1%), ‘PI89772’ (FI 2%) ‘PI548316’ (FI 8%) and ‘PI548402’ (FI 3%). The soil collected from Yichun in China (YC3) contained Hg Type 0 (SCN race 3) was ‘Peking’ (FI 0%), ‘PI 88788’ (FI 0%), ‘PI 90763’ (FI 6%) and ‘Pickett’ (FI 9%). Essex was the susceptible genotype used across assays to determine 100% FI. Therefore, the standard differentials showed these HG Types to all be variants on Hg Type 0 (Niblack et al. 2003) corresponding to race 3 (Riggs and Schmitt 1988). Nematode population PA3 was described previously [16, 23].
Assays of root development in the absence of infestations
Root development was assayed in NILs and transgenic plants using conditions identical to the SCN assays or SDS assays (as noted in text) except that the roots were not infested with either pathogen. Roots were weighed separately from shoots during destructive sampling. These plant organs were used for RNA and protein extractions.
SDS measures following F. virguliforme infestations
Greenhouse assays of the effects of F. virguliforme infestations followed the methods previously described . A culture of F. virguliforme virulent strain ‘Mont-1’ (NRRL 22292; MAFF 238545) was provided by Dr. A. Fakhoury (SIUC). Briefly, the strain was grown on ten potato dexrose agar plates. Hyphae and spores were washed with distilled water and 10 μl of this is used for spore count on a hemocytometer under a microscope. Spore counts of 104 spores/cm3 of sand and soil mix were used. Seed were sown in sterilized 1:1 (v/v) of sand and soil inoculated with F. virguliforme virulent strain Mont-1 in 10 cm square pots. Pots were kept in trays filled to the 5 cm level with water to keep the lower half of the soil saturated. The water was inoculated with 104 spores/cm3 to avoid dilution of the inocula in the pots. Transgenic plants, non transgenic X5, selected ExF RILs and a set of breeders advanced lines were grown in the green house at the Southern Illinois University Horticulture Research Center in Carbondale, IL. Experiments were conducted from October 2009 to May 2012. Plants were grown with a 14 h photoperiod under supplemental lights. The air temperature ranged from 20 ± 2°C at night to 27 ± 2°C during the day in the green house. Leaf symptoms were rated every 7 days from infestation to senescence. The standard sudden death syndrome DS score comparable across studies was rated at 21 days after inoculation, determined on the basis of the degree of leaf damage (chlorosis/necrosis) on each plant, and was rated on a scale of 1 to 9 (1 = 0–10%/1–5%, 2 = 10–20%/6–10%, 3 =20–40%/10–20%, 4 =40–60%/20–40%, 5 = 60%/_40% of leaf surface chlorosis/necrosis, respectively, 6 = up to 33% premature defoliation, 7 = up to 66% premature defoliation 8 = 66% premature defoliation, and 9 = premature death of plant). At 28 days after infection (dai) roots were washed, photographed and a root sample (1 g) taken. Root rot severity (RS) was scored on a scale of 1–5 where 1 was unaffected; 2 was discolored; 3 was discolored and partly rotted; 4 was discolored and heavily rotted; and 5 was discolored and necrotic. Plants were weighed, the root to shoot ratio estimated visually and then repotted into the media of the reciprocal genotype to test for pot effects. The experiments were repeated on 3 occasions using 3–5 plants of each genotype. Highly resistant plants of RIL EF23, and highly susceptible plants of RIL EF85 were assayed in parallel with each test.
DNA and RNA for genotype analysis
DNA was isolated following . Concentrations of DNAs were calculated by measuring absorbances at 260 and 280 nm. Total RNA was isolated with TrizolTM (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s instructions. First strand cDNA was synthesis carried out using oligo dT primers using a cDNA synthesis kit, according to manufacturer (Invitrogen). Presence of the Rhg1/Rfs2 resistance alleles was confirmed by PCR analysis using TMD1 an indel marker in the RLK intron. Several designs of TMD1 primers have been reported [11, 16, 23]. Used here were the primers pair; forward 5′- CAC CTG CAT CAA GAT GAA CA -3′ and reverse 5′- GCC TAT TAC TTG GGA CCC AA -3′ (Additional file 4: Table S2). Genotyping by markers linked to Rhg1/Rfs2 used about 50 ng of DNA for microsatellite analysis on PAGE after  and on agarose gels after .
BAC sequence and allele comparisons
Southern hybridizations were performed following the standard procedure described in  to identify paralagous BACs. BACs B73p06 and H38f23 were sequenced at TIGR (nee JCVI). Briefly, the entire BAC was sheared by nebulization to provide fragments in the 3–5 kbp or 9–11 kbp range. The fragments were ligated into pHOS2 and used for Sanger DNA sequencing. BACs were sequenced to 8–12 fold redundancy and assembled. Assembly quality was judged by BLAST comparisions for sequences from A3244 BACs [11, 16, 50] and Williams 82 genome sequence .
Allelic discrimination at the RLK within the Rhg1/Rfs2 locus
The SNP genotyping assay within the gene encoding the RLK was performed using a custom TaqmanTM Kit. Three probes were designed for the SNPs at 1,486 bp, 506 bp and 2,040 bp (relative to the translation start site) to distinguish the 8 commonest alleles of the RLK (Additional file 4: Table S2). Six process were designed to the non-synonymous SNPs in the GmRLK18-1 preotein to distinguish the seven allotypes. The PCR reactions were carried out using a 3 step PCR protocol with one hold at 95°C for 10 minutes followed by 35 cycles that included a denaturation cycle of 95°C for 30 sec, annealing at 58°C for 10 seconds and an extension at 68°C for 20 sec.
Primers for SNPs within the Rhg1/Rfs2 locus on Lg G were used in fine melt curve assays as described previously  with the following modifications. Briefly, genomic DNA was used; multiple amplicon sizes were detected on PAGE gels; an ABI7900 with HTM software was used; melt curve data were normalized by both local and global metrics. Primers used were to the SNP at 2090–1 for transgenic cDNA between X5 and X5 transgenics. The target was SNP G to T (G for resistant and T for susceptible plants). SNP2090-1 forward primer was 5′- GTT GGT TGA TCC AGA AGG GTT -3′ and the reverse was 5′- CTA AGC TTC CTG AGG CCT TG -3′.
For detection of mRNA products of GmRLK18-1 (Figure 3) the qRT-PCR methods described previously [45–47] were followed. Briefly, mRNA concentrations were estimated from Cot curve analyses using cDNAs.Actin1 was the control used for comparisons.
Allelic discrimination across the region flanking the Rfs2/Rhg1 locus
Screens for new recombination events used the flanking markers Satt309 and Satt038. Recombinants identified with those markers were screened with TMD1 and three SNP primers for TaqmanTM assays used were set as described in [45–47] (Afzal, 2008b; 2012) with the following modifications. The SNPs used were; AX196295 10893 between the laccase and RLK; AX196295 37583 CR-G in the laccase; and AX196295 37581 CR-G in the hypothetical antiporter gene.
Total root protein extraction, SDS–PAGE and Western hybridization
Protein from root material was isolated from infested and non-infested roots was extracted after . Total protein concentration was determined using a non-interfering protein assay. For the Western hybridizations, a custom made antibody generated against a peptide [C]TL SRL KTL DIS NNA LNG NLP ATL SNL S from the LRR domain of RLK at RHG1/RFS2 was used (Alpha diagnostics, San Antonio, Texas).
For soybean transformation, the cassette included pSBHB94 that was a 9.772 kbp insert sub-cloned from BAC B21d09 by nebulization, size fractionated to 9–11 kbp and ligated into pHOS2. The plasmid pSBHB94 encompassed the sequences found from 30,423–40,194 bp in BAC B73p06. Transformation, selection and plant regeneration were conducted after . Briefly, proliferative embryogenic cultures of soybean cv. X5 (AAFC breeding line X2650-7-2-3) or Westag 97 were co-bombarded with the pHOS_SBHB94 and Hyg R  constructs; transgenic events were selected and maintained on 55 mg L−1 hygromycin; embryos were matured on antibiotic-free medium, air desiccated and converted on B5 medium ; tissue cultures and regenerating plantlets were maintained at 20 C and 20 h photoperiod. The plantlets were transferred to soil and plants were regenerated under controlled conditions as in . Primary transgenic (T0) plants were tested for the presence of the pHOS::SBHB94 transgene using PCR with the TMD1 primers. Fourteen primary lines were obtained and grown for seed. T2 seed from T1 plants was tested for transgene segregation to identify homozygous T1 individuals. Event 6B3-7D2 and 6B3-7D3 provided seed for the X5 experiments described here. Repeats of the experiments were made on seed of events 8B2-7D2 and 8B2-15D1 in Westag 97.
The genotype of X5 was rhg1-e/rfs2, rhg1-e/rfs2, rhg4rhg4. Purified stable transgenics were of genotype rhg1-e/rfs2rhg1-e/rfs2, rhg4rhg4:: GmRLK18-1, GmRLK18-1 (so Rhg1- a/Rfs2Rhg1- a/Rfs2) as shown by markers TMD1 and A2D8. Expression of the transgene was established by qRT-PCR from cDNA with allele specific Taqman probes and HRM of amplicons. Protein allotypes were identified by two dimensional PAGE  followed by Western hybridization [43, 44].
Field trials were conducted at the ARC in Carbondale during 2010 and 2011 using conditions described in . Plants were arrayed in 12′ plots arrayed in a randomized complete block with 4 replications. Insect herbivory was measured as described in . Briefly, the pest incidence was calculated as the number of individual plants within a given line that were affected by herbivorous insects; the pest severity was the percent defoliation; and both were measured once a week from the R1 to R6 growth stages. The major defoliating pest was the Japanese Beetle (Popillia japonica, Newman). Plant biomass was measured for 5 plant per plot. The seed yield of each plot was measured after harvest.
TMD1 marker, FJ520231; Corrected RLK at Rhg1/Rfs2 gene AF506516 and mRNA AF506517; SIUC-Satt122, bankit1155667; BAC pB73P06 complete sequence JN597009; pSBHB94 complete sequence HQ008939; BAC pH38f23 complete sequence TBD.
Three supplementary figures and three supplementary tables are published online and links to 7 sequences at Genbank are provided.
Receptor like kinase
soybean cyst nematode
sudden death syndrome
bacterial artificial chromosome
Genebank System National Institute of Agrobiological Sciences, Tsukuba, Japan
The Agriculture Research Service Culture Collection National Center for Agricultural Utilization Research, USDA/ARS, Peoria, IL USA.
The physical map was supported by the NSF under Grant No. 872635 and proteins analysis under Grant No. 0487654. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the NSF. We dedicate this manuscript to the memory of Dr. Oval Myers Jr.
- Roy KW, Rupe JCD, Hershman DE, Abney TS: Sudden death syndrome of soybean. Plant Dis. 1997, 81: 1100-1111. 10.1094/PDIS.1922.214.171.1240.View ArticleGoogle Scholar
- Wrather JA, Anderson TR, Arsyad DM, Tan Y, Ploper LD, Porta-Puglia A, Ram HH, Yorinori JT: Soybean disease loss estimates for the top ten soybean-producing countries in 1998. Can J Plant Path. 2001, 23: 115-121. 10.1080/07060660109506918.View ArticleGoogle Scholar
- Li S, Hartman GL GL, Chen Y: Evaluation of aggressiveness of Fusarium virguliforme isolates that cause soybean sudden death syndrome. Journal of Plant Pathology. 2009, 91: 77-86.Google Scholar
- Kolander TM: The host range of Fusarium virguliforme on rotational crops and common plant species and its survival and growth on crop residue. 2010, University of Minnesota, MS Thesis, 196-http://purl.umn.edu/60700,Google Scholar
- Njiti VN, Johnson JE, Torto GA, Gray L, Luo Y, Gibson PT, Lightfoot DA: Inoculum rate influences selection for field resistance to sudden death syndrome in the greenhouse. Crop Sci. 2001, 41: 1726-1733. 10.2135/cropsci2001.1726.View ArticleGoogle Scholar
- Covert SF, Aoki T, O’Donnell K, Starkey D, Holliday A, Geiser DM, Cheung F, Town CD, Strom A, Juba J, Scandiani M, Yang XB: Sexual reproduction in the soybean sudden death syndrome pathogen Fusarium tucumaniae. Fungal Genetics and Biology. 2007, 44: 799-807. 10.1016/j.fgb.2006.12.009.View ArticlePubMedGoogle Scholar
- Njiti VN, Doubler TW, Suttner RJ, Gray LE, Gibson PT, Lightfoot DA: Resistance to soybean sudden death syndrome and root colonization by Fusarium solani f. sp. glycine in near-isogenic lines. Crop Sci. 1998, 38: 472-477. 10.2135/cropsci1998.0011183X003800020033x.View ArticleGoogle Scholar
- Kazi S, Shultz JL, Bashir R, Afzal AJ, Njiti VN, Lightfoot DA: Separate loci underlie resistance to soybean sudden death syndrome in ‘Hartwig’ by ‘Flyer’. Theor Appl Genet. 2008, 116: 967-977. 10.1007/s00122-008-0728-0.View ArticlePubMedGoogle Scholar
- Brar HK, Swaminathan S, Bhattacharyya MK: The Fusarium virguliforme toxin FvTox1 causes foliar sudden death syndrome-like symptoms in soybean. Mol Plant Microbe Interact. 2011, 24 (10): 1179-118. [Epub ahead of print]View ArticlePubMedGoogle Scholar
- Iqbal MJ, Yaegashi S, Ahsan R, Shopinski KL, Lightfoot DA: Root response to Fusarium solani f. sp. glycines: Temporal accumulation of transcripts in partially resistant and susceptible soybean. Theor Appl Genet. 2005, 110: 1429-1438. 10.1007/s00122-005-1969-9.View ArticlePubMedGoogle Scholar
- Triwitayakorn K, Njiti VN, Iqbal MJ, Yaegashi S, Town CD, Lightfoot DA: Genomic analysis of a region encompassing QRfs1 and QRfs2: genes that underlie soybean resistance to sudden death syndrome. Genome. 2005, 48: 125-138. 10.1139/g04-103.View ArticlePubMedGoogle Scholar
- Radwan O, Liu Y, Clough SJ: Transcriptional analysis of soybean root response to Fusarium virguliforme, the causal agent of sudden death syndrome. Mol Plant Microbe Interact. 2011, 24: 958-972. 10.1094/MPMI-11-10-0271.View ArticlePubMedGoogle Scholar
- Arelli PR: Inheritance of resistance to Heterodera glycines race 3 in soybean accessions. Plant Dis. 1994, 78: 898-900. 10.1094/PD-78-0898.View ArticleGoogle Scholar
- Gelin JR, Arelli PR, Rojas-Cifuentes PR: Using independent culling to screen plant introductions for combined resistance to soybean cyst nematode and sudden death syndrome. Crop Sci. 2006, 46: 2081-2083. 10.2135/cropsci2005.12.0505.View ArticleGoogle Scholar
- Prabhu RR, Njiti V, Johnson JE, Schmidt ME, Klein RJ, Lightfoot DA: Selecting soybean cultivars for dual resistance to cyst nematode sudden death syndrome with two DNA markers. Crop Sci. 1999, 39: 982-987. 10.2135/cropsci1999.0011183X003900040005x.View ArticleGoogle Scholar
- Ruben E, Jamai A, Afzal AJ, Njiti VN, Triwitayakorn K, Iqbal MJ, Yaegashi S, Bashir R, Kazi S, Arelli PR, Town CD, Ishihara H, Meksem K, Lightfoot DA: Genomic analysis of the Rhg1 locus: candidate genes that underlie soybean resistance to the cyst nematode. Mol Genet Genom. 2006, 276: 503-516. 10.1007/s00438-006-0150-8.View ArticleGoogle Scholar
- Kazi S, Shultz JL, Afzal J, Hashmi R, Jasim M, Bond J, Arelli PR, Lightfoot DA: Iso-lines and inbred-lines confirmed loci that underlie resistance from cultivar 'Hartwig' to three soybean cyst nematode populations. Theor Appl Genet. 2009, 120: 633-640.View ArticlePubMedGoogle Scholar
- Hussey RS, Grundler FMW: Nematode parasitism of plants. The physiology and biochemistry of free-living and plant parasitic nematodes. Edited by: Perry RN, Wright DJ. 1998, CAB International Wallingford, USA, 213-243.Google Scholar
- Caldwell BE, Brim CA, Ross JP: Inheritance of resistance of soybeans to the cyst nematode Heterodera glycines. Agron J. 1960, 52: 635-636. 10.2134/agronj1960.00021962005200110007x.View ArticleGoogle Scholar
- Niblack TL, Noel GR, Lambert KL: The Illinois SCN type test: practical application of the Hg Type classification system. J Nematol. 2003, 35: 355-345.Google Scholar
- Webb DM, Baltazar BM, Arelli PR, Schupp J, Clayton K, Keim P, Beavis WD: Genetic-mapping of soybean cyst-nematode race-3 resistance loci in the soybean PI 437654. Theor Appl Genet. 1995, 91: 574-581.View ArticlePubMedGoogle Scholar
- Vierling RA, Faghihi J, Ferris VR, Ferris JM: Association of RFLP markers with loci conferring broad-based resistance to the soybean cyst nematode:Heterodera glycines. Theor Appl Genet. 1996, 92: 83-86. 10.1007/BF00222955.View ArticlePubMedGoogle Scholar
- Meksem K, Pantazopoulos P, Njiti VN, Hyten DL, Arelli PR, Lightfoot DA: 'Forrest' resistance to the soybean cyst nematode is bigenic: saturation mapping of the Rhg1 and Rhg4 loci. Theor Appl Genet. 2001, 103: 710-717. 10.1007/s001220100597.View ArticleGoogle Scholar
- Concibido VC, Diers BW, Arelli PR: A decade of QTL mapping for cyst nematode resistance in soybean. Crop Sci. 2004, 44: 1121-1131. 10.2135/cropsci2004.1121.View ArticleGoogle Scholar
- Li Y, Zhang C, Gao Z-S, Johannes M, Smulders M, Ma Z, Liu Z-X, Nan H-Y, Chang R-Z Qiu L-J: Development of SNP markers and haplotype analysis of the candidate gene for rhg1, which confers resistance to soybean cyst nematode in soybean. Mol Breed. 2009, 24: 63-76. 10.1007/s11032-009-9272-0.View ArticleGoogle Scholar
- Bekal S, Niblack TL, Lambert KN: A chorismate mutase from the soybean cyst nematode Heterodera glycines shows polymorphisms that correlate with virulence. Mol Plant Microb Interact. 2003, 16: 439-446. 10.1094/MPMI.2003.16.5.439.View ArticleGoogle Scholar
- Palmateer AJ, Schmidt ME, Stetina SR, Russin JS: Temperature effects on race determination in Heterodera glycines. Journal of Nematology. 2000, 32: 349-355.PubMed CentralPubMedGoogle Scholar
- Mahalingam R, Skorupska HT: Cytological expression of early response to infection by Heterodera glycines Ichinohe in resistant PI 437654 soybean. Genome. 1996, 39: 986-998. 10.1139/g96-123.View ArticlePubMedGoogle Scholar
- Davis EL, Hussey RS, Baum TJ: Getting to the roots of parasitism by nematodes. Trends Parasitol. 2004, 20: 134-41. 10.1016/j.pt.2004.01.005.View ArticlePubMedGoogle Scholar
- Gao BR, Allen G, Davis EL, Baum TJ, Hussey RS: Molecular characterization and developmental expression of a cellulose-binding protein gene in the soybean cyst nematode Heterodera glycines. Int J Parasitol. 2004, 34: 1377-1383. 10.1016/j.ijpara.2004.09.001.View ArticlePubMedGoogle Scholar
- Lightfoot DA: Soybean genomics: Developments through the use of cultivar Forrest. Internat J of Plant Genom. 2008, 2008: 1-22.View ArticleGoogle Scholar
- Hartwig EE, Epps JM: Registration of Forrest soybeans. Crop Sci. 1973, 13: 287-View ArticleGoogle Scholar
- Lightfoot DA, Njiti VN, Gibson PT, Kassem MA, Iqbal MJ, Meksem K: Registration of the Essex X Forrest recombinant inbred line mapping population. Crop Sci. 2005, 45: 1678-1681. 10.2135/cropsci2004.0279.View ArticleGoogle Scholar
- Njiti VN, Meksem K, Iqbal MJ, Johnson JE, Zobrist KF, Kilo VY, Lightfoot DA: Common loci underlie field resistance to soybean sudden death syndrome in Forrest, Pyramid, Essex, and Douglas. Theor Appl Genet. 2002, 104: 294-300. 10.1007/s001220100682.View ArticlePubMedGoogle Scholar
- Wang D, Arelli PR, Shoemaker RC, Diers BW: Loci underlying resistance to Race 3 of soybean cyst nematode in Glycine soja plant introduction 468916. Theor Appl Genet. 2001, 103: 561-566. 10.1007/PL00002910.View ArticleGoogle Scholar
- Guo B, Sleper DA, Arelli PR, Shannon JG, Nguyen HT: Identification of QTL associated with resistance to soybean cyst nematode races 2, 3 and 5 in soybean PI 90763. Theor Appl Genet. 2005, 111: 965-971. 10.1007/s00122-005-0031-2.View ArticlePubMedGoogle Scholar
- Yue P, Arelli PR, Sleper DA: Molecular characterization of resistance to Heterodera glycines in soybean PI438489B. Theor Appl Genet. 2001, 102: 921-92. 10.1007/s001220000453.View ArticleGoogle Scholar
- Wu X, Blake S, Sleper DA, Shannon JG, Cregan P, Nguyen HT: QTL, additive and epistatic effects for SCN resistance in PI 437654. Theor Appl Genet. 2009, 118: 1093-1105. 10.1007/s00122-009-0965-x.View ArticlePubMedGoogle Scholar
- Yuan J, Njiti VN, Meksem K, Iqbal MJ, Triwitayakorn K, Kassem MA, Davis GT, Schmidt ME, Lightfoot DA: Quantitative trait loci in two soybean recombinant inbred line populations segregating for yield and disease resistance. Crop Sci. 2002, 42: 271-277. 10.2135/cropsci2002.0271.View ArticlePubMedGoogle Scholar
- Kopisch-Obuch FJ, McBroom KL, Diers BW: Association between SCN resistance loci and yield in soybean. Crop Sci. 2005, 45: 956-965. 10.2135/cropsci2004.0441.View ArticleGoogle Scholar
- Brucker E, Niblack T, Kopisch-Obuch FJ, Diers BW: The effect of rhg1 on reproduction of Heterodera glycines in the field and greenhouse and associated effects on agronomic traits. Crop Sci. 2005, 45: 1721-1727. 10.2135/cropsci2004.0477.View ArticleGoogle Scholar
- Kopisch-Obuch FW, Diers BW: Segregation at the SCN resistance locus rhg1 in soybean is distorted by an association between the resistance allele and reduced field emergence. Theor Appl Genet. 2006, 112: 199-207. 10.1007/s00122-005-0104-2.View ArticlePubMedGoogle Scholar
- Afzal AJ: Structure-function analysis of a candidate receptor like kinase protein in soybean cyst nematode resistance and identification of accessory proteins involved in plant defense. 2007, PhD Dissertation. Southern Illinois University, Carbondale, 1-234.Google Scholar
- Afzal AJ, Wood AJ, Lightfoot DA: Plant receptor-like serine threonine kinases: roles in signaling and plant defense. Molec Plant Microb Interact. 2008, 21: 507-517. 10.1094/MPMI-21-5-0507.View ArticleGoogle Scholar
- Afzal AJ, Srour A, Saini N, El Shemy H, Lightfoot DA: The multigeneic Rhg1 Locus: A model for the effects on root development, nematode resistance and recombination suppression. Nature Preceedings. 2008, hdl:10101/npre.2008.2726.1Google Scholar
- Afzal AJ, Natarajan A, Saini N, Iqbal MJ, Geisler MA: El Shemy H, Mungur R, Willmitzer L, Lightfoot DA: The nematode resistance allele at the Rhg1 locus alters the proteome and metabolome of soybean roots. Plant Physiology. 2009, 151: 1264-1280. 10.1104/pp.109.138149.PubMed CentralView ArticlePubMedGoogle Scholar
- Afzal AJ, Saini N, Srour A, Hemmati N, El Shemy H, Lightfoot DA: Recombination suppression at the Rhg1 locus underlying nematode resistance. Theor Appl Genet. 2012, 124: 1027-1039. 10.1007/s00122-011-1766-6.View ArticlePubMedGoogle Scholar
- Iqbal MJ, Ahsan R, Afzal AJ, Jamai A, Meksem K, El Shemy H, Lightfoot DA: Analysis of the activity of the soybean laccase encoded within the Rhg1/Rfs2 locus. Curr Iss Mol Biol. 2009, 11: i11-19.Google Scholar
- Severin A, Woody J, Bolon Y-T, Joseph B, Diers B, Farmer A, Muehlbauer G, Nelson R, Grant D, Specht J, Graham M, Cannon S, May G, Vance C, Shoemaker R: RNA-Seq atlas of Glycine max: A guide to the soybean transcriptome. BMC Plant Biol. 2010, 10: 160-176. 10.1186/1471-2229-10-160.PubMed CentralView ArticlePubMedGoogle Scholar
- Hauge BM, Wang ML, Parsons JD, Parnell LD: Methods of introgressing nucleic acid molecules associated with soybean cyst nematode resistance into soybean. US Patent. 2006, 7: 154,021-Google Scholar
- Lightfoot DA, Meksem K: Isolated soybean cyst nematode and sudden death syndrome polypeptides. US Patent. 2011, 7: 902,337-Google Scholar
- Replogle A, Wang J, Bleckmann A, Hussey RS, Baum TJ, Sawa S, Davis EL, Wang X, Simon R, Mitchum MG: Nematode CLE signaling in Arabidopsis requires CLAVATA2 and CORYNE. Plant J. 2011, 65: 430-440. 10.1111/j.1365-313X.2010.04433.x.View ArticlePubMedGoogle Scholar
- Melito S, Heuberger AL, Cook D, Diers BW, Macguidwin AE, Bent AF: A nematode demographics assay in transgenic roots reveals no significant impacts of the Rhg1 locus LRR-Kinase on soybean cyst nematode resistance. BMC Plant Biol. 2010, 10: 104-112. 10.1186/1471-2229-10-104.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, et al: Genome sequence of the palaeopolyploid soybean. Nature. 2010, 463: 178-183. 10.1038/nature08670.View ArticlePubMedGoogle Scholar
- Li Y, Zhang C, Gao Z-S, et al, et al: Genome sequence for 31 soybean genotypes, development of SNP markers and haplotype analysis. Nat Genet. 2010, 124: 63-76.Google Scholar
- Xiang Y, Cao Y, Xu C, Li X, Wang S: Xa3 conferring resistance for rice bacterial blight and encoding a receptor kinase-like protein is the same as Xa26. Theor Appl Genet. 2006, 113: 1347-1355. 10.1007/s00122-006-0388-x.View ArticlePubMedGoogle Scholar
- Nieto C, Piron F, Dalmais M, Marco CF, Moriones E, Gómez-Guillamón ML, Truniger V, Gómez P, Garcia-Mas J, Aranda MA, Bendahmane A: EcoTILLING for the identification of allelic variants of melon eIF4E a factor that controls virus susceptibility. BMC Plant Biol. 2007, 7: 34-44. 10.1186/1471-2229-7-34.PubMed CentralView ArticlePubMedGoogle Scholar
- Sessa G, D'Ascenzo M, Martin GB: Thr38 and Ser198 are Pto auto-phosphorylation sites required for the AvrPto–Pto-mediated hypersensitive response. EMBO J. 2000, 19: 2257-2269. 10.1093/emboj/19.10.2257.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang X, Goshe MB, Soderblom EJ, Phinney BS, Kuchar JA, Li J, Asami T, Yoshida S, Huber SC, Clouse SD: Identification and functional analysis of in vivo phosphorylation sites of the Arabidopsis BRASSINOSTEROID -INSENSITIVE1 receptor kinase. Plant Cell. 2005, 17: 1685-1703. 10.1105/tpc.105.031393.PubMed CentralView ArticlePubMedGoogle Scholar
- Lease KA, Lau NY, Schuster RA, Torii KU, Walker JC: Receptor serine/threonine protein kinases in signalling: Analysis of the erecta receptor-like kinase of Arabidopsis thaliana. New Phytol. 2001, 151: 133-143. 10.1046/j.1469-8137.2001.00150.x.View ArticleGoogle Scholar
- Liu X, Liu S, Jamai A, Bendahmane A, Lightfoot DA, Mitchum MG, Meksem K: Soybean cyst nematode resistance in soybean is independent of the Rhg4 locus LRR RLK gene. Funct Integ Genomics. 2011, DOI:10.1007/s10142-011-0225-4.Google Scholar
- Yesudas CR, Sharma H, Lightfoot DA: Identification of QTL in soybean underlying resistance to herbivory by Japanese beetles: Popillia japonica, Newman. Theor Appl Genet. 2010, 121: 353-362. 10.1007/s00122-010-1314-9.View ArticlePubMedGoogle Scholar
- Makandar R, Nalam V, Chaturvedi R, Jeannotte R, Sparks AA, Shah J: Involvement of salicylate and jasmonate signaling pathways in Arabidopsis interaction with Fusarium graminearum. Molec Plant Microb Interact. 2010, 23: 861-870. 10.1094/MPMI-23-7-0861.View ArticleGoogle Scholar
- Hyten DL, Choi IY, Song Q, et al: Highly variable patterns of linkage disequilibrium in multiple soybean populations. Genetics. 2007, 175: 1937-1944. 10.1534/genetics.106.069740.PubMed CentralView ArticlePubMedGoogle Scholar
- Afzal AJ, Lightfoot DA: Soybean disease resistance protein RHG1-LRR domain expressed, purified and refolded from Escherichia coli inclusion bodies: preparation for a functional analysis. Protein Expr Purif. 2007, 53: 346-355. 10.1016/j.pep.2006.12.017.View ArticlePubMedGoogle Scholar
- Afzal AJ, Srour A, Liu T, Vasudaven S, Samudrala R, Lightfoot DA: A model structure, ligand binding and dimerization of allotypes of a leucine rich repeat domain associated with soybean resistance to cyst nematodes and Fusarium virguliforme. BMC Plant Biol. 2012, (in review)Google Scholar
- Simmonds DH: Genetic transformation of soybean with biolistics. Molecular Methods of Plant Analysis. Vol. 23. Edited by: Jackson JF, Linskens HF. 2003, Springer-Verlag Berlin, Heidelberg Germany, 159-174.Google Scholar
- Shopinski K, Iqbal J, Afzal J, Shultz J, Jayaraman D, Lightfoot DA: Development of a pooled probe method for locating small gene families in a physical map of soybean using stress related paralogues and a BAC minimum tile path. Plant Methods. 2006, 2: 20-28. 10.1186/1746-4811-2-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuan J, Haroon M, Lightfoot DA, Pelletier Y, Liu Q, Li X-Q: A high-resolution melting approach for analyzing allelic expression dynamics. Curr Issue Mol Biol. 2008, 11: i1-9.Google Scholar
- Gritz L, Davies J: Plasmid-encoded hygromycin B resistance : The sequence of hygromycin B phosphotransferase gene and its expression in Escherichia coli and Saccharomyces cerevisiae. Gene. 1983, 25: 179-188. 10.1016/0378-1119(83)90223-8.View ArticlePubMedGoogle Scholar
- Gamborg OL, Miller RA, Ojima K: Nutrient requirements of suspension cultures of soybean root cells. Exp Cell Res. 1968, 50: 151-158. 10.1016/0014-4827(68)90403-5.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.