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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1531-3) contains supplementary material, which is available to authorized users.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean

Copy number variation (CNV) is implicated in important traits in multiple crop plants, but can be challenging to genotype using conventional methods. The Rhg1 locus of soybean, which confers resistance to soybean cyst nematode (SCN), is a CNV of multiple 31.2‐kb genomic units each containing four genes. Reliable, high‐throughput methods to quantify Rhg1 and other CNVs for selective breeding were developed. The CNV genotyping assay described here uses a homeologous gene copy within the paleopolyploid soybean genome to provide the internal control for a single‐tube TaqMan copy number assay. Using this assay, CNV in breeding populations can be tracked with high precision. We also show that extensive CNV exists within Fayette, a released, inbred SCN‐resistant soybean cultivar with a high copy number at Rhg1 derived from a single donor parent. Copy number at Rhg1 is therefore unstable within a released variety over a relatively small number of generations. Using this assay to select for individuals with altered copy number, plants were obtained with both increased copy number and increased SCN resistance relative to control plants. Thus, CNV genotyping technologies can be used as a new type of marker‐assisted selection to select for desirable traits in breeding populations, and to control for undesirable variation within cultivars.     

Cloning and sequence diversity analysis of <Emphasis Type="Italic">GmHs1</Emphasis><Superscript><Emphasis Type="Italic">pro-1</Emphasis></Superscript> in Chinese domesticated and wild soybeans

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is the most important pathogen in soybean production worldwide and causes substantial yield losses. An apparent narrow genetic base of SCN resistance was observed in current elite soybean cultivars, and searching for novel SCN resistance genes as well as novel resistance sources rather than focusing on the two important genes rhg1 and Rhg4 has become another major objective in soybean research. In the present paper we report a 1,477 bp Hs1 ^pro-1 homolog, named GmHs1 ^pro-1 . This gene was cloned from soybean variety Wenfeng 7 based on information for Hs1 ^pro-1 , a beet cyst nematode resistance gene in sugar beet. It has two domains, Hs1pro-1_N and Hs1pro-1_C, both of which are believed to confer resistance to nematodes. Of the 1,477 bp sequence in GmHs1 ^pro-1 , an open reading frame of 1,314 bp, encoding a protein with 437 amino acids, was flanked by a 5′-untranslated region of 27 bp and a 3′-untranslated region of 135 bp. Fourteen single-nucleotide polymorphisms (SNPs) were observed in 44 soybean accessions including 23 wild soybeans, 8 landraces, and 13 soybean varieties (or lines), among which 5 in wild soybeans and 3 in landrace accessions were unique. Sequence diversity analysis on the 44 soybean accessions showed π = 0.00168 and θ = 0.00218 for GmHs1 ^pro-1 ; landraces had the highest diversity, followed by wild soybeans, with varieties (or lines) having the lowest. Although we did not detect a significant effect of selection on GmHs1 ^pro-1 in the three populations, sequence diversity, unique SNPs, and phylogenetic analysis indicated a slight domestication bottleneck and an intensive selection bottleneck. High sequence diversity, more unique SNPs, and broader representation across the phylogenetic tree in wild soybeans and landraces indicated that wild collections and landrace accessions are invaluable germplasm for broadening the genetic base of elite soybean varieties resistant to SCN. C. Yuan and G. Zhou contributed to this paper equally.     

Transgenic soybean overexpressing GmSAMT1 exhibits resistance to multiple‐HG types of soybean cyst nematode Heterodera glycines

Soybean (Glycine max (L.) Merr.) salicylic acid methyl transferase (GmSAMT1) catalyses the conversion of salicylic acid to methyl salicylate. Prior results showed that when GmSAMT1 was overexpressed in transgenic soybean hairy roots, resistance is conferred against soybean cyst nematode (SCN), Heterodera glycines Ichinohe. In this study, we produced transgenic soybean overexpressing GmSAMT1 and characterized their response to various SCN races. Transgenic plants conferred a significant reduction in the development of SCN HG type 1.2.5.7 (race 2), HG type 0 (race 3) and HG type 2.5.7 (race 5). Among transgenic lines, GmSAMT1 expression in roots was positively associated with SCN resistance. In some transgenic lines, there was a significant decrease in salicylic acid titer relative to control plants. No significant seed yield differences were observed between transgenics and control soybean plants grown in one greenhouse with 22 °C day/night temperature, whereas transgenic soybean had higher yield than controls grown a warmer greenhouse (27 °C day/23 °C night) temperature. In a 1‐year field experiment in Knoxville, TN, there was no significant difference in seed yield between the transgenic and nontransgenic soybean under conditions with negligible SCN infection. We hypothesize that GmSAMT1 expression affects salicylic acid biosynthesis, which, in turn, attenuates SCN development, without negative consequences to soybean yield or other morphological traits. Thus, we conclude that GmSAMT1 overexpression confers broad resistance to multiple SCN races, which would be potentially applicable to commercial production.     

Soybean cyst nematode resistance in soybean is independent of the <Emphasis Type="Italic">Rhg4</Emphasis> locus <Emphasis Type="Italic">LRR-RLK</Emphasis> gene

To test the function of candidate genes in soybean for resistance to the soybean cyst nematode (SCN), a large collection of EMS-mutants from the SCN-resistant soybean cultivar “Forrest” was developed for Targeting Induced Local Lesions IN Genomes (TILLING). Additionally, due to the complexity of the soybean genome, an integrated set of genomic and genetic analysis tools was employed to complement the TILLING approach. The efficiency of this integrated set of tools was tested using a candidate soybean gene for resistance to SCN, encoding a leucine-rich repeat receptor-like kinase (LRR-RLK) that was identified by map-based cloning at the Rhg4 locus. The Rhg4 locus is one of the major quantitative trait loci controlling soybean resistance against SCN race 3 (HG type 0) in cv. Forrest, but the gene(s) sequence for resistance remains to be determined. Using TILLING, a Forrest mutant containing a nonsense mutation in the LRR domain of the candidate resistance protein was identified and confirmed; however, the SCN-resistant phenotype of the mutant was not altered. Haplotyping and EcoTILLING of recombinant inbred lines along with complementation analysis corroborated the TILLING result and ruled out the possibility of functional redundancy by a second copy of the LRR-RLK gene identified in the soybean genome. This study validates the use of TILLING, in combination with an integrated set of genomic tools, as an efficient means of testing candidate genes for SCN resistance in soybean.     

Recombination suppression at the dominant <Emphasis Type="Italic">Rhg1/Rfs2</Emphasis> locus underlying soybean resistance to the cyst nematode

Host resistance to “yellow dwarf” or “moonlight” disease cause by any population (Hg type) of Heterodera glycines I., the soybean cyst nematode (SCN), requires a functional allele at rhg1. The host resistance encoded appears to mimic an apoptotic response in the giant cells formed at the nematode feeding site about 24–48 h after nematode feeding commences. Little is known about how the host response to infection is mediated but a linked set of 3 genes has been identified within the rhg1 locus. This study aimed to identify the role of the genes within the locus that includes a receptor-like kinase (RLK), a laccase and an ion antiporter. Used were near isogeneic lines (NILs) that contrasted at their rhg1 alleles, gene-based markers, and a new Hg type 0 and new recombination events. A syntenic gene cluster on Lg B1 was found. The effectiveness of SNP probes from the RLK for distinguishing homolog sequence variants on LgB1 from alleles at the rhg1 locus on LgG was shown. The resistant allele of the rhg1 locus was shown to be dominant in NILs. None of the recombination events were within the cluster of the three candidate genes. Finally, rhg1 was shown to reduce the plant root development. A model for rhg1 as a dominant multi-gene resistance locus based on the developmental control was inferred.     

Rhg1 alleles from soybean PI 437654 and PI 88788 respond differentially to isolates of Heterodera glycines in the greenhouse

The production of resistant soybean [Glycine max (L.) Merr.] cultivars is the most effective means for controlling losses from soybean cyst nematode (SCN) (Heterodera glycines Ichinohe). The major resistance gene in most SCN resistance sources is rhg1, which has been mapped as a quantitative trait locus onto linkage group G. Our objective was to determine whether the SCN resistance sources PI 437654 and PI 88788 have different functional alleles at rhg1 based on resistance phenotypes. Populations segregating for resistance alleles at rhg1 from both PI 88788 and PI 437654 and at Rhg4, a second SCN resistance gene from PI 437654, were developed. These populations were screened for resistance to the H. glycines inbred isolates PA3 (HG type 7) and TN14 (HG type 1.2.5.7) in the greenhouse and evaluated with molecular markers linked to both rhg1 and Rhg4. Each isolate test was repeated, and the evaluations were done on a single-plant and a line-mean basis in Test 1, and solely on a single-plant basis in Test 2. Across two tests with the TN14 isolate, plants with the PI 437654 allele for a marker linked to rhg1 had significantly (P<0.0001) less SCN reproduction than plants carrying the PI 88788 allele. A marker linked to Rhg4, however, was not significantly associated with resistance to TN14. Across two tests with the PA3 isolate, alleles of rhg1 from both sources gave a resistant reaction, although plants homozygous for the PI 88788 allele had significantly (P<0.05) greater resistance than plants with the PI 437654 allele. The marker allele from PI 437654 linked to Rhg4 was significantly (P<0.0005) associated with greater resistance than the PI 88788 allele in both PA3 tests, and resistance was dominant. There was a significant interaction between alleles at rhg1 and Rhg4 in both PA3 tests. These results suggest that PI 437654 and PI 88788 each have a different functional SCN resistance allele at or close to rhg1. These allelic differences have implications that breeders should consider before incorporation into cultivars.     

Overexpression of a soybean salicylic acid methyltransferase gene confers resistance to soybean cyst nematode

Salicylic acid plays a critical role in activating plant defence responses after pathogen attack. Salicylic acid methyltransferase (SAMT) modulates the level of salicylic acid by converting salicylic acid to methyl salicylate. Here, we report that a SAMT gene from soybean (GmSAMT1) plays a role in soybean defence against soybean cyst nematode (Heterodera glycines Ichinohe, SCN). GmSAMT1 was identified as a candidate SCN defence‐related gene in our previous analysis of soybean defence against SCN using GeneChip microarray experiments. The current study started with the isolation of the full‐length cDNAs of GmSAMT1 from a SCN‐resistant soybean line and from a SCN‐susceptible soybean line. The two cDNAs encode proteins of identical sequences. The GmSAMT1 cDNA was expressed in Escherichia coli. Using in vitro enzyme assays, E. coli‐expressed GmSAMT1 was confirmed to function as salicylic acid methyltransferase. The apparent Km value of GmSAMT1 for salicylic acid was approximately 46 μm . To determine the role of GmSAMT1 in soybean defence against SCN, transgenic hairy roots overexpressing GmSAMT1 were produced and tested for SCN resistance. Overexpression of GmSAMT1 in SCN‐susceptible backgrounds significantly reduced the development of SCN, indicating that overexpression of GmSAMT1 in the transgenic hairy root system could confer resistance to SCN. Overexpression of GmSAMT1 in transgenic hairy roots was also found to affect the expression of selected genes involved in salicylic acid biosynthesis and salicylic acid signal transduction.     

The soybean-Phytophthora resistance locus <Emphasis Type="Italic">Rps1</Emphasis>-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences

Background  

A series of Rps (resistance to P ytophthora s ojae) genes have been protecting soybean from the root and stem rot disease caused by the Oomycete pathogen, Phytophthora sojae. Five Rps genes were mapped to the Rps1 locus located near the 28 cM map position on molecular linkage group N of the composite genetic soybean map. Among these five genes, Rps1-k was introgressed from the cultivar, Kingwa. Rps1-k has been providing stable and broad-spectrum Phytophthora resistance in the major soybean-producing regions of the United States. Rps1-k has been mapped and isolated. More than one functional Rps1-k gene was identified from the Rps1-k locus. The clustering feature at the Rps1-k locus might have facilitated the expansion of Rps1-k gene numbers and the generation of new recognition specificities. The Rps1-k region was sequenced to understand the possible evolutionary steps that shaped the generation of Phytophthora resistance genes in soybean.     

Genetic mapping of soybean cyst nematode race-3 resistance loci in the soybean PI 437.654

Resistance to the soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is difficult to evaluate in soybean [Glycine max (L.) Merr.] breeding. PI 437.654 has resistance to more SCN race isolates than any other known soybean. We screened 298 F₆₇ recombinant-inbred lines from a cross between PI 437.654 and BSR101 for SCN race-3 resistance, genetically mapped 355 RFLP markers and the I locus, and tested these markers for association with resistance loci. The Rhg ₄ resistance locus was within 1 cM of the I locus on linkage group A. Two additional QTLs associated with SCN resistance were located within 3cM of markers on groups G and M. These two loci were not independent because 91 of 96 lines that had a resistant-parent marker type on group G also had a resistant-parent marker type on group M. Rhg ₄ and the QTL on G showed a significant interaction by together providing complete resistance to SCN race-3. Individually, the QTL on G had greater effect on resistance than did Rhg ₄, but neither locus alone provided a degree of resistance much different from the susceptible parent. The nearest markers to the mapped QTLs on groups A and G had allele frequencies from the resistant parent indicating 52 resistant lines in this population, a number not significantly different from the 55 resistant lines found. Therefore, no QTLs from PI 437.654 other than those mapped here are expected to be required for resistance to SCN race-3. All 50 lines that had the PI 437.654 marker type at the nearest marker to each of the QTLs on groups A and G were resistant to SCN race-3. We believe markers near to these QTLs can be used effectively to select for SCN race-3 resistance, thereby improving the ability to breed SCN-resistant soybean varieties.     

High-throughput genotyping for a polymorphism linked to soybean cyst nematode resistance gene Rhg4 by using TaqmanTM probes

An individual soybean breeder can generate over one hundred thousand new genotypes each year. The efficiency of selection in these populations could be improved if these genotypes were effectively screened with one DNA marker that identified an important gene, and if laboratory throughput was high and costs were low. Our aim was to develop a rapid genotyping procedure for resistance to the soybean cyst nematode. A high-throughput genotyping method was developed with fluorogenic probes to distinguish between two insertion polymorphisms in alleles of an AFLP marker that is located about 50 kbp from the Rhg4 gene candidate. The assay uses the 5 exonuclease activity of Taq polymerase in conjunction with fluorogenic probes for each allele. The method can be used for scoring the polymorphism in a recombinant inbred line population and for screening parent lines in a breeding program. The Taqman^TM method of determining genotype was accurate in 90% of scores in the RIL population compared to 95% accuracy with electrophoresis. Among 94 cultivars that are parents in our breeding program allele 2 that is derived from the sources of resistance to SCN was common in resistant cultivars (30 of 56) but rare in susceptible cultivars (3 of 38). Therefore, this method can be applied to automated large-scale genotyping for soybean breeding programs.     

Genome-wide association study for soybean cyst nematode resistance in Chinese elite soybean cultivars

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a highly recalcitrant endoparasite of soybean roots, causing more yield loss than any other pest. To identify quantitative trait loci (QTL) controlling resistance to SCN (HG type 2.5.7, race 1), a genome-wide association study (GWAS) was performed. The association panel, consisting of 120 Chinese soybean cultivars, was genotyped with 7189 single nucleotide polymorphism (SNPs). A total of 6204 SNPs with minor allele frequency >0.05 were used to estimate linkage disequilibrium (LD) and population structure. The mean level of LD measured by r ² declined very rapidly to half its maximum value (0.51) at 220 kb. The overall population structure was approximately coincident with geographic origin. The GWAS results identified 13 SNPs in 7 different genomic regions significantly associated with SCN resistance. Of these, three SNPs were localized in previously mapped QTL intervals, including rhg1 and Rhg4. The GWAS results also detected 10 SNPs in 5 different genomic regions associated with SCN resistance. The identified loci explained an average of 95.5% of the phenotypic variance. The proportion of phenotypic variance was due to additive genetic variance of the validated SNPs. The present study identified multiple new loci and refined chromosomal regions of known loci associated with SCN resistance. The loci and trait-associated SNPs identified in this study can be used for developing soybean cultivars with durable resistance against SCN.     

Genomic analysis of the rhg1 locus: candidate genes that underlie soybean resistance to the cyst nematode

The rhg1 gene or genes lie at a recessive or co-dominant locus, necessary for resistance to all Hg types of the soybean (Glycine max (L.) Merr.) cyst nematode (Heterodera glycines I.). The aim here was to identify nucleotide changes within a candidate gene found at the rhg1 locus that were capable of altering resistance to Hg types 0 (race 3). A 1.5 ± 0.25 cM region of chromosome 18 (linkage group G) was shown to encompass rhg1 using recombination events from four near isogenic line populations and nine DNA markers. The DNA markers anchored two bacterial artificial chromosome (BAC) clones 21d9 and 73p6. A single receptor like kinase (RLK; leucine rich repeat-transmembrane-protein kinase) candidate resistance gene was amplified from both BACs using redundant primers. The DNA sequence showed nine alleles of the RLK at Rhg1 in the soybean germplasm. Markers designed to detect alleles showed perfect association between allele 1 and resistance to soybean cyst nematode Hg types 0 in three segregating populations, fifteen additional selected recombination events and twenty-two Plant Introductions. A quantitative trait nucleotide in the RLK at rhg1 was inferred that alters A47 to V47 in the context of H297 rather than N297. Contiguous DNA sequence of 315 kbp of chromosome 18 (about 2 cM) contained additional gene candidates that may modulate resistance to other Hg-types including a variant laccase, a hydrogen-sodium ion antiport and two proteins of unknown function. A molecular basis for recessive and co-dominant resistance that involves interactions among paralagous disease-resistance genes was inferred that would improve methods for developing new nematode-resistant soybean cultivars.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Linkage disequilibrium at the APA insecticidal seed protein locus of common bean (<Emphasis Type="Italic">Phaseolus vulgaris</Emphasis> L.)

Background  

An interesting seed protein family with a role in preventing insect herbivory is the multi-gene, APA family encoding the α-amylase inhibitor, phytohemagglutinin and arcelin proteins of common bean (Phaseolus vulgaris). Variability for this gene family exists and has been exploited to breed for insect resistance. For example, the arcelin locus has been successfully transferred from wild to cultivated common bean genotypes to provide resistance against the bruchid species Zabrotes subfasciatus although the process has been hampered by a lack of genetic tools for and understanding about the locus. In this study, we analyzed linkage disequilibrium (LD) between microsatellite markers at the APA locus and bruchid resistance in a germplasm survey of 105 resistant and susceptible genotypes and compared this with LD in other parts of the genome.