Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

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

Ali Srour12, Ahmed J Afzal123, Laureen Blahut-Beatty4, Naghmeh Hemmati1, Daina H Simmonds4, Wenbin Li5, Miao Liu5, Christopher D Town6, Hemlata Sharma17, Prakash Arelli8 and David A Lightfoot1259*

Author affiliations

1 Department of Molecular Biology, Microbiology and Biochemistry, Southern Illinois University at Carbondale, Carbondale, IL 62901, USA

2 Department of Plant Soil and Agricultural Systems, Southern Illinois University at Carbondale, Carbondale, IL 62901-4415, USA

3 Department of Horticulture and Crop Science, Ohio State University, 2021 Coffey Rd, Columbus, OH 43210, USA

4 Agriculture and Agri-Food Canada, Building 21, 960 Carling Ave, Ottawa, ON K1A 0C6, USA

5 Key Laboratory of Soybean Biology in the Chinese Ministry of Education, Harbin University, Harbin, China

6 JCVI, Rockville, MD, USA

7 Department of Plant Breeding & Genetics, Rajasthan College of Agriculture, MPUAT, Udaipur, India

8 USDA, Crop Genetics Research Unit, Jackson, TN, USA

9 Genomics Core Facility; Center for Excellence the Illinois Soybean Center, Southern Illinois University at Carbondale, Carbondale, IL 62901-4415, USA

For all author emails, please log on.

Citation and License

BMC Genomics 2012, 13:368  doi:10.1186/1471-2164-13-368

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/368


Received:17 November 2011
Accepted:12 June 2012
Published:2 August 2012

© 2012 Srour et al.; licensee 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.

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords:
Segregation; Pleiotropy; Rhg1/Rfs2; Soybean; Resistance; Soybean cyst nematode (SCN); Sudden death syndrome (SDS); Insect herbivory

Background

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 [1]. 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 [2]. 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 [5]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 [9]; 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) [13]. Interestingly most of these PIs were also partially resistant to SDS [14]. 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. virguliformeH. 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 [27] 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’ [13]. Three recessive loci (rhg1rhg3) 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 ([17] 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) [22], 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 [42]. 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 [21]. 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 [49]. 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 [53]. 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 [53]. 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.

Results

Allelic variations at the Rhg1/Rfs2 locus

The Forrest (resistance) allele of the Rhg1/Rfs2 locus was analyzed by sequencing the entire BAC B73p06 to 8 fold redundancy (Figure 1; GenBank HQ008938). The BAC sequence was compared to the sequences of (susceptible) Williams 82 whole genome shotgun sequence [54] and a BAC contig developed from (susceptible) ‘Asgrow 3244’ [50] from the same region. The BAC insert encompassed 82,157 bp that was predicted to encode nine genes comprised of 57 exons (2–8 per gene; Table 1). The GmRLK18-1 was gene 5, the laccase gene 6 and the antiporter gene 7. Sequences of the resistant and susceptible alleles were colinear with no large (more than 100 bp) insertions, deletions or inversions. Across the entire BAC there were 800 SNPs between Forrest and the sequences of alleles for susceptibility. SNPs were found in the promoter and enhancer regions of all 9 genes. However there were only 31 SNPs within genes. Only 11 of the 31 caused amino acid changes. There were just 6 of the 9 proteins changed by those 11 SNPs.

thumbnailFigure 1. Marker map of the genomic region aroundRhg1/Rfs2and the most identical paralog (homeolog) ofRhg1/Rfs2with locus, BAC and gene ideograms. Panel A shows the marker map of the genomic region around Rhg1/Rfs2 (Lg G; chromosome 18) with locus ideograms. Sequence coordinates were from the susceptible cultivar Asgrow 3244 [16,49]. The GmRLK18-1 gene encoding the RLK was shown as a black block arrow. The genes encoding the laccase and antiporter were shown as opposite white block arrows. All other genes were shown as grey block arrows. Locations of overlapping BAC clones B73p06 and B21d09 that both encoded the GmRLK18-1 at the Rhg1/Rfs2 locus (Lg G; chromosome 18) were shown below the ideogram. Panel B shows the BAC clone B73p06 that encoded the Rhg1/Rfs2-a locus. The gene encoding the RLK was shown as a black block arrow. The genes encoding the laccase and antiporter were shown as opposite white block arrows. All other genes were shown as grey block arrows. The extent of the pSBHB94 (HQ008939), the 9.772 kbp subclone from BAC B21d09 that was used for soybean transformations was shown as a blue arrow. The plasmid pSBHB94 encompassed from 30,423- 40–194 bp. Sequence coordinates were from the complete sequence of the BAC derived from resistant cultivar Forrest (HQ008938). Panel C showed a syntenic homeolog of Rhg1-a /Rfs2/ found in the sequence of BAC H38F23 from Lg B1 (chromosome 11). The homeolog of the gene encoding the RLK was show as a black block arrow. The homeologs of the genes encoding the laccase and antiporter were shown as opposite white block arrows. All other syntenic genes were shown as grey block arrows. The marker TMD1 amplified a fragment from Rhg1-a /Rfs2/Rhg1-a/Rfs2 of 303+ 15 bp (resistant allele was the smaller) and of 362 bp from a syntenic homeolog of Rhg1-a /Rfs2Rfs2/ found in the sequence of BAC H38f23 from Lg B1 (chromosome 11). Sequence coordinates were from the complete sequence of the BAC derived from resistant cultivar Forrest (HQ008940). Panel D shows an ideogram of markers, the GmRLK18-1 gene and gene features found in the insert of plasmid pSBHB94 that encompassed from 30,423- 40,194 bp of B73p06. Panel E shows the position of markers in the transcribed region of the GmRLK18-1 gene. Panel F shows the position of amino acid substitutions in the protein encoded by the cDNA.

Table 1. Genes predicted from DNA sequence within theRhg1-a/Rfs2locus encoded by BAC 73p06 (81, 157 bp) from Forrest (resistant) and sequence contig from A3244 (susceptible; AX196295)

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’) [49]. 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.

Additional file 1:. Figure S1. Paralogs ofRhg1in the soybean genome. Panel (A) shows LRR probe (200bp) hybridized to Forrest MTP. (B) Southern hybridization of LRR probe (200bp) to the MTP positives. Five out of the 7 MTP clones hybridized after BAC clone purification and restriction digestion with HindIII. The lower panel (C) shows the kinase domain probe (200bp) hybridized to Forrest MTP. Panel (D) shows Southern hybridization of the same kinase probe to the MTP positives. Three out of 5 MTP clones hybridized after BAC clone purification and restriction digestion with HindIII. Panel F shows ideograms of the genes predicted from the genome sequences centered on GmRLK18-1 and GmRLK11-1 Panel F shows an alignment of the genome sequences of 70kbp centered on GmRLK18-1 and GmRLK11-1 showing the extent of synteny.

Format: PPT Size: 1.2MB Download file

This file can be viewed with: Microsoft PowerPoint ViewerOpen Data

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 [16] and [25], 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 2: 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.

Additional file 2:. Table S1. Comparisons of sequence identity between Forrest alleles of GmRLK18-1 and GmRLK11-1 the most similar and syntenic RLK like protein. The amino acid identity was 78% in the signal peptide (residues 1–61); 94% in the ten LRRs (141–471), 93% in the transmembrane domain (485–507) and 97% in the kinase domain (569–840). Residues that differ in alloproteins of GmRLK18-1 are in bold. Four of the six are identical in the homeoprotein the other 2 are identical to the susceptible allele. The neighboring laccase and the antiporter also showed 85–96% amino acid identity.

Format: DOC Size: 34KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

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 [47].

Allele and paralog discrimination in NILs and transgenic plants

Since paralogs with homeolog sequence variants (HSVs) appeared to exist for each gene in the cluster, it was important to distinguish alleles precisely and separately from their most similar HSVs. In the NIL population RLK alleles and HSVs were distinguished using both the SIUC-TMD1 marker (Figure 2) and SNP Ala87Val from the LRR region by Taqman (Additional file 3: Figure S2; Additional file 4: Table S2). In genomic DNA TMD1 distinguished both the Peking and Peking-like allele of GmRLK18-1 from the other 7 alleles and the homeologous alleles of GmRLK11-1. To distinguish Forrest from both Essex and X5 the SNP probe to Ala87Val was again used for both genomic DNA and cDNA.

Additional file 3:. Figure S2. Detection of the SNP polymorphism at position 1486 in the LRR region ofRhg1–a and –e using an allelic discriminatory assay. A Famlabeled probe was used for the detection of resistant haplotypes1 and 2 (red) and Hex labeled probe for the detection of susceptible haplotypes2, 3 and 4 (blue). A total of 16 individuals from the 110 PIs were selected for the analysis. The Panel shows relative fluorescent signal intensity for each of the 16 plant introductions. The two groups form separate clusters.

Format: PDF Size: 250KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 4:. Table S2. The sequence of the microsatellite primers that were used forRhg1-afine map development (from Triwitayakorn et al., 2005).

Format: DOC Size: 161KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

thumbnailFigure 2. TheRhg1/Rfs2genomic region altered root development. Panel A; soybean NILs at 2 weeks pre-SCN inoculation show different root morphologies, Rhg1/Rfs2 inhibited germination and early root growth. Panel B; the root morphologies co-segregated with the allele of GmRLK18-1 at the Rhg1/Rfs2 locus as shown by the 3 bands generated by intron located marker TMD1 (arrowed). R indicated the resistant allele from Forrest; S the susceptible allele from Essex; and H indicated the heterozygous state for the Rhg1/Rfs2 locus. Panel C showed the distribution of root mass (rfw was root fresh weight) among 20 heterozygous plants. Panel D showed the distribution of root mass among 20 susceptible plants. Panel E showed the distribution of root mass among 20 resistant 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% [23]. 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 [21].

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 [16] 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 [44]. 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 [44] and so protein phosphorylation.

Transcripts of the GmRLK18-1 alleles at Rhg1/Rfs2 were found in soybean roots

Expression judged by examination of EST libraries in silico, cDNA libraries by hybridization, RNA sequencing and mRNA populations assayed by qRT-PCR showed both transcripts encoded by GmRLK18-1-a and GmRLK18-1-e alleles (Figure 3) were present in both non-infested roots and SCN-infested roots, as reported in [16]. The transcript and protein abundances were not increased by more than 2 fold in either H. glycines or F. virguliforme infested plants compared to non-infested plants (Figure 3). Low transcript and protein abundances appeared to result from largely constitutive expression under the conditions and stages of development tested to date.

thumbnailFigure 3. Expression of Forrest (−a) and Essex (−e) alleles of theGmRLK18-1atRhg1/Rfs2in non-transgenic soybean NIL roots. Panel (A); Western hybridization using an anti-RHG1/RFS2 (GmRLK18-1) antibody (Afzal et al. 2007) from roots of; Forrest allele in ExF34-23 (a); Essex allele in ExF34-3(b); E.coli expressed proteins RHG1/RFS2-LRR-Shrt (c); RHG1/RFS2-LRR-Long (d); and RHG 4 (e). Panel (B); agarose gel electrophoresis of cDNA amplified using Rhg1/Rfs2 LRR flanking primers from; non-infested NIL 34–23 (1); SCN infested NIL 34–23 (2); SCN infested NIL 34–3 (3); and Forrest genomic DNA (4). Negative control without template is shown in lane 5.

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.

In several transgenic lines from (T1-T3) used for SDS and SCN assays, lines were identified with the GmRLK18-1-a transgene in the homozygous state. In these plants the transgenes were both transcribed and translated. For example, in the progeny of fertile T0 line 6B3 the progeny T1 lines like 7D2 and sixteen T3 derived plants expressed the GmRLK18-1-a transgene as both mRNA and protein (Figure 4). Expression of the transgene alleles was equal to the endogenous gene alleles as judged by protein abundance. GmRLK1 transcripts and proteins were detected in about 50% of the progenies of the 6B3 primary transformants. In selfed lines expressing GmRLK18-1-a the transgene allele was detected as distinct mRNAs and alloproteins. Expression was directed by the native promoter and the enhancer elements contained on the proximal 6 kbp of the 9.772 kbp fragment of pSBHB94. That plasmid was sub-cloned from BAC B21d09 (Figure 1; HQ008939) during BAC sequencing.

thumbnailFigure 4. Soybean transgenic plants expressed the mRNA and protein from the Forrest allele ofGmRLK18-1, the RLK at theRfs2/Rhg1locus. Panel A shows PCR from leaf samples of progeny plants derived from a primary transgenic event 6B3-7D2(1) with TMD1 primers. Lanes contain transgenic plants 1 to 13. The arrow shows the double band for Gm18RLK-1-a positive sample at 314 bp for lines 1,3–5, 7,8,10 and 12. M was the marker; H was the no DNA (water) control; P was the Rhg1 plasmid pSBHB94; X5 was the control plant. Panel B shows PCR from cDNA leaf samples of sixteen transgenic lines derived from event 6B3-7D2(1) with HRM primers. Green lines are from transgenic plants. Red melt curve was a resistant control blue line was a susceptible control. Panel C shows a Western of a 2D gel from roots of a transgenic plant probed with the anti-RLK peptide antibody. An alloprotein at pI 8.42 and 92.41 kDa was found in the non transgenic cv X5 but the presence of the Forrest alloprotein at pI 8.44 and and 92.39 kDa was found in transgenic plants derived from event 6B3-7D2(1) expressing GmRLK18-1-a. GmRLK18-1 was shown to be a very low abundance protein impossible to visualize without immune-staining.

In transgenic soybean plants GmRLK18-1-a allele provided resistance to both root infection and root rot by Fusarium virguliforme. That root resistance underlay a significant reduction in leaf symptoms and delay of the senescence caused by SDS (Figure 5; Table 2B). Resistance to SDS was effective throughout the life of the plants which flowered and set seed. The non-transgenic X5 plants proved to be highly susceptible to SDS and showed all the expected phenotypes of root rot and leaf scorch. The phenotypes among the susceptible plants included a gradual worsening of leaf scorch symptoms from 3.0 at 21 days after infection (dai) to 8.5 by 56 dai. Only at near maturity did the susceptible plants showed symptoms characteristic of SDS. For example, leaflet abscission occurred from the top of the petiole instead of the base of the petiole. However, the transgenic plants did not show early senescence or reduced pod set. Equally, the characteristic root symptoms of SDS were very reduced in the RLK transgenic plants. Symptoms reduced included the degree of root rot and browning of the root cortex. Senescence was delayed by 14 days in GmRLK18-1 transgenics compared to controls and more pods were set in each of the three repeats of the experiment. Therefore, the GmRLK18-1-a provided a very high degree of resistance to the transgenic plants in roots and leaves and pods.

thumbnailFigure 5. Resistance toFusarium virguliformeand partial resistance toHeterodera glycinescaused by the Forrest allele of a receptor like kinase (GmRLK18-1-a) found at theRfs2/Rhg1locus as transgene in primary transgenic lines (cv ‘X5’). Panels A-H show the SDS assays. F. virguliforme was used at 104 cfu per cm3 of soil. The experiment was carried out on 3 separate occasions. Leaf scorch was scored as DS at 7, 14, 21, 28, 35, 42, 49 and 56 dai (days after infestation). Derived crossed lines (X5::GmRLK18-1-a x WENIL35 and cultivars X5:: GmRLK18-1-a x WENIL35xEF2) were included in runs 2 and 3. Panel A &B shows stable soybean transgenics with and without the 10 kbp GmRLK18-1 (Rhg1/Rfs2) subclone at 21 dai. Panel C shows the GmRLK18-1-a transgene reduced root rot at 28 dai. Panel D shows leaf symptoms at 28 dai. Panel E shows plants at 56 dai where X5 is senescent with abscission of leaflets from erect petioles and X5::GmRLK18-1-a is still green and filling pods. Panel F shows selected leaflets at 28 dai with a 1–9 range in DS scores arranged in order of severity from bottom left to top right. Panels G-J show the SCN assays. Panel G shows SCN development by the RLK in X5 transgenics. Panel H shows SCN development by the Rhg1-a allele in resistant NIL 34–23. Panel I shows normal SCN development in the susceptible X5. Panel J shows normal SCN development in susceptible NIL 34–3. Panel K shows a GmRLK18-1-a transgenic plant of cultivar X5 that was defoliated by insect herbivory in the field. Panel L shows an non-transgenic X5 plant with much less leaf area loss.

Table 2. 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

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.

Discussion

The GmRLK18-1 at the Rhg1/Rfs2 locus was shown to underlie resistance to root infection by F. virguliforme and the subsequent leaf scorch, as predicted by [11]. Some resistance to SCN was found but this was partial as expected in the absence of a resistance allele at Rhg4 [23,61]. The RLK simultaneously contributed resistance to both pathogens, establishing pleiotropy (Figure 6). Reduced seedling root growth was part of the resistance mechanism and this may underlie the reduced yield of resistant cultivars [39-41]. In Arabidopsis CLAVATA1, the RLK that regulated meristem development, also had an effect on nematode resistance [52].

thumbnailFigure 6. Genetic model for the function of theRhg1/Rfs2locus in association with other known resistance loci. Black arrows show positive interactions, blue arrows show inhibitions. In this model four phenotypic events are controlled by the 3 genes at Rhg1/Rfs2 and four unlinked genes Rhg2-4 and Rzd1 (sup-Rhg1/Rfs2). Root growth inhibition occurs despite suppression of the RLK by Rzd1. Pleiotropy or close linkage between loci underlying resistance was found at Rhg1/Rfs2 (this work) and Rhg5/Rfs3 [8,17].

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 [62]. Arabidopsis resistance to F. graminearum was recently shown to require operation of both SA and JA pathways [63]. 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 [46]. 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.

Conclusions

The Rhg1/Rfs2 locus was shown to include the GmRLK18-1. Previously [53] 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 [25] were confirmed by BAC fragment sequencing here. Fortunately the amino acid changes made to [16] 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 [52]. 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 [2]. 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 [53]. 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 [53] 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.

Methods

Plant materials

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 [67]. 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.

SCN inoculations

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 [16]. 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 [5]. 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 [11]. 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 [39] and on agarose gels after [8].

BAC sequence and allele comparisons

Southern hybridizations were performed following the standard procedure described in [68] 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 [54].

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 [69] 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 [46]. 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).

Transgenic Plants

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 [67]. 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 HygR[70] 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 [71]; 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 [67]. 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 [46] followed by Western hybridization [43,44].

Field trials were conducted at the ARC in Carbondale during 2010 and 2011 using conditions described in [62]. Plants were arrayed in 12′ plots arrayed in a randomized complete block with 4 replications. Insect herbivory was measured as described in [62]. 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.

Abbreviations

RLK, Receptor like kinase; SCN, soybean cyst nematode; SDS, sudden death syndrome; BAC, bacterial artificial chromosome; MAFF, Genebank System National Institute of Agrobiological Sciences, Tsukuba, Japan; NRRL, The Agriculture Research Service Culture Collection National Center for Agricultural Utilization Research, USDA/ARS, Peoria, IL USA.

Competing interest

Authors declare that they have no competing interest.

Authors’ contributions

AS and AJA contributed equally to the NILs analyses. AS was responsible for the transgenic plant analyses in vitro, DALightfoot for the assays in vivo. LB-B and DHS generated all the transgenic plants described here. WL, ML and PA provided SCN assays with YC3 and PA3. NH, and DAL analyszed BAC sequences and generated Figure 1 and Additional file 1: Table S1. CDT provide the BAC sequences. HS carried out the SDS assays on NILs. DAL wrote the manuscript and edited the final version. All authors read and approved the final manuscript.

Data depositions

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.

Acknowledgements

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.

References

  1. Roy KW, Rupe JCD, Hershman DE, Abney TS: Sudden death syndrome of soybean.

    Plant Dis 1997, 81:1100-1111. Publisher Full Text OpenURL

  2. 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. Publisher Full Text OpenURL

  3. 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. OpenURL

  4. Kolander TM: The host range of Fusarium virguliforme on rotational crops and common plant species and its survival and growth on crop residue. University of Minnesota, MS Thesis; 2010:196.

    http://purl.umn.edu/60700 webcite

    OpenURL

  5. 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. Publisher Full Text OpenURL

  6. 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. PubMed Abstract | Publisher Full Text OpenURL

  7. 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. Publisher Full Text OpenURL

  8. 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. PubMed Abstract | Publisher Full Text OpenURL

  9. 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]

    PubMed Abstract | Publisher Full Text OpenURL

  10. 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. PubMed Abstract | Publisher Full Text OpenURL

  11. 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. PubMed Abstract | Publisher Full Text OpenURL

  12. 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. PubMed Abstract | Publisher Full Text OpenURL

  13. Arelli PR: Inheritance of resistance to Heterodera glycines race 3 in soybean accessions.

    Plant Dis 1994, 78:898-900. Publisher Full Text OpenURL

  14. 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. Publisher Full Text OpenURL

  15. 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. Publisher Full Text OpenURL

  16. 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. Publisher Full Text OpenURL

  17. 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. PubMed Abstract | Publisher Full Text OpenURL

  18. Hussey RS, Grundler FMW: Nematode parasitism of plants. In The physiology and biochemistry of free-living and plant parasitic nematodes. Edited by Perry RN, Wright DJ. CAB International Wallingford, USA; 1998:213-243. OpenURL

  19. Caldwell BE, Brim CA, Ross JP: Inheritance of resistance of soybeans to the cyst nematode Heterodera glycines.

    Agron J 1960, 52:635-636. Publisher Full Text OpenURL

  20. 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. OpenURL

  21. 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. OpenURL

  22. 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. Publisher Full Text OpenURL

  23. 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. Publisher Full Text OpenURL

  24. Concibido VC, Diers BW, Arelli PR: A decade of QTL mapping for cyst nematode resistance in soybean.

    Crop Sci 2004, 44:1121-1131. Publisher Full Text OpenURL

  25. 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. Publisher Full Text OpenURL

  26. 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. Publisher Full Text OpenURL

  27. Palmateer AJ, Schmidt ME, Stetina SR, Russin JS: Temperature effects on race determination in Heterodera glycines.

    Journal of Nematology 2000, 32:349-355. PubMed Abstract | PubMed Central Full Text OpenURL

  28. 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. PubMed Abstract | Publisher Full Text OpenURL

  29. Davis EL, Hussey RS, Baum TJ: Getting to the roots of parasitism by nematodes.

    Trends Parasitol 2004, 20:134-41. PubMed Abstract | Publisher Full Text OpenURL

  30. 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. PubMed Abstract | Publisher Full Text OpenURL

  31. Lightfoot DA: Soybean genomics: Developments through the use of cultivar Forrest.

    Internat J of Plant Genom 2008, 2008:1-22. OpenURL

  32. Hartwig EE, Epps JM: Registration of Forrest soybeans.

    Crop Sci 1973, 13:287. OpenURL

  33. 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. Publisher Full Text OpenURL

  34. 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. PubMed Abstract | Publisher Full Text OpenURL

  35. 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. Publisher Full Text OpenURL

  36. 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. PubMed Abstract | Publisher Full Text OpenURL

  37. Yue P, Arelli PR, Sleper DA: Molecular characterization of resistance to Heterodera glycines in soybean PI438489B.

    Theor Appl Genet 2001, 102:921-92. Publisher Full Text OpenURL

  38. 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. PubMed Abstract | Publisher Full Text OpenURL

  39. 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. PubMed Abstract | Publisher Full Text OpenURL

  40. Kopisch-Obuch FJ, McBroom KL, Diers BW: Association between SCN resistance loci and yield in soybean.

    Crop Sci 2005, 45:956-965. Publisher Full Text OpenURL

  41. 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. Publisher Full Text OpenURL

  42. 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. PubMed Abstract | Publisher Full Text OpenURL

  43. 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. PhD Dissertation. Southern Illinois University, Carbondale; 2007:1-234. OpenURL

  44. 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. Publisher Full Text OpenURL

  45. 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.1

    OpenURL

  46. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. 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. PubMed Abstract | Publisher Full Text OpenURL

  48. 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. OpenURL

  49. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  50. 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. OpenURL

  51. Lightfoot DA, Meksem K: Isolated soybean cyst nematode and sudden death syndrome polypeptides.

    US Patent 2011, 7:902,337. OpenURL

  52. 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. PubMed Abstract | Publisher Full Text OpenURL

  53. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  54. Schmutz J, Cannon SB, Schlueter J, et al.: Genome sequence of the palaeopolyploid soybean.

    Nature 2010, 463:178-183. PubMed Abstract | Publisher Full Text OpenURL

  55. Li Y, Zhang C, Gao Z-S, et al.: Genome sequence for 31 soybean genotypes, development of SNP markers and haplotype analysis.

    Nat Genet 2010, 124:63-76. OpenURL

  56. 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. PubMed Abstract | Publisher Full Text OpenURL

  57. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  58. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. 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. Publisher Full Text OpenURL

  61. 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. Publisher Full Text OpenURL

  62. 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. PubMed Abstract | Publisher Full Text OpenURL

  63. 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. Publisher Full Text OpenURL

  64. Hyten DL, Choi IY, Song Q, et al.: Highly variable patterns of linkage disequilibrium in multiple soybean populations.

    Genetics 2007, 175:1937-1944. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. 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. PubMed Abstract | Publisher Full Text OpenURL

  66. 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)

    OpenURL

  67. Simmonds DH: Genetic transformation of soybean with biolistics. In Molecular Methods of Plant Analysis. Vol. 23. Edited by Jackson JF, Linskens HF. Springer-Verlag Berlin, Heidelberg Germany; 2003:159-174. OpenURL

  68. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  69. 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. OpenURL

  70. 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. PubMed Abstract | Publisher Full Text OpenURL

  71. Gamborg OL, Miller RA, Ojima K: Nutrient requirements of suspension cultures of soybean root cells.

    Exp Cell Res 1968, 50:151-158. PubMed Abstract | Publisher Full Text OpenURL